The E3 ubiquitin-protein ligase MDM2 is a novel interactor of the von Hippel–Lindau tumor suppressor

Mutations of the von Hippel–Lindau (pVHL) tumor suppressor are causative of a familiar predisposition to develop different types of cancer. pVHL is mainly known for its role in regulating hypoxia-inducible factor 1 α (HIF-1α) degradation, thus modulating the hypoxia response. There are different pVHL isoforms, including pVHL30 and pVHL19. However, little is known about isoform-specific functions and protein–protein interactions. Integrating in silico predictions with in vitro and in vivo assays, we describe a novel interaction between pVHL and mouse double minute 2 homolog (MDM2). We found that pVHL30, and not pVHL19, forms a complex with MDM2, and that the N-terminal acidic tail of pVHL30 is required for its association with MDM2. Further, we demonstrate that an intrinsically disordered region upstream of the tetramerization domain of MDM2 is responsible for its isoform-specific association with pVHL30. This region is highly conserved in higher mammals, including primates, similarly to what has been already shown for the N-terminal tail of pVHL30. Finally, we show that overexpression of pVHL30 and MDM2 together reduces cell metabolic activity and necrosis, suggesting a synergistic effect of these E3 ubiquitin ligases. Collectively, our data show an isoform-specific interaction of pVHL with MDM2, suggesting an interplay between these two E3 ubiquitin ligases.

functions 9,22 . Conversely, rather than playing a role as tumor suppressor, pVHL172 contributes to renal carcinoma by up-regulating a subset of pro-tumorigenic genes 23 , i.e. Tgfb1, Mmp1, and Mmp13. pVHL30 was proposed to specifically interact with p14ARF 22 , suggesting a further connection between oxygen sensing and apoptosis. Although these observations suggest that the different pVHL isoforms have redundant as well as isoform-specific functions in the cells, little is known about pVHL isoform-specific protein-protein structural and functional interactions. Here, we report a novel direct interaction between mouse double minute 2 homolog (MDM2) and pVHL30. By undertaking a complementary approach, spanning from in silico analyses to yeast two-hybrid and mammalian cell assays, we defined the molecular details driving this novel interaction. Finally, we showed that a short intrinsically disordered segment of MDM2 is responsible for its association with pVHL in an isoformspecific fashion restricted to pVHL30.

Results
In silico protein-protein interacting network analysis identifies MDM2 as a novel interactor of pVHL. We previously documented a binary interaction between pVHL30 and p14ARF sustained by the acidic N-terminal tail of pVHL30, suggesting a functional asymmetry among pVHL isoforms 22 . p14ARF tumor suppressor is well known for its role in regulating p53 stability 24 , sequestering MDM2 into a subnuclear compartment 25 , while inhibiting p53 degradation 24 . p53 is also stabilized through pVHL 26 . In Zebrafish MDM2 degradation was regulated upon pVHL interaction with programmed cell death 5 protein (PDCD5), a negative regulator of MDM2. This evidence suggests a direct functional effect of pVHL loss of function on the p53/ MDM2 pathway. To address this question, we used STRING to generate an interaction network centered on pVHL, HIF-1α, p14ARF, p53, and MDM2 27 . The resulting network was composed of twenty nodes connected by 52 edges, each of them representing an experimentally validated interaction (Fig. 1). The network is significantly more connected than expected by chance (p-value < 1.0e −16 ), implying that these proteins are functionally interconnected. We asked whether non-trivial functional subnetworks can be linked to specific cell functions. We identified three distinct clusters formed by proteins linked to apoptosis, protein ubiquitination and hypoxia response. The first cluster was formed by six proteins, including p53, MDM2, and p14ARF. Gene ontology terms enrichment predicts this cluster to be mainly localized in the nucleus and participating in the regulation of transcription in response to hypoxia and cellular stress. The second cluster was composed of seven proteins, including pVHL, and, as expected, it was linked to protein degradation. The last cluster was the most heterogeneous and included proteins associated with both hypoxia-and stress-response, e.g. nutrient deprivation and viral aggression. Of note, we found a number of experimentally validated interactions connecting all of these three clusters, suggesting that the hypoxia response is functionally interconnected with cell cycle regulation. Manually curated enrichment of protein-protein interactions across network nodes further extended the connection between clusters I and II. Considering the concept of "interactors of interactors", a well-established property of protein-protein interaction networks 28,29 and the already validated physical interaction between pVHL, p53, and p14ARF, our in silico analysis predicts that MDM2 is a novel direct interactor of pVHL.
pVHL30 specifically interacts with MDM2. To validate the in silico-resolved pVHL association with MDM2, we used yeast two-hybrid (Y2H) assays. We tested whether pVHL binds to MDM2 in an isoform-specific fashion. To address this question, we generated vectors to express human MDM2 fused with Gal4-binding domain (Gal4-BD) as bait, and human pVHL30 and pVHL19 fused with Gal4-activation domain (Gal4-AD) as prey. Protein expression in yeast cells was verified by Western blotting ( Figure S1). Transient expression of pVHL30, pVHL19, and MDM2 alone did not result in yeast growth in selective medium. On the other hand, coexpression of pVHL30 and MDM2 together resulted in increased yeast cell growth in selective medium, thereby . Although Y2H data cannot be considered quantitative per se, these results suggest a weak interaction, as the yeast growth rate was less sustained compared to the positive control. Importantly, co-expression of pVHL19 and MDM2 did not result in cell growth in selective medium, indicating that pVHL19 is unable to associate with MDM2 under these experimental conditions. Next, we asked whether pVHL binds to MDM2 in an isoform-specific manner in mammalian cells. To test this hypothesis, we used human embryo kidney (HEK) 293 T cells transiently transfected with vectors expressing either HA-tagged pVHL30 (pVHL30-HA) or pVHL19-HA with Flag-tagged MDM2. Using immunoprecipitation and Western blotting analysis, we found that pVHL30 forms a complex with MDM2 ( Fig. 2B and Figure S3). To determine whether pVHL30 binds endogenous MDM2 in an isoform-specific fashion, we generated vectors to express green fluorescent protein (GFP)-tagged pVHL30 (pVHL30-GFP) and pVHL19-GFP. Overexpression of pVHL30-GFP, and not pVHL19-GFP, in motor neuron-neuroblastoma hybrid MN-1 cells resulted in the pull down of endogenous MDM2 ( Fig. 2C and Figure S4). Finally, we asked whether endogenous pVHL forms a complex with endogenous MDM2 (Fig. 2D and Figure S5). Immunoprecipitation of endogenous pVHL resulted in the pull down of endogenous MDM2. Notably, reciprocal immunoprecipitation of MDM2 specifically pulled down the pVHL isoform with the higher molecular weight, which corresponds to pVHL30. These observations indicate that pVHL30 specifically forms a complex with MDM2 in yeast as well as human-and mouse-derived mammalian cells. Moreover, our results suggest a functional role of the pVHL30 N-terminal tail in driving the interaction with MDM2.
MDM2 C-terminal domain binds to pVHL30. Next, we sought to determine which region of MDM2 specifically binds pVHL30. We generated vectors expressing MDM2 N-terminus (1-150), central acidic domain (151-350), and C-terminus (351-491) (Fig. 3A). By Western blotting we confirmed the expression of truncated MDM2 proteins in yeast ( Figure S1). Yeast cells co-expressing pVHL30 together with the N-terminus or central acidic domain of MDM2 did not grow in selective medium, indicating no binding ( Fig. 3B and Figure S6). Rather, yeast cells co-expressing pVHL30 and the MDM2 C-terminus grew in selective medium. These results indicate that the last 140 C-terminal residues of MDM2 (MDM2-CT, aa 351-491) are key for the pVHL30/ Analysis of the MDM2 primary structure with FELLS revealed hallmarks of an intrinsically disordered region spanning residues 329-433, terminating before the C3HC4-type zinc finger domain (RING, aa 434-491), which mediates MDM2 protein-protein and protein-RNA interactions [30][31][32] (Fig. 3C). By a Y2H assay, cells coexpressing both RING domain and pVHL30 did not grow in selective medium, indicating that the MDM2 RING domain does not associate with pVHL30 (Fig. 3D,E and Figure S7). Yeast cells co-expressing pVHL30 and the MDM2 351-433 fragment also failed to grow in selective medium. Elongation of the C-terminal end of this fragment to include the first part of the RING domain (MDM2 351-451) restored the interaction with pVHL30, despite the overall disruption of the RING domain structure. Collectively, these findings suggest that at least a portion of the RING domain mediates or modulates MDM2-pVHL30 association. We thus wondered whether the MDM2 disordered region (aa 329-433) mediates the interaction with pVHL30, possibly by an additional and/or cooperative mechanism within the MDM2 C-terminus. Considering that the pVHL30 N-terminal tail, which is specifically involved in this interaction, is mostly acidic and characterized by the presence of eight GxEEx repetitions 22 , we investigated whether clusters of positively charged residues within the MDM2 disordered region account for pVHL30 binding. In silico characterization with FELLS identified a region spanning residues 329-351 as a good candidate for pVHL30 binding (Fig. 3C). This region is predicted as fully disordered, with a typical alternation of charged residues. Of note, the position of a stretch of five lysine residues suggests these to form a positively charged cluster, which we hypothesize to form electrostatic interactions with the acidic pVHL30 N-terminus. In Y2H assay, co-expression of pVHL30 and the MDM2 329-433 fragment resulted in cell growth in selective medium, indicating that residues 329-433 The pVHL30 N-terminus and β-domain are required for interaction with MDM2. pVHL has multiple binding surfaces 33 , each responsible for specific protein-protein interactions 34 . To identify the pVHL30 surface required for MDM2 binding, we generated yeast vectors expressing the pVHL N-terminus (aa 1-53), β-region (aa 54-157), and C-terminal α-domain (aa 158-213), and we performed Y2H assays ( Figure S1). The N-terminus and β-and α-regions did not associate with MDM2 full length ( Figure S8) and MDM2 329-433 fragment (Fig. 4A,B and Figure S9). These findings show that the pVHL30 N-terminal tail, although necessary ( Fig. 2B), is not sufficient for MDM2 binding. Rather, our data suggest that pVHL30/MDM2 interaction is cooperatively driven by concurrent participation of the intrinsically disordered N-terminal tail and at least another pVHL30 domain. Consistent with this hypothesis, yeast cells expressing the pVHL30-N-β fragment (1-157) containing both the acid tail and ß-domain together with MDM2 329-433 fragment grew in selective medium, indicating that these portions of pVHL30 and MDM2 are critical binding surfaces required for such interaction (Fig. 4C,D). Taken together, this experimental evidence supports the idea that pVHL30/MDM2 interaction occurs through the MDM2 C-terminal region (aa 329-491) combining the binding of multiple interaction motifs within the 350-433 region to the pVHL-β domain. We propose that the pVHL30/MDM2 interaction relies on two factors, on one side the electrostatic interactions of the positively charged MDM2 residues (localized in the 329-350 stretch) with the negatively charged pVHL-N-terminus, and on the other side the positive contribution of (or part of) the RING domain.
Ab initio structure prediction and intrinsically disorder flavors characterization of MDM2 fragments. Limited structural information is available about the MDM2 region mediating the association with pVHL30 as no crystal structure is currently available for this region. Since intrinsically disordered regions, albeit lacking a fixed tertiary structure, are often associated with protein cellular function(s) depending on chemical property intrinsic to their primary structure [35][36][37] , we investigated whether structural elements in MDM2 explain the different binding behaviors of MDM2 regions. To this aim, we performed ab initio structure prediction and disorder flavors characterization of the different MDM2 fragments tested in yeast. All predictions   5A and Figure S10). The MDM2 351-451 fragment, containing the minimal sequences required for association with pVHL30, is characterized by a relatively high content of secondary structure. A small β-hairpin followed by a variable number of alpha bulges seems to be shared among different predictions. The β-hairpin is lost in predictions of MDM2 329-433 (which binds pVHL30), while alpha bulges are conserved. MDM2 329-451 is characterized by unfolded-like conformations albeit presenting few conserved secondary structure elements. On the other hand, the shorter MDM2 351-433 fragment, which lacks both flanking regions promoting MDM2-pVHL30 binding, is predicted as the most compact. Average secondary structure content, calculated upon 5000 3D structure predictions for each fragment, indicated that alpha bulges are conserved in both interacting and not interacting MDM2 fragments. These bulges mostly localize in the stretch 351-433, with flaking regions more frequently predicted as coil (Fig. 5B). Detailed characterization of residue-residue interactions for each model performed with RING 38 shows that the MDM2 351-433 (unable to bind pVHL30) is characterized by the highest hydrogen bonds, Van der Waals and ionic interactions (Fig. 5C). This finding supports the concept that this region, when isolated from flanking regions, acquires a more globular-like conformation. Conversely, on average the other fragments have comparable hydrogen bonds paired with a lower number of Van der Waals and ionic interactions, thus www.nature.com/scientificreports/ being relatively rich in secondary structure and more prone to adopt extended conformations. Collectively, all of these findings suggest that shortening MDM2 329-451, by removing both N-and C-flanking regions, stabilizes the polypeptide, thus lowering the intrinsically disordered behavior and inducing a conformational switch to order. Of note, multiple sequence alignment shows that the positively charged MDM2 329-351 region is maintained throughout evolution and is mostly conserved in Artiodactyla and Primates ( Figure S11). As the number of acidic repeats forming pVHL30 N-terminal tail increases from rodents to primates 22 , data suggest that the isoform-dependent MDM2/pVHL30 association could have specifically evolved in higher mammals. We further observed that the cluster of positive charges in region 329-351 appears counterbalanced by an opposed concentration of negative charges in region 410-433 ( Figure S11), supporting the hypothesis that the two flanking sequences interact with the pVHL30 acidic N-terminus and β-domain, respectively, as pVHL β-domain presents a number of conserved positive residues mediating multiple protein-protein interactions 10,28,33 . Notably, MDM2 residues spanning 351-433 fragment includes six ATM (ataxia-telangiectasia mutated kinase) phosphorylation sites (i.e. S386, S395, S407, T419, S425, S429, Figure S11), strongly suggesting that post-translational modifications may operate as regulator of MDM2/pVHL30 association, e.g. introducing electrostatic repulsion against the pVHL30 acidic tail upon phosphorylation.
pVHL30 and MDM2 reduce cell metabolic activity and necrosis. To gain insights into the biological relevance of pVHL30/MDM2 interaction, we performed an MTT assay. Cell metabolic activity correlates with the rate of reduction of the tetrazolium salt to a purple formazan precipitate that occurs in metabolically active cells. We transiently transfected MN-1 cells with vectors expressing HA-tagged pVHL30 and Flag-tagged MDM2, and we verified exogenous protein expression by Western blotting (Fig. 6A). By MTT assay, we found that overexpression of either pVHL30 or MDM2 did not modify cell metabolism (Fig. 6B). Rather, co-expression of both E3 ubiquitin-protein ligases significantly reduced cell metabolic activity by 19%. Then, to evaluate whether the overexpression of HA-tagged pVHL30 and Flag-tagged MDM2 affects cell death, we analyzed apoptosis by flow cytometry using Annexin V and propidium iodide (Fig. 6C). We found no differences in early and late apoptosis as well as necrosis in the cells grown in normal conditions (data not shown). Nutrient deprivation for 24 h significantly reduced cell necrosis upon overexpression of both E3 ubiquitin ligases, suggesting a cooperative and synergistic effect of these E3 ubiquitin-protein ligases on cell metabolism and survival (Fig. 6C).
In summary, our results show an isoform-specific pVHL protein-protein structural and functional interaction. Together with our previous finding of pVHL30-p14ARF specific interaction 9,22 , these observations support the idea that pVHL30 and the other isoforms have non-redundant functions in the cells.

Discussion
Functional intersections between oxygen sensing, apoptosis and cell cycle regulation are well described in the literature 26,28,39,40 , although their functional inter-connection and mutual regulation are still debated. In particular, under prolonged hypoxia, p53 has been reported to accumulate in the cell, resulting in the direct repression of HIF-1α transcriptional activity 41 . A ternary complex between p53, HIF-1α, and MDM2 was proposed to promote pVHL-independent degradation of HIF-1α 39 . pVHL has been shown to directly associate with p53, thereby inducing both suppression of MDM2-mediated ubiquitination and nuclear export of p53 26 . Degradation of p53 is also regulated by p14ARF tumor suppressor 24 , which inhibits the MDM2 E3 ubiquitin-protein ligase activity, leading to p53 accumulation and apoptosis 42 . Here, we describe a novel direct interaction between pVHL and MDM2. This association is pVHL30-isoform-specific and was observed in yeast and mammalian cells. Importantly, we provide evidence that this interaction occurs through the intrinsically disordered N-terminal tail of pVHL30, which indeed is not present in pVHL19. Notably, this acidic tail is necessary, but not sufficient for binding to MDM2, and it requires the pVHL β domain. This domain is known to engage a number of different protein-protein associations 34 , and consistent with this idea it harbors the HIF-1α binding site. HIF-1α associates www.nature.com/scientificreports/ with MDM2 39 , an interaction thought to play a role in modulating hypoxia-induced p53 activation. Our data, while pointing out a novel interaction between MDM2 and pVHL30, suggest that pVHL30 and pVHL19 are characterized by different functional specializations. We speculate that transient macromolecular complexes composed of proteins participating in the regulation of both hypoxia sensing and cell cycle regulation may constantly be formed and destroyed in different cell compartments. Ideally, under mild hypoxia conditions the increased amount of free pVHL may trigger a pVHL30 functional relocation and MDM2 association. Of note, the MDM2 region that we found to participate in this association is immediately downstream of two overlapping binding sites for p14ARF (aa 210-304) and USP7 (aa 242-331), a deubiquitinating enzyme involved in p53 stabilization 43 . Within the MDM2 351-433 region, multiple phosphorylation sites for the ATM kinase are found to be relevant for activation of the DNA damage response 44 . Our data show that the pVHL30/MDM2 interaction is sustained by at least two linear interacting peptides (LIPs) 45 present on the MDM2 ab initio models and putatively modulated by either indirect (e.g. variation of oxygen concentration) and/or direct (e.g. sitespecific phosphorylation) mechanisms. We show that overexpression of both pVHL30 and MDM2 reduces cell metabolic activity and, upon nutrient deprivation cell necrosis, suggesting activation of a specific metabolic reprogramming to escape stress condition. During the early cancer development, when the growing mass is still characterized by insufficient vascularization and diffused hypoxia, a similar behavior may lead to the selection of aggressive clones from heterogeneous tumor cells and promote a malignant progression 46 . Notably, during the revision of this manuscript, pVHL and MDM2 were shown to interact, leading to pVHL neddylation within a region crucial for the p53-pVHL interaction, preventing this association and induction of antiangiogenic factors 47 . Consistent with these observations, here we show for the first time that MDM2/pVHL interaction is isoform-specific. Thus, while pVHL30 is neddylated by MDM2, pVHL19 may not be a target of MDM2, further supporting isoform-specific functions, as we have previously proposed 22,28,48 . In a growing cancer environment, frequently characterized by heterogeneous clone population, the impact of hypoxic and nutritional stresses is amplified by both chaotic and insufficient neovasculature 49 . Considered this complex scenario, it is possible that different pVHL isoforms drive multiple modulation of the MDM2/p53/HIF-1α axis. On one hand, data about pVHL30/MDM2 interaction point to a starvation-induced adaptative role and modulation of angiogenesis. On the other hand, during sustained hypoxia pVHL19 can continue suppressing MDM2-mediated ubiquitination of p53, thus inducing cell cycle arrest and apoptosis 26 . In the context of transient complexes, the observations that a MDM2 motif driving its association with pVHL30 is immediately adjacent to other binding sites and can be modified by PTM is intriguing. It has been proposed that the interaction of pVHL30 with p14ARF induces its nuclear-to-nucleolar translocation, with concurrent dissociation of pVHL30 from the pVHL-elongin C-elongin B (VCB) complex 9 . A similar mechanism also triggers p14ARF-dependent functional inhibition of MDM2 42 . VCB-unbound pVHL30 can interact and stabilize p53 through its α-domain, by suppressing Mdm2-mediated ubiquitination as well as nuclear export 26,34 . Here, we demonstrated that both the N-terminal tail and β-domain are necessary for interaction with MDM2. Such kind of association involves different regions of the protein and is common among linear motif-mediated protein-protein interactions 50,51 . When bound to MDM2, pVHL30 may conceal ATM-specific phosphorylation sites yielding a further layer of regulation of MDM2 activity. Indeed, it has been already demonstrated that ATM induces p53 accumulation and activation by regulating MDM2 oligomerization 44 , while pVHL associates with ATM impairing MDM2-mediated ubiquitination and p53 nuclear export 26 . Acting as a protein-hub, pVHL30 may concurrently form alternative ternary complexes under different conditions. The presence of an USP7 binding site in the same MDM2 region suggests that deubiquitination of several target proteins, such as MDM2 and p53, promotes their stabilization 43,52 . Consistent with this idea, other proteins can mediate USP7 association with p53 53 . In our recent investigation 22 , we found two putative USP7 binding sites on the pVHL30 N-terminus (residues [35][36][37][38][39][40][41][42][43][44]. Collectively, we propose that in addition to suggesting a novel HIF-1α independent function, the pVHL30 association with MDM2 may play a novel role in MDM2 stabilization and mutual regulation between oxygen sensing and apoptosis.

Methods
Bioinformatics analysis. A network of interacting proteins was built with Cytoscape 54 retrieving interaction data from STRING 27 (medium confidence 0.400), imposing no more than 5 interactors for the first shell and 10 for second shell. To reduce the number of false positive, text mining derived data were manually removed. The final network is composed of 20 nodes connected by 52 edges. The network shows an average node degree of 5.2 and an average local clustering coefficient of 0.721, with expected number of edges of 25. Human MDM2 and pVHL sequences were retrieved from UniProt 55 and visualized with Jalview 56 . Orthologous sequences were retrieved from OMA browser and aligned with a T-Coffee 57 (default parameters). Sequence features and secondary structure content were investigated with FELLS 58 , functional domains mapped on the sequence using InterPro 59 and Pfam 60 . Linear motifs were predicted with ELM 61 and intrinsic disorder with MobiDB 3.0 62 . Ab initio protein structures prediction was performed with Rosetta 3.8 63 using specific protocol to model intrinsically disordered regions 64 . For each fragment 5000 decoys were generated and clustered using Rosetta clustering module. DSSP 65 was used to calculate secondary structure content for each 3D structure model, while networks of interacting residues were predicted with RING 2.0 38 .
Plasmid constructs. The recombinant plasmids pGADVHL30, pGADVHL-Nt, pGADVHL-β and pcDNA3.1 HA-VHL30 were constructed as previously described 28 . Genes encoding the human MDM2 and pVHL30 proteins were obtained from commercial plasmid provided by GenScript (GenEZ plasmid OHu28568 and OHu23297) and cDNA transferred into pGBKT7 and pGADT7 plasmids (Clontech) to perform yeast twohybrid assays (Y2H). PCR was used to amplify cDNAs using primers showed in www.nature.com/scientificreports/ vectors. PCR obtained were cloned into pGADT7 or pGBKT7 linearized using EcoRI enzyme, following the In-Fusion HD Cloning Kit (Clontech) protocol. All Y2H recombinant plasmids constructed (Table S1) express the proteins in fusion with the Gal4 DNA binding domain (DBD) or the Gal4 DNA activation domain (AD) expressed by pGBKT7 and pGADT7, respectively. Chimeric proteins in yeast cells were revealed using anti HA and cMyc antibodies. GeneScript plasmid OHu28568 was directly used to overexpress MDM2 in mammalian cells, as the cDNA 5′-end is associated with the FLAG epitope. For Co-immunoprecipitation experiments, the pVHL30 and pVHL19 cDNAs sequences were cloned in the BamHI/EcoRI sites of the pcDNA3.1 vector, plus a N-terminal HA epitope sequence. The recombinant plasmids were used to overexpress pVHL30 and pVHL19 as N-terminal HA-tagged proteins. Sanger sequencing was used to verify all cloned sequences.
Yeast two-hybrid (Y2H) assays. The Matchmaker Gold Two Hybrid System (Clontech) was used to detect direct interaction between MDM2 and pVHL (default conditions). Positive interactions were revealed on multiple conditions by using selective medium lacking histidine and supplemented with 30, 60 or 90 mM 3-AT. Y190 cells were co-transformed with pGBKT7 and pGADT7 empty plasmids as negative control, whereas positive control expresses pGADT7-T (Gal4 AD-SV40 large T-antigen) and pGBKT7-53 (Gal4 DBD-murine p53, fragment 72-390). To exclude erroneous interpretation (i.e. auto-activation), Yeast Y190 cells co-expressing both fusion protein with Gal4 DBD and AD domain alone (and vice-versa) were tested. Expression of the Gal4fusion proteins in yeast cells has been checked by Western blot analysis ( Figure S1). For each transformation three independent clones were serially diluted (tenfold) and spotted either on solid selective medium lacking histidine (30 mM 3-AT and 60 mM 3-AT) or permissive medium, to check cell growth. Yeast strain growth was constantly monitored for 3-8 days at 30 °C. Total yeast proteins have been obtained by TCA-based solubilization of yeast cells 66 , followed by standard Western Blot analysis 67 , using either anti-HA (Abcam, ab16918), or anti-Myc (Abcam, ab127421) antibodies, to reveal the Gal4AD-pVHL or the Gal4BD-MDM2, respectively. Binding of two specific proteins is evaluated by co-expressing the proteins fused with either DNA-binding domain (bait), or the activation domain (prey), of the Gal4 transcription factor. Physical association between bait and prey allows activation of the reporter gene leading to the yeast cell growth on selective medium lacking histidine.
MTT assay. Motor neuron-neuroblastoma hybrid MN-1 cells were seeded on 24-well plate and transfected at 70-80% confluence using Lipofectamine 2000 following manufacturer protocols (Invitrogen). After 48 h, 50 µl of MTT (Thiazolyl Blue Tetrazolium Bromide, code M5655 Sigma-5 mg/ml in PBS1X) was added to each well and the plate incubated for 10′ at 37 °C. Medium was eliminated, and DMSO was added to solubilize the water-insoluble MTT formazan crystals. The absorbance was recorded at 570 nm with spectrophotometer Nan-odrop2000 (Thermo Fisher). www.nature.com/scientificreports/

Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.