An allosteric intramolecular PDZ-PDZ interaction modulates PTP-BL PDZ2 binding specificity

The Human PDZome 2.0: Characterization of a New Resource to Test for PDZ Interactions by Yeast Two-Hybrid

PSD95-disc large-zonula occludens (PDZ) domains are globular modules of 80-90 amino acids that co-evolved with multicellularity. They commonly bind to carboxy-terminal sequences of a plethora of membrane-associated proteins and influence their trafficking and signaling. We previously built a PDZ resource (PDZome) allowing us to unveil human PDZ interactions by Yeast two-hybrid. Yet, this resource is incomplete according to the current knowledge on the human PDZ proteome. Here we built the PDZome 2.0 library for Yeast two-hybrid, based on a PDZ library manually curated from online resources. The PDZome2.0 contains 305 individual clones (266 PDZ domains in isolation and 39 tandems), for which all boundaries were designed based on available PDZ structures. Using as bait the E6 oncoprotein from HPV16, a known promiscuous PDZ interactor, we show that PDZome 2.0 outperforms the previous resource.

Investigating the allosteric response of the PICK1 PDZ domain to different ligands with all-atom simulations

The PDZ family is comprised of small modular domains that play critical roles in the allosteric modulation of many cellular signaling processes by binding to the C-terminal tail of different proteins. As dominant modular proteins that interact with a diverse set of peptides, it is of particular interest to explore how different binding partners induce different allosteric effects on the same PDZ domain. Because the PICK1 PDZ domain can bind different types of ligands, it is an ideal test case to answer this question and explore the network of interactions that give rise to dynamic allostery. Here, we use all-atom molecular dynamics simulations to explore dynamic allostery in the PICK1 PDZ domain by modeling two PICK1 PDZ systems: PICK1 PDZ-DAT and PICK1 PDZ-GluR2. Our results suggest that ligand binding to the PICK1 PDZ domain induces dynamic allostery at the αA helix that is similar to what has been observed in other PDZ domains. We found that the PICK1 PDZ-ligand distance is directly correlated with both dynamic changes of the αA helix and the distance between the αA helix and βB strand. Furthermore, our work identifies a hydrophobic core between DAT/GluR2 and I35 as a key interaction in inducing such dynamic allostery. Finally, the unique interaction patterns between different binding partners and the PICK1 PDZ domain can induce unique dynamic changes to the PICK1 PDZ domain. We suspect that unique allosteric coupling patterns with different ligands may play a critical role in how PICK1 performs its biological functions in various signaling networks.

Atomic resolution protein allostery from the multi-state structure of a PDZ domain

Recent methodological advances in solution NMR allow the determination of multi-state protein structures and provide insights into structurally and dynamically correlated protein sites at atomic resolution. This is demonstrated in the present work for the well-studied PDZ2 domain of protein human tyrosine phosphatase 1E for which protein allostery had been predicted. Two-state protein structures were calculated for both the free form and in complex with the RA-GEF2 peptide using the exact nuclear Overhauser effect (eNOE) method. In the apo protein, an allosteric conformational selection step comprising almost 60% of the domain was detected with an "open" ligand welcoming state and a "closed" state that obstructs the binding site by changing the distance between the β-sheet 2, α-helix 2, and sidechains of residues Lys38 and Lys72. The observed induced fit-type apo-holo structural rearrangements are in line with the previously published evolution-based analysis covering ~25% of the domain with only a partial overlap with the protein allostery of the open form. These presented structural studies highlight the presence of a dedicated highly optimized and complex dynamic interplay of the PDZ2 domain owed by the structure-dynamics landscape.

SenseNet, a tool for analysis of protein structure networks obtained from molecular dynamics simulations

Computational methods play a key role for investigating allosteric mechanisms in proteins, with the potential of generating valuable insights for innovative drug design. Here we present the SenseNet (“Structure ENSEmble NETworks”) framework for analysis of protein structure networks, which differs from established network models by focusing on interaction timelines obtained by molecular dynamics simulations. This approach is evaluated by predicting allosteric residues reported by NMR experiments in the PDZ2 domain of hPTP1e, a reference system for which previous computational predictions have shown considerable variance. We applied two models based on the mutual information between interaction timelines to estimate the conformational influence of each residue on its local environment. In terms of accuracy our prediction model is comparable to the top performing model published for this system, but by contrast benefits from its independence from NMR structures. Our results are complementary to experimental data and the consensus of previous predictions, demonstrating the potential of our new analysis tool SenseNet. Biochemical interpretation of our model suggests that allosteric residues in the PDZ2 domain form two distinct clusters of contiguous sidechain surfaces. SenseNet is provided as a plugin for the network analysis software Cytoscape, allowing for ease of future application and contributing to a system of compatible tools bridging the fields of system and structural biology.

Allosterism in the PDZ Family

Dynamic allosterism allows the propagation of signal throughout a protein. The PDZ (PSD-95/Dlg1/ZO-1) family has been named as a classic example of dynamic allostery in small modular domains. While the PDZ family consists of more than 200 domains, previous efforts have primarily focused on a few well-studied PDZ domains, including PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ. Taken together, experimental and computational studies have identified regions of these domains that are dynamically coupled to ligand binding. These regions include the αA helix, the αB lower-loop, and the αC helix. In this review, we summarize the specific residues on the αA helix, the αB lower-loop, and the αC helix of PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ that have been identified as participants in dynamic allostery by either experimental or computational approaches. This review can serve as an index for researchers to look back on the previously identified allostery in the PDZ family. Interestingly, our summary of previous work reveals clear consistencies between the domains. While the PDZ family has a low sequence identity, we show that some of the most consistently identified allosteric residues within PTP-BL PDZ2 and PSD-95 PDZ3 domains are evolutionarily conserved. These residues include A46/A347, V61/V362, and L66/L367 on PTP-BL PDZ2 and PSD-95 PDZ3, respectively. Finally, we expose a need for future work to explore dynamic allostery within (1) PDZ domains with multiple binding partners and (2) multidomain constructs containing a PDZ domain.

Supertertiary protein structure affects an allosteric network

Protein function can be allosterically regulated by changes in structure or dynamics. PDZ domains are classic examples for studies of allostery in single protein domains. However, PDZ domains are often found in multidomain proteins; in particular, PDZ3 is located in a supramodule containing three domains. The allosteric network in PDZ3 has never been studied in the presence of the adjacent domains. Here we map the allosteric network for a PDZ3:ligand complex, both in isolation and in the context of a supramodule. We demonstrate that the allosteric network is highly dependent on this supertertiary structure, with broad implications for studies of allostery in single domains.

The notion that protein function is allosterically regulated by structural or dynamic changes in proteins has been extensively investigated in several protein domains in isolation. In particular, PDZ domains have represented a paradigm for these studies, despite providing conflicting results. Furthermore, it is still unknown how the association between protein domains in supramodules, consitituting so-called supertertiary structures, affects allosteric networks. Here, we experimentally mapped the allosteric network in a PDZ:ligand complex, both in isolation and in the context of a supramodular structure, and show that allosteric networks in a PDZ domain are highly dependent on the supertertiary structure in which they are present. This striking sensitivity of allosteric networks to the presence of adjacent protein domains is likely a common property of supertertiary structures in proteins. Our findings have general implications for prediction of allosteric networks from primary and tertiary structures and for quantitative descriptions of allostery.

Conformational flexibility of GRASPs and their constituent PDZ subdomains reveals structural basis of their promiscuous interactome

The Golgi complex is a central component of the secretory pathway, responsible for several critical cellular functions in eukaryotes. The complex is organized by the Golgi matrix that includes the Golgi reassembly and stacking protein (GRASP), which was shown to be involved in cisternae stacking and lateral linkage in metazoan. GRASPs also have critical roles in other processes, with an unusual ability to interact with several different binding partners. The conserved N terminus of the GRASP family includes two PSD-95, DLG, and ZO-1 (PDZ) domains. Previous crystallographic studies of orthologues suggest that PDZ1 and PDZ2 have similar conformations and secondary structure content. However, PDZ1 alone mediates nearly all interactions between GRASPs and their partners. In this work, NMR, synchrotron radiation CD, and molecular dynamics (MD) were used to examine the structure, flexibility, and stability of the two constituent PDZ domains. GRASP PDZs are structured in an unusual β₃ α₁ β₄ β₅ α₂ β₆ β₁ β₂ secondary structural arrangement and NMR data indicate that the PDZ1 binding pocket is formed by a stable β₂ -strand and a more flexible and unstable α₂ -helix, suggesting an explanation for the higher PDZ1 promiscuity. The conformational free energy profiles of the two PDZ domains were calculated using MD simulations. The data suggest that, after binding, the protein partner significantly reduces the conformational space that GRASPs can access by stabilizing one particular conformation, in a partner-dependent fashion. The structural flexibility of PDZ1, modulated by PDZ2, and the coupled, coordinated movement between the two PDZs enable GRASPs to interact with multiple partners, allowing them to function as promiscuous, multitasking proteins.

The binding affinity of PTPN13's tandem PDZ2/3 domain is allosterically modulated

Background

Protein tyrosine phosphatase PTPN13, also known as PTP-BL in mice, is a large multi-domain non-transmembrane scaffolding protein with a molecular mass of 270 kDa. It is involved in the regulation of several cellular processes such as cytokinesis and actin-cytoskeletal rearrangement. The modular structure of PTPN13 consists of an N-terminal KIND domain, a FERM domain, and five PDZ domains, followed by a C-terminal protein tyrosine phosphatase domain. PDZ domains are among the most abundant protein modules and they play a crucial role in signal transduction of protein networks.

Results

Here, we have analysed the binding characteristics of the isolated PDZ domains 2 and 3 from PTPN13 and compared them to the tandem domain PDZ2/3, which interacts with 12 C-terminal residues of the tumour suppressor protein of APC, using heteronuclear multidimensional NMR spectroscopy. Furthermore, we could show for the first time that PRK2 is a weak binding partner of PDZ2 and we demonstrate that the presence of PDZ3 alters the binding affinity of PDZ2 for APC, suggesting an allosteric effect and thereby modulating the binding characteristics of PDZ2. A HADDOCK-based molecular model of the PDZ2/3 tandem domain from PTPN13 supports these results.

Conclusions

Our study of tandem PDZ2/3 in complex with APC suggests that the interaction of PDZ3 with PDZ2 induces an allosteric modulation within PDZ2 emanating from the back of the domain to the ligand binding site. Thus, the modified binding preference of PDZ2 for APC could be explained by an allosteric effect and provides further evidence for the pivotal function of PDZ2 in the PDZ123 domain triplet within PTPN13.

PDZ Domains as Drug Targets

Protein-protein interactions within protein networks shape the human interactome, which often is promoted by specialized protein interaction modules, such as the postsynaptic density-95 (PSD-95), discs-large, zona occludens 1 (ZO-1) (PDZ) domains. PDZ domains play a role in several cellular functions, from cell-cell communication and polarization, to regulation of protein transport and protein metabolism. PDZ domain proteins are also crucial in the formation and stability of protein complexes, establishing an important bridge between extracellular stimuli detected by transmembrane receptors and intracellular responses. PDZ domains have been suggested as promising drug targets in several diseases, ranging from neurological and oncological disorders to viral infections. In this review, the authors describe structural and genetic aspects of PDZ-containing proteins and discuss the current status of the development of small-molecule and peptide modulators of PDZ domains. An overview of potential new therapeutic interventions in PDZ-mediated protein networks is also provided.

Structural and functional characterization of the PDZ domain of the human phosphatase PTPN3 and its interaction with the human papillomavirus E6 oncoprotein

The human protein tyrosine phosphatase non-receptor type 3 (PTPN3) is a PDZ (PSD-95/Dlg/ZO-1) domain-containing phosphatase with a tumor-suppressive or a tumor-promoting role in many cancers. Interestingly, the high-risk genital human papillomavirus (HPV) types 16 and 18 target the PDZ domain of PTPN3. The presence of a PDZ binding motif (PBM) on E6 confers interaction with a number of different cellular PDZ domain-containing proteins and is a marker of high oncogenic potential. Here, we report the molecular basis of interaction between the PDZ domain of PTPN3 and the PBM of the HPV E6 protein. We combined biophysical, NMR and X-ray experiments to investigate the structural and functional properties of the PDZ domain of PTPN3. We showed that the C-terminal sequences from viral proteins encompassing a PBM interact with PTPN3-PDZ with similar affinities to the endogenous PTPN3 ligand MAP kinase p38γ. PBM binding stabilizes the PDZ domain of PTPN3. We solved the X-ray structure of the PDZ domain of PTPN3 in complex with the PBM of the HPV E6 protein. The crystal structure and the NMR chemical shift mapping of the PTPN3-PDZ/peptide complex allowed us to pinpoint the main structural determinants of recognition of the C-terminal sequence of the E6 protein and the long-range perturbations induced upon PBM binding.

Emerging Themes in PDZ Domain Signaling: Structure, Function, and Inhibition

Post-synaptic density-95, disks-large and zonula occludens-1 (PDZ) domains are small globular protein-protein interaction domains widely conserved from yeast to humans. They are composed of ∼90 amino acids and form a classical two α-helical/six β-strand structure. The prototypical ligand is the C-terminus of partner proteins; however, they also bind internal peptide sequences. Recent findings indicate that PDZ domains also bind phosphatidylinositides and cholesterol. Through their ligand interactions, PDZ domain proteins are critical for cellular trafficking and the surface retention of various ion channels. In addition, PDZ proteins are essential for neuronal signaling, memory, and learning. PDZ proteins also contribute to cytoskeletal dynamics by mediating interactions critical for maintaining cell-cell junctions, cell polarity, and cell migration. Given their important biological roles, it is not surprising that their dysfunction can lead to multiple disease states. As such, PDZ domain-containing proteins have emerged as potential targets for the development of small molecular inhibitors as therapeutic agents. Recent data suggest that the critical binding function of PDZ domains in cell signaling is more than just glue, and their binding function can be regulated by phosphorylation or allosterically by other binding partners. These studies also provide a wealth of structural and biophysical data that are beginning to reveal the physical features that endow this small modular domain with a central role in cell signaling.

Molecular Basis of Class III Ligand Recognition by PDZ3 in Murine Protein Tyrosine Phosphatase PTPN13

Protein tyrosine phosphatase PTPN13, also known as PTP-BL in mice, represents a large multi-domain non-transmembrane scaffolding protein that contains five consecutive PDZ domains. Here, we report the solution structures of the extended murine PTPN13 PDZ3 domain in its apo form and in complex with its physiological ligand, the carboxy-terminus of protein kinase C-related kinase-2 (PRK2), determined by multidimensional NMR spectroscopy. Both in its ligand-free state and when complexed to PRK2, PDZ3 of PTPN13 adopts the classical compact, globular D/E fold. PDZ3 of PTPN13 binds five carboxy-terminal amino acids of PRK2 via a groove located between the EB-strand and the DB-helix. The PRK2 peptide resides in the canonical PDZ3 binding cleft in an elongated manner and the amino acid side chains in position P0 and P-2, cysteine and aspartate, of the ligand face the groove between EB-strand and DB-helix, whereas the PRK2 side chains of tryptophan and alanine located in position P-1 and P-3 point away from the binding cleft. These structures are rare examples of selective class III ligand recognition by a PDZ domain and now provide a basis for the detailed structural investigation of the promiscuous interaction between the PDZ domains of PTPN13 and their ligands. They will also lead to a better understanding of the proposed scaffolding function of these domains in multi-protein complexes assembled by PTPN13 and could ultimately contribute to low molecular weight antagonists that might even act on the PRK2 signaling pathway to modulate rearrangements of the actin cytoskeleton.

Seeking allosteric networks in PDZ domains

Ever since Ranganathan and coworkers subjected the covariation of amino acid residues in the postsynaptic density-95/Discs large/Zonula occludens 1 (PDZ) domain family to a statistical correlation analysis, PDZ domains have represented a paradigmatic family to explore single domain protein allostery. Nevertheless, several theoretical and experimental studies in the past two decades have contributed contradicting results with regard to structural localization of the allosteric networks, or even questioned their actual existence in PDZ domains. In this review, we first describe theoretical and experimental approaches that were used to probe the energetic network(s) in PDZ domains. We then compare the proposed networks for two well-studied PDZ domains namely the third PDZ domain from PSD-95 and the second PDZ domain from PTP-BL. Our analysis highlights the contradiction between the different methods and calls for additional work to better understand these allosteric phenomena.

Role of Terahertz (THz) Fluctuations in the Allosteric Properties of the PDZ Domains

With the aim of investigating the relationship between the fast fluctuations of proteins and their allosteric behavior, we perform molecular dynamics simulations of two model PDZ domains with differential allosteric responses. We focus on protein dynamics in the THz regime (0.1–3 THz) as opposed to lower frequencies. By characterizing the dynamic modulation of the protein backbone induced by ligand binding in terms of single residue and pairwise distance fluctuations, we identify a response nucleus modulated by the ligand that is visible only at THz frequencies. The residues of this nucleus undergo a significant stiffening and an increase in mutual coordination upon binding. Additionally, we find that the dynamic modulation is significantly more intense for the side chains, where it is also redistributed to distal regions not immediately in contact with the ligand allowing us to better define the response nucleus at THz frequencies. The overlap between the known allosterically responding residues of the investigated PDZ domains and the modulated region highlighted here suggests that fast THz dynamics could play a role in allosteric mechanisms.

Electric-field-stimulated protein mechanics

The internal mechanics of proteins-the coordinated motions of amino acids and the pattern of forces constraining these motions-connects protein structure to function. Here we describe a new method combining the application of strong electric field pulses to protein crystals with time-resolved X-ray crystallography to observe conformational changes in spatial and temporal detail. Using a human PDZ domain (LNX2^PDZ2) as a model system, we show that protein crystals tolerate electric field pulses strong enough to drive concerted motions on the sub-microsecond timescale. The induced motions are subtle, involve diverse physical mechanisms, and occur throughout the protein structure. The global pattern of electric-field-induced motions is consistent with both local and allosteric conformational changes naturally induced by ligand binding, including at conserved functional sites in the PDZ domain family. This work lays the foundation for comprehensive experimental study of the mechanical basis of protein function.

From Binding-Induced Dynamic Effects in SH3 Structures to Evolutionary Conserved Sectors

Src Homology 3 domains are ubiquitous small interaction modules known to act as docking sites and regulatory elements in a wide range of proteins. Prior experimental NMR work on the SH3 domain of Src showed that ligand binding induces long-range dynamic changes consistent with an induced fit mechanism. The identification of the residues that participate in this mechanism produces a chart that allows for the exploration of the regulatory role of such domains in the activity of the encompassing protein. Here we show that a computational approach focusing on the changes in side chain dynamics through ligand binding identifies equivalent long-range effects in the Src SH3 domain. Mutation of a subset of the predicted residues elicits long-range effects on the binding energetics, emphasizing the relevance of these positions in the definition of intramolecular cooperative networks of signal transduction in this domain. We find further support for this mechanism through the analysis of seven other publically available SH3 domain structures of which the sequences represent diverse SH3 classes. By comparing the eight predictions, we find that, in addition to a dynamic pathway that is relatively conserved throughout all SH3 domains, there are dynamic aspects specific to each domain and homologous subgroups. Our work shows for the first time from a structural perspective, which transduction mechanisms are common between a subset of closely related and distal SH3 domains, while at the same time highlighting the differences in signal transduction that make each family member unique. These results resolve the missing link between structural predictions of dynamic changes and the domain sectors recently identified for SH3 domains through sequence analysis.

Small protein domains as Src Homology 3 often act as docking sites and serve as regulatory elements. To understand their role in the regulation of a protein’s activity, one needs to understand how their backbone and sidechain dynamics are affected when binding to peptides. We have therefore computationally analyzed eight different SH3 domain structures, predicting dynamical effects induced by binding through our MCIT approach that has been shown to correlate well with experimental data. We show first that binding the Src SH3 domain triggers a particular cascade of dynamic effects, which are compatible with an induced fit mechanism reported before. We then combined the predictions for the eight SH3 domains into different consensus models, with the aim of analyzing, for the first time from a structural perspective, commonalities and differences in the transduction mechanisms among these SH3 domains. These consensus results are, on one hand, in agreement with the domain sectors recently identified for the entire family of SH3 domains. On the other hand, they reveal also that differences exist between the different subgroups that were studied here, requiring extensive experimental investigations of the importance of these differences for the proteins wherein these SH3 domains can be found.

Ligand binding to the PDZ domains of postsynaptic density protein 95

Cellular scaffolding and signalling is generally governed by multidomain proteins, where each domain has a particular function. Postsynaptic density protein 95 (PSD-95) is involved in synapse formation and is a typical example of such a multidomain protein. Protein-protein interactions of PSD-95 are well studied and include the following three protein ligands: (i)N-methyl-d-aspartate-type ionotropic glutamate receptor subunit GluN2B, (ii) neuronal nitric oxide synthase and (iii) cysteine-rich protein (CRIPT), all of which bind to one or more of the three PDZ domains in PSD-95. While interactions for individual PDZ domains of PSD-95 have been well studied, less is known about the influence of neighbouring domains on the function of the respective individual domain. We therefore performed a systematic study on the ligand-binding kinetics of PSD-95 using constructs of different size for PSD-95 and its ligands. Regarding the canonical peptide-binding pocket and relatively short peptides (up to 15-mer), the PDZ domains in PSD-95 by and large work as individual binding modules. However, in agreement with previous studies, residues outside of the canonical binding pocket modulate the affinity of the ligands. In particular, the dissociation of the 101 amino acid CRIPT from PSD-95 is slowed down at least 10-fold for full-length PSD-95 when compared with the individual PDZ3 domain.

Mechanisms of Differential Allosteric Modulation in Homologous Proteins: Insights from the Analysis of Internal Dynamics and Energetics of PDZ Domains

Allostery is a general phenomenon in proteins whereby a perturbation at one site reverberates into a functional change at another one, through modulation of its conformational dynamics. Herein, we address the problem of how the molecular signal encoded by a ligand is differentially transmitted through the structures of two homologous PDZ proteins: PDZ2, which responds to binding with structural and dynamical changes in regions distal from the ligand site, and PDZ3, which is characterized by less-intense dynamical variations. We use novel methods of analysis of MD simulations in the unbound and bound states to investigate the determinants of the differential allosteric behavior of the two proteins. The analysis of the correlations between the redistribution of stabilization energy and local fluctuation patterns highlights the nucleus of residues responsible for the stabilization of the 3D fold, the stability core, as the substructure that defines the difference in the allosteric response: in PDZ2, it undergoes a consistent dynamic and energetic reorganization, whereas in PDZ3, it remains largely unperturbed. Specifically, we observe for PDZ2 a significant anticorrelation between the motions of distal loops and residues of the stability core and differences in the correlation patterns between the bound and unbound states. Such variation is not observed in PDZ3, indicating that its energetics and internal dynamics are less affected by the presence/absence of the ligand. Finally, we propose a model with a direct link between the modulation of the structural, energetic and dynamic properties of a protein, and its allosteric response to a perturbation.

Deciphering the unconventional peptide binding to the PDZ domain of MAST2

Peptide binding on to microtubule-associated serine threonine kinase 2 (MAST2)—PDZ (PSD-95, Dlg1, Zo-1) prevents dimerization of the domain without directly competing with the monomer interface. Peptide binding affects positions distal from the binding groove through a network of residues undergoing subtle changes of conformation and dynamics.

Structures of the human Pals1 PDZ domain with and without ligand suggest gated access of Crb to the PDZ peptide-binding groove

Many components of epithelial polarity protein complexes possess PDZ domains that are required for protein interaction and recruitment to the apical plasma membrane. Apical localization of the Crumbs (Crb) transmembrane protein requires a PDZ-mediated interaction with Pals1 (protein-associated with Lin7, Stardust, MPP5), a member of the p55 family of membrane-associated guanylate kinases (MAGUKs). This study describes the molecular interaction between the Crb carboxy-terminal motif (ERLI), which is required for Drosophila cell polarity, and the Pals1 PDZ domain using crystallography and fluorescence polarization. Only the last four Crb residues contribute to Pals1 PDZ-domain binding affinity, with specificity contributed by conserved charged interactions. Comparison of the Crb-bound Pals1 PDZ structure with an apo Pals1 structure reveals a key Phe side chain that gates access to the PDZ peptide-binding groove. Removal of this side chain enhances the binding affinity by more than fivefold, suggesting that access of Crb to Pals1 may be regulated by intradomain contacts or by protein-protein interaction.

Disorder-to-order transition of MAGI-1 PDZ1 C-terminal extension upon peptide binding: thermodynamic and dynamic insights

PDZ domains are highly abundant protein-protein interaction modules commonly found in multidomain scaffold proteins. The PDZ1 domain of MAGI-1, a protein present at cellular tight junctions that contains six PDZ domains, is targeted by the E6 oncoprotein of the high-risk human papilloma virus. Thermodynamic and dynamic studies using complementary isothermal titration calorimetry and nuclear magnetic resonance (NMR) (15)N heteronuclear relaxation measurements were conducted at different temperatures to decipher the molecular mechanism of this interaction. Binding of E6 peptides to the MAGI-1 PDZ1 domain is accompanied by an unusually large and negative change in heat capacity (ΔC(p)) that is attributed to a disorder-to-order transition of the C-terminal extension of the PDZ1 domain upon E6 binding. Analysis of temperature-dependent thermodynamic parameters and (15)N NMR relaxation data of a PDZ1 mutant in which this disorder-to-order transition was abolished allows the unusual thermodynamic signature of E6 binding to be correlated to local folding of the PDZ1 C-terminal extension. Comparison of the exchange contributions observed for wild-type and mutant proteins explains how variation in the solvent-exposed area may compensate for the loss of conformational entropy and further designates a distinct set of a few residues that mediate this local folding phenomena.

Post-translational modifications modulate ligand recognition by the third PDZ domain of the MAGUK protein PSD-95

The relative promiscuity of hub proteins such as postsynaptic density protein-95 (PSD-95) can be achieved by alternative splicing, allosteric regulation, and post-translational modifications, the latter of which is the most efficient method of accelerating cellular responses to environmental changes in vivo. Here, a mutational approach was used to determine the impact of phosphorylation and succinimidation post-translational modifications on the binding affinity of the postsynaptic density protein-95/discs large/zonula occludens-1 (PDZ3) domain of PSD-95. Molecular dynamics simulations revealed that the binding affinity of this domain is influenced by an interplay between salt-bridges linking the α3 helix, the β2–β3 loop and the positively charged Lys residues in its high-affinity hexapeptide ligand KKETAV. The α3 helix is an extra structural element that is not present in other PDZ domains, which links PDZ3 with the following SH3 domain in the PSD-95 protein. This regulatory mechanism was confirmed experimentally via thermodynamic and NMR chemical shift perturbation analyses, discarding intra-domain long-range effects. Taken together, the results presented here reveal the molecular basis of the regulatory role of the α3 extra-element and the effects of post-translational modifications of PDZ3 on its binding affinity, both energetically and dynamically.

Kinetic response of a photoperturbed allosteric protein

By covalently linking an azobenzene photoswitch across the binding groove of a PDZ domain, a conformational transition, similar to the one occurring upon ligand binding to the unmodified domain, can be initiated on a picosecond timescale by a laser pulse. The protein structures have been characterized in the two photoswitch states through NMR spectroscopy and the transition between them through ultrafast IR spectroscopy and molecular dynamics simulations. The binding groove opens on a 100-ns timescale in a highly nonexponential manner, and the molecular dynamics simulations suggest that the process is governed by the rearrangement of the water network on the protein surface. We propose this rearrangement of the water network to be another possible mechanism of allostery.

The human PDZome: a gateway to PSD95-Disc large-zonula occludens (PDZ)-mediated functions

Protein–protein interactions organize the localization, clustering, signal transduction, and degradation of cellular proteins and are therefore implicated in numerous biological functions. These interactions are mediated by specialized domains able to bind to modified or unmodified peptides present in binding partners. Among the most broadly distributed protein interaction domains, PSD95-disc large-zonula occludens (PDZ) domains are usually able to bind carboxy-terminal sequences of their partners. In an effort to accelerate the discovery of PDZ domain interactions, we have constructed an array displaying 96% of the human PDZ domains that is amenable to rapid two-hybrid screens in yeast. We have demonstrated that this array can efficiently identify interactions using carboxy-terminal sequences of PDZ domain binders such as the E6 oncoviral protein and protein kinases (PDGFRβ, BRSK2, PCTK1, ACVR2B, and HER4); this has been validated via mass spectrometry analysis. Taking advantage of this array, we show that PDZ domains of Scrib and SNX27 bind to the carboxy-terminal region of the planar cell polarity receptor Vangl2. We also have demonstrated the requirement of Scrib for the promigratory function of Vangl2 and described the morphogenetic function of SNX27 in the early Xenopus embryo. The resource presented here is thus adapted for the screen of PDZ interactors and, furthermore, should facilitate the understanding of PDZ-mediated functions.

Calcineurin A versus NS5A-TP2/HD domain containing 2: a case study of site-directed low-frequency random mutagenesis for dissecting target specificity of peptide aptamers

We previously identified a peptide aptamer (named R5G42) via functional selection for its capacity to slow cell proliferation. A yeast two-hybrid screen of human cDNA libraries, using R5G42 as “bait,” allowed the identification of two binding proteins with very different functions: calcineurin A (CnA) (PP2B/PPP3CA), a protein phosphatase well characterized for its role in the immune response, and NS5A-TP2/HD domain containing 2, a much less studied protein induced subsequent to hepatitis C virus non-structural protein 5A expression in HepG2 hepatocellular carcinoma cells, with no known activity. Our objective in the present study was to dissect the dual target specificity of R5G42 in order to have tools with which to better characterize the actions of the peptide aptamers toward their individual targets. This was achieved through the selection of random mutants of the variable loop, derived from R5G42, evaluating their specificity toward CnA and NS5A-TP2 and analyzing their sequence. An interdisciplinary approach involving biomolecular computer simulations with integration of the sequence data and yeast two-hybrid binding phenotypes of these mutants yielded two structurally distinct conformers affording the potential molecular basis of the binding diversity of R5G42. Evaluation of the biological impact of CnA- versus NS5A-TP2-specific peptide aptamers indicated that although both contributed to the anti-proliferative effect of R5G42, CnA-binding was essential to stimulate the nuclear translocation of nuclear factor of activated T cells, indicative of the activation of endogenous CnA. By dissecting the target specificity of R5G42, we have generated novel tools with which to study each target individually. Apta-C8 is capable of directly activating CnA independent of binding to NS5A-TP2 and will be an important tool in studying the role of CnA activation in the regulation of different signaling pathways, whereas Apta-E1 will allow dissection of the function of NS5A-TP2, serving as an example of the usefulness of peptide aptamer technology for investigating signaling pathways.

Accurate prediction of the dynamical changes within the second PDZ domain of PTP1e

Experimental NMR relaxation studies have shown that peptide binding induces dynamical changes at the side-chain level throughout the second PDZ domain of PTP1e, identifying as such the collection of residues involved in long-range communication. Even though different computational approaches have identified subsets of residues that were qualitatively comparable, no quantitative analysis of the accuracy of these predictions was thus far determined. Here, we show that our information theoretical method produces quantitatively better results with respect to the experimental data than some of these earlier methods. Moreover, it provides a global network perspective on the effect experienced by the different residues involved in the process. We also show that these predictions are consistent within both the human and mouse variants of this domain. Together, these results improve the understanding of intra-protein communication and allostery in PDZ domains, underlining at the same time the necessity of producing similar data sets for further validation of thses kinds of methods.

Intra-protein communication has recently attracted an increasing interest from the scientific community, because of its important functional consequences: allostery and signalling. Unravelling how information is processed and transferred within a protein structure requires the study of the dynamical effects of, for instance, binding events, which may be captured experimentally by NMR relaxation experiments. Given the complexity of this experimental analysis, computational approaches, often based on molecular dynamics simulations, have been proposed for predicting these dynamical effects, using protein structural information as input. We examine here the accuracy of these predictors in the context of a well-studied domain, i.e. the second PSD95/Disc-large/ZO-1 domain (or PDZ domain) of PTP1e, and compare it to our approach that combines Monte-Carlo sampling of the conformational space of the side-chains and an information theoretical analysis. The results we discuss in this manuscript show clearly that the latter method provides very accurate predictions when compared to the experimental results, and has a better predictive quality compared to other computational approaches. The predictions, which are consistent between closely related structures, and the global network perspective provided by this approach, improve our understanding of intra-protein communication and allostery in these domains.

Ligand binding studied by 2D IR spectroscopy using the azidohomoalanine label

We explore the capability of the azidohomoalanine (Aha) as a vibrational label for 2D IR spectroscopy to study the binding of the target peptide to the PDZ2 domain. The Aha label responds sensitively to its local environment and its peak extinction coefficient of 350-400 M(-1) cm(-1) is high enough to routinely measure it in the low millimolar concentration regime. The central frequency, inhomogeneous width and spectral diffusion times deduced from the 2D IR line shapes of the Aha label at various positions in the peptide sequence is discussed in relationship to the known X-ray structure of the peptide bound to the PDZ2 domain. The results suggest that the Aha label introduces only a small perturbation to the overall structure of the peptide in the binding pocket. Finally, Aha is a methionine analog that can be incorporated also into larger proteins at essentially any position using protein expression. Altogether, Aha thus fulfills the requirements a versatile label should have for studies of protein structure and dynamics by 2D IR spectroscopy.

Ligand binding by PDZ domains

The postsynaptic density protein-95/disks large/zonula occludens-1 (PDZ) protein domain family is one of the most common protein-protein interaction modules in mammalian cells, with paralogs present in several hundred human proteins. PDZ domains are found in most cell types, but neuronal proteins, for example, are particularly rich in these domains. The general function of PDZ domains is to bring proteins together within the appropriate cellular compartment, thereby facilitating scaffolding, signaling, and trafficking events. The many functions of PDZ domains under normal physiological as well as pathological conditions have been reviewed recently. In this review, we focus on the molecular details of how PDZ domains bind their protein ligands and their potential as drug targets in this context.

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.

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.

Biochemical and computational analysis of LNX1 interacting proteins

PDZ (Post-synaptic density, 95 kDa, Discs large, Zona Occludens-1) domains are protein interaction domains that bind to the carboxy-terminal amino acids of binding partners, heterodimerize with other PDZ domains, and also bind phosphoinositides. PDZ domain containing proteins are frequently involved in the assembly of multi-protein complexes and clustering of transmembrane proteins. LNX1 (Ligand of Numb, protein X 1) is a RING (Really Interesting New Gene) domain-containing E3 ubiquitin ligase that also includes four PDZ domains suggesting it functions as a scaffold for a multi-protein complex. Here we use a human protein array to identify direct LNX1 PDZ domain binding partners. Screening of 8,000 human proteins with isolated PDZ domains identified 53 potential LNX1 binding partners. We combined this set with LNX1 interacting proteins identified by other methods to assemble a list of 220 LNX1 interacting proteins. Bioinformatic analysis of this protein list was used to select interactions of interest for future studies. Using this approach we identify and confirm six novel LNX1 binding partners: KCNA4, PAK6, PLEKHG5, PKC-alpha1, TYK2 and PBK, and suggest that LNX1 functions as a signalling scaffold.

Change in allosteric network affects binding affinities of PDZ domains: analysis through perturbation response scanning

The allosteric mechanism plays a key role in cellular functions of several PDZ domain proteins (PDZs) and is directly linked to pharmaceutical applications; however, it is a challenge to elaborate the nature and extent of these allosteric interactions. One solution to this problem is to explore the dynamics of PDZs, which may provide insights about how intramolecular communication occurs within a single domain. Here, we develop an advancement of perturbation response scanning (PRS) that couples elastic network models with linear response theory (LRT) to predict key residues in allosteric transitions of the two most studied PDZs (PSD-95 PDZ3 domain and hPTP1E PDZ2 domain). With PRS, we first identify the residues that give the highest mean square fluctuation response upon perturbing the binding sites. Strikingly, we observe that the residues with the highest mean square fluctuation response agree with experimentally determined residues involved in allosteric transitions. Second, we construct the allosteric pathways by linking the residues giving the same directional response upon perturbation of the binding sites. The predicted intramolecular communication pathways reveal that PSD-95 and hPTP1E have different pathways through the dynamic coupling of different residue pairs. Moreover, our analysis provides a molecular understanding of experimentally observed hidden allostery of PSD-95. We show that removing the distal third alpha helix from the binding site alters the allosteric pathway and decreases the binding affinity. Overall, these results indicate that (i) dynamics plays a key role in allosteric regulations of PDZs, (ii) the local changes in the residue interactions can lead to significant changes in the dynamics of allosteric regulations, and (iii) this might be the mechanism that each PDZ uses to tailor their binding specificities regulation.

PDZ domain proteins (PDZs) act as adapters in organizing functional protein complexes. Through dynamic interactions, PDZs play a key role in mediating key cellular functions in the cell, and they are linked to currently challenging diseases including Alzheimer's, Parkinson's and cancer. Moreover, they are associated with allosteric regulations in mediating signaling. Therefore, it is critical to have knowledge of how the allosteric transition occurs in PDZs. We investigate the allosteric response of the two most studied PDZs, PSD-95 and hPTP1E, using the perturbation response scanning (PRS) approach. The method treats the protein as an elastic network and uses linear response theory (LRT) to obtain residue fluctuations upon exerting directed random forces on selected residues. With this efficient and fast approach, we identify the key residues that mediate long-range communication and find the allosteric pathways. Although the structures of PSD-95 and hPTP1E are very similar, our analysis predicts that their allosteric pathways are different. We also observe a significant change in allosteric pathways and a decrease in binding affinity upon removal of the distal α3 helix of PSD-95. This approach enables us to understand how dynamic interactions play an important role in allosteric regulations.

Domain orientation in the N-Terminal PDZ tandem from PSD-95 is maintained in the full-length protein

Tandem PDZ domains have been suggested to form structurally independent supramodules. However, dissimilarity between crystallography and NMR models emphasize their malleable conformation. Studies in full-length scaffold proteins are needed to examine the effect of tertiary interactions within their native context. Using single-molecule fluorescence to characterize the N-terminal PDZ tandem in PSD-95, we provide the first direct evidence that PDZ tandems can be structurally independent within a full-length scaffold protein. Molecular refinement using our data converged on a single structure with an antiparallel alignment of the ligand-binding sites. Devoid of interaction partners, single-molecule conditions captured PSD-95 in its unbound, ground state. Interactions between PDZ domains could not be detected while fluctuation correlation spectroscopy showed that other conformations are dynamically sampled. We conclude that ultra-weak interactions stabilize the conformation providing a "low-relief" energy landscape that allows the domain orientation to be flipped by environmental interactions.

Crystallographic and nuclear magnetic resonance evaluation of the impact of peptide binding to the second PDZ domain of protein tyrosine phosphatase 1E

PDZ (PSD95/Discs large/ZO-1) domains are ubiquitous protein interaction motifs found in scaffolding proteins involved in signal transduction. Despite the fact that many PDZ domains show a limited tendency to undergo structural change, the PDZ family has been associated with long-range communication and allostery. One of the PDZ domains studied most in terms of structure and biophysical properties is the second PDZ ("PDZ2") domain from protein tyrosine phosphatase 1E (PTP1E, also known as PTPL1). Previously, we showed through NMR relaxation studies that binding of the RA-GEF2 C-terminal peptide substrate results in long-range propagation of side-chain dynamic changes in human PDZ2 [Fuentes, E. J., et al. (2004) J. Mol. Biol. 335, 1105-1115]. Here, we present the first X-ray crystal structures of PDZ2 in the absence and presence of RA-GEF2 ligand, determined to resolutions of 1.65 and 1.3 Å, respectively. These structures deviate somewhat from previously determined NMR structures and indicate that very minor structural changes in PDZ2 accompany peptide binding. NMR residual dipolar couplings confirm the crystal structures to be accurate models of the time-averaged atomic coordinates of PDZ2. The impact on side-chain dynamics was further tested with a C-terminal peptide from APC, which showed results nearly identical to those of RA-GEF2. Thus, allosteric transmission in PDZ2 induced by peptide binding is conveyed purely and robustly by dynamics. (15)N relaxation dispersion measurements did not detect appreciable populations of a kinetic structural intermediate. Collectively, for ligand binding to PDZ2, these data support a lock-and-key binding model from a structural perspective and an allosteric model from a dynamical perspective, which together suggest a complex energy landscape for functional transitions within the ensemble.

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.

Conserved tertiary couplings stabilize elements in the PDZ fold, leading to characteristic patterns of domain conformational flexibility

Single-domain allostery has been postulated to occur through intramolecular pathways of signaling within a protein structure. We had previously investigated these pathways by introducing a local thermal perturbation and analyzed the anisotropic propagation of structural changes throughout the protein. Here, we develop an improved approach, the Rotamerically Induced Perturbation (RIP), that identifies strong couplings between residues by analyzing the pathways of heat-flow resulting from thermal excitation of rotameric rotations at individual residues. To explore the nature of these couplings, we calculate the complete coupling maps of 5 different PDZ domains. Although the PDZ domain is a well conserved structural fold that serves as a scaffold in many protein-protein complexes, different PDZ domains display unique patterns of conformational flexibility in response to ligand binding: some show a significant shift in a set of alpha-helices, while others do not. Analysis of the coupling maps suggests a simple relationship between the computed couplings and observed conformational flexibility. In domains where the alpha-helices are rigid, we find couplings of the alpha-helices to the body of the protein, whereas in domains having ligand-responsive alpha-helices, no couplings are found. This leads to a model where the alpha-helices are intrinsically dynamic but can be damped if sidechains interact at key tertiary contacts. These tertiary contacts correlate to high covariation contacts as identified by the statistical coupling analysis method. As these dynamic modules are exploited by various allosteric mechanisms, these tertiary contacts have been conserved by evolution.

A regulated complex of the scaffolding proteins PDZK1 and EBP50 with ezrin contribute to microvillar organization

We characterize a ternary complex of PDZK1, EBP50, and ezrin that is regulated by their individual inter- and intramolecular interactions. PDZK1 is shown to undergo cell confluence-dependent nucleocytoplasmic shuttling that regulates the formation of this complex. A functional redundancy between PDZK1 and EBP50 in microvilli maintenance is shown.

PDZK1 and ezrin, radixin, moesin binding phosphoprotein 50 kDa (EBP50) are postsynaptic density 95/disc-large/zona occludens (PDZ)-domain–containing scaffolding proteins found in the apical microvilli of polarized epithelial cells. Binary interactions have been shown between the tail of PDZK1 and the PDZ domains of EBP50, as well as between EBP50 and the membrane–cytoskeletal linking protein ezrin. Here, we show that these molecules form a regulated ternary complex in vitro and in vivo. Complex formation is cooperative because ezrin positively influences the PDZK1/EBP50 interaction. Moreover, the interaction of PDZK1 with EBP50 is enhanced by the occupancy of EBP50's adjacent PDZ domain. The complex is further regulated by location, because PDZK1 shuttles from the nucleus in low confluence cells to microvilli in high confluence cells, and this regulates the formation of the PDZK1/EBP50/ezrin complex in vivo. Knockdown of EBP50 decreases the presence of microvilli, a phenotype that can be rescued by EBP50 re-expression or expression of a PDZK1 chimera that is directly targeted to ezrin. Thus, when appropriately located, PDZK1 can provide a function necessary for microvilli formation normally provided by EBP50. By entering into the ternary complex, PDZK1 can both enhance the scaffolding at the apical membrane as well as augment EBP50's role in microvilli formation.

Hidden dynamic allostery in a PDZ domain

Structure-function relationships in proteins are predicated on the spatial proximity of noncovalently interacting groups of atoms. Thus, structural elements located away from a protein's active site are typically presumed to serve a stabilizing or scaffolding role for the larger structure. Here we report a functional role for a distal structural element in a PDZ domain, even though it is not required to maintain PDZ structure. The third PDZ domain from PSD-95/SAP90 (PDZ3) has an unusual additional third alpha helix (alpha3) that packs in contiguous fashion against the globular domain. Although alpha3 lies outside the active site and does not make direct contact with C-terminal peptide ligand, removal of alpha3 reduces ligand affinity by 21-fold. Further investigation revealed that the difference in binding free energies between the full-length and truncated constructs is predominantly entropic in nature and that without alpha3, picosecond-nanosecond side-chain dynamics are enhanced throughout the domain, as determined by (2)H methyl NMR relaxation. Thus, the distal modulation of binding function appears to occur via a delocalized conformational entropy mechanism. Without removal of alpha3 and characterization of side-chain dynamics, this dynamic allostery would have gone unnoticed. Moreover, what appeared at first to be an artificial modification of PDZ3 has been corroborated by experimentally verified phosphorylation of alpha3, revealing a tangible biological mechanism for this novel regulatory scheme. This hidden dynamic allostery raises the possibility of as-yet unidentified or untapped allosteric regulation in this PDZ domain and is a very clear example of function arising from dynamics rather than from structure.

Sending signals dynamically

Proteins mediate transmission of signals along intercellular and intracellular pathways and between the exterior and the interior of a cell. The dynamic properties of signaling proteins are crucial to their functions. We discuss emerging paradigms for the role of protein dynamics in signaling. A central tenet is that proteins fluctuate among many states on evolutionarily selected energy landscapes. Upstream signals remodel this landscape, causing signaling proteins to transmit information to downstream partners. New methods provide insight into the dynamic properties of signaling proteins at the atomic scale. The next stages in the signaling hierarchy-how multiple signals are integrated and how cellular signaling pathways are organized in space and time-present exciting challenges for the future, requiring bold multidisciplinary approaches.

Channel-interacting PDZ protein, ‘CIPP’, interacts with proteins involved in cytoskeletal dynamics

Neuronal CIPP (channel-interacting PDZ protein) is a multivalent PDZ protein that interacts with specific channels and receptors highly expressed in the brain. It is composed of four PDZ domains that behave as a scaffold to clusterize functionally connected proteins. In the present study, we selected a set of potential CIPP interactors that are involved directly or indirectly in mechanisms of cytoskeletal remodelling and membrane protrusion formation. For some of these, we first proved the direct binding to specific CIPP PDZ domains considered as autonomous elements, and then confirmed the interaction with the whole protein. In particular, the small G-protein effector IRSp53 (insulin receptor tyrosine kinase substrate protein p53) specifically interacts with the second PDZ domain of CIPP and, when co-transfected in cultured mammalian cells with a tagged full-length CIPP, it induces a marked reorganization of CIPP cytoplasmic localization. Large punctate structures are generated as a consequence of CIPP binding to the IRSp53 C-terminus. Analysis of the puncta nature, using various endocytic markers, revealed that they are not related to cytoplasmic vesicles, but rather represent multi-protein assemblies, where CIPP can tether other potential interactors.