PDZ domains: folding and binding

The PDZ domain is one of the most common protein-protein interaction domains in humans, and it is found in all kingdoms of life. We will review recent progress in the understanding of biophysical aspects of PDZ domains with emphasis on the folding and binding reactions. Finally, we discuss an intriguing correlation between stability and binding of peptide for PDZ2 from PTP-BL.     

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.     

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.     

Functional dynamics of PDZ binding domains: a normal-mode analysis

Postsynaptic density-95/disks large/zonula occludens-1 (PDZ) domains are relatively small (80-120 residues) protein binding modules central in the organization of receptor clusters and in the association of cellular proteins. Their main function is to bind C-terminals of selected proteins that are recognized through specific amino acids in their carboxyl end. Binding is associated with a deformation of the PDZ native structure and is responsible for dynamical changes in regions not in direct contact with the target. We investigate how this deformation is related to the harmonic dynamics of the PDZ structure and show that one low-frequency collective normal mode, characterized by the concerted movements of different secondary structures, is involved in the binding process. Our results suggest that even minimal structural changes are responsible for communication between distant regions of the protein, in agreement with recent NMR experiments. Thus, PDZ domains are a very clear example of how collective normal modes are able to characterize the relation between function and dynamics of proteins, and to provide indications on the precursors of binding/unbinding events.     

Supramodular structure and synergistic target binding of the N-terminal tandem PDZ domains of PSD-95

PDZ domain proteins play critical roles in binding, clustering and subcellular targeting of membrane receptors and ion channels. PDZ domains in multi-PDZ proteins often are arranged in groups with highly conserved spacing and intervening sequences; however, the functional significance of such tandem arrangements of PDZs is unclear. We have solved the three-dimensional structure of the first two PDZ domains of postsynaptic density protein-95 (PSD-95 PDZ1 and PDZ2), which are closely linked to each other in the PSD-95 family of scaffold proteins. The two PDZs have limited freedom of rotation and their C-terminal peptide-binding grooves are aligned with each other with an orientation preference for binding to pairs of C termini extending in the same direction. Increasing the spacing between PDZ1 and PDZ2 resulted in decreased binding between PDZ12 and its dimeric targets. The same mutation impaired the functional ability of PSD-95 to cluster Kv1.4 potassium channels in heterologous cells. The data presented provide a molecular basis for preferential binding of PSD-95 to multimeric membrane proteins with appropriate C-terminal sequences.     

Distinct ligand specificity of the Tiam1 and Tiam2 PDZ domains

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

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

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

Non-carbohydrate binding partners/domains of animal lectins

• 1.1. Protein-carbohydrate interactions are involved in a large number of biologically important recognition processes.
• 2.2. Among the participating classes of proteins lectins are defined as carbohydrate-binding proteins other than an antibody or an enzyme.
• 3.3. In addition to the essential carbohydrate-binding domain other functionally and/or structurally important sites, defined by sequence comparison or by experimental demonstration of protein-protein interactions, can be present within the lectin molecule and may be relevant for its physiological significance.
• 4.4. Sequence motifs of lectins for protein-protein interactions include amino acid structures designed for cell adhesion, growth regulatory biosignalling, intracellular routing and enzymatic activity.
• 5.5. Elucidation of the complete functional role(s) of a lectin requires accurate delineation of its carbohydrate and, if present, of its protein ligands.
• 6.6. Presence of more than one carbohydrate-binding domain in a single lectin, potential ligand properties of the glycopart of a lectin, regulatory interplay between different sites and possible interaction of complementarily shaped peptide sequences to the sugar-recognizing site should all be assessed in the quest to comprehensively explain the physiological role(s) of a lectin.

     

Degenerate specificity of PDZ domains from RhoA-specific nucleotide exchange factors PDZRhoGEF and LARG

PDZ domains are ubiquitous protein-protein interaction modules which bind short, usually carboxyterminal fragments of receptors, other integral or membrane-associated proteins, and occasionally cytosolic proteins. Their role in organizing multiprotein complexes at the cellular membrane is crucial for many signaling pathways, but the rules defining their binding specificity are still poorly understood and do not readily explain the observed diversity of their known binding partners. Two homologous RhoA-specific, multidomain nucleotide exchange factors PDZRhoGEF and LARG contain PDZ domains which show a particularly broad recognition profile, as suggested by the identification of five diverse biological targets. To investigate the molecular roots of this phenomenon, we constructed a phage display library of random carboxyterminal hexapeptides. Peptide variants corresponding to the sequences identified in library selection were synthesized and their affinities for both PDZ domains were measured and compared with those of peptides derived from sequences of natural partners. Based on the analysis of the binding sequences identified for PDZRhoGEF, we propose a sequence for an 'optimal' binding partner. Our results support the hypothesis that PDZ-peptide interactions may be best understood when one considers the sum of entropic and dynamic effects for each peptide as a whole entity, rather than preferences for specific residues at a given position.     

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.     

Interacting partners for kringle domains of plasminogen: common binding with K1 and K5 domains

We have identified MAZR and Rgl2 as specific interacting partners for kringle domains in angiostatin (K1-4) and K5 using yeast two hybrid screening. Both K1 and K1-4 have strong interaction with MAZR and Rgl2 whereas K5 only binds with Rgl2. No interaction of K2, K3, and K4 with either of these binding proteins was detected. We suggest that a common binding motif may exist near LBS-4 that is required for binding with Rgl2 but not with MAZR.     

Elucidation of the binding preferences of peptide recognition modules: SH3 and PDZ domains

Peptide-binding domains play a critical role in regulation of cellular processes by mediating protein interactions involved in signalling. In recent years, the development of large-scale technologies has enabled exhaustive studies on the peptide recognition preferences for a number of peptide-binding domain families. These efforts have provided significant insights into the binding specificities of these modular domains. Many research groups have taken advantage of this unprecedented volume of specificity data and have developed a variety of new algorithms for the prediction of binding specificities of peptide-binding domains and for the prediction of their natural binding targets. This knowledge has also been applied to the design of synthetic peptide-binding domains in order to rewire protein-protein interaction networks. Here, we describe how these experimental technologies have impacted on our understanding of peptide-binding domain specificities and on the elucidation of their natural ligands. We discuss SH3 and PDZ domains as well characterized examples, and we explore the feasibility of expanding high-throughput experiments to other peptide-binding domains.     

Characterization of molecular recognition features, MoRFs, and their binding partners

Molecular Recognition Features (MoRFs) are short, interaction-prone segments of protein disorder that undergo disorder-to-order transitions upon specific binding, representing a specific class of intrinsically disordered regions that exhibit molecular recognition and binding functions. MoRFs are common in various proteomes and occupy a unique structural and functional niche in which function is a direct consequence of intrinsic disorder. Example MoRFs collected from the Protein Data Bank (PDB) have been divided into three subtypes according to their structures in the bound state: alpha-MoRFs form alpha-helices, beta-MoRFs form beta-strands, and iota-MoRFs form structures without a regular pattern of backbone hydrogen bonds. These example MoRFs were indicated to be intrinsically disordered in the absence of their binding partners by several criteria. In this study, we used several geometric and physiochemical criteria to examine the properties of 62 alpha-, 20 beta-, and 176 iota-MoRF complex structures. Interface residues were examined by calculating differences in accessible surface area between the complex and isolated monomers. The compositions and physiochemical properties of MoRF and MoRF partner interface residues were compared to the interface residues of homodimers, heterodimers, and antigen-antibody complexes. Our analysis indicates that there are significant differences in residue composition and several geometric and physicochemical properties that can be used to discriminate, with a high degree of accuracy, between various interfaces in protein interaction data sets. Implications of these findings for the development of MoRF-partner interaction predictors are discussed. In addition, structural changes upon MoRF-to-partner complex formation were examined for several illustrative examples.     

Identification of functional domains and novel binding partners of STIM proteins

With a signal trap method, we previously identified stromal interaction molecule (STIM: originally named as SIM) as a protein, which has a signal peptide in 1996. However, recent works have accumulated evidences that STIM1 and STIM2 reside in endoplasmic reticulum (ER) and that both mainly sense ER Ca(2+) depletion, which plays an essential role in store operated calcium entry. In the present study, we extensively analyzed the domain functions and associated molecules of STIMs. A STIM1 mutant lacking the coiled-coil domains was massively expressed on the cell surface while mutants with the coiled-coil domains localized in ER. In addition, STIM1 mutants with the coiled-coil domains showed a longer half-life of proteins than those without them. These results are likely to indicate that the coiled-coil domains of STIM1 are essential for its ER-retention and its stability. Furthermore, we tried to comprehensively identify STIM1-associated molecules with mass spectrometry analysis of co-immunoprecipitated proteins for STIM1. This screening clarified that both STIM1 and STIM2 have a capacity to bind to a chaperone, calnexin as well as two protein-transporters, exportin1 and transportin1. Of importance, our result that glycosylation on STIM1 was not required for the association between STIM1 and calnexin seems to indicate that calnexin might function on STIM1 beyond a chaperone protein. Further information concerning regulatory mechanisms for STIM proteins including the data shown here will provide a model of Ca(2+) control as well as a useful strategy to develop therapeutic drugs for intracellular Ca(2+)-related diseases including inflammation and allergy.     

Syntenin-syndecan binding requires syndecan-synteny and the co-operation of both PDZ domains of syntenin

Syntenin is an adaptor-like molecule that binds to the cytoplasmic domains of all four vertebrate syndecans. Syntenin-syndecan binding involves the C-terminal part of syntenin that contains a tandem of PDZ domains. Here we provide evidence that each PDZ domain of syntenin can interact with a syndecan. Isolated or combined mutations of the carboxylate binding lysines in the inter-betaAbetaB loops and of the alphaB1 residues in either one or both the PDZ domains of syntenin all reduce syntenin-syndecan binding in yeast two-hybrid, blot-overlay, and surface plasmon resonance assays. PDZ2 mutations have more pronounced effects on binding than PDZ1 mutations, but complete abrogation of syntenin-syndecan binding requires the combination of both the lysine and the alphaB1 mutations in both the PDZ domains of syntenin. Isothermal calorimetric titration of syntenin with syndecan peptide reveals the presence of two binding sites in syntenin. Yet, unlike a tandem of two PDZ2 domains and a reconstituted PDZ1+PDZ2 tandem, a tandem of two PDZ1 domains and isolated PDZ1 or PDZ2 domains do not interact with syndecan bait. We conclude to a co-operative binding mode whereby neither of these two PDZ domains is sufficient by itself but where PDZ2 functions as a "major" or "high affinity" syndecan binding domain, and PDZ1 functions as an "accessory" or "low affinity" syndecan binding domain. The paired, but not the isolated PDZ domains of syntenin bind also strongly to the immobilized cytoplasmic domains of neurexin and B-class ephrins. By inference, these data suggest a model whereby recruitment of syntenin to membrane surfaces requires two compatible types of bait that are in "synteny" (occurring together in location) and engages both PDZ domains of syntenin. The synteny of compatible bait may result from the assemblies and co-assemblies of syndecans and other similarly suited partners in larger supramolecular complexes. In general, an intramolecular combination of PDZ domains that are weak, taken individually, would appear to be designed to detect rather than drive the formation of specific molecular assemblies.     

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.     

Rat RAMP domains involved in adrenomedullin binding specificity

When coexpressed with receptor activity-modifying protein (RAMP)2 or -3, calcitonin receptor-like receptor (CRLR) functions as an adrenomedullin (AM) receptor (CRLR/RAMP2 or -3). Coexpression of rat (r)CRLR with rRAMP deletion mutants in HEK293T cells revealed that deletion of residues 93-99 from rRAMP2 or residues 58-64 from rRAMP3 significantly inhibits high-affinity [125I]AM binding and AM-evoked cAMP production, despite full cell surface expression of the receptor heterodimer. Apparently, these two seven-residue segments are key determinants of high-affinity agonist binding to rAM receptors and of receptor functionality. Consequently, their deletion yields peptides that are able to serve as negative regulators of AM receptor function.