Biochemical Large-Scale Interaction Analysis of Murine Olfactory Receptors and Associated Signaling Proteins with Post-Synaptic Density 95, Drosophila Discs Large, Zona-Occludens 1 (PDZ) Domains*

G protein-coupled receptors (GPCRs) constitute the largest family among mammalian membrane proteins and are capable of initiating numerous essential signaling cascades. Various GPCR-mediated pathways are organized into protein microdomains that can be orchestrated and regulated through scaffolding proteins, such as PSD-95/discs-large/ZO1 (PDZ) domain proteins. However, detailed binding characteristics of PDZ–GPCR interactions remain elusive because these interactions seem to be more complex than previously thought. To address this issue, we analyzed binding modalities using our established model system. This system includes the 13 individual PDZ domains of the multiple PDZ domain protein 1 (MUPP1; the largest PDZ protein), a broad range of murine olfactory receptors (a multifaceted gene cluster within the family of GPCRs), and associated olfactory signaling proteins. These proteins were analyzed in a large-scale peptide microarray approach and continuative interaction studies. As a result, we demonstrate that canonical binding motifs were not overrepresented among the interaction partners of MUPP1. Furthermore, C-terminal phosphorylation and distinct amino acid replacements abolished PDZ binding promiscuity. In addition to the described in vitro experiments, we identified new interaction partners within the murine olfactory epithelium using pull-down-based interactomics and could verify the partners through co-immunoprecipitation. In summary, the present study provides important insight into the complexity of the binding characteristics of PDZ–GPCR interactions based on olfactory signaling proteins, which could identify novel clinical targets for GPCR-associated diseases in the future.

G protein-coupled receptors (GPCRs) constitute the largest family among mammalian membrane proteins and are capable of initiating numerous essential signaling cascades. Various GPCR-mediated pathways are organized into protein microdomains that can be orchestrated and regulated through scaffolding proteins, such as PSD-95/ discs-large/ZO1 (PDZ) domain proteins. However, detailed binding characteristics of PDZ-GPCR interactions remain elusive because these interactions seem to be more complex than previously thought. To address this issue, we analyzed binding modalities using our established model system. This system includes the 13 individual PDZ domains of the multiple PDZ domain protein 1 (MUPP1; the largest PDZ protein), a broad range of murine olfactory receptors (a multifaceted gene cluster within the family of GPCRs), and associated olfactory signaling proteins. These proteins were analyzed in a large-scale peptide microarray approach and continuative interaction studies. As a result, we demonstrate that canonical binding motifs were not overrepresented among the interaction partners of MUPP1. Furthermore, C-terminal phosphorylation and distinct amino acid replacements abolished PDZ binding promiscuity. In addition to the described in vitro experiments, we identified new interaction partners within the murine olfactory epithelium using pulldown-based interactomics and could verify the partners through co-immunoprecipitation. In summary, the present study provides important insight into the complexity of the PDZ domain proteins comprise one of the largest families among interaction domain scaffolding proteins and are highly abundant in various multicellular eukaryotic species. These proteins fulfill important physiological functions in a broad range of different tissues and cells as they orchestrate complex protein networks. Among putative PDZ interaction partners, one important protein family is the group of GPCRs 1 , constituting the largest family of membrane proteins in mammals (1). Here, signal efficiency, speed, desensitization, and internalization can be modulated by PDZ proteins (2)(3)(4)(5). Olfactory receptors (ORs) represent a multigene family within this group of seven-transmembrane domain proteins and encompass 2% of the mammalian genome (6). Belonging to class I GPCRs, ORs share many general features of this receptor family, making them an interesting target for interactions involving PDZ proteins. Until recently, an organizing complex builder, such as the inactivation no afterpotential D (InaD) protein in the visual system of Drosophila melanogaster (7,8), could not be clearly identified for olfactory signaling.
The multiple PDZ domain protein 1, with 13 individual PDZ domains, represents the largest of the described PDZ pro-teins to date (9) and interacts with different GPCRs (10 -12). One well-described example is its interaction with GABA B receptors, leading to enhanced receptor stability at the plasma membrane and prolonged signaling duration (2). In previous studies, we demonstrated that PDZ domains 1 ϩ 2 can interact with a selected subset of ORs (13). Furthermore, we showed that MUPP1 binds to a specific OR and that most of the described proteins are involved in mammalian olfactory signal transduction in the native system, making MUPP1 a promising candidate for orchestrating the olfactory system (14).
Many PDZ-ligand interactions depend on classical binding motifs at the ligand's C-terminal end, thereby building weak transient protein complexes (15,16). However, an increasing number of PDZ interactions have emerged that seem to provide more complex binding modalities, differing from the canonical interactions (17,18). Ligand binding seems not to be exclusively restricted to C-terminal sites, and PDZ domains cannot be distinctly classified but are evenly distributed throughout a selective space (17, 19 -21). Therefore, it is of great interest to analyze OR-PDZ interactions to characterize the putative binding requirements and to further investigate the role of MUPP1 in olfactory signaling.
In the present study, we characterized the binding modalities between the 13 individual PDZ domains of MUPP1 and a broad range of murine olfactory receptors in a large-scale approach, indicating that classical binding motifs were not overrepresented among the evaluated binding partners. In addition, we identified new binding partners from the murine olfactory epithelium using pull-down-based interactomics.

EXPERIMENTAL PROCEDURES
DNA Constructs-GST-fusion constructs of the PDZ domains 1-13 of MUPP1 were generated by cloning the single domains into the pGEX-3X expression vector (Amersham Biosciences, Buckinghamshire, UK), and pCDNA3_MUPP1 served as the template. All of the constructs were verified by sequencing. For a detailed primer list, see the previous publication (14).
Immunohistochemistry-Heterozygous OMP-GFP transgenic mice (22) were raised and maintained according to institutional and governmental instructions. Horizontal cryosections were performed on P18 mice with a layer thickness of 14 m. The sections were blocked for 1 h (PBS -, 1% gelatin) at room temperature and incubated with the corresponding antibody solution overnight at 4°C. After several washing steps, the slices were coated with the appropriate fluorophore-coupled secondary antibody for 1 h at room temperature. After embedding the slices with Prolong Antifade Gold (Life Technologies, Carlsbad, CA), fluorescence images were obtained using a confocal microscope (Zeiss LSM 510 Meta) with a 40ϫ oil immersion objective. Control experiments without primary antibodies showed a very low level of background staining (data not shown).
Mass Spectrometry-To reduce unspecific binding events, 50 l of glutathione-covered beads coupled with the GST-tagged PDZ domain were blocked in 100 l of TBS ϩ/ϩ (100 mM Tris-HCl, 150 mM NaCl) containing 1% BSA for 30 min under ice-cold conditions. The entire suspension was centrifuged for 3 min at 1,000 g, the supernatant was subsequently decanted, and the beads were resuspended in ice-cold interaction buffer (50 mM ammonium bicarbonate, 480 M MgCl 2 , 100 mM KCl, 500 M EDTA, 1 mM DTT, 20 mM HEPES, 1% BSA). The interaction was initiated by the incubation of GST-tagged PDZ domain with 4 l of total protein lysate of murine olfactory epithelium overnight at 4°C. As a background control, we performed the analysis three times using purified GST in place of the fusion proteins. After this incubation procedure, the protein-coupled beads were washed four times in 100 l of ice-cold wash buffer (50 mM ammonium bicarbonate, 480 M MgCl 2 , 100 mM KCl, 500 M EDTA, 1 mM DTT, 20 mM HEPES) followed by centrifugation for 3 min at 1,000 g at 4°C. The beads were resuspended in 100 of l 50 mM ammonium bicarbonate buffer (Fluka, Buchs, Switzerland) and 20% methanol (Avantor, Center Valley, PA). Thereafter, 2 g of sequencing-grade trypsin (Promega, Fitchburg, WI) was added for an overnight proteolytic digestion at 37°C. Subsequently, the extracted peptides were diluted in 0.1% trifluoroacetic acid and 2% acetonitrile and were loaded directly onto a fused silica capillary (inner diameter: 100 M, Polymicro, Lisle, IL) that was packed with a reversed phase (RP, C18 Luna, 3 m, 100 Å, Phenomenex, Torrance, CA) and strong cation exchanger (SCX, PolySULPHOETHYL A, 5 m, 200 Å, PolyLC, Columbia, SC) chromatography resin (10 -12 cm RP/4 cm SCX/3 cm RP). A multidimensional protein identification technology analysis was performed using a quaternary U-HPLC pump that was directly connected by a polyether ether ketone (PEEK) microcross to a Thermo LTQ Orbitrap or an LTQ XL ion trap mass spectrometer (Thermo Fisher Scientific, Waltham, MA). For optimal peptide separation, an effective flow rate of 200 (Ϯ50) nl/min, a spray voltage of 1.8 kV, and a transfer capillary temperature of 180°C were applied. For an eightstep multidimensional protein identification technology experiment, each step was represented by one instrument method file with identical settings for the HPLC gradient program. The mass spectrometer was set to detect a full MS spectrum between 400 -2,000 m/z for the precursor ion followed by MS/MS scans of the top 10 precursor ions from the preceding MS scan to detect fragment ions by the ion trap. MS/MS data interpretation was accomplished with Proteome Discoverer software version 1.2.0.208. The data were searched against the UniProt murine database (release 10.0, 20,752 entries) with semitryptic peptides, a mass accuracy of 10 ppm or 2.5 Da, a fragment ion tolerance of 1 Da, and the oxidation of methionine as a variable modification, allowing two missed cleavage sites. All of the accepted results had an XCorr of 1.8 for singly charged peptides, 2.5 for doubly charged peptides, and 3.5 for triply charged peptides (23,24). At least two different peptides must be identified per protein. The FDR value was greater than 2.7% based on a decoy database search. The analysis of the spectral count data was performed using the SAINTexpress software (25,26) and an AvgP cutoff Ͼ 0.66 was used to determine positive hits for further testing.
Co-Immunoprecipitations and Western Blotting-For co-immunoprecipitations, olfactory epithelium lysate from BL6 mice (ϾP20) was ground in RIPA buffer (150 nM NaCl, 50 nM Tris-HCl, 1% Nonidet, 0.5% deoxycholate (w/v), 0.1% SDS (w/v), and proteinase inhibitor) and processed using the Catch & Release system (Catch & Release v 2.0, Millipore, Darmstadt, Germany). After centrifugation at 1,000 g for 10 min, the lysate (4 mg of total protein) was incubated with primary antibodies (each 0.8 g) overnight at 4°C in a spin column. After several washing steps, the proteins were eluted using a denaturing elution buffer. Additionally, an input control and a negative control using a normal rabbit IgG-antibody were prepared in the same manner. The samples were loaded onto an SDS gel to separate them according to size and immunoblotted on a nitrocellulose membrane (Porablot NCL, Machery-Nagel, Dü ren, Germany). The interactions were detected through specific antibodies (1:50 -1:250) using the ECL Western blotting detection reagent (GE Healthcare, Little Chalfont, UK).
Phosphoinositide (4,5)-bisphosphate (PIP 2 )/Phosphatidylinositol Agarose Beads Interaction Assay-Phosphoinositide(4,5)-bisphosphate (PIP 2 ) and phosphatidylinositol-coupled agarose beads (Molecular Probes) were added to cell/tissue lysates in PIP 2 -binding buffer and rocked at 4°C for 3 h. The beads were washed four times in PIP 2 -binding buffer and then eluted using Laemmli buffer (1:20 DTT) and boiled for 5 min at 95°C. Then, 15 l of the supernatant was then loaded onto an SDS gel, along with a sample of the cell/tissue lysate as an input control and a negative control using uncoated agarose beads. Western blotting was performed as described, and anti-MUPP1 or anti-GFP antibodies were used to detect any interaction between MUPP1 and PIP 2 .
PIP 2 Strip TM Protein-Lipid Overlay Assay-PIP strips are nitrocellulose membranes containing 100 pmol samples of 15 different phospholipids and a blank control sample. Hana3A (27) cells were transfected with MUPP1-GFP and harvested 24 h after transfection. The cell pellet was homogenized in 800 l of RingerЈs solution containing protease inhibitors. Then, the lysate was diluted in 3% fatty acid-free BSA in TBS-T and added to the PIP strips for 2.5 h at room temperature. After several washing steps, the blot was analyzed using anti-GFP and anti-rabbit secondary antibody coupled to horseradish peroxidase.
Peptide Microarrays-CelluSpots TM Peptide Arrays (Intavis AG, Cologne, Germany) were blocked for 2 h at room temperature with 3% protein mixture (3K Eiweiss Shake, Layenberger, Rodenbach, Germany) in NaCl/Tris/Tween. The arrays were incubated with single PDZ domains (ϳ2 g of protein) of MUPP1 (1-13) that were fused to HA overnight at 4°C. Several washing steps were performed using NaCl/Tris/Tween. Then, the microarrays were incubated with anti-HA antibodies (1:250) for 4 h at room temperature and detected with HRP-coupled secondary antibodies using ECL plus a Western blotting detection reagent (GE Healthcare). Hits were considered positive when a signal was detectable in three out of three experiments.
GST Fusion Peptides-Single PDZ domains of MUPP1 (1-13) were produced in Escherichia coli BL21 DE3pLysS. The culture was induced using 0.5-1 mM isopropyl ␤-D-thiogalactoside (IPTG) and purified via glutathione-coupled Sepharose beads (Becton-Dickinson Biosciences, Franklin Lakes, NJ). Accurate protein sizes were verified using SDS gelelectrophoresis and Coomassie staining. The protein amount was measured using a Bradford protein assay (Coomassie Plus Protein Assay, Thermo Scientific).
Data Analysis-Microarray data were obtained using VilberLourmat Fusion Software. The spot intensities were analyzed with TIGR Spotfinder 3.2.1. The intensities (0 -255) were calculated as mean values for spots where a positive signal was observed three out of three times. The individual highest intensity of one array was set as 100%, while all of the other intensity values were normalized to this intensity. To depict the data, heat maps were designed using the R Program (3.13.0, www.r-project.org). 3-D scatter plots (mass, isoelectric point, and hydropathicity (Grand average of hydropathicity (GRAVY))) were obtained with ProtoParam (www.expasy.org). The plots of bound and unbound interaction partners were prepared with Sigma Plot 8.02. Further mass spectrometry analysis was performed with SAINTex-press (25, 26) (http://saint-apms.sourceforge.net/). PSPP (0.8.4) was used for binary logistic regression analysis (http://www.gnu.org/ software/pspp/).

Proteomic Analysis of Interactions between Mouse
Olfactory Receptor C Termini and MUPP1 PDZ Domains-To unravel the interaction characteristics of GPCRs and PDZ domains more precisely, we selected the protein family of murine ORs and the largest described PDZ protein to date, MUPP1, as model systems. We analyzed the binding of the 13 single PDZ domains of MUPP1 to OR C termini using a large-scale in vitro approach. We designed peptide microarrays consisting of peptides corresponding to the 20 C-terminal amino acids of 304 murine ORs. Based on a phylogenetic tree of all of the mouse ORs, we randomly chose target sequences from every mouse OR subfamily and immobilized the corresponding C-terminal peptides onto a cellulose layer. We also produced the individual PDZ domains 1-13 of MUPP1 in E. coli and purified them via a GST-tag using glutathione-coupled beads. For interaction studies, we incubated the microarrays with 2 g of purified HA-tagged PDZ domains; bound PDZ domains were detected by HA antibody incubation and chemiluminescence and purified GST was used as a negative control (Fig. 1A). Peptide interaction assays were performed three times for each PDZ domain, and the interactions were scored as positive when the signal was detected in all of the experiments. We analyzed the relative strength of the observed interactions between the 13 PDZ domains and their interacting receptors by normalizing the intensity values to the individual highest intensity of the spots on the microarrays and calculating the mean value of three replicates. The resulting intensity values are given in Fig. 1A as a color-coded heat map (for a detailed list of OR positions and intensity values, see Supplemental Table 1). With this approach, we could confirm known interactions for OR2AG1 (13) and mGluR7 (14,28) (14). Most of the ORs bound to more than one of the domains that were used for the interaction analysis (Fig. 1B). A total of 62% of the selected OR subsets showed interaction with different MUPP1 PDZ domains (Fig. 1C), while only 28% of the ORs bore classical C-terminal PDZ interaction motifs (15). Remarkably, 81% of the ORs that bound to PDZ domains did not contain classical C-terminal PDZ interaction motifs (Fig. 1D). In addition, we quantified the relative binding intensities of these PDZ domains to mutant OR C termini peptides of 15 amino acids in length (Fig. 1E), showing that X-W-X-W mutations abolished PDZ-GPCR binding, while X-A-X-A mutations could even enhance the relative interaction strength.  Table 4. (B) Bar charts comparing the C termini, while the third intracellular loop (IC3) seemed to be a promising candidate for PDZ interaction (Fig. 1E). Together, the 13 MUPP1 domains interacted with 189 of the 304 OR C-terminal peptides that were immobilized on the microarrays. Many ORs were attached to more than one specific domain, resulting in a total of 1,229 interactions.
Physicochemical Properties of MUPP1Јs PDZ Domains-We further analyzed the physicochemical properties of each PDZ domain of MUPP1 compared with the PDZ domains of other scaffolding proteins (Supplemental Table 7) and calculated 3-D scatter plots displaying the hydropathicity (GRAVY), pI, and mass of the PDZ domains ( Fig. 2A). The PDZ domains of MUPP1 are very similar regarding their size and their hydropathicity (-0.153 to 0.327) and cluster in relatively narrow ranges compared with PDZ domains in general (-1.499 to 0.651 for our selection).
In addition, we analyzed the physicochemical properties of GPCR interactors and noninteractors and designed 3-D scatter plots representing the results. We did not observe any prominent differences between interacting and noninteracting OR C termini (Fig. 2B) concerning their chemical properties; both had relatively high pI values and medium hydropathicity. Box whisker plots show a nearly equal distribution of all of the properties between interacting and noninteracting peptides ( Fig 2C). Furthermore, a binary logistic regression analysis showed only a 68.09% correct prediction based on the three factors (Fig 2D). A general binding affinity regarding these chemical properties could therefore not be determined for GPCR-PDZ interactions.
Phosphorylation of PDZ Binding Partners Alters PDZ Interactions-To account for the significance of phosphorylated PDZ binding partners regarding the modulation of interactions, we designed peptide microarrays of 20 amino acids in length with unmodified and phosphorylated C-terminal and IC3 (intracellular loop 3) peptides of deorphanized ORs and the olfactory-relevant Cl Ϫ channel family TMEM16 (29 -32) and performed interaction assays with single 13 PDZ domains of MUPP1 as described above (Fig. 3A). Remarkably, the phosphorylation of both the C-terminal region as well as IC3 seems to strongly influence PDZ binding. All of the observed interactions (32 in total) between the 13 PDZ domains of MUPP1 and the tested IC3 regions of different ORs and ion channels were completely abolished after phosphorylation (Fig. 3A). The interactions of the tested PDZ domains with the C-terminal regions of diverse ORs were altered to both higher (8) and lower (2) relative binding affinities after phosphorylation. However, specifying a general tendency, the number of interactions was reduced to nearly half (71/36) after phosphorylation (Fig. 3B).

Interaction Analysis of the PDZ Domains of MUPP1 and GPCR-Associated Signaling
Components-In addition to the interactions of PDZ domains with GPCRs, we were also interested in binding events to GPCR-related and other signaling proteins. Therefore, we probed a signaling peptide-containing microarray with purified 13 PDZ domains of MUPP1 and detected putative interactions as described above (Fig.  4A). The C-terminal peptide sequences of these microarrays were of 20 amino acids in length. We detected interactions for MUPP1 PDZ domains with phospholipase C ␥1 (PLC-␥1), phosphoinositide-3-kinase, TRPC6, Whirlin, SANS, PACS-1, GluR2, and G olf .
Because phosphoinositide signaling is often cross-linked to GPCR-mediated signal transduction cascades, we were interested in the importance of the PDZ domain binding for these interactions. We found that recombinantly expressed MUPP1 as well as MUPP1 isolated from murine olfactory epithelium can be precipitated with PIP 2 -coated beads, while a pull down with unphosphorylated phosphatidylinositol beads did not show any interaction (Fig. 4B). Additionally, we detected interactions of recombinantly expressed MUPP1 with various phosphoinositides using a protein-lipid overlay assay (Fig. 4C).
Pull-Down-Based Mass Spectrometry of PDZ (MUPP1) Interaction Partners-In addition to the characterization of PDZ-GPCR interactions in the heterologous system, we wanted to examine the detailed binding events of our PDZ domains and proteins. Therefore, we performed pull-down experiments with nine single PDZ domains of MUPP1 and proteins from murine OE lysate and analyzed PDZ interaction partners by mass spectrometry. We used the PDZ domains that yielded the highest and purest peptide amounts for further examination. Analyses were performed three times using individual OE and PDZ samples. The background estimation was performed three times by pulling down OE proteins with purified GST protein and acted as a control for the mass spectrometry data. For the probabilistical analysis of these protein-protein interaction datasets, we used the SAINTexpress Software (25,26) to evaluate the quality of the identified interactions. The raw data from the mass spectrometry analysis are listed in Supplemental Table 5, while the SAINTexpress evaluation of  Table 6. A list of the most interesting interactions regarding the olfactory system with an AvgP cutoff at 0.66 is shown in Fig. 5. The tested PDZ domains of MUPP1 interacted with various olfactory signaling proteins, such as G olf , G s ␤-subunit, calmodulin (CaM), olfactory marker protein (OMP), and clathrin as well as with other promising signaling proteins and ion channels, e.g. KCNK3 or VLGR1.
To validate the mass spectrometry data, we performed co-immunoprecipitations using OE lysate and antibodies against MUPP1 and other signaling proteins to pull out protein complexes and detect specific interaction partners (Fig. 6A). We used a normal rabbit IgG-antibody for control pull downs and OE lysate as an input control. Thereby, we could validate the interactions for MUPP1 with CaM, OMP, VDAC1, Golf, and VLGR1, indicating the accuracy of the mass-spectrometry-based approach. To observe the detailed localization of MUPP1 interaction partners, we performed immunohistological stainings using a transgenic heterozygous OMP-GFP mouse strain (22) and specific antibodies (Fig. 6B). Thereby, we could detect the colocalization of MUPP1 with G olf , CaM, AC3, OMP, and VDAC1 within the apical region of the murine OE, indicating a close proximity of these proteins within this tissue. DISCUSSION PDZ domains are one of the most frequently encountered domains that are specifically responsible for protein-protein  Table 5. interactions. These domains consist of ϳ90 amino acids and can occur in multiple copies inside a protein (33,34). PDZ proteins can fulfill various functions, build large protein networks, and are involved in the modulation of many different cellular processes. Especially in neuronal signal transduction processes, the spatial and temporal organization of its components is of utter importance. Therefore, PDZ domains are key players in these systems, regulating the interactions of all of the involved signaling proteins and enhancing the efficiency and speed of signal processing (35,36). Among the multitude of PDZ interactions partners, one important group is the diverse family of GPCRs (37). These membrane proteins are essential for various cellular processes and can be involved in different diseases, such as cancer, obesity, diabetes, or central nervous system disorders (38 -40), and therefore make excellent therapeutic targets. Thus far, the detailed binding characteristics and regulatory mechanisms of PDZ-GPCR interactions are still under debate and must be further investigated. Here, we carried out a comprehensive interaction study of the 13 single PDZ domains of MUPP1 and 304 murine olfactory receptors, representing the protein family of GPCRs, using peptide microarrays. The analyzed PDZ domains of MUPP1 show very conserved chemical characteristics compared with the PDZ domains of other proteins, sug-FIG. 6. MUPP1 interacts and colocalizes with different signaling proteins. (A) Co-immunoprecipitations of MUPP1 with different potential binding partners as identified through mass spectrometric analysis. Rows 1-3: Co-immunoprecipitations of endogenous MUPP1 with calmodulin, olfactory marker protein (OMP) and voltage-dependent anion channel 1 (VDAC1), including normal rabbit IgG-antibody controls and input controls. Row 4: Co-Immunoprecipitations of endogenous G olf and VLGR1 with MUPP1, including normal rabbit IgG-antibody control and input control. Murine OE lysate was immunoprecipitated with the protein mentioned above and detected with antibodies against the interacting protein. (B) Confocal images of horizontal cryosections using heterozygous OMP-GFP mice stained with specific antibodies against MUPP1 (red) and G olf , calmodulin, AC3, OMP, VDAC1 (blue). Scale bar ϭ 10 m. Arrows in overlay images indicate colocalization of MUPP1 and G olf , calmodulin, AC3, OMP, VDAC1 (pink/purple). gesting a specialized scaffolding complex. Interestingly, the chemical properties of PDZ-interacting and noninteracting peptides seem to be relatively similar to each other, demonstrating the complexity of PDZ interactions. In contrast to studies in which classical hieratic binding motifs were presented (15), we demonstrated a disagreement between those motifs and our practical findings. The proportion of interacting and noninteracting peptides to classical binding motifs was nearly vice versa. However, PDZ domains do not fall in discrete classes but are distributed throughout selectivity space, thereby avoiding cross reactivity (17,41). Furthermore, extensive interaction studies in Caenorhabditis elegans revealed frequent noncanonical interactions between PDZ domains and their binding partners, indicating the extreme fine-tuning of PDZ-ligand interactions (18). The individual domains bind with different affinities to various immobilized ORs, suggesting the presence of a highly dynamic and flexible interactome (42)(43)(44).
Additionally, we analyzed the interactions of PDZ domains and various signaling peptides. Here, we could validate a published MUPP1 interactor (G olf ) (14) but also identify new PDZ binding partners, such as PLC-␥1, PIP 2 , TRPC6, SANS, or PACS-1. This identification indicates that the putative interactions of PDZ proteins are not exclusively with peptides but also with phospholipids. Different studies have demonstrated the relevance of PDZ-phosphoinositide interactions and have claimed the necessity to further characterize such binding events (45). Interestingly, only phosphorylated versions of phosphatidylinositol were able to interact with MUPP1 in either of our approaches, further indicating the importance of such modifications for PDZ binding. PDZ-lipid interactions are important for cell signaling, receptor trafficking, and receptor cargo recycling (46). Furthermore, an interaction of MUPP1 and tandem-PH-domain-containing protein-1 (TAPP1) is important for PIP 2 -mediated signaling (47), which further strengthens the idea of an interactome consisting of PDZ proteins, GPCRs, signaling peptides, and lipids. Interestingly, associated enzymes, such as PLC-␥1, were also able to interact with the PDZ domains of MUPP1. In this context, different PLC subtypes bear classical PDZ binding motifs, but detailed binding modalities must be further determined (48).
We were interested in the mechanisms regulating and influencing interactions of PDZ proteins with GPCRs. Because posttranslational modifications, such as phosphorylation, are important regulatory instruments, we focused on PDZ binding to GPCR peptides with phosphorylated C termini and IC3. We found that PDZ binding promiscuity was diminished for most interactions after introducing a phospho-group to both IC3 and C termini. This finding agrees with other studies demonstrating that the phosphorylation of PDZ ligands modifies binding properties (49). The phosphorylation of AMPA receptor subunit GluR2 decreases binding to the PDZ protein GRIP1 but not to PICK1, indicating the differential dynamics of such a regulatory mechanism (50). Similar to our results, in some cases, the binding affinities of PDZ-ligand interactions are even enhanced after phosphorylation (51). Additionally, distinct C-terminal mutations (X-A-X-A and X-W-X-W) of the PDZ ligands led to altered relative binding affinities. Alanine substitutions caused even higher relative intensity values, whereas sterically complex tryptophan mutations abolished PDZ-ligand binding in all of the tested scenarios. Tryptophans at C-terminal positions 0 and -2 can rarely be found within PDZ interactors, and these mutations are capable of inhibiting interactions (13,14,17). These observations indicate the importance of C-terminal positions 0 and -2 and suggest that underlying amino acids can be variable within a certain range.
In addition to in vitro interaction analyses, we focused on PDZ-ligand binding using murine olfactory epithelium and pull-down-associated mass spectrometry. We could thereby validate the published interactions of MUPP1Јs PDZ domains with G olf and CaM (14) and identify unknown binding partners. Here, OMP was able to bind to several PDZ domains and has been reported to play important roles in olfactory signaling and development of the olfactory system (52)(53)(54). VLGR1, the largest known GPCR, which is part of a PDZ-organized protein network that is associated with a type of sensineuronal hearing loss called Usher syndrome (55,56) is another new, interesting interaction partner of MUPP1. To validate the accuracy of our mass-spectrometry-based interaction analyses, we performed additional co-immunoprecipitations and colocalization studies that confirmed our findings.
In conclusion, we provide here a comprehensive overview of the binding characteristics of diverse PDZ-GPCR interactions based on the murine olfactory system, indicating the complexity und variability of these binding events, and found that phosphorylation of both the receptorЈs C terminus and IC3 alters PDZ binding promiscuity. Thus, we anticipate that our study will greatly enhance the understanding of these important protein-protein and protein-lipid interactions and might identify novel clinical targets for GPCR-associated diseases.