Interactome Mapping of the Phosphatidylinositol 3-Kinase-Mammalian Target of Rapamycin Pathway Identifies Deformed Epidermal Autoregulatory Factor-1 as a New Glycogen Synthase Kinase-3 Interactor*

The phosphatidylinositol 3-kinase-mammalian target of rapamycin (PI3K-mTOR) pathway plays pivotal roles in cell survival, growth, and proliferation downstream of growth factors. Its perturbations are associated with cancer progression, type 2 diabetes, and neurological disorders. To better understand the mechanisms of action and regulation of this pathway, we initiated a large scale yeast two-hybrid screen for 33 components of the PI3K-mTOR pathway. Identification of 67 new interactions was followed by validation by co-affinity purification and exhaustive literature curation of existing information. We provide a nearly complete, functionally annotated interactome of 802 interactions for the PI3K-mTOR pathway. Our screen revealed a predominant place for glycogen synthase kinase-3 (GSK3) A and B and the AMP-activated protein kinase. In particular, we identified the deformed epidermal autoregulatory factor-1 (DEAF1) transcription factor as an interactor and in vitro substrate of GSK3A and GSK3B. Moreover, GSK3 inhibitors increased DEAF1 transcriptional activity on the 5-HT1A serotonin receptor promoter. We propose that DEAF1 may represent a therapeutic target of lithium and other GSK3 inhibitors used in bipolar disease and depression.

The phosphatidylinositol 3-kinases (PI3Ks) 1 are a conserved family of lipid kinases that phosphorylate phosphatidylinositol on phosphoinositide 3Ј-hydroxyl group, forming second messengers. Class IA PI3K (hereafter "PI3K") is activated by growth factor receptor tyrosine kinases, and upon stimulation by insulin or IGF1, PI3K triggers the formation of phosphatidylinositol 3,4,5-trisphosphate, which recruits and activates kinases such as AKT and PDK1, mediating most effects of insulin and IGF1 on cell metabolism, growth, proliferation, and differentiation (1,2). Downstream of AKT, the mammalian target of rapamycin (mTOR) kinase is an essential activator of protein synthesis, promoting cell growth and proliferation (1,3,4). mTOR is regulated by growth factors through AKT, by energy availability through the AMP-activated kinase (AMPK), and by amino acid content through class III PI3K. Glycogen synthase kinase-3 (GSK3) is another major target of the PI3K pathway, and its inhibitory phosphorylation by AKT relieves its negative impact on cell cycle progression and cell growth (5). The PI3K-mTOR pathway is central for cell metabolism and proliferation, and its perturbation is implicated in many human diseases (1,3,4,6). Mutations leading to PI3K-mTOR pathway activation are important steps in the initiation and progression of tumors and are frequently encountered in human cancers. In contrast, downregulation of the PI3K pathway impairs cell responses to insulin, leading to type 2 diabetes. Perturbations of the PI3K-mTOR pathway are also linked to muscle atrophy and autoimmune and cardiovascular diseases. Moreover, GSK3 dysregulation is associated with mood disorders and Alzheimer disease (7). Several components of the PI3K-mTOR pathway are promising targets for antitumoral, metabolic, and neurological therapies (3,4,7).
The complexity of the PI3K-mTOR pathway necessitates innovative strategies to identify its exact involvement in physiology and pathology and to predict the consequences of its manipulation in therapy. A better comprehension of cell responses to PI3K-mTOR pathway activation may come from the identification of new regulators or effectors of this pathway, and this goal can now be reached via high throughput approaches (8 -11). We conducted a large scale yeast twohybrid screen of 33 components of the PI3K-mTOR pathway. The resulting interactions were supplemented with a manually curated set of literature interactions, providing a comprehensive and annotated interactome for the PI3K-mTOR pathway. Our screen revealed a predominant place for GSK3A, GSK3B, and AMPK and highlights their role in cancer, metabolic diseases, immune response, and neurological disorders. In particular, we characterized a functional interaction of GSK3A and GSK3B with deformed epidermal autoregulatory factor-1 (DEAF1) transcription factor in the serotonergic pathway.

Cloning, Yeast Two-hybrid (Y2H) Screens, and Co-affinity Purification (Co-AP) Experiments
Detailed descriptions of cloning, Y2H screens, and co-AP experiments are available in the supplemental material and methods. Briefly, full-length ORFs for our baits were cloned, using the Gateway technology, as DNA binding domain (DB) and activation domain (AD) expression vectors and transformed in MaV203 and MaV103 yeast strains. DB expression vectors were used for screening with an AD-cDNA library from E10.5 mouse embryo, and both DB and AD expression vectors were used for screening with the hORFeome1.1 library as described previously (9,12,13). The activation of three reporter genes (HIS3, URA3, and lacZ) was assessed in comparison with controls, and an interaction was considered positive if at least two reporter genes were activated. For the cDNA screen, after selection of in-frame cDNAs, interactions were retested using the gap repair technique (14). The gap repair technique was also used to test interactions with GSK3A/GSK3B and RHEB/RHEBL1 constructs. A subset of 40 interactions was further tested in HEK293T cells using full-length ORFs cloned as GST-bait and Myc-prey vectors using the Gateway technology. Experiments were essentially performed as described (9).

Literature-completed Interactome
For each selected PI3K-mTOR component, the Biomolecular Interaction Network Database (BIND) (15), Molecular INTeraction database (MINT) (16), human protein reference database (HPRD) (17), Agile Protein Interaction DataAnalyzer (APID) (18), and PubMed databases were investigated for protein-protein interactions (PPIs). Each PPI was verified in the corresponding study to check the exact reference of the protein and the technique. Only binary interactions were retained, corresponding to the followings techniques: yeast or mammalian two-hybrid experiments, in vitro kinase or other enzymatic assays with relevant controls, in vitro binding of recombinant proteins purified from mammalian cell-free systems, binding of a protein purified from mammalian cell-free systems to a membraneimmobilized protein, crystallography, and surface plasmon resonance analysis. Careful attention was paid to check that the interaction could not be indirect due to a third component. Co-immunoprecipitation and GST pulldown performed in vivo were not retained. However, in vitro kinase assays involving mTOR or STK11/LKB1 kinases were often performed with kinases isolated from cells as a complex (mammalian target of rapamycin complexes 1 and 2 for mTOR and in complex with MO25 and STE20-related adapter protein for STK11). If controls for these specific kinase activities were suitable, interactions were included with the corresponding indication. Interactions involving a purified protein for which the exact isoform could not be determined were not retained. Described interactions mostly involve human proteins but sometimes involve proteins from mouse or other mammals. For each interaction, the PMID reference referring to the study describing the binary interaction is indicated.
Functional annotation of each interactor for its molecular and subcellular functions and its pathway involvement was deduced from NCBI-Gene and PubMed data banks. Proteins were classified according to a limited number of chosen terms rather than to gene ontology annotations that were either much too detailed or not accurate enough according to the considered protein. As a protein can be engaged in several PPIs with different components of the PI3K-mTOR pathway, interactors may be redundant, whereas interactions are not. Therefore, the annotation of interactors is formally referred to as the annotation of interactions. Figs. 2, 3, and 4 were realized using Cytoscape software (19).

Kinase Assay with GSK3, MLK3, and DEAF1 Constructs
Protein Purification-GST-GSK3A, GST-GSK3B, GST-GSK3A-KA, and GST-GSK3B-KA Gateway vectors were transfected in HEK293T cells as described for co-AP experiments. Cells were lysed with co-AP lysis buffer for 30 min on ice. Lysates were sonicated and precleared by centrifugation for 10 min at 14,000 rpm at 4°C. Precleared lysates were incubated with immobilized glutathione beads (Perbio) for 1 h at 4°C. Beads were washed extensively three times with lysis buffer. GST fusion proteins were eluted with PBS, 50 mM glutathione (Perbio) for 15 min at room temperature. Aliquots of eluates were run on a polyacrylamide-SDS gel in parallel to a BSA scale. The gel was stained with a Coomassie solution, and the quantity of GST fusion proteins was estimated. Adjusted amounts of proteins were then analyzed by Western blot using anti-GST antibody (see Fig. 6B, left gel). Myc-DEAF1, Myc-DEAF1-1, Myc-DEAF1-2, Myc-DEAF1-3 Gateway and pCDN3-M2-MLK3 (kindly provided by K. A. Gallo) vectors were transfected in HEK293T cells as described for co-AP experiments. Each well of 6-well plates was lysed with 250 l of IP buffer 1 (300 mM NaCl, 50 mM Tris, pH 7.4, 0.5 mM EDTA, 1% Triton X-100, Complete protease inhibitor mixture (Roche Applied Science)) for 30 min on ice. Lysates were sonicated and precleared by centrifugation for 10 min at 14,000 rpm at 4°C. 250 l of IP buffer 2 (50 mM Tris, pH 7.4, 0.5 mM EDTA, 1% Triton X-100, Complete protease inhibitor mixture (Roche Applied Science)) and mouse anti-Myc 9E10 antibody (Covance) or anti-M2 FLAG (Sigma) were added to the lysates for overnight agitation at 4°C. Protein A-and G-Sepharose beads were added for 1 h at 4°C and then extensively washed with IP buffer 3 (150 mM NaCl, 50 mM Tris, pH 7.4, 0.5 mM EDTA, 1% Triton X-100, Complete protease inhibitor mixture (Roche Applied Science)). All three IP buffers were supplemented with phosphatase inhibitors (200 M NaF, 1 mM sodium orthovanadate, 100 M ␤-glycerophosphate). Finally, beads were resuspended in 1ϫ kinase buffer (25 mM Hepes, pH 7.4, 10 mM MgCl 2 ), and 10 l of beads suspensions were run on a polyacrylamide-SDS gel, transferred to PVDF membranes, and analyzed by Western blotting with anti-Myc (see Fig. 6C) and -M2 (see Fig. 6B, right gel) antibodies to estimate the quantity of Myc-DEAF1 and M2-MLK3 proteins.
Kinase Assay-Similar amounts of GST-GSK3 proteins (around 70 ng) were added to Myc-DEAF1 or M2-MLK3 recombinant proteins (around 50 ng). The mixture was adjusted to final concentrations of 25 mM Hepes, pH 7.4 and 10 mM MgCl 2 . 30Ci of [␥-32 P]ATP were added to the mixture. After incubation for 30 min at 30°C, the reaction was stopped by addition of 5ϫ loading buffer and denaturation at 95°C for 5 min. Aliquots run on polyacrylamide-SDS gels and transferred to PVDF membranes were revealed using a phosphorimaging system (Fuji) for radioactivity measures.
Quantification of Myc-DEAF1 Constructs-The membranes were then subjected to Western blotting with mouse anti-Myc (9E10) and goat anti-mouse IRDye 800 secondary antibody (LI-COR) to reveal and quantify Myc-DEAF1 proteins on an Odyssey LI-COR scanner (ScienceTec). Adjusted radioactive signal from Myc-DEAF1 proteins was calculated by the ratio radioactive signal/protein quantity.
Transfections and Drug Treatment-HEK293T cells were grown at 37°C in 5% CO 2 in Dulbecco's modified Eagle's medium supplemented with 10% fetal bovine serum. Cells were seeded in 12-well plates and transfected using calcium chloride. Each well was transfected with 200 ng of 5-HT1A(C) vector; 50 ng of pCMV-SPORT-␤gal; 150 ng of GST, Deaf1-Myc, or Myc vector; and 400 ng of pcDNA3 to complete the DNA mixture. Drugs were added at transfection. Lithium and sodium chloride were added at a final concentration of 5 mM. Azakenpaullone (Calbiochem) was diluted in DMSO and added at a final concentration of 1 M or 500 nM. An equivalent volume of DMSO (0.25 l) was added in control wells. Medium was changed 12 h after transfection, and drugs were replaced at the same concentrations.
RN46A cells (kindly provided by Dr. Scott Whittemore's laboratory) were grown at 33°C in 5% CO 2 in Neurobasal medium (Invitrogen) supplemented with 10% fetal bovine serum and glutamine. Cells were seeded in 12-well plates. The transfection mixture was composed of 100 l of Opti-MEM (Invitrogen); 4 l of FuGENE HD (Roche Applied Science); 200 ng of 5-HT1A(C) vector; 50 ng of CMV-GFP vector; 150 ng of GST, Deaf1-Myc, or Myc vector; and 400 ng of pcDNA3. Drugs were added in the same conditions as for HEK293T cells except that that the medium was not changed after transfection.
Luciferase Assay-Cells were harvested 36 h after transfection. HEK293T cells were lysed with 250 l of Passive Lysis Buffer (Promega). 10 l of lysate were used for luciferase activity detection using the Luciferase Assay System (Promega). 50 l of lysate were used for ␤-galactosidase activity detection and were added to 50 l of detection buffer containing 1 mM MgCl 2 , 50 M ␤-mercaptoethanol, and 40 g of ortho-nitrophenyl-␤-galactoside. After incubation at 37°C, the plate was read at 405 nm.
RN46A cells were dissociated with Versene solution (Invitrogen) and resuspended in 1 ml of PBS. 200 l were used for fluorescenceactivated cell sorting analysis (BD Biosciences) of GFP-transfected cells. The remaining cells were lysed with 100 l of Reporter Lysis Buffer (Promega). 10 l of lysate were used for luciferase activity detection using the Luciferase Assay System (Promega). 10 l of lysate were used for the protein assay using Bio-Rad Protein Assay.
Statistical Analysis-Independent experiments were performed at least six times. Adjusted luciferase activity was corrected for transfection efficiency and cell density by calculating the ratio of luciferase activity/␤-galactosidase activity for HEK293T cells. As ␤-galactosidase activity detection was not sensitive enough for normalization in RN46A cells (data not shown), adjusted luciferase activity was calculated as the ratio of luciferase activity/GFP transfection rate/protein concentration. Statistical analysis was performed using Student's t test.
Selected components were cloned by Gateway recombinational cloning into Y2H destination vectors to generate GAL4 DB or AD fusion proteins (9,22). Insulin receptor, IRS2, IRS4, and SHIP1 failed to clone, leaving 33 constructs for Y2H screening.
Yeast Two-hybrid Screens and Co-affinity Purification Experiments-As the PI3K-mTOR pathway is conserved and ubiquitously expressed, it was important to screen using the broadest possible libraries. Two stringent Y2H screens were carried out, one with an E10.5 mouse embryo library and the other with the human ORFeome v1.1 ORF library (9,22). Both screens followed optimized protocols in which high quality parameters are assessed (14,23,24).
In total, the cDNA and the ORFeome screens identified 68 and eight PPIs, respectively, involving 15 PI3K-mTOR pathway components ( Fig. 2 and supplemental Table I). Two PPIs were found in both screens, showing the complementarity of these screens. Three PPIs have been described to occur in the literature, and four other interactions found with GSK3A had been previously characterized with its close paralog GSK3B, giving 67 of 74 (91%) new potential interactions (supplemental Table II). There were 35 PPIs (47%) found at FIG. 1. Description of PI3K-mTOR pathway. A, functional relationships between the proteins of the PI3K-mTOR pathway and major outputs of the pathway activation. Modules are color-coded: orange, green, pink, and blue correspond to the receptor, the PI3K, the mTOR, and the related protein modules, respectively. Functional paralogs are depicted as joined rectangles. PI, phosphatidylinositol; PIP, phosphatidylinositide phosphate. B, current designations of the PI3K-mTOR pathway components are presented in parallel to their official symbols and to their gene ID. Genes indicated in italic could not be included in our yeast-two hybrid screen but were integrated in the literature-completed interactome. G␤L, G protein ␤ subunit-like; IR, insulin receptor; PIKE, phosphoinositide 3-kinase enhancer; PKB, protein kinase B; PTEN, phosphatase and tensin homolog; IGF1R, insulin-like growth factor 1 receptor. least twice independently ( Fig. 2 and supplemental Table  I). GSK3A and PRKAA1, the AMPK ␣1 catalytic subunit, were the most connected proteins with 17 and 27 PPIs, respectively.
Paralogs in the screens had sharply different numbers of PPIs. We found four PPIs for RHEBL1 but none for RHEB and 17 PPIs for GSK3A but two for GSK3B. Using the gap repair technique (14), we tested the interaction between RHEBL1 interactors and RHEB. Y2H interactions were retrieved with RHEB and RHEBL1 only when their last amino acids that constitute a prenylation motif (25) were removed (supplemental Fig. 1). In this case, RHEB interacted with three of four interactors of RHEBL1 (supplemental Fig. 1, Fig. 2, and supplemental Table I).
To validate Y2H interactions, we tested half (n ϭ 40) of the PPIs by co-AP assays in human HEK293T cells (9). Transfected proteins were expressed in 33 co-APs, and 19 co-APs (58%) tested positive (supplemental Fig. 2 and supplemental Table I), which is in line with expectation from previous such validation efforts (26,27). The protein interactions reported FIG. 2. Yeast-two hybrid interaction network for PI3K-mTOR pathway. Y2H baits are depicted are triangles, and the discovered interactors are depicted as circles. Baits for which no interaction was found are not represented. The color of the links is indicative of the Y2H method that led to its discovery. The thickness of the links is indicative of the number of times an interaction was found in the Y2H cDNA or high throughput (HTP) screen. The molecular function of each protein is described by a color code. GAP, GTPase-activating protein.
Literature-completed Interactome-We manually curated from the literature all binary interactions for each PI3K-mTOR pathway component. Verified interactions were added to our new interactions in a resulting "literature-completed interactome" containing 802 distinct interactions (supplemental Table II). The interactors were annotated for their molecular function and their subcellular functions based on validated data from the literature (supplemental Table II and Fig. 3). Accordingly, 648 interactors were annotated for their participation in particular signaling pathways (supplemental Table II and Fig. 4).
To test whether the kinase activity of GSK3A and GSK3B was necessary for these Y2H interactions, we cloned the kinase-dead forms, GSK3A-KA and GSK3B-KA, in which Lys-148 and Lys-85, respectively, were replaced by an alanine (32,33). To prevent a deleterious phosphorylation on Ser-21 and Ser-9, we also constructed GSK3A-SA-KA and GSK3B-SA-KA forms, respectively, combining both described mutations. The KA mutations abolished most interactions for both kinases independently of the presence of the inactivation site. There was only one interactor for GSK3B-KA and two for GSK3A-KA (Fig. 5B). We tested seven of these interactions by co-AP in HEK293T cells (Fig. 5C and supplemental Table I). GSK3A-SA and GSK3B-SA interacted with all tested GSK3A⌬ interactors; GSK3A-wt and GSK3B-wt interacted with six and four GSK3A⌬ interactors, respectively; and GSK3A-KA and GSK3B-KA interacted with three and four GSK3A⌬ interactors, respectively. Seemingly, the co-AP assay was less sensitive to GSK3 kinase activity than Y2H.
In conclusion, all GSK3A⌬ tested interactors were shared by GSK3B and were favored by constitutive activation of the kinases, whereas kinase inactivation abolished most interactions. Different sets of interactors were observed for wildtype, constitutively active, and inactive forms of GSK3A and GSK3B (Fig. 5, B and C), suggesting specificities for these two kinases.
DEAF1 Is a New GSK3 Substrate-Our Y2H screens identified the DEAF1 (also called NUDR) transcription factor as a new interactor for GSK3A and GSK3B, and these interactions were confirmed by co-AP assay in human cells (Figs. 2 and 5, supplemental Fig. 2, and supplemental Table I).

PI3K-mTOR Interactome
kinase activity affected DEAF1 transcriptional activity. A major transcriptional target for DEAF1 is the 5-HT1A serotonin receptor gene (20,36). We examined DEAF1 transcriptional activity upon GSK3 inhibition with a luciferase reporter under the control of 5-HT1A promoter (20). Transfection of the fulllength Myc-DEAF1 construct into HEK293T cells caused luciferase activity to increase ( Fig. 7A; p Ͻ 0.005). Both lithium and azakenpaullone, two inhibitors of GSK3 acting through different mechanisms (37,38), significantly increased Myc-DEAF1 activity as reflected by increased luciferase activity in HEK293T cells (Fig. 7A; p Ͻ 0.001). Similarly, Myc-DEAF1 transfection into rat RN46A serotonergic cells markedly increased luciferase activity ( Fig. 7B; p Ͻ 0.005), and treatment by lithium or azakenpaullone stimulated Myc-DEAF1 activa- FIG. 6. DEAF1 in vitro phosphorylation by GSK3A and GSK3B kinases. A, DEAF1-FL, DEAF1-1, DEAF1-2, and DEAF1-3 constructs used for kinase assay are represented in parallel to the cDNA clone that was isolated from the Y2H screen. Numbers in italic indicate the first and last amino acids of the constructs in comparison with DEAF1-FL. Conserved residues between human and mouse that are potential sites of phosphorylation by GSK3A and GSK3B are shown on the DEAF1-FL construct. B, quantification of purified wild-type (WT) and kinase-dead (KA) GST-GSK3A and GST-GSK3B as well as M2-MLK3 by Western blot with anti-GST and -M2 antibodies, respectively. C, upper parts, 32 P radioactive signal from kinase assay reveals MLK3 (positive control), DEAF1-FL, DEAF1-1, DEAF1-2, and DEAF1-3 protein phosphorylation by wild-type GSK3A and GSK3B kinases. Lower parts, quantitative fluorescent Western blot for Myc on the kinase assay membranes reveals the respective amount of Myc-DEAF1 proteins. D, adjusted quantification of DEAF1 construct phosphorylation by GSK3A and GSK3B was calculated as the ratio ( 32 P radioactive signal/Myc signal quantification) for each DEAF1 construct. a.u., arbitrary units. tion of luciferase activity. Whereas lithium stimulated luciferase activity in RN46A cells as strongly as in HEK293T cells (p Ͻ 0.001), azakenpaullone only moderately stimulated luciferase activity (p Ͻ 0.05 at 1 M and p Ͻ 0.025 at 500 nM). DEAF1 is broadly expressed (39), and RT-PCR experiments confirmed that DEAF1 was expressed in HEK293T and RN46A cells (data not shown; see "Experimental Procedures"). We tested the effect of GSK3 inhibition upon endogenous DEAF1 activity. In both cell lines, luciferase activity under the control of 5-HT1A promoter was significantly increased by GSK3 inhibitors with the strongest stimulation by lithium in RN46A cells (Fig. 7, C and D; p Ͻ 0.001). In summary, DEAF1 activation of 5-HT1A promoter activity increased upon GSK3 inhibition in HEK293T and RN46A cells, suggesting an inhibitory role of GSK3 phosphorylation on DEAF1.

Protein-Protein Interaction
Mapping for PI3K-mTOR Pathway-Y2H mapping for 33 components of the PI3K-mTOR pathway against two different libraries identified 74 interactions, 67 of which have not been described previously. This screen revealed two highly connected proteins: PRKAA1, the AMPK ␣1 catalytic subunit (27 interactors) and GSK3A (17 interactors). Literature curation provided a set of 802 functionally annotated PPIs for the PI3K-mTOR pathway (literature-completed interactome) that offers a comprehensive picture of the connectivity and the predominant cellular processes involving the PI3K-mTOR pathway (supplemental Data 2). The low overlap between our screen and the interactions from the literature (seven of 74) may be explained by two parameters: sampling sensitivity and assay sensitivity. Sampling sensitivity refers to the interactions that can be identified in a single trial of an assay. Sampling sensitivity was estimated at 45% in a comparable Y2H screen (24), and six screens were needed to reach 90% saturation. Successive iterations of our screen would likely achieve a higher degree of saturation. In parallel, assay sensitivity designates the fraction of all interactions that can be identified by an assay with its specific experimental conditions. In parallel to the high strin- FIG. 7. Activation of 5-HT1A promoter by DEAF1 and GSK3 inhibitors. A and B, results of luciferase assays on HEK293T cells (A) and RN46A cells (B) show the influence of DEAF1 on 5-HT1A promoter activity alone or in addition to GSK3 inhibitors, i. e. lithium chloride and azakenpaullone, in comparison with control conditions (Myc vector, sodium chloride, and dimethyl sulfoxide, respectively). C and D, results of luciferase assays on HEK293T cells (C) and RN46A cells (D) show the influence of GSK3 inhibitors on 5-HT1A promoter activity in comparison with control conditions. Independent experiments were performed at least six times. The luciferase level was adjusted to the transfection conditions (see "Experimental Procedures"). The histograms show means Ϯ S.E. Statistical analysis was performed using Student's t test. *, p Ͻ 0.05; **, p Ͻ 0.025; ***, p Ͻ 0.01; ****, p Ͻ 0.005; *****, p Ͻ 0.001. a.u., arbitrary units.

PI3K-mTOR Interactome
gency of the Y2H criteria used in this study (low copy plasmids and two reporters activated), its intrinsic limitations may include the cDNA and ORF isoforms and the level of representation in the libraries used, the activation status, and the post-translational modifications of the baits in yeast that may be different from the ones encountered in previous screens. Moreover, full-length baits, instead of domains, may have prevented some PPIs (27), and lipid-interacting domains present in several components of the PI3K-mTOR pathway may have impaired the nuclear translocation required in Y2H screening. Interestingly, assay sensitivity for different methods of PPI identification was estimated to a similar value around 20 -35%, and our screen probably lies in the same range (24). Altogether, sampling and assay sensitivities may explain the complementarities of the different approaches used in PPI identification for PI3K-mTOR components. The newly discovered interactions from this study provide insights on GSK3A and GSK3B specificities and suggest new paths for the comprehension of PI3K-mTOR pathway involvement in cancer, metabolism, diabetes, immune response, and neurobiology.
Implications for GSK3A and GSK3B Specificities-GSK3A and GSK3B kinases are highly related paralog kinases that share a large redundancy, although recent studies revealed specific functions for each of them (29, 32, 40 -43). Understanding isoform redundancies and specificities represents an important stake for cell biology and drug development (44). However, GSK3A received little attention compared with GSK3B. Our screens identified 17 common interactors for GSK3A and GSK3B, reinforcing a functional redundancy between them. Y2H interactions were nearly abrogated when kinase-dead forms of GSK3A and GSK3B were used, indicating that these interactors may be substrates for both kinases.
However, all interactors did not interact with the same strength and with the equivalent constructs for each kinase. For example, C14orf129 interacted with C-GSK3A and constitutively active GSK3A but with kinase-dead GSK3B in Y2H. These different affinities for GSK3A and GSK3B, depending on their activation status, may lead in vivo to a preferential interaction with one of these two kinases or to a sequestration by one of them. Beyond an apparent redundancy, the status of GSK3A and GSK3B for their kinase activity may represent a determinant element of their specificities. Characterization of these respective interactions could help understanding the coordinated activities of these kinases and developing focused strategies for targeting specific interactions.
PI3K-mTOR Pathway and Cancer-Many components of the PI3K-mTOR pathway act as oncogenes or tumor suppressor genes, and mutations leading to an activation of the PI3K-mTOR pathway may account for 30% of human cancers (45). Many drugs targeting PI3K, AKT, or mTOR are currently tested for cancer therapy (44). Several interactions discovered in our Y2H may participate in specific aspects of tumor initiation or progression.
Our screen identified a new interaction between TSC1 and MAPK14 (p38␣), one of the four members of the p38 mitogenactivated protein kinase family. p38 mitogen-activated protein kinases regulate inflammatory and stress responses and behave as tumor suppressors (46). The TSC1-TSC2 complex inhibits mTOR, a translational activator promoting cell growth and proliferation, and mutations in TSC1 or TSC2 lead to a tumor syndrome (47). p38 may regulate mTOR activation through TSC1 phosphorylation.
New interactions may link the PI3K-mTOR pathway to the Wnt-␤-catenin pathway, which is also involved in tumorigenesis (48). GSK3 is at the crossroad of these two signaling pathways. Upon activation of Wnt signaling pathway, GSK3, in complex with Axin, phosphorylates ␤-catenin and targets it to degradation (34). In our Y2H screens, GSK3A and GSK3B interacted with LZTS2 and ZBED3. LZTS2 tumor suppressor promotes ␤-catenin nuclear export, impairing its transcriptional activity (49). ZBED3, an Axin-binding protein, inhibits ␤-catenin phosphorylation by GSK3B (50). The interactions among GSK3, LZTS2, and ZBED3 suggest new coordination for ␤-catenin inhibition, ensuring a precise control of this potent oncogene.
AMPK, the "cell energy sensor," inhibits mTOR through the activation of the TSC1-TSC2 complex (51), and AMPK-activating drugs are candidate for cancer therapies (52). Interestingly, LZTS2 was also found as a new interactor for PRKAA1, suggesting a putative coordination between GSK3 and AMPK in ␤-catenin regulation. A similar link between these two kinases already exists in sequential phosphorylation and activation of TSC2 (53). This interaction between AMPK and LZTS2 may reveal another role for AMPK in cancer.
Moreover, AMPK is involved in the control of epithelial polarity, actin organization, and cell proliferation downstream of LKB1 (54 -56). We found a new interaction between PRKAA1 and Abelson interactor protein-1 (ABI1). ABI1 regulates actin polymerization, and its expression correlates with cell migration and invasiveness of cancer cell lines (57,58). ABI1 may thus mediate part of AMPK control of cell polarity and tumor progression.
Finally, we describe a new interaction between PRKAA1 and KIF1C, a plus-end microtubule motor that promotes podosome dynamics (59). KIF1C expression is highly predictive of brain metastasis in lung cancer (60), and immunotherapy against KIF1C appears promising in glioma (61). We suggest that KIF1C interaction with AMPK may modulate its activity in cancer progression.
PI3K-mTOR Pathway, Metabolism, and Diabetes-AMPK is activated upon a high AMP/ATP ratio and turns on ATPgenerating processes while switching off ATP-consuming activities (62). AMPK-activating drugs, such as aminoimidazole carboxamide ribonucleotide and metformin, are used in type 2 diabetes treatment. It was recently shown that AMPK regulates the circadian clock of peripheral organs (63). We describe a new interaction between PRKAA1 and BHLHB2 (or DEC1), a transcription factor regulating the expression of metabolic and clock genes that are required for circadian rhythm (64,65). This interaction may participate in the control of circadian rhythm by AMPK, likely adapting metabolism to animal activity.
Carboxypeptidase E (CPE) was another new interactor of PRKAA1. Mutation of CPE is associated with type 2 diabetes (66 -68). The regulation of CPE activity by AMPK could contribute to the therapeutic benefits of AMPK activation in type 2 diabetes.
GSK3 is another prime target for type 2 diabetes management (48). In our screen, GSK3A and GSK3B interacted with ARHGEF11 (or PDZ-RhoGEF), a small GTPase activator. ARHGEF11 polymorphism is associated with insulin resistance and type 2 diabetes (69 -71). Several small GTPases are implicated in the trafficking of glucose transporter 4 (GLUT4) following insulin stimulation (72). We propose that GSK3 interaction with ARHGEF11 may regulate glucose uptake in response to insulin through the modulation of GLUT4 transport.
Finally, we describe a new interaction between RHEB/ RHEBL1 and APPL1, an adaptor protein that mediates adiponectin-induced sensitization of the insulin pathway (73). APPL1 activates AMPK, leading to S6K inhibition and subsequent release of IRS1 inhibition (see Fig. 1) (74). This likely occurs through RHEB inhibition as RHEB overexpression impairs the APPL1 effect (74). APPL1 direct interaction with RHEB and RHEBL1 may thus reinforce APPL1-mediated sensitization of insulin signaling.
PI3K-mTOR Pathway and Immune Response-The PI3K-mTOR pathway participates in many aspects of the immune response, and rapamycin, an mTOR inhibitor, is used as an immunosuppressive drug (75,76). We found a new interaction between RHEB/RHEBL1, two mTOR activators, and TRAF2, an adaptor protein that lies downstream of TNFR1 and CD40 receptors and mediates the activation of the MAPK and nuclear factor B pathways (77). TRAF2 participates in various aspects of the immune response (78,79), and RHEB negatively regulates MAPK activation through B-Raf inhibition (80). TRAF2 interaction with RHEB and RHEBL1 may represent another possible level of control of immune responses by PI3K downstream of tumor necrosis factor ␣ and CD40 ligand.
In addition, we identified an interaction between PRKAA1 and RC3H1 (or Roquin). In mouse, Roquin mutation impairs the degradation of particular mRNAs, leading to the development of severe autoimmunity (81,82). AMPK may influence Roquin activity and modulate certain aspects of autoimmunity.
PI3K-mTOR Pathway and Neurobiology-During development, GSK3 regulates neurite extension and neuronal architecture (5). In the adult, GSK3 has been implicated in specific behavioral responses (83). In our Y2H screens, GSK3A and GSK3B interacted with ARHGEF11, which increases EAAT4 glutamate transporter activity in brain and participates in neuronal morphogenesis (84,85). GSK3 interaction with ARHGEF11 may explain part of GSK3 involvement in synaptic plasticity. More broadly, we suggest that this interaction may participate in regulated trafficking in different contexts, such as the insulin pathway (see above) and neurogenesis.
Our Y2H screens also found SPG21 (or Maspardin) interacting with GSK3A. SPG21 mutation leads to hereditary spastic paraplegia, a neurodegenerative disease associated with dementia (86). Spastic disorders are also encountered in hereditary cases of Alzheimer disease (87). Importantly, GSK3 is involved in the progression of Alzheimer disease, and chemical inhibition of GSK3 alleviates ␤-amyloid accumulation (7). The interaction between GSK3A and SPG21 may reveal a mechanistic proximity between hereditary spastic paraplegia and Alzheimer disease.
Possible Involvement of DEAF1-GSK3 Interaction in Mood Disorders-Finally, we identified DEAF1 transcription factor as an interactor and an in vitro substrate for GSK3A and GSK3B. GSK3 inhibitors increased DEAF1 transcriptional activity on the 5-HT1A promoter in both HEK293T and RN46A cells, suggesting an inhibitory regulation of DEAF1 by GSK3.
Variations in 5-hydroxytryptamine (serotonin) (5-HT) availability and signaling cascade are associated with psychiatric disorders such as major depression disease. Among several receptors, the 5-HT1A receptor for serotonin is a major target of antidepressant treatments (88 -91). Presynaptic 5-HT1A receptors act as inhibitory autoreceptors in serotonergic neurons located in midbrain raphe nuclei. On the contrary, postsynaptic 5-HT1A receptors are found in 5-HT-responsive areas that are limbic regions and specific cortical layers.
In previous reports, DEAF1 repressed 5-HT1A expression in HEK293T cells and in the serotonergic RN46A cells that express 5-HT1A but enhanced 5-HT1A expression in non-serotonergic 5-HT1A-expressing cells such as SN48, NG108-15, and SKN-SH cells (20,36). However, DEAF1 enhanced 5-HT1A expression in several hippocampal and septal cells (20). These data likely reflect an opposite regulation on presynaptic versus postsynaptic 5-HT1A expression (36) and perhaps a phenotypic derivation of the HEK293T cells and RN46A cells we used. The RN46A cells we used did not express 5-HT1A (data not shown) and may behave as nonserotonergic cells. The stimulation of 5-HT1A expression by DEAF1 that we observed in both cell lines is in accordance with neither of them expressing 5-HT1A.

PI3K-mTOR Interactome
In parallel, lithium has long been used in the treatment of bipolar disease and as an adjuvant to antidepressant drugs (38,83,88). The preeminent mechanism of action of lithium in psychiatric disorders is likely represented by GSK3 inhibition (7,38). Accordingly, lithium behavioral responses are increased by GSK3B haploinsufficiency and reproduced by another GSK3 inhibitor (98), whereas GSK3B haploinsufficiency and GSK3B inhibition alleviates aberrant behaviors due to 5-HT deficiency (99). Other mood-stabilizing agents such as clozapine and valproate also inhibit GSK3 activity (83,100). Conversely, an increased GSK3B activity is observed in the prefrontal cortex of major depression disease subjects (101), whereas both 5-HT1A and DEAF1 were decreased in the prefrontal cortex of depressed women, suggesting a positive regulation of 5-HT1A transcription by DEAF1 in these areas (95). Importantly, treatment with lithium or divalproex increased 5-HT1A expression in bipolar disorder patients (102). Our finding that GSK3 phosphorylates DEAF1 and impairs DEAF1-driven expression of 5-HT1A suggests that abnormal GSK3B activation reduces 5-HT1A expression, contributing to depressive symptoms, whereas GSK3 inhibition by lithium or another drug would increase DEAF1-driven 5-HT1A expression in 5-HT-responsive areas and restore 5-HT transmission. GSK3 and DEAF1 interaction might explain part of GSK3 involvement in bipolar disease and depression, and DEAF1 may represent a therapeutic target of lithium and other GSK3 inhibitors used in these disorders.