Like other protein-protein interaction domains, PDZ domains are involved in many key cellular processes. These processes often require that specific multiprotein complexes be assembled, a task that PDZ domains accomplish by binding to specific peptide motifs in target proteins. However, a growing number of experimental studies show that PDZ domains (like other protein-protein interaction domains) can engage in a variety of interactions and bind distinct peptide motifs. Such promiscuity in ligand recognition raises intriguing questions about the molecular and thermodynamic mechanisms that can sustain it. To identify possible sources of promiscuity and selectivity underlying PDZ domain interactions, we performed molecular dynamics simulations of 20 to 25 ns on a set of 12 different PDZ domain complexes (for the proteins PSD-95, Syntenin, Erbin, GRIP, NHERF, Inad, Dishevelled, and Shank). The electrostatic, nonpolar, and configurational entropy binding contributions were evaluated using the MM/PBSA method combined with a quasi-harmonic analysis. The results revealed that PDZ domain interactions are characterized by overwhelmingly favorable nonpolar contributions and almost negligible electrostatic components, a mix that may readily sustain promiscuity. In addition, despite the structural similarity in fold and in recognition modes, the entropic and other dynamical aspects of binding were remarkably variable not only across PDZ domains but also for the same PDZ domain bound to distinct ligands. This variability suggests that entropic and dynamical components can play a role in determining selectivity either of PDZ domain interactions with peptide ligands or of PDZ domain complexes with downstream effectors.

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.

PDZ domains are protein adapter modules present in a few hundred human proteins. They play important roles in scaffolding and signal transduction. PDZ domains usually bind to the C termini of their target proteins. To assess the binding mechanism of this interaction we have performed the first in-solution kinetic study for PDZ domains and peptides corresponding to target ligands. Both PDZ3 from postsynaptic density protein 95 and PDZ2 from protein tyrosine phosphatase L1 bind their respective target peptides through an apparent A + B --> A.B mechanism without rate-limiting conformational changes. But a mutant with a fluorescent probe (Trp) outside of the binding pocket suggests that slight changes in the structure take place upon binding in protein tyrosine phosphatase-L1 PDZ2. For PDZ3 from postsynaptic density protein 95 the pH dependence of the binding reaction is consistent with a one-step mechanism with one titratable group. The salt dependence of the interaction shows that the formation of electrostatic interactions is rate-limiting for the association reaction but not for dissociation of the complex.

PDZ domains are important scaffolding modules that typically bind to the C-termini of their interaction partners. Several structures of such complexes have been solved, revealing a conserved binding site in the PDZ domain and an extended conformation of the bound peptide. A compendium of information regarding PDZ complexes demonstrates that dissimilar C-terminal peptides bind to the same PDZ domain, and different PDZ domains can bind the same peptides. A detailed understanding of the PDZ-peptide recognition is needed to elucidate this complexity. To this end, we have designed a family of docking protocols for PDZ domains (termed PDZ-DocScheme) that is based on simulated annealing molecular dynamics and rotamer optimization, and is applicable to the docking of long peptides (20-40 rotatable bonds) to both known PDZ structures and to the more complicated problem of homology models of these domains. The resulting protocol reproduces the structures of PDZ complexes with peptides 4-8 amino acids long within 1-2 A from the experimental structure when the docking is performed to the original structure. If the structure of the target PDZ domain is an apo structure or a homology model, the docking protocol yields structures within 3 A in 9 out of 12 test cases. The automated docking procedure PDZ-DocScheme can serve in the generation of a structural context for validation of PDZ domain specificity from mutagenesis and ligand binding data.

We develop a procedure for exploring the free energy landscape of protein-peptide binding at atomic detail and apply it to PDZ domain-peptide interactions. The procedure involves soft constraints on receptor proteins providing limited chain flexibility, including backbone motions. Peptide chains are left fully flexible and kept in spatial proximity of the protein through periodic boundary conditions. By extensive Monte Carlo simulations, full representative conformational ensembles at temperatures where bound and unbound states coexist are obtained. To make this approach computationally feasible, we develop an effective all-atom energy function centering on hydrophobicity, hydrogen bonding, and electrostatic interactions. Our initial focus is a set of 11 PDZ domain-peptide pairs with experimentally determined complex structures. Minimum-energy conformations are found to be highly similar to the respective native structures in eight of the cases (all-atom peptide RMSDs <6 A). Having achieved that, we turn to a more complete characterization of the bound peptide state through a clustering scheme applied on the full ensembles of peptide structures. We find a significant diversity among bound peptide conformations for several PDZ domains, in particular involving the N terminal side of the peptide chains. Our computational model is then tested further on a set of nine PDZ domain-peptide pairs where the peptides are not originally present in the experimentally determined structures. We find a similar success rate in terms of the nativeness of minimum-energy conformations. Finally, we investigate the ability of our approach to capture variations in binding affinities for different peptide sequences. This is done in particular for a set of related sequences binding to the third PDZ domain of PSD-95 with encouraging results.

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 betaA-betaB 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. 2009 Wiley-Liss, Inc.