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.     

Energetics of peptide recognition by the second PDZ domain of human protein tyrosine phosphatase 1E

Formation of protein-protein assemblies is essential in maintaining cell structure and function. Conservation of structural motifs and binding sites is the result of evolutionary pressure for solutions compatible with both molecular economy and regulation. PDZ domains are a typical example: A conserved fold governs specificity toward recognition of C-terminal protein sequences by small sequential and/or structural deviations within a canonical binding mode. The energetic principles underlying the strength and specificity of PDZ-protein interactions are practically unknown. We use the second PDZ domain (PDZ2) of the human protein tyrosine phosphatase (hPTP1E) as a model to study the energetics of peptide binding to a class I PDZ domain. Calorimetric experiments reveal the enthalpy, entropy, and heat capacity changes accompanying PDZ2 binding to the C-terminal pentadecapeptide derived from the guanine nucleotide exchange factor RA-GEF2. Association is driven by favorable enthalpy and entropy changes below 18 degrees C. Above that temperature the entropy change opposes complex formation. Structure-based predictions poorly reproduce the observed thermodynamic profile of the PDZ-peptide complex. On the basis of MD simulations and experimental findings by others we suggest that changes in the dynamics of the PDZ domain upon peptide binding make a large contribution to the observed thermodynamic parameters. Possible impacts of subtle, ligand-induced structural "stiffening" of PDZ domains are discussed. In our hands, the C-terminal segment of the tumor suppressor APC binds much less tightly to PDZ2 than what has been proposed earlier from surface plasmon resonance experiments.     

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.     

Role of electrostatic interactions in PDZ domain ligand recognition

PDZ domains are protein-protein interaction modules that normally recognize short C-terminal peptides. The apparent requirement for a ligand with a free terminal carboxylate group has led to the proposal that electrostatic interactions with the terminus play a significant role in recognition. However, this model has been called into question by the more recent finding that PDZ domains can recognize some internal peptide motifs that occur within a specific secondary structure context. Although these motifs bind at the same interface, they lack a terminal charge. Here we have investigated the role of electrostatics in PDZ-mediated recognition in the mouse alpha1-syntrophin PDZ domain by examining the salt dependence of binding to both terminal and internal ligands and the effects of mutating a conserved basic residue previously proposed to play a role in electrostatic recognition. These studies indicate that direct electrostatic interactions with the peptide terminus do not play a significant energetic role in binding. Additional chemical modification studies of the peptide terminus support a model in which steric and hydrogen bonding complementarity play a primary role in recognition specificity. Peptides with a free carboxy terminus, or presented within a specific structural context, can satisfy these requirements.     

Signaling pathways of PDZ2 domain: a molecular dynamics interaction correlation analysis

PDZ domains are found in many signaling proteins. One of their functions is to provide scaffolds for forming membrane-associated protein complexes by binding to the carboxyl termini of their partners. PDZ domains are thought also to play a signal transduction role by propagating the information that binding has occurred to remote sites. In this study, a molecular dynamics (MD) simulation-based approach, referred to as an interaction correlation analysis, is applied to the PDZ2 domain to identify the possible signal transduction pathways. A residue correlation matrix is constructed from the interaction energy correlations between all residue pairs obtained from the MD simulations. Two continuous interaction pathways, starting at the ligand binding pocket, are identified by a hierarchical clustering analysis of the residue correlation matrix. One pathway is mainly localized at the N-terminal side of helix alpha1 and the adjacent C-terminus of loop beta1-beta2. The other pathway is perpendicular to the central beta-sheet and extends toward the side of PDZ2 domain opposite to the ligand binding pocket. The results complement previous studies based on multiple sequence analysis, NMR, and MD simulations. Importantly, they reveal the energetic origin of the long-range coupling. The PDZ2 results, as well as the earlier rhodopsin analysis, show that the interaction correlation analysis is a robust approach for determining pathways of intramolecular signal transduction.     

Characterization of domain-peptide interaction interface: a case study on the amphiphysin-1 SH3 domain

Many important protein-protein interactions are mediated by peptide recognition modular domains, such as the Src homology 3 (SH3), SH2, PDZ, and WW domains. Characterizing the interaction interface of domain-peptide complexes and predicting binding specificity for modular domains are critical for deciphering protein-protein interaction networks. Here, we propose the use of an energetic decomposition analysis to characterize domain-peptide interactions and the molecular interaction energy components (MIECs), including van der Waals, electrostatic, and desolvation energy between residue pairs on the binding interface. We show a proof-of-concept study on the amphiphysin-1 SH3 domain interacting with its peptide ligands. The structures of the human amphiphysin-1 SH3 domain complexed with 884 peptides were first modeled using virtual mutagenesis and optimized by molecular mechanics (MM) minimization. Next, the MIECs between domain and peptide residues were computed using the MM/generalized Born decomposition analysis. We conducted two types of statistical analyses on the MIECs to demonstrate their usefulness for predicting binding affinities of peptides and for classifying peptides into binder and non-binder categories. First, combining partial least squares analysis and genetic algorithm, we fitted linear regression models between the MIECs and the peptide binding affinities on the training data set. These models were then used to predict binding affinities for peptides in the test data set; the predicted values have a correlation coefficient of 0.81 and an unsigned mean error of 0.39 compared with the experimentally measured ones. The partial least squares-genetic algorithm analysis on the MIECs revealed the critical interactions for the binding specificity of the amphiphysin-1 SH3 domain. Next, a support vector machine (SVM) was employed to build classification models based on the MIECs of peptides in the training set. A rigorous training-validation procedure was used to assess the performances of different kernel functions in SVM and different combinations of the MIECs. The best SVM classifier gave satisfactory predictions for the test set, indicated by average prediction accuracy rates of 78% and 91% for the binding and non-binding peptides, respectively. We also showed that the performance of our approach on both binding affinity prediction and binder/non-binder classification was superior to the performances of the conventional MM/Poisson-Boltzmann solvent-accessible surface area and MM/generalized Born solvent-accessible surface area calculations. Our study demonstrates that the analysis of the MIECs between peptides and the SH3 domain can successfully characterize the binding interface, and it provides a framework to derive integrated prediction models for different domain-peptide systems.     

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.     

An In Silico Analysis of the Binding Modes and Binding Affinities of Small Molecule Modulators of PDZ-Peptide Interactions

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known K_i or K_d values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high K_i or K_d moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.     

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.     

Characterization of domain-peptide interaction interface: prediction of SH3 domain-mediated protein-protein interaction network in yeast by generic structure-based models

Determination of the binding specificity of SH3 domain, a peptide recognition module (PRM), is important to understand their biological functions and reconstruct the SH3-mediated protein-protein interaction network. In the present study, the SH3-peptide interactions for both class I and II SH3 domains were characterized by the intermolecular residue-residue interaction network. We developed generic MIEC-SVM models to infer SH3 domain-peptide recognition specificity that achieved satisfactory prediction accuracy. By investigating the domain-peptide recognition mechanisms at the residue level, we found that the class-I and class-II binding peptides have different binding modes even though they occupy the same binding site of SH3. Furthermore, we predicted the potential binding partners of SH3 domains in the yeast proteome and constructed the SH3-mediated protein-protein interaction network. Comparison with the experimentally determined interactions confirmed the effectiveness of our approach. This study showed that our sophisticated computational approach not only provides a powerful platform to decipher protein recognition code at the molecular level but also allows identification of peptide-mediated protein interactions at a proteomic scale. We believe that such an approach is general to be applicable to other domain-peptide 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.     

Convergent and divergent ligand specificity among PDZ domains of the LAP and zonula occludens (ZO) families

We present a detailed comparative analysis of the PDZ domains of the human LAP proteins Erbin, Densin-180, and Scribble and the MAGUK ZO-1. Phage-displayed peptide libraries and in vitro affinity assays were used to define ligand binding profiles for each domain. The analysis reveals the importance of interactions with all four C-terminal residues of the ligand, which constitute a core recognition motif, and also the role of interactions with more upstream ligand residues that support and modulate the core binding interaction. In particular, the results highlight the importance of site(-1), which interacts with the penultimate residue of ligand C termini. Site(-1) was found to be monospecific in the Erbin PDZ domain (accepts tryptophan only), bispecific in the first PDZ domain of ZO-1 (accepts tryptophan or tyrosine), and promiscuous in the Scribble PDZ domains. Furthermore, it appears that the level of promiscuity within site(-1) greatly influences the range of potential biological partners and functions that can be associated with each protein. These findings show that subtle changes in binding specificity can significantly alter the range of biological partners for PDZ domains, and the insights enhance our understanding of this diverse family of peptide-binding modules.     

Structure of PICK1 and other PDZ domains obtained with the help of self-binding C-terminal extensions

PDZ domains are protein-protein interaction modules that generally bind to the C termini of their target proteins. The C-terminal four amino acids of a prospective binding partner of a PDZ domain are typically the determinants of binding specificity. In an effort to determine the structures of a number of PDZ domains we have included appropriate four residue extensions on the C termini of PDZ domain truncation mutants, designed for self-binding. Multiple truncations of each PDZ domain were generated. The four residue extensions, which represent known specificity sequences of the target PDZ domains and cover both class I and II motifs, form intermolecular contacts in the expected manner for the interactions of PDZ domains with protein C termini for both classes. We present the structures of eight unique PDZ domains crystallized using this approach and focus on four which provide information on selectivity (PICK1 and the third PDZ domain of DLG2), binding site flexibility (the third PDZ domain of MPDZ), and peptide-domain interactions (MPDZ 12th PDZ domain). Analysis of our results shows a clear improvement in the chances of obtaining PDZ domain crystals by using this approach compared to similar truncations of the PDZ domains without the C-terminal four residue extensions.     

Comprehensive Assessment of the Relationship Between Site−2 Specificity and Helix α2 in the Erbin PDZ Domain

PDZ domains are key players in signalling pathways. These modular domains generally recognize short linear C-terminal stretches of sequences in proteins that organize the formation of complex multi-component assemblies. The development of new methodologies for the characterization of the molecular principles governing these interactions is critical to fully understand the functional diversity of the family and to elucidate biological functions for family members. Here, we applied an in vitro evolution strategy to explore comprehensively the capacity of PDZ domains for specific recognition of different amino acids at a key position in C-terminal peptide ligands. We constructed a phage-displayed library of the Erbin PDZ domain by randomizing the binding site⁻² and adjacent residues, which are all contained in helix α2, and we selected for variants binding to a panel of peptides representing all possible position⁻² residues. This approach generated insights into the basis for the common natural class I and II specificities, demonstrated an alternative basis for a rare natural class III specificity for Asp⁻², and revealed a novel specificity for Arg⁻² that has not been reported in natural PDZ domains. A structure of a PDZ-peptide complex explained the minimum requirement for switching specificity from class I ligands containing Thr/Ser⁻² to class II ligands containing hydrophobic residues at position⁻². A second structure explained the molecular basis for the specificity for ligands containing Arg⁻². Overall, the evolved PDZ variants greatly expand our understanding of site⁻² specificities and the variants themselves may prove useful as building blocks for synthetic biology.     

Energetic determinants of animal cell polarity regulator Par-3 interaction with the Par complex

The animal cell polarity regulator Par-3 recruits the Par complex (consisting of Par-6 and atypical PKC, aPKC) to specific sites on the cell membrane. Although numerous physical interactions have been reported between Par-3 and the Par complex, it is unclear how each of these interactions contributes to the overall binding. Using a purified, intact Par complex and a quantitative binding assay, here, we found that the energy required for this interaction is provided by the second and third PDZ protein interaction domains of Par-3. We show that both Par-3 PDZ domains bind to the PDZ-binding motif of aPKC in the Par complex, with additional binding energy contributed from the adjacent catalytic domain of aPKC. In addition to highlighting the role of Par-3 PDZ domain interactions with the aPKC kinase domain and PDZ-binding motif in stabilizing Par-3–Par complex assembly, our results indicate that each Par-3 molecule can potentially recruit two Par complexes to the membrane during cell polarization. These results provide new insights into the energetic determinants and structural stoichiometry of the Par-3–Par complex assembly.     

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.     

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.     

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.     

Redox-regulated affinity of the third PDZ domain in the phosphotyrosine phosphatase PTP-BL for cysteine-containing target peptides

PDZ domains are protein-protein interaction modules that are crucial for the assembly of structural and signalling complexes. They specifically bind to short C-terminal peptides and occasionally to internal sequences that structurally resemble such peptide termini. The binding of PDZ domains is dominated by the residues at the P(0) and P(-2) position within these C-terminal targets, but other residues are also important in determining specificity. In this study, we analysed the binding specificity of the third PDZ domain of protein tyrosine phosphatase BAS-like (PTP-BL) using a C-terminal combinatorial peptide phage library. Binding of PDZ3 to C-termini is preferentially governed by two cysteine residues at the P(-1) and P(-4) position and a valine residue at the P(0) position. Interestingly, we found that this binding is lost upon addition of the reducing agent dithiothrietol, indicating that the interaction is disulfide-bridge-dependent. Site-directed mutagenesis of the single cysteine residue in PDZ3 revealed that this bridge formation does not occur intermolecularly, between peptide and PDZ3 domain, but rather is intramolecular. These data point to a preference of PTP-BL PDZ3 for cyclic C-terminal targets, which may suggest a redox state-sensing role at the cell cortex.     

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.