A specificity map for the PDZ domain family

PDZ domains are protein-protein interaction modules that recognize specific C-terminal sequences to assemble protein complexes in multicellular organisms. By scanning billions of random peptides, we accurately map binding specificity for approximately half of the over 330 PDZ domains in the human and Caenorhabditis elegans proteomes. The domains recognize features of the last seven ligand positions, and we find 16 distinct specificity classes conserved from worm to human, significantly extending the canonical two-class system based on position -2. Thus, most PDZ domains are not promiscuous, but rather are fine-tuned for specific interactions. Specificity profiling of 91 point mutants of a model PDZ domain reveals that the binding site is highly robust, as all mutants were able to recognize C-terminal peptides. However, many mutations altered specificity for ligand positions both close and far from the mutated position, suggesting that binding specificity can evolve rapidly under mutational pressure. Our specificity map enables the prediction and prioritization of natural protein interactions, which can be used to guide PDZ domain cell biology experiments. Using this approach, we predicted and validated several viral ligands for the PDZ domains of the SCRIB polarity protein. These findings indicate that many viruses produce PDZ ligands that disrupt host protein complexes for their own benefit, and that highly pathogenic strains target PDZ domains involved in cell polarity and growth.     

Classification of PDZ domains.

A diverse family of PDZ domains has been identified, but the rules that govern their ligand specificity are not clear. Here we propose a novel classification of PDZ domains based on the nature of amino acids in the two critical positions in the PDZ domain fold. Using these principles, we classified PDZ domains present in the SMART database. Using yeast two-hybrid, in vitro pull-down and plasmon surface resonance assays, we demonstrated that in agreement with their position in the proposed classification the Mint1-1, hINADL-5, and PAR6 PDZ domains display similar dual ligand specificity. The proposed classification helps to organize PDZ domain containing proteins.     

Distinct binding specificity of the multiple PDZ domains of INADL, a human protein with homology to INAD from Drosophila melanogaster

PDZ domains are protein-protein interaction modules that typically bind to short peptide sequences at the carboxyl terminus of target proteins. Proteins containing multiple PDZ domains often bind to different trans-membrane and intracellular proteins, playing a central role as organizers of multimeric complexes. To characterize the rules underlying the binding specificity of different PDZ domains, we have assembled a novel repertoire of random peptides that are displayed at high density at the carboxyl terminus of the capsid D protein of bacteriophage lambda. We have exploited this combinatorial library to determine the peptide binding preference of the seven PDZ domains of human INADL, a multi-PDZ protein that is homologous to the INAD protein of Drosophila melanogaster. This approach has permitted the determination of the consensus ligand for each PDZ domain and the assignment to class I, class II, and to a new specificity class, class IV, characterized by the presence of an acidic residue at the carboxyl-terminal position. Homology modeling and site-directed mutagenesis experiments confirmed the involvement of specific residues at contact positions in determining the domain binding preference. However, these experiments failed to reveal simple rules that would permit the association of the chemical characteristics of any given residue in the peptide binding pocket to the preference for specific amino acid sequences in the ligand peptide. Rather, they suggested that to infer the binding preference of any PDZ domain, it is necessary to simultaneously take into account all contact positions by using computational procedures. For this purpose we extended the SPOT algorithm, originally developed for SH3 domains, to evaluate the probability that any peptide would bind to any given PDZ domain.     

The multiple‐specificity landscape of modular peptide recognition domains

Modular protein interaction domains form the building blocks of eukaryotic signaling pathways. Many of them, known as peptide recognition domains, mediate protein interactions by recognizing short, linear amino acid stretches on the surface of their cognate partners with high specificity. Residues in these stretches are usually assumed to contribute independently to binding, which has led to a simplified understanding of protein interactions. Conversely, we observe in large binding peptide data sets that different residue positions display highly significant correlations for many domains in three distinct families (PDZ, SH3 and WW). These correlation patterns reveal a widespread occurrence of multiple binding specificities and give novel structural insights into protein interactions. For example, we predict a new binding mode of PDZ domains and structurally rationalize it for DLG1 PDZ1. We show that multiple specificity more accurately predicts protein interactions and experimentally validate some of the predictions for the human proteins DLG1 and SCRIB. Overall, our results reveal a rich specificity landscape in peptide recognition domains, suggesting new ways of encoding specificity in protein interaction networks.     

Genome-wide analysis of PDZ domain binding reveals inherent functional overlap within the PDZ interaction network

Binding selectivity and cross-reactivity within one of the largest and most abundant interaction domain families, the PDZ family, has long been enigmatic. The complete human PDZ domain complement (the PDZome) consists of 267 domains and we applied here a Bayesian selectivity model to predict hundreds of human PDZ domain interactions, using target sequences of 22,997 non-redundant proteins. Subsequent analysis of these binding scores shows that PDZs can be divided into two genome-wide clusters that coincide well with the division between canonical class 1 and 2 PDZs. Within the class 1 PDZs we observed binding overlap at unprecedented levels, mediated by two residues at positions 1 and 5 of the second α-helix of the binding pocket. Eight PDZ domains were subsequently selected for experimental binding studies and to verify the basics of our predictions. Overall, the PDZ domain class 1 cross-reactivity identified here implies that auxiliary mechanisms must be in place to overcome this inherent functional overlap and to minimize cross-selectivity within the living cell. Indeed, when we superimpose PDZ domain binding affinities with gene ontologies, network topology data and the domain position within a PDZ superfamily protein, functional overlap is minimized and PDZ domains position optimally in the binding space. We therefore propose that PDZ domain selectivity is achieved through cellular context rather than inherent binding specificity.     

Alteration of the C-Terminal Ligand Specificity of the Erbin PDZ Domain by Allosteric Mutational Effects

Modulation of protein binding specificity is important for basic biology and for applied science. Here we explore how binding specificity is conveyed in PDZ (postsynaptic density protein-95/discs large/zonula occludens-1) domains, small interaction modules that recognize various proteins by binding to an extended C terminus. Our goal was to engineer variants of the Erbin PDZ domain with altered specificity for the most C-terminal position (position 0) where a Val is strongly preferred by the wild-type domain. We constructed a library of PDZ domains by randomizing residues in direct contact with position 0 and in a loop that is close to but does not contact position 0. We used phage display to select for PDZ variants that bind to 19 peptide ligands differing only at position 0. To verify that each obtained PDZ domain exhibited the correct binding specificity, we selected peptide ligands for each domain. Despite intensive efforts, we were only able to evolve Erbin PDZ domain variants with selectivity for the aliphatic C-terminal side chains Val, Ile and Leu. Interestingly, many PDZ domains with these three distinct specificities contained identical amino acids at positions that directly contact position 0 but differed in the loop that does not contact position 0. Computational modeling of the selected PDZ domains shows how slight conformational changes in the loop region propagate to the binding site and result in different binding specificities. Our results demonstrate that second-sphere residues could be crucial in determining protein binding specificity.     

Structure-Based Prediction of the Peptide Sequence Space Recognized by Natural and Synthetic PDZ Domains

Protein-protein recognition, frequently mediated by members of large families of interaction domains, is one of the cornerstones of biological function. Here, we present a computational, structure-based method to predict the sequence space of peptides recognized by PDZ domains, one of the largest families of recognition proteins. As a test set, we use a considerable amount of recent phage display data that describe the peptide recognition preferences for 169 naturally occurring and engineered PDZ domains. For both wild-type PDZ domains and single point mutants, we find that 70-80% of the most frequently observed amino acids by phage display are predicted within the top five ranked amino acids. Phage display frequently identified recognition preferences for amino acids different from those present in the original crystal structure. Notably, in about half of these cases, our algorithm correctly captures these preferences, indicating that it can predict mutations that increase binding affinity relative to the starting structure. We also find that we can computationally recapitulate specificity changes upon mutation, a key test for successful forward design of protein-protein interface specificity. Across all evaluated data sets, we find that incorporation backbone sampling improves accuracy substantially, irrespective of using a crystal or NMR structure as the starting conformation. Finally, we report successful prediction of several amino acid specificity changes from blind tests in the DREAM4 peptide recognition domain specificity prediction challenge. Because the foundational methods developed here are structure based, these results suggest that the approach can be more generally applied to specificity prediction and redesign of other protein-protein interfaces that have structural information but lack phage display data.     

The Tiam1 PDZ Domain Couples to Syndecan1 and Promotes Cell-Matrix Adhesion

The T-cell lymphoma invasion and metastasis gene 1 (Tiam1) is a guanine exchange factor (GEF) for the Rho-family GTPase Rac1 that is crucial for the integrity of adherens junctions, tight junctions, and cell-matrix interactions. This GEF contains several protein-protein interaction domains, including a PDZ domain. Earlier studies identified a consensus PDZ-binding motif and a synthetic peptide capable of binding to the Tiam1 PDZ domain, but little is known about its ligand specificity and physiological role in cells. Here, we investigated the structure, specificity, and function of the Tiam1 PDZ domain. We determined the crystal structures of the Tiam1 PDZ domain free and in complex with a “model” peptide, which revealed the structural basis for ligand specificity. Protein database searches using the consensus PDZ-binding motif identified two eukaryotic cell adhesion proteins, Syndecan1 and Caspr4, as potential Tiam1 PDZ domain binding proteins. Equilibrium binding experiments confirmed that C-terminal peptides derived from Syndecan1 and Caspr4 bound the Tiam1 PDZ domain. NMR chemical shift perturbation experiments indicated that the Tiam1 PDZ/Syndecan1 and PDZ/Caspr4 complexes were structurally distinct and identified key residues likely to be responsible for ligand selectivity. Moreover, cell biological analysis established that Syndecan1 is a physiological binding partner of Tiam1 and that the PDZ domain has a function in cell-matrix adhesion and cell migration. Collectively, our data provide insight into the structure, specificity, and function of the Tiam1 PDZ domain. Importantly, our data report on a physiological role for the Tiam1 PDZ domain and establish a novel link between two previously unrelated signal transduction pathways, both of which are implicated in cancer.     

Comprehensive analysis of all evolutionary paths between two divergent PDZ domain specificities

To understand the molecular evolution of functional diversity in protein families, we comprehensively investigated the consequences of all possible mutation combinations separating two peptide‐binding domains with highly divergent specificities. We analyzed the Erbin PDZ domain (Erbin‐PDZ), which exhibits canonical type I specificity, and a synthetic Erbin‐PDZ variant (E‐14) that differs at six positions and exhibits an atypical specificity that closely resembles that of the natural Pdlim4 PDZ domain (Pdlim4‐PDZ). We constructed a panel of 64 PDZ domains covering all possible transitions between Erbin‐PDZ and E‐14 (i.e., the panel contained variants with all possible combinations of either the Erbin‐PDZ or E‐14 sequence at the six differing positions). We assessed the specificity profiles of the 64 PDZ domains using a C‐terminal phage‐displayed peptide library containing all possible genetically encoded heptapeptides. The specificity profiles clustered into six distinct groups, showing that intermediate domains can be nodes for the evolution of divergent functions. Remarkably, three substitutions were sufficient to convert the specificity of Erbin‐PDZ to that of Pdlim4‐PDZ, whereas Pdlim4‐PDZ contains 71 differences relative to Erbin‐PDZ. X‐ray crystallography revealed the structural basis for specificity transition: a single substitution in the center of the binding site, supported by contributions from auxiliary substitutions, altered the main chain conformation of the peptide ligand to resemble that of ligands bound to Pdlim4‐PDZ. Our results show that a very small set of mutations can dramatically alter protein specificity, and these findings support the hypothesis whereby complex protein functions evolve by gene duplication followed by cumulative mutations.     

Distinct ligand specificity of the Tiam1 and Tiam2 PDZ domains

Guanine nucleotide exchange factor proteins of the Tiam family are activators of the Rho GTPase Rac1 and critical for cell morphology, adhesion, migration, and polarity. These proteins are modular and contain a variety of interaction domains, including a single post-synaptic density-95/discs large/zonula occludens-1 (PDZ) domain. Previous studies suggest that the specificities of the Tiam1 and Tiam2 PDZ domains are distinct. Here, we sought to conclusively define these specificities and determine their molecular origin. Using a combinatorial peptide library, we identified a consensus binding sequence for each PDZ domain. Analysis of these consensus sequences and binding assays with peptides derived from native proteins indicated that these two PDZ domains have overlapping but distinct specificities. We also identified residues in two regions (S(0) and S(-2) pockets) of the Tiam1 PDZ domain that are important determinants of ligand specificity. Site-directed mutagenesis of four nonconserved residues in these two regions along with peptide binding analyses confirmed that these residues are crucial for ligand affinity and specificity. Furthermore, double mutant cycle analysis of each region revealed energetic couplings that were dependent on the ligand being investigated. Remarkably, a Tiam1 PDZ domain quadruple mutant had the same specificity as the Tiam2 PDZ domain. Finally, analysis of Tiam family PDZ domain sequences indicated that the PDZ domains segregate into four distinct families based on the residues studied here. Collectively, our data suggest that Tiam family proteins have highly evolved PDZ domain-ligand interfaces with distinct specificities and that they have disparate PDZ domain-dependent biological functions.     

Structural and functional analysis of the PDZ domains of human HtrA1 and HtrA3

High-temperature requirement A (HtrA) and its homologs contain a serine protease domain followed by one or two PDZ domains. Bacterial HtrA proteins and the mitochondrial protein HtrA2/Omi maintain cell function by acting as both molecular chaperones and proteases to manage misfolded proteins. The biological roles of the mammalian family members HtrA1 and HtrA3 are less clear. We report a detailed structural and functional analysis of the PDZ domains of human HtrA1 and HtrA3 using peptide libraries and affinity assays to define specificity, structural studies to view the molecular details of ligand recognition, and alanine scanning mutagenesis to investigate the energetic contributions of individual residues to ligand binding. In common with HtrA2/Omi, we show that the PDZ domains of HtrA1 and HtrA3 recognize hydrophobic polypeptides, and while C-terminal sequences are preferred, internal sequences are also recognized. However, the details of the interactions differ, as different domains rely on interactions with different residues within the ligand to achieve high affinity binding. The results suggest that mammalian HtrA PDZ domains interact with a broad range of hydrophobic binding partners. This promiscuous specificity resembles that of bacterial HtrA family members and suggests a similar function for recognizing misfolded polypeptides with exposed hydrophobic sequences. Our results support a common activation mechanism for the HtrA family, whereby hydrophobic peptides bind to the PDZ domain and induce conformational changes that activate the protease. Such a mechanism is well suited to proteases evolved for the recognition and degradation of misfolded proteins.     

An allosteric intramolecular PDZ-PDZ interaction modulates PTP-BL PDZ2 binding specificity

PDZ (acronym of the synapse-associated protein PSD-95/SAP90, the septate junction protein Discs-large, and the tight junction protein ZO-1) domains are abundant small globular protein interaction domains that mainly recognize the carboxyl termini of their target proteins. Detailed knowledge on PDZ domain binding specificity is a prerequisite for understanding the interaction networks they establish. We determined the binding preference of the five PDZ domains in the protein tyrosine phosphatase PTP-BL by screening a random C-terminal peptide lambda phage display library. Interestingly, the potential of PDZ2 to interact with class III-type ligands was found to be modulated by the presence of PDZ1. Structural studies revealed a direct and specific interaction of PDZ1 with a surface on PDZ2 that is opposite the peptide binding groove. Long-range allosteric effects that cause structural changes in the PDZ2 peptide binding groove thus explain the altered PDZ2 binding preference. Our results experimentally corroborate that the molecular embedding of PDZ domains is an important determinant of their ligand binding specificity.     

Membrane localization is critical for activation of the PICK1 BAR domain

The PSD-95/Discs-large/ZO-1 homology (PDZ) domain protein, protein interacting with C kinase 1 (PICK1) contains a C-terminal Bin/amphiphysin/Rvs (BAR) domain mediating recognition of curved membranes; however, the molecular mechanisms controlling the activity of this domain are poorly understood. In agreement with negative regulation of the BAR domain by the N-terminal PDZ domain, PICK1 distributed evenly in the cytoplasm, whereas truncation of the PDZ domain caused BAR domain-dependent redistribution to clusters colocalizing with markers of recycling endosomal compartments. A similar clustering was observed both upon truncation of a short putative α-helical segment in the linker between the PDZ and the BAR domains and upon coexpression of PICK1 with a transmembrane PDZ ligand, including the alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor GluR2 subunit, the GluR2 C-terminus transferred to the single transmembrane protein Tac or the dopamine transporter C-terminus transferred to Tac. In contrast, transfer of the GluR2 C-terminus to cyan fluorescent protein, a cytosolic protein, did not elicit BAR domain-dependent clustering. Instead, localizing PICK1 to the membrane by introducing an N-terminal myristoylation site produced BAR domain-dependent, but ligand-independent, PICK1 clustering. The data support that in the absence of PDZ ligand, the PICK1 BAR domain is inhibited through a PDZ domain-dependent and linker-dependent mechanism. Moreover, they suggest that unmasking of the BAR domain's membrane-binding capacity is not a consequence of ligand binding to the PDZ domain per se but results from, and coincides with, recruitment of PICK1 to a membrane compartment.     

Energetic determinants of internal motif recognition by PDZ domains

PDZ domains are protein-protein interaction modules that organize intracellular signaling complexes. Most PDZ domains recognize specific peptide motifs followed by a required COOH-terminus. However, several PDZ domains have been found which recognize specific internal peptide motifs. The best characterized example is the syntrophin PDZ domain which, in addition to binding peptide ligands with the consensus sequence -E-S/T-X-V-COOH, also binds the neuronal nitric oxide synthase (nNOS) PDZ domain in a manner that does not depend on its precise COOH-terminal sequence. In the structure of the syntrophin-nNOS PDZ heterodimer complex, the two PDZ domains interact in a head-to-tail fashion, with an internal sequence from the nNOS PDZ domain binding precisely at the peptide binding groove of the syntrophin PDZ domain. To understand the energetic basis of this alternative mode of PDZ recognition, we have undertaken an extensive mutagenic and biophysical analysis of the nNOS PDZ domain and its interaction with the syntrophin PDZ domain. Our data indicate that the presentation of the nNOS internal motif within the context of a rigid beta-hairpin conformation is absolutely essential to binding; amino acids crucial to the structural integrity of the hairpin are as important or more important than residues that make direct contacts. The results reveal the general rules of PDZ recognition of diverse ligand types.     

Sequence-specific long range networks in PSD-95/discs large/ZO-1 (PDZ) domains tune their binding selectivity

Protein-protein interactions mediated by modular protein domains are critical for cell scaffolding, differentiation, signaling, and ultimately, evolution. Given the vast number of ligands competing for binding to a limited number of domain families, it is often puzzling how specificity can be achieved. Selectivity may be modulated by intradomain allostery, whereby a remote residue is energetically connected to the functional binding site via side chain or backbone interactions. Whereas several energetic pathways, which could mediate intradomain allostery, have been predicted in modular protein domains, there is a paucity of experimental data to validate their existence and roles. Here, we have identified such functional energetic networks in one of the most common protein-protein interaction modules, the PDZ domain. We used double mutant cycles involving site-directed mutagenesis of both the PDZ domain and the peptide ligand, in conjunction with kinetics to capture the fine energetic details of the networks involved in peptide recognition. We performed the analysis on two homologous PDZ-ligand complexes and found that the energetically coupled residues differ for these two complexes. This result demonstrates that amino acid sequence rather than topology dictates the allosteric pathways. Furthermore, our data support a mechanism whereby the whole domain and not only the binding pocket is optimized for a specific ligand. Such cross-talk between binding sites and remote residues may be used to fine tune target selectivity.     

Molecular Basis of Class III Ligand Recognition by PDZ3 in Murine Protein Tyrosine Phosphatase PTPN13

Protein tyrosine phosphatase PTPN13, also known as PTP-BL in mice, represents a large multi-domain non-transmembrane scaffolding protein that contains five consecutive PDZ domains. Here, we report the solution structures of the extended murine PTPN13 PDZ3 domain in its apo form and in complex with its physiological ligand, the carboxy-terminus of protein kinase C-related kinase-2 (PRK2), determined by multidimensional NMR spectroscopy. Both in its ligand-free state and when complexed to PRK2, PDZ3 of PTPN13 adopts the classical compact, globular D/E fold. PDZ3 of PTPN13 binds five carboxy-terminal amino acids of PRK2 via a groove located between the EB-strand and the DB-helix. The PRK2 peptide resides in the canonical PDZ3 binding cleft in an elongated manner and the amino acid side chains in position P0 and P-2, cysteine and aspartate, of the ligand face the groove between EB-strand and DB-helix, whereas the PRK2 side chains of tryptophan and alanine located in position P-1 and P-3 point away from the binding cleft. These structures are rare examples of selective class III ligand recognition by a PDZ domain and now provide a basis for the detailed structural investigation of the promiscuous interaction between the PDZ domains of PTPN13 and their ligands. They will also lead to a better understanding of the proposed scaffolding function of these domains in multi-protein complexes assembled by PTPN13 and could ultimately contribute to low molecular weight antagonists that might even act on the PRK2 signaling pathway to modulate rearrangements of the actin cytoskeleton.     

Structural Basis for Exquisite Specificity of Affinity Clamps, Synthetic Binding Proteins Generated through Directed Domain-interface Evolution

We have established a new protein-engineering strategy termed “directed domain-interface evolution” that generates a binding site by linking two protein domains and then optimizing the interface between them. Using this strategy, we have generated synthetic two-domain “affinity clamps” using PDZ and fibronectin type III (FN3) domains as the building blocks. While these affinity clamps all had significantly higher affinity toward a target peptide than the underlying PDZ domain, two distinct types of affinity clamps were found in terms of target specificity. One type conserved the specificity of the parent PDZ domain, and the other increased the specificity dramatically. Here, we characterized their specificity profiles using peptide phage-display libraries and scanning mutagenesis, which suggested a significantly enlarged recognition site of the high-specificity affinity clamps. The crystal structure of a high-specificity affinity clamp showed extensive contacts with a portion of the peptide ligand that is not recognized by the parent PDZ domain, thus rationalizing the improvement of the specificity of the affinity clamp. A comparison with another affinity clamp structure showed that, although both had extensive contacts between PDZ and FN3 domains, they exhibited a large offset in the relative position of the two domains. Our results indicate that linked domains could rapidly fuse and evolve as a single functional module, and that the inherent plasticity of domain interfaces allows for the generation of diverse active-site topography. These attributes of directed domain-interface evolution provide facile means to generate synthetic proteins with a broad range of functions.     

Quantification of PDZ domain specificity, prediction of ligand affinity and rational design of super-binding peptides

Transient macromolecular complexes are often formed by protein-protein interaction domains (e.g. PDZ, SH2, SH3, WW) which recognize linear sequence motifs with in vitro affinities typically in the micromolar range. The analysis of the resulting interaction networks requires a quantification of domain specificity and selectivity towards all possible ligands with physiologically relevant affinity. As representative examples, we determined specificity as a function of ligand sequence-dependent affinity contributions by statistical analysis of peptide library screens for the AF6, ERBIN and SNA1 (alpha-1-syntrophin) PDZ domains. For this purpose, the three PDZ domains were first screened for binding with a peptide library comprising 6223 human C termini created by SPOT synthesis. Based on the detected ligand preferences, we designed focused peptide libraries (profile libraries). These libraries were used to quantify the affinity contributions of the four C-terminal ligand residues by means of ANOVA models (analysis of variance) relating the C-terminal ligand sequences to the corresponding dissociation constants. Our models agreed well with experimentally determined dissociation constants and allowed us to design super binding peptides. The latter were shown experimentally to bind to their cognate PDZ domains with the highest affinity. In addition, we determined structure-activity relationships and thereby rationalized the position-specific affinity contributions. Furthermore, we used the statistical models to predict the dissociation constants for the complete ligand sequence space and thus determined the specificity overlap for the three investigated PDZ domains (). Altogether, we present an efficient method for profiling protein-protein interaction domains that provides a biophysical picture of specificity and selectivity. This approach allows the rational design of functional experiments and provides a basis for simulating interaction networks in the field of systems biology.     

Computer-aided design of a PDZ domain to recognize new target sequences

PDZ domains are small globular domains that recognize the last 4-7 amino acids at the C-terminus of target proteins. The specificity of the PDZ-ligand recognition is due to side chain-side chain interactions, as well as the positioning of an alpha-helix involved in ligand binding. We have used computer-aided protein design to produce mutant versions of a Class I PDZ domain that bind to novel Class I and Class II target sequences both in vitro and in vivo, thus providing an alternative to primary antibodies in western blotting, affinity chromatography and pull-down experiments. Our results suggest that by combining different backbone templates with computer-aided protein design, PDZ domains could be engineered to specifically recognize a large number of proteins.     

Canonical and Noncanonical Sites Determine NPT2A Binding Selectivity to NHERF1 PDZ1

Na⁺/H⁺ Exchanger Regulatory Factor-1 (NHERF1) is a scaffolding protein containing 2 PDZ domains that coordinates the assembly and trafficking of transmembrane receptors and ion channels. Most target proteins harboring a C-terminus recognition motif bind more-or-less equivalently to the either PDZ domain, which contain identical core-binding motifs. However some substrates such as the type II sodium-dependent phosphate co-transporter (NPT2A), uniquely bind only one PDZ domain. We sought to define the structural determinants responsible for the specificity of interaction between NHERF1 PDZ domains and NPT2A. By performing all-atom/explicit-solvent molecular dynamics (MD) simulations in combination with biological mutagenesis, fluorescent polarization (FP) binding assays, and isothermal titration calorimetry (ITC), we found that in addition to canonical interactions of residues at 0 and -2 positions, Arg at the -1 position of NPT2A plays a critical role in association with Glu43 and His27 of PDZ1 that are absent in PDZ2. Experimentally introduced mutation in PDZ1 (Glu43Asp and His27Asn) decreased binding to NPT2A. Conversely, introduction of Asp183Glu and Asn167His mutations in PDZ2 promoted the formation of favorable interactions yielding micromolar K _Ds. The results describe novel determinants within both the PDZ domain and outside the canonical PDZ-recognition motif that are responsible for discrimination of NPT2A between two PDZ domains. The results challenge general paradigms for PDZ recognition and suggest new targets for drug development.