Identification of specificity and promiscuity of PDZ domain interactions through their dynamic behavior

PDZ domains (PDZs), the most common interaction domain proteins, play critical roles in many cellular processes. PDZs perform their job by binding specific protein partners. However, they are very promiscuous, binding to more than one protein, yet selective at the same time. We examined the binding related dynamics of various PDZs to have insight about their specificity and promiscuity. We used full atomic normal mode analysis and a modified coarse‐grained elastic network model to compute the binding related dynamics. In the latter model, we introduced specificity for each single parameter constant and included the solvation effect implicitly. The modified model, referred to as specific‐Gaussian Network Model (s‐GNM), highlights some interesting differences in the conformational changes of PDZs upon binding to Class I or Class II type peptides. By clustering the residue fluctuation profiles of PDZs, we have shown: (i) binding selectivities can be discriminated from their dynamics, and (ii) the dynamics of different structural regions play critical roles for Class I and Class II specificity. s‐GNM is further tested on a dual‐specific PDZ which showed only Class I specificity when a point mutation exists on the βA‐βB loop. We observe that the binding dynamics change consistently in the mutated and wild type structures. In addition, we found that the binding induced fluctuation profiles can be used to discriminate the binding selectivity of homolog structures. These results indicate that s‐GNM can be a powerful method to study the changes in binding selectivities for mutant or homolog PDZs. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Lock and chop: A novel method for the generation of a PICK1 PDZ domain and piperidine‐based inhibitor co‐crystal structure

The membrane protein interacting with kinase C1 (PICK1) plays a trafficking role in the internalization of neuron receptors such as the amino‐3‐hydroxyl‐5‐methyl‐4‐isoxazole‐propionate (AMPA) receptor. Reduction of surface AMPA type receptors on neurons reduces synaptic communication leading to cognitive impairment in progressive neurodegenerative diseases such as Alzheimer disease. The internalization of AMPA receptors is mediated by the PDZ domain of PICK1 which binds to the GluA2 subunit of AMPA receptors and targets the receptor for internalization through endocytosis, reducing synaptic communication. We planned to block the PICK1‐GluA2 protein–protein interaction with a small molecule inhibitor to stabilize surface AMPA receptors as a therapeutic possibility for neurodegenerative diseases. Using a fluorescence polarization assay, we identified compound BIO124 as a modest inhibitor of the PICK1‐GluA2 interaction. We further tried to improve the binding affinity of BIO124 using structure‐aided drug design but were unsuccessful in producing a co‐crystal structure using previously reported crystallography methods for PICK1. Here, we present a novel method through which we generated a co‐crystal structure of the PDZ domain of PICK1 bound to BIO124.     

Rapid characterization of the binding property of HtrA2/Omi PDZ domain by validation screening of PDZ ligand library

HtrA2/Omi is a mammalian mitochondrial serine protease, and was found to have dual roles in mammalian cells, not only acting as an apoptosis-inducing protein but also maintaining mitochondrial homeostasis. PDZ domain is one of the most important protein-protein interaction modules and is involved in a variety of important cellular functions, such as signal transduction, degradation of proteins, and formation of cytoskeleton. Recently, it was reported that the PDZ domain of HtrA2/Omi might regulate proteolytic activity through its interactions with ligand proteins. In this study, we rapidly characterized the binding properties of HtrA2/Omi PDZ domain by validation screening of the PDZ ligand library with yeast two-hybrid approach. Then, we predicted its novel ligand proteins in human proteome and reconfirmed them in the yeast two-hybrid system. Finally, we analyzed the smallest networks bordered by the shortest path length between the protein pairs of novel interactions to evaluate the confidence of the identified interactions. The results revealed some novel binding properties of HtrA2/Omi PDZ domain. Besides the reported Class II PDZ motif, it also binds to Class I and Class III motifs, and exhibits restricted variability at P⁻³, which means that the P⁻³ residue is selected according to the composition of the last three residues. Seven novel ligand proteins of HtrA2/Omi PDZ domain were discovered, providing significant clues for further clarifying the roles of HtrA2/Omi. Moreover, this study proves the high efficiency and practicability of the newly developed validation screening of candidate ligand library method for binding property characterization of peptide-binding domains.     

Rapid characterization of the binding property of HtrA2/Omi PDZ domain by validation screening of PDZ ligand library

HtrA2/Omi is a mammalian mitochondrial serine protease, and was found to have dual roles in mammalian cells, not only acting as an apoptosis-inducing protein but also maintaining mitochondrial homeostasis. PDZ domain is one of the most important protein-protein interaction modules and is involved in a variety of important cellular functions, such as signal transduction, degradation of proteins,and formation of cytoskeleton. Recently, it was reported that the PDZ domain of HtrA2/Omi might regulate proteolytic activity through its interactions with ligand proteins. In this study, we rapidly characterized the binding properties of HtrA2/Omi PDZ domain by validation screening of the PDZ ligand library with yeast two-hybrid approach. Then, we predicted its novel ligand proteins in human proteome and reconfirmed them in the yeast two-hybrid system. Finally, we analyzed the smallest networks bordered by the shortest path length between the protein pairs of novel interactions to evaluate the confidence of the identified interactions. The results revealed some novel binding properties of HtrA2/Omi PDZ domain. Besides the reported Class Ⅱ PDZ motif, it also binds to Class Ⅰ and Class Ⅲ motifs, and exhibits restricted variability at P-3, which means that the P-3 residue is selected according to the composition of the last three residues. Seven novel ligand proteins of HtrA2/Omi PDZ domain were discovered, providing significant clues for further clarifying the roles of HtrA2/Omi.Moreover, this study proves the high efficiency and practicability of the newly developed validation screening of candidate ligand library method for binding property characterization of peptide-binding domains.     

A flexible docking scheme to explore the binding selectivity of PDZ domains

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTA LIGAND , we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density‐95/Dlg/ZO‐1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately.     

Zinc binding site in PICK1 is dominantly located at the CPC motif of its PDZ domain

PICK1 ( p rotein i nteracting with C k inase 1) containing a PDZ domain, a BAR domain, and two short acidic regions is as an adaptor protein that plays an important role in α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid receptor trafficking, cell morphology and migration, as well as in some diseases such as cancer, schizophrenia and pain. To better understand the physiological function of PICK1, we expressed the recombinant PICK1 and its truncated mutants in E.coli, and measured their zinc binding properties by fluorescence and competition assay. It is shown that PICK1 has one Zn²⁺-binding site. The Zn²⁺-binding properties of PICK1 are not appreciably affected after the removal of BARC domain (involving BAR domain and C-terminal acidic region). Deleting the N-terminal acidic region of NPDZ domain (involving PDZ domain and N-terminal acidic region) in PICK1 impairs its Zn²⁺-binding capacity.The mutation of the CPC (Cys-Pro-Cys) motif in the PDZ domain of PICK1 abolishes the ability of Zn²⁺-binding. In addition, Zn²⁺ can enhance the lipid-binding ability of PDZ domain as observed in both protein-lipid overlay assay and fluorescence analysis. The results presented in this report suggested that Zn²⁺ plays a regulatory role in the trafficking of PICK1 from the cytoplasm to cell membrane.     

Peptide binding to the PDZ3 domain by conformational selection

The PDZ domains, a large family of peptide recognition proteins, bind to the C‐terminal segment of membrane ion channels and receptors thereby mediating their localization. The peptide binding process is not known in detail and seems to differ among different PDZ domains. For the third PDZ domain of the synaptic protein PSD‐95 (PDZ3), a lock‐and‐key mechanism was postulated on the basis of the almost perfect overlap of the crystal structures in the presence and absence of its peptide ligand. Here, peptide binding to PDZ3 is investigated by explicit solvent molecular dynamics (MD) simulations (for a total of 1.3 μs) and the cut‐based free energy profile method for determining free energy barriers and basins. The free energy landscape of apo PDZ3 indicates that there are multiple basins within the native state. These basins differ by the relative orientation of the α2 helix and β2 strand, the two secondary structure elements that make up the peptide binding site. Only the structure with the smallest aperture of the binding site is populated in the MD simulations of the complex whose analysis reveals that the peptide ligand binds to PDZ3 by selecting one of three conformations. Thus, the dynamical information obtained by the atomistic simulations increment the static, that is, partial, picture of the PDZ3 binding mechanism based on the X‐ray crystallography data. Importantly, the simulation results show for the first time that conformational selection is a possible mechanism of peptide binding by PDZ domains in general. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

NMR evaluation of interactions between substituted-indole and PDZ1 domain of PSD-95

We synthesized small organic molecules designed as PDZ ligands. These indole-based compounds were evaluated for their interaction with the PDZ1 domain of the post-synaptic density 95 (PSD-95) protein. Three molecules were found to interact with the targeted PDZ protein by NMR. One of them showed chemical shift perturbations closely related to the natural ligands.     

Solution structure of the PDZ2 domain from cytosolic human phosphatase hPTP1E complexed with a peptide reveals contribution of the beta2-beta3 loop to PDZ domain-ligand interactions

The solution structure of the second PDZ domain from human phosphatase hPTP1E in complex with a C-terminal peptide from the guanine nucleotide exchange factor RA-GEF-2 has been determined using 2D and 3D heteronuclear NMR experiments. Compared to previously solved structures, the hPTP1E complex shows an enlarged interaction surface with the C terminus of the bound peptide. Novel contacts were found between the long structured beta2/beta3 loop of the PDZ domain and the sixth amino acid residue from the C terminus of the peptide. This work underlines the importance of the beta2/beta3 loop for ligand selection by PDZ domains.     

Solution structure of the RIM1alpha PDZ domain in complex with an ELKS1b C-terminal peptide

PDZ domains are widespread protein modules that commonly recognize C-terminal sequences of target proteins and help to organize macromolecular signaling complexes. These sequences usually bind in an extended conformation to relatively shallow grooves formed between a beta-strand and an alpha-helix in the corresponding PDZ domains. Because of this binding mode, many PDZ domains recognize primarily the C-terminal and the antepenultimate side-chains of the target protein, which commonly conform to motifs that have been categorized into different classes. However, an increasing number of PDZ domains have been found to exhibit unusual specificities. These include the PDZ domain of RIMs, which are large multidomain proteins that regulate neurotransmitter release and help to organize presynaptic active zones. The RIM PDZ domain binds to the C-terminal sequence of ELKS with a unique specificity that involves each of the four ELKS C-terminal residues. To elucidate the structural basis for this specificity, we have determined the 3D structure in solution of an RIM/ELKS C-terminal peptide complex using NMR spectroscopy. The structure shows that the RIM PDZ domain contains an unusually deep and narrow peptide-binding groove with an exquisite shape complementarity to the four ELKS C-terminal residues in their bound conformation. This groove is formed, in part, by a set of side-chains that is conserved selectively in RIM PDZ domains and that hence determines, at least in part, their unique specificity.     

Structure determination and ligand interactions of the PDZ2b domain of PTP-Bas (hPTP1E): splicing-induced modulation of ligand specificity

Two versions of the PDZ2 domain of the protein tyrosine phosphatase PTP-Bas/human PTP-BL are generated by alternative splicing. The domains differ by the insertion of five amino acid residues and their affinity to the tumour suppressor protein APC. Whereas PDZ2a is able to bind APC in the nanomolar range, PDZ2b shows no apparent interaction with APC. Here the solution structure of the splicing variant of PDZ2 with the insertion has been determined using 2D and 3D heteronuclear NMR experiments. The structural reason for the changed binding specificity is the reorientation of the loop with extra five amino acid residues, which folds back onto beta-strands two and three. In addition the side-chain of Lys32 closes the binding site of the APC binding protein and the two helices, especially alpha-helix 2, change their relative position to the protein core. Consecutively, the binding site is sterically no longer fully accessible. From the NMR-titration studies with a C-terminal APC-peptide the affinity of the peptide with the protein can be estimated as 540(+/-40)microM. The binding site encompasses part of the analogous binding site of PDZ2a as already described previously, yet specific interaction sites are abolished by the insertion of amino acids in PDZ2b. As shown by high-affinity chromatography, GST-PDZ2b and GST-PDZ2a bind to phosphatidylinositol 4,5-bisphosphate (PIP(2)) micelles with a dissociation constant K(D) of 21 microM and 55 microM, respectively. In line with these data PDZ2b binds isolated, dissolved PIP(2) and PIP(3) (phosphatidylinositol 3,4,5-trisphosphate) molecules specifically with a lower K(D) of 230(+/-20)microM as detected by NMR spectroscopy. The binding site could be located by our studies and involves the residues Ile24, Val26, Val70, Asn71, Gly77, Ala78, Glu85, Arg88, Gly91 and Gln92. PIP(2) and PIP(3) binding takes place in the groove of the PDZ domain that is normally part of the APC binding site.     

Ligand-dependent dynamics and intramolecular signaling in a PDZ domain

Allosteric communication is a fundamental process that proteins use to propagate signals from one site to functionally important distal sites. Although allostery is usually associated with multimeric proteins and enzymes, “long-range” communication may be a fundamental property of proteins. In some cases, communication occurs with minimal structural change. PDZ (post-synaptic density-95/discs large/zonula occludens-1) domains are small, protein-protein binding modules that can use multiple surfaces for docking diverse molecules. Furthermore, these domains have long-range energetic couplings that link the ligand-binding site to distal regions of the structure. Here, we show that allosteric behavior in a representative member of the PDZ domain family may be directly detected using side-chain methyl dynamics measurements. The changes in side-chain dynamics parameters in the second PDZ domain from the human tyrosine phosphatase 1E (hPTP1E) were determined upon binding a peptide target. Long-range dynamic effects were detected that correspond to previously observed pair-wise energetic couplings. These results provide one of the first experimental examples for the potential role of ps-ns timescale dynamics in propagating long-range signals within a protein, and reinforce the idea that dynamic fluctuations in proteins contribute to allosteric signal transduction.     

Proteome‐wide analysis of phospho‐regulated PDZ domain interactions

A key function of reversible protein phosphorylation is to regulate protein–protein interactions, many of which involve short linear motifs (3–12 amino acids). Motif‐based interactions are difficult to capture because of their often low‐to‐moderate affinities. Here, we describe phosphomimetic proteomic peptide‐phage display, a powerful method for simultaneously finding motif‐based interaction and pinpointing phosphorylation switches. We computationally designed an oligonucleotide library encoding human C‐terminal peptides containing known or predicted Ser/Thr phosphosites and phosphomimetic variants thereof. We incorporated these oligonucleotides into a phage library and screened the PDZ (PSD‐95/Dlg/ZO‐1) domains of Scribble and DLG1 for interactions potentially enabled or disabled by ligand phosphorylation. We identified known and novel binders and characterized selected interactions through microscale thermophoresis, isothermal titration calorimetry, and NMR. We uncover site‐specific phospho‐regulation of PDZ domain interactions, provide a structural framework for how PDZ domains accomplish phosphopeptide binding, and discuss ligand phosphorylation as a switching mechanism of PDZ domain interactions. The approach is readily scalable and can be used to explore the potential phospho‐regulation of motif‐based interactions on a large scale.     

Alternate modes of binding in two crystal structures of alkaline phosphatase-inhibitor complexes

Two high resolution crystal structures of Escherichia coli alkaline phosphatase (AP) in the presence of phosphonate inhibitors are reported. The phosphonate compounds, phosphonoacetic acid (PAA) and mercaptomethylphosphonic acid (MMP), bind competitively to AP with dissociation constants of 5.5 and 0.6 mM, respectively. The structures of the complexes of AP with PAA and MMP were refined at high resolution to crystallographic R-values of 19.0 and 17.5%, respectively. Refinement of the AP-inhibitor complexes was carried out using X-PLOR. The final round of refinement was done using SHELXL-97. Crystallographic analyses of the inhibitor complexes reveal different binding modes for the two phosphonate compounds. The significant difference in binding constants can be attributed to these alternative binding modes observed in the high resolution X-ray structures. The phosphinyl group of PAA coordinates to the active site zinc ions in a manner similar to the competitive inhibitor and product inorganic phosphate. In contrast, MMP binds with its phosphonate moiety directed toward solvent. Both enzyme-inhibitor complexes exhibit close contacts, one of which has the chemical and geometrical potential to be considered an unconventional hydrogen bond of the type C-H...X.     

High‐energy water sites determine peptide binding affinity and specificity of PDZ domains

PDZ domains have well known binding preferences for distinct C‐terminal peptide motifs. For most PDZ domains, these motifs are of the form [S/T]‐W‐[I/L/V]. Although the preference for S/T has been explained by a specific hydrogen bond interaction with a histidine in the PDZ domain and the (I/L/V) is buried in a hydrophobic pocket, the mechanism for Trp specificity at the second to last position has thus far remained unknown. Here, we apply a method to compute the free energies of explicit water molecules and predict that potency gained by Trp binding is due to a favorable release of high‐energy water molecules into bulk. The affinities of a series of peptides for both wild‐type and mutant forms of the PDZ domain of Erbin correlate very well with the computed free energy of binding of displaced waters, suggesting a direct relationship between water displacement and peptide affinity. Finally, we show a correlation between the magnitude of the displaced water free energy and the degree of Trp‐sensitivity among subtypes of the HTRA PDZ family, indicating a water‐mediated mechanism for specificity of peptide binding.     

Clustering and synaptic targeting of PICK1 requires direct interaction between the PDZ domain and lipid membranes

Protein interacting with c kinase 1 (PICK1) regulates the trafficking of receptors and ion-channels such as AMPA receptors. Traditionally, the PICK1 PDZ domain is regarded as an adaptor capable of binding to receptors trafficked by PICK1, and the lipid-binding BAR domain functions to tether PICK1 directly to membranes. Here, we show that the PICK1 PDZ domain can directly interact with lipid membranes. The PDZ domain and lipid membrane interaction is mediated by both a polybasic amino-acid cluster and a conserved 'Cys-Pro-Cys' motif located away from the peptide ligand-binding groove. Disruption of the PDZ and lipid membrane interaction totally abolished synaptic targeting of PICK1. Although mutation of the CPC motif did not affect the interaction between PICK1 and AMPA receptors, the mutant PICK1 was unable to cluster the GluR2 subunit of the receptor. In neurons, PICK1 containing the same mutation displayed dramatically compromised capacity in the trafficking of AMPA receptors. Taken together, our findings not only uncovered the novel lipid membrane-binding property of the PICK1 PDZ domain, but also provided direct evidence supporting the functional relevance of the PDZ-lipid interaction.     

Solution structure and backbone dynamics of the AF-6 PDZ domain/Bcr peptide complex

The human AF-6, a scaffold protein between cell membrane-associated proteins and the actin cytoskeleton, plays an important role in special cell-cell junctions and signal transduction. It can be phosphorylated by the protein kinase Bcr, which allows efficient binding of the C terminus of Bcr to the PDZ domain of AF-6 and consequently enhances the binding affinity of AF-6 to Ras. Formation of the AF-6, Bcr, and Ras ternary complex results in down-regulation of the Ras-mediated signal transduction pathway. To better understand the molecular basis for the recognition of the AF-6 PDZ domain and Bcr, we solve the solution structure of the AF-6 PDZ domain complexed with the C-terminal peptide of Bcr and explore the interactions between them in detail. Compared with previously reported structures, the complex exhibits a noncanonical binding mode of PDZ/peptide. Owing to the distinct residues involved in the AF-6 PDZ domain and Bcr peptide interaction, the interaction mode does not adapt to the existing classification rules that have been put forward, based on the ligand or the PDZ domain specificity. Furthermore, the PDZ domain of AF-6 can bind to the C terminus of Bcr efficiently after phosphorylation of AF-6 by the Bcr kinase. The phosphorylation may induce a conformational change of AF-6, which makes the binding surface on the PDZ domain accessible to Bcr for efficient binding. This study not only characterizes the structural details of the AF-6 PDZ/Bcr peptide complex, but also provides a potential target for future drug design and disease therapy.