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.     

Functional genomics of intracellular peptide recognition domains with combinatorial biology methods

Phage-displayed peptide libraries have been used to identify specific ligands for peptide-binding domains that mediate intracellular protein-protein interactions. These studies have provided significant insights into the specificities of particular domains. For PDZ domains that recognize C-terminal sequences, the information has proven useful in identifying natural binding partners from genomic databases. For SH3 domains that recognize internal proline-rich motifs, the results of database searches with phage-derived ligands have been compared with the results of yeast-two-hybrid experiments to produce overlap networks that reliably predict natural protein-protein interactions. In addition, libraries of phage-displayed PDZ and SH3 domains have been used to identify the residues responsible for ligand recognition, and also to engineer domains with altered specificities.     

Plasticity of PDZ domains in ligand recognition and signaling

The PDZ domain is a protein-protein interacting module that plays an important role in the organization of signaling complexes. The recognition of short intrinsically disordered C-terminal peptide motifs is the archetypical PDZ function, but the functional repertoire of this versatile module also includes recognition of internal peptide sequences, dimerization and phospholipid binding. The PDZ function can be tuned by various means such as allosteric effects, changes of physiological buffer conditions and phosphorylation of PDZ domains and/or ligands, which poses PDZ domains as dynamic regulators of cell signaling. This review is focused on the plasticity of the PDZ interactions.     

Binding free energy landscape of domain-peptide interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step.     

Analysis of PDZ domain-ligand interactions using carboxyl-terminal phage display

PDZ domains mediate protein-protein interactions at specialized subcellular sites, such as epithelial cell tight junctions and neuronal post-synaptic densities. Because most PDZ domains bind extreme carboxyl-terminal sequences, the phage display method has not been amenable to the study of PDZ domain binding specificities. For the first time, we demonstrate the functional display of a peptide library fused to the carboxyl terminus of the M13 major coat protein. We used this library to analyze carboxyl-terminal peptide recognition by two PDZ domains. For each PDZ domain, the library provided specific ligands with sub-micromolar binding affinities. Synthetic peptides and homology modeling were used to dissect and rationalize the binding interactions. Our results establish carboxyl-terminal phage display as a powerful new method for mapping PDZ domain binding specificity.     

Characterization of PDZ domain‐peptide interaction interface based on energetic patterns

PDZ domain is one of the abundant modular domains that recognize short peptide sequences to mediate protein–protein interactions. To decipher the binding specificity of PDZ domain, we analyzed the interactions between 11 mouse PDZ domains and 217 peptides using a method called MIECSVM, which energetically characterizes the domain‐peptide interaction using molecular interaction energy components (MIECs) and predicts binding specificity using support vector machine (SVM). Cross‐validation and leave‐one‐domain‐out test showed that the MIEC‐SVM using all 44 PDZ‐peptide residue pairs at the interaction interface outperformed the sequence‐based methods in the literature. A further feature (residue pair) selection procedure illustrated that 16 residue pairs were uninformative to the binding specificity, even though they contributed significantly (～50%) to the binding energy. If only using the 28 informative residue pairs, the performance of the MIEC‐SVM on predicting the PDZ binding specificity was significantly improved. This analysis suggests that the informative and uninformative residue interactions between the PDZ domain and the peptide may represent those contributing to binding specificity and affinity, respectively. We performed additional structural and energetic analyses to shed light on understanding how the PDZ‐peptide recognition is established. The success of the MIEC‐SVM method on PDZ domains in this study and SH3 domains in our previous studies illustrates its generality on characterizing protein‐ peptide interactions and understanding protein recognition from a structural and energetic viewpoint.     

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.     

Phage display can select over-hydrophobic sequences that may impair prediction of natural domain-peptide interactions

MOTIVATION: The phage display peptide selection approach is widely used for defining binding specificities of globular domains. PDZ domains recognize partner proteins via C-terminal motifs and are often used as a model for interaction predictions. Here, we investigated to which extent phage display data that were recently published for 54 human PDZ domains can be applied to the prediction of human PDZ-peptide interactions. RESULTS: Promising predictions were obtained for one-third of the 54 PDZ domains. For the other two-thirds, we detected in the phage display peptides an important bias for hydrophobic amino acids that seemed to impair correct predictions. Therefore, phage display-selected peptides may be over-hydrophobic and of high affinity, while natural interaction motifs are rather hydrophilic and mostly combine low affinity with high specificity. We suggest that potential amino acid composition bias should systematically be investigated when applying phage display data to the prediction of specific natural domain-linear motif interactions.     

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.     

Synthesis and screening of support-bound combinatorial peptide libraries with free C-termini: determination of the sequence specificity of PDZ domains

Preparation of support-bound combinatorial peptide libraries with free C-termini has been challenging in the past because solid-phase peptide synthesis usually starts from the C-terminus, which must be covalently attached to the solid support. In this work, we have developed a general methodology to synthesize and screen one-bead-one-compound peptide libraries containing free C-termini. TentaGel microbeads (90 mum) were spatially segregated into outer and inner layers, and peptides were synthesized on the beads in the conventional C --> N manner, with their C-termini attached to the support through an ester linkage on the bead surface but through an amide bond in the bead interior. The surface peptides were cyclized between their N-terminal amine and a carboxyl group installed at a C-terminal linker sequence, while the internal peptides were kept in the linear form. Base hydrolysis of the ester linkage in the cyclic peptides regenerated linear peptides that contained a free alpha-carboxyl group at their C-termini but remained covalently attached to the resin via the N-termini ("inverted" peptides). An inverted peptide library containing five random residues (theoretical diversity of 3.2 x 10 (6)) was synthesized and screened for binding to four postsynaptic density-95/discs large/zona occluden-1 (PDZ) domains of sodium-hydrogen exchanger regulatory factor-1 (NHERF1) and channel-interacting PDZ domain protein (CIPP). The identity of the binding peptides was determined by sequencing the linear encoding peptides inside the bead by partial Edman degradation/mass spectrometry. Consensus recognition motifs were identified for the PDZ domains, and representative peptides were resynthesized and confirmed for binding to their cognate PDZ domains. This method should be generally applicable to all PDZ domains as well as other protein domains and enzymes that recognize the C-terminus of their target proteins.     

Characterization of PDZ domain-peptide interaction interface based on energetic patterns

PDZ domain is one of the abundant modular domains that recognize short peptide sequences to mediate protein-protein interactions. To decipher the binding specificity of PDZ domain, we analyzed the interactions between 11 mouse PDZ domains and 2387 peptides using a method called MIEC-SVM, which energetically characterizes the domain-peptide interaction using molecular interaction energy components (MIECs) and predicts binding specificity using support vector machine (SVM). Cross-validation and leave-one-domain-out test showed that the MIEC-SVM using all 44 PDZ-peptide residue pairs at the interaction interface outperformed the sequence-based methods in the literature. A further feature (residue pair) selection procedure illustrated that 16 residue pairs were uninformative to the binding specificity, even though they contributed significantly (~50%) to the binding energy. If only using the 28 informative residue pairs, the performance of the MIEC-SVM on predicting the PDZ binding specificity was significantly improved. This analysis suggests that the informative and uninformative residue interactions between the PDZ domain and the peptide may represent those contributing to binding specificity and affinity, respectively. We performed additional structural and energetic analyses to shed light on understanding how the PDZ-peptide recognition is established. The success of the MIEC-SVM method on PDZ domains in this study and SH3 domains in our previous studies illustrates its generality on characterizing protein-peptide interactions and understanding protein recognition from a structural and energetic viewpoint.     

PDZ domains facilitate binding of high temperature requirement protease A (HtrA) and tail-specific protease (Tsp) to heterologous substrates through recognition of the small stable RNA A (ssrA)-encoded peptide

The Escherichia coli protease HtrA has two PDZ domains, and sequence alignments predict that the E. coli protease Tsp has a single PDZ domain. PDZ domains are composed of short sequences (80-100 amino acids) that have been implicated in a range of protein:protein interactions. The PDZ-like domain of Tsp may be involved in binding to the extreme COOH-terminal sequence of its substrate, whereas the HtrA PDZ domains are involved in subunit assembly and are predicted to be responsible for substrate binding and subsequent translocation into the active site. E. coli has a system of protein quality control surveillance mediated by the ssrA-encoded peptide tagging system. This system tags misfolded proteins or protein fragments with an 11-amino acid peptide that is recognized by a battery of cytoplasmic and periplasmic proteases as a degradation signal. Here we show that both HtrA and Tsp are able to recognize the ssrA-encoded peptide tag with apparent K(D) values of approximately 5 and 390 nm, respectively, and that their PDZ-like domains mediate this recognition. Fusion of the ssrA-encoded peptide tag to the COOH terminus of a heterologous protein (glutathione S-transferase) renders it sensitive to digestion by Tsp but not HtrA. These observations support the prediction that the HtrA PDZ domains facilitate substrate binding and the differential proteolytic responses of HtrA and Tsp to SsrA-tagged glutathione S-transferase are interpreted in terms of the structure of HtrA.     

Domain swapping within PDZ2 is responsible for dimerization of ZO proteins

ZO-1 is a multidomain protein involved in cell-cell junctions and contains three PDZ domains, which are necessary for its function in vivo. PDZ domains play a central role in assembling diverse protein complexes through their ability to recognize short peptide motifs on other proteins. We determined the structure of the second of the three PDZ domains of ZO-1, which is known to promote dimerization as well as bind to C-terminal sequences on connexins. The dimer is stabilized by extensive symmetrical domain swapping of beta-strands, which is unlike any other known mechanism of PDZ dimerization. The canonical peptide-binding groove remains intact in both subunits of the PDZ2 dimer and is created by elements contributed from both monomers. This unique structure reveals an additional example of how PDZ domains dimerize and has multiple implications for both peptide binding and oligomerization in vivo.     

Structures of a human papillomavirus (HPV) E6 polypeptide bound to MAGUK proteins: mechanisms of targeting tumor suppressors by a high-risk HPV oncoprotein

Human papillomavirus (HPV) E6 oncoprotein targets certain tumor suppressors such as MAGI-1 and SAP97/hDlg for degradation. A short peptide at the C terminus of E6 interacts specifically with the PDZ domains of these tumor suppressors, which is a property unique to high-risk HPVs that are associated with cervical cancer. The detailed recognition mechanisms between HPV E6 and PDZ proteins are unclear. To understand the specific binding of cellular PDZ substrates by HPV E6, we have solved the crystal structures of the complexes containing a peptide from HPV18 E6 bound to three PDZ domains from MAGI-1 and SAP97/Dlg. The complex crystal structures reveal novel features of PDZ peptide recognition that explain why high-risk HPV E6 can specifically target these cellular tumor suppressors for destruction. Moreover, a new peptide-binding loop on these PDZs is identified as interacting with the E6 peptide. Furthermore, we have identified an arginine residue, unique to high-risk HPV E6 but outside the canonical core PDZ recognition motif, that plays an important role in the binding of the PDZs of both MAGI-I and SAP97/Dlg, the mutation of which abolishes E6's ability to degrade the two proteins. Finally, we have identified a dimer form of MAGI-1 PDZ domain 1 in the cocrystal structure with E6 peptide, which may have functional relevance for MAGI-1 activity. In addition to its novel insights into the biochemistry of PDZ interactions, this study is important for understanding HPV-induced oncogenesis; this could provide a basis for developing antiviral and anticancer compounds.     

PDZ domains and their binding partners: structure, specificity, and modification

PDZ domains are abundant protein interaction modules that often recognize short amino acid motifs at the C-termini of target proteins. They regulate multiple biological processes such as transport, ion channel signaling, and other signal transduction systems. This review discusses the structural characterization of PDZ domains and the use of recently emerging technologies such as proteomic arrays and peptide libraries to study the binding properties of PDZ-mediated interactions. Regulatory mechanisms responsible for PDZ-mediated interactions, such as phosphorylation in the PDZ ligands or PDZ domains, are also discussed. A better understanding of PDZ protein-protein interaction networks and regulatory mechanisms will improve our knowledge of many cellular and biological processes.     

Proteome scanning to predict PDZ domain interactions using support vector machines

Background  

PDZ domains mediate protein-protein interactions involved in important biological processes through the recognition of short linear motifs in their target proteins. Two recent independent studies have used protein microarray or phage display technology to detect PDZ domain interactions with peptide ligands on a large scale. Several computational predictors of PDZ domain interactions have been developed, however they are trained using only protein microarray data and focus on limited subsets of PDZ domains. An accurate predictor of genomic PDZ domain interactions would allow the proteomes of organisms to be scanned for potential binders. Such an application would require an accurate and precise predictor to avoid generating too many false positive hits given the large amount of possible interactors in a given proteome. Once validated these predictions will help to increase the coverage of current PDZ domain interaction networks and further our understanding of the roles that PDZ domains play in a variety of biological processes.     

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.     

Single-amino acid substitutions alter the specificity and affinity of PDZ domains for their ligands

PDZ domains are modular protein-protein interaction domains that bind to specific C-terminal sequences of membrane proteins and/or to other PDZ domains. Certain PDZ domains in PSD-95 and syntrophins interact with C-terminal peptide ligands and heterodimerize with the extended nNOS PDZ domain. The capacity to interact with nNOS correlates with the presence of a Lys residue in the carboxylate- binding loop of these PDZ domains. Here, we report that substitution of an Arg for Lys-165 in PSD-95 PDZ2 disrupted its interaction with nNOS, but not with the C terminus of the Shaker-type K(+) channel Kv1.4. The same mutation affected nNOS binding to alpha1- and beta1-syntrophin PDZ domains to a lesser extent, due in part to the stabilizing effect of tertiary interactions with the canonical nNOS PDZ domain. PDZ domains with an Arg in the carboxylate-binding loop do not bind nNOS; however, substitution with Lys or Ala was able to confer nNOS binding. Our results indicate that the carboxylate-binding loop Lys or Arg is a critical determinant of nNOS binding and that the identity of this residue can profoundly alter one mode of PDZ recognition without affecting another. We also analyzed the effects of mutating Asp-143, a residue in the alphaB helix of alpha1-syntrophin that forms a tertiary contact with the nNOS PDZ domain. This residue is important for both nNOS and C-terminal peptide binding and confers a preference for peptides with a positively charged residue at position -4. On this basis, we have identified the C terminus of the Kir2.1 channel as a possible binding partner for syntrophin PDZ domains. Together, our results demonstrate that single-amino acid substitutions alter the specificity and affinity of PDZ domains for their ligands.     

The emerging contribution of sequence context to the specificity of protein interactions mediated by PDZ domains

The canonical binding mode of PDZ domains to target motifs involves a small interface, unlikely to fully account for PDZ-target interaction specificities. Here, we review recent work on sequence context, defined as the regions surrounding not only the PDZ domains but also their target motifs. We also address the theoretical problem of defining the core of PDZ domains and the practical issue of designing PDZ constructs. Sequence context is found to introduce structural diversity, to impact the stability and solubility of constructs, and to deeply influence binding affinity and specificity, thereby increasing the difficulty of predicting PDZ-motif interactions. We expect that sequence context will have similar importance for other protein interactions mediated by globular domains binding to short linear motifs.     

Binding of Vac14 to neuronal nitric oxide synthase: Characterisation of a new internal PDZ-recognition motif

PDZ domains mediate protein interactions primarily through either classical recognition of carboxyl-terminal motifs or PDZ/PDZ domain associations. Several studies have also described internal modes of PDZ recognition, most of which depend on β-finger structures. Here, we describe a novel interaction between the PDZ domain of nNOS and Vac14, the activator of the PtdIns(3)P 5-kinase PIKfyve. Binding assays using various Vac14 deletion constructs revealed a β-finger independent interaction that is based on a novel internal motif. Mutational analyses reveal essential residues within the motif allowing us to define a new type of PDZ domain interaction.