The structural and dynamic response of MAGI-1 PDZ1 with noncanonical domain boundaries to the binding of human papillomavirus E6

PDZ domains are protein interaction domains that are found in cytoplasmic proteins involved in signaling pathways and subcellular transport. Their roles in the control of cell growth, cell polarity, and cell adhesion in response to cell contact render this family of proteins targets during the development of cancer. Targeting of these network hubs by the oncoprotein E6 of “high-risk” human papillomaviruses (HPVs) serves to effect the efficient disruption of cellular processes. Using NMR, we have solved the three-dimensional solution structure of an extended construct of the second PDZ domain of MAGI-1 (MAGI-1 PDZ1) alone and bound to a peptide derived from the C-terminus of HPV16 E6, and we have characterized the changes in backbone dynamics and hydrogen bonding that occur upon binding. The binding event induces quenching of high-frequency motions in the C-terminal tail of the PDZ domain, which contacts the peptide upstream of the canonical X-[T/S]-X-[L/V] binding motif. Mutations designed in the C-terminal flanking region of the PDZ domain resulted in a significant decrease in binding affinity for E6 peptides. This detailed analysis supports the notion of a global response of the PDZ domain to the binding event, with effects propagated to distal sites, and reveals unexpected roles for the sequences flanking the canonical PDZ domain boundaries.     

The Tiam1 PDZ Domain Couples to Syndecan1 and Promotes Cell-Matrix Adhesion

The T-cell lymphoma invasion and metastasis gene 1 (Tiam1) is a guanine exchange factor (GEF) for the Rho-family GTPase Rac1 that is crucial for the integrity of adherens junctions, tight junctions, and cell-matrix interactions. This GEF contains several protein-protein interaction domains, including a PDZ domain. Earlier studies identified a consensus PDZ-binding motif and a synthetic peptide capable of binding to the Tiam1 PDZ domain, but little is known about its ligand specificity and physiological role in cells. Here, we investigated the structure, specificity, and function of the Tiam1 PDZ domain. We determined the crystal structures of the Tiam1 PDZ domain free and in complex with a “model” peptide, which revealed the structural basis for ligand specificity. Protein database searches using the consensus PDZ-binding motif identified two eukaryotic cell adhesion proteins, Syndecan1 and Caspr4, as potential Tiam1 PDZ domain binding proteins. Equilibrium binding experiments confirmed that C-terminal peptides derived from Syndecan1 and Caspr4 bound the Tiam1 PDZ domain. NMR chemical shift perturbation experiments indicated that the Tiam1 PDZ/Syndecan1 and PDZ/Caspr4 complexes were structurally distinct and identified key residues likely to be responsible for ligand selectivity. Moreover, cell biological analysis established that Syndecan1 is a physiological binding partner of Tiam1 and that the PDZ domain has a function in cell-matrix adhesion and cell migration. Collectively, our data provide insight into the structure, specificity, and function of the Tiam1 PDZ domain. Importantly, our data report on a physiological role for the Tiam1 PDZ domain and establish a novel link between two previously unrelated signal transduction pathways, both of which are implicated in cancer.     

Design, synthesis, structure and binding properties of PDZ binding, cyclic β-finger peptides

Protein interaction domains (PIDs) play a critical role in signal transduction. One PID of great interest is the PDZ domain, a 100 amino-acid-residue domain. Most PDZ domains recognize short, C-terminal peptide motives. In the heterodimer of the nNOS-PDZ domain and the α-syntrophin-PDZ domain, however, one PDZ domain forms a β-finger that binds to the other PDZ domain. We show here that cyclic peptides derived from the β-finger of the nNOS-PDZ domain can bind the syntrophin-PDZ domain in the same manner as the whole domain. The structure of three “finger-peptides” of different size has been determined and the binding investigated using calorimetry and NMR-titration experiments.     

Energetic determinants of internal motif recognition by PDZ domains

PDZ domains are protein-protein interaction modules that organize intracellular signaling complexes. Most PDZ domains recognize specific peptide motifs followed by a required COOH-terminus. However, several PDZ domains have been found which recognize specific internal peptide motifs. The best characterized example is the syntrophin PDZ domain which, in addition to binding peptide ligands with the consensus sequence -E-S/T-X-V-COOH, also binds the neuronal nitric oxide synthase (nNOS) PDZ domain in a manner that does not depend on its precise COOH-terminal sequence. In the structure of the syntrophin-nNOS PDZ heterodimer complex, the two PDZ domains interact in a head-to-tail fashion, with an internal sequence from the nNOS PDZ domain binding precisely at the peptide binding groove of the syntrophin PDZ domain. To understand the energetic basis of this alternative mode of PDZ recognition, we have undertaken an extensive mutagenic and biophysical analysis of the nNOS PDZ domain and its interaction with the syntrophin PDZ domain. Our data indicate that the presentation of the nNOS internal motif within the context of a rigid beta-hairpin conformation is absolutely essential to binding; amino acids crucial to the structural integrity of the hairpin are as important or more important than residues that make direct contacts. The results reveal the general rules of PDZ recognition of diverse ligand types.     

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

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

Structural basis of the Na+/H+ exchanger regulatory factor PDZ1 interaction with the carboxyl-terminal region of the cystic fibrosis transmembrane conductance regulator

The PDZ1 domain of the Na(+)/H(+) exchanger regulatory factor (NHERF) binds with nanomolar affinity to the carboxyl-terminal sequence QDTRL of the cystic fibrosis transmembrane conductance regulator (CFTR) and plays a central role in the cellular localization and physiological regulation of this chloride channel. The crystal structure of human NHERF PDZ1 bound to the carboxyl-terminal peptide QDTRL has been determined at 1.7-A resolution. The structure reveals the specificity and affinity determinants of the PDZ1-CFTR interaction and provides insights into carboxyl-terminal leucine recognition by class I PDZ domains. The peptide ligand inserts into the PDZ1 binding pocket forming an additional antiparallel beta-strand to the PDZ1 beta-sheet, and an extensive network of hydrogen bonds and hydrophobic interactions stabilize the complex. Remarkably, the guanido group of arginine at position -1 of the CFTR peptide forms two salt bridges and two hydrogen bonds with PDZ1 residues Glu(43) and Asn(22), respectively, providing the structural basis for the contribution of the penultimate amino acid of the peptide ligand to the affinity of the interaction.     

Beyond the Binding Site: The Role of the β2 – β3 Loop and Extra-Domain Structures in PDZ Domains

A general paradigm to understand protein function is to look at properties of isolated well conserved domains, such as SH3 or PDZ domains. While common features of domain families are well understood, the role of subtle differences among members of these families is less clear. Here, molecular dynamics simulations indicate that the binding mechanism in PSD95-PDZ3 is critically regulated via interactions outside the canonical binding site, involving both the poorly conserved loop and an extra-domain helix. Using the CRIPT peptide as a prototypical ligand, our simulations suggest that a network of salt-bridges between the ligand and this loop is necessary for binding. These contacts interconvert between each other on a time scale of a few tens of nanoseconds, making them elusive to X-ray crystallography. The loop is stabilized by an extra-domain helix. The latter influences the global dynamics of the domain, considerably increasing binding affinity. We found that two key contacts between the helix and the domain, one involving the loop, provide an atomistic interpretation of the increased affinity. Our analysis indicates that both extra-domain segments and loosely conserved regions play critical roles in PDZ binding affinity and specificity.     

Alteration of the C-Terminal Ligand Specificity of the Erbin PDZ Domain by Allosteric Mutational Effects

Modulation of protein binding specificity is important for basic biology and for applied science. Here we explore how binding specificity is conveyed in PDZ (postsynaptic density protein-95/discs large/zonula occludens-1) domains, small interaction modules that recognize various proteins by binding to an extended C terminus. Our goal was to engineer variants of the Erbin PDZ domain with altered specificity for the most C-terminal position (position 0) where a Val is strongly preferred by the wild-type domain. We constructed a library of PDZ domains by randomizing residues in direct contact with position 0 and in a loop that is close to but does not contact position 0. We used phage display to select for PDZ variants that bind to 19 peptide ligands differing only at position 0. To verify that each obtained PDZ domain exhibited the correct binding specificity, we selected peptide ligands for each domain. Despite intensive efforts, we were only able to evolve Erbin PDZ domain variants with selectivity for the aliphatic C-terminal side chains Val, Ile and Leu. Interestingly, many PDZ domains with these three distinct specificities contained identical amino acids at positions that directly contact position 0 but differed in the loop that does not contact position 0. Computational modeling of the selected PDZ domains shows how slight conformational changes in the loop region propagate to the binding site and result in different binding specificities. Our results demonstrate that second-sphere residues could be crucial in determining protein binding specificity.     

Beyond the binding site: the role of the β₂-β₃ loop and extra-domain structures in PDZ domains

A general paradigm to understand protein function is to look at properties of isolated well conserved domains, such as SH3 or PDZ domains. While common features of domain families are well understood, the role of subtle differences among members of these families is less clear. Here, molecular dynamics simulations indicate that the binding mechanism in PSD95-PDZ3 is critically regulated via interactions outside the canonical binding site, involving both the poorly conserved β?-β? loop and an extra-domain helix. Using the CRIPT peptide as a prototypical ligand, our simulations suggest that a network of salt-bridges between the ligand and this loop is necessary for binding. These contacts interconvert between each other on a time scale of a few tens of nanoseconds, making them elusive to X-ray crystallography. The loop is stabilized by an extra-domain helix. The latter influences the global dynamics of the domain, considerably increasing binding affinity. We found that two key contacts between the helix and the domain, one involving the β?-β? loop, provide an atomistic interpretation of the increased affinity. Our analysis indicates that both extra-domain segments and loosely conserved regions play critical roles in PDZ binding affinity and specificity.     

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.     

The Erbin PDZ domain binds with high affinity and specificity to the carboxyl termini of delta-catenin and ARVCF

Erbin is a recently described member of the LAP (leucine-rich repeat and PDZ domain) protein family. We used a C-terminally displayed phage peptide library to identify optimal ligands for the Erbin PDZ domain. Phage-selected peptides were type 1 PDZ ligands that bound with high affinity and specificity to the Erbin PDZ domain in vitro. These peptides most closely resembled the C-terminal PDZ domain-binding motifs of three p120-related catenins: delta-catenin, ARVCF, and p0071 (DSWV-COOH). Analysis of the interactions of the Erbin PDZ domain with synthetic peptides matching the C termini of ARVCF or delta-catenin also demonstrated specific high affinity binding. We characterized the interactions between the Erbin PDZ domain and both ARVCF and delta-catenin in vitro and in vivo. The Erbin PDZ domain co-localized and coprecipitated with ARVCF or delta-catenin complexed with beta-catenin and E/N-cadherin. Mutagenesis and peptide competition experiments showed that the association of Erbin with the cadherin-catenin complex was mediated by the interaction of its PDZ domain with the C-terminal PDZ domain-binding motifs (DSWV-COOH) of ARVCF and delta-catenin. Finally, we showed that endogenous delta-catenin and Erbin co-localized in and co-immunoprecipitated from neurons. These results suggest that delta-catenin and ARVCF may function to mediate the association of Erbin with the junctional cadherin-catenin complex. They also demonstrate that C-terminal phage-display technology can be used to predict physiologically relevant ligands for PDZ domains.     

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

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

A Structural Portrait of the PDZ Domain Family

PDZ (PSD-95/Discs-large/ZO1) domains are interaction modules that typically bind to specific C-terminal sequences of partner proteins and assemble signaling complexes in multicellular organisms. We have analyzed the existing database of PDZ domain structures in the context of a specificity tree based on binding specificities defined by peptide-phage binding selections. We have identified 16 structures of PDZ domains in complex with high-affinity ligands and have elucidated four additional structures to assemble a structural database that covers most of the branches of the PDZ specificity tree. A detailed comparison of the structures reveals features that are responsible for the diverse specificities across the PDZ domain family. Specificity differences can be explained by differences in PDZ residues that are in contact with the peptide ligands, but these contacts involve both side-chain and main-chain interactions. Most PDZ domains bind peptides in a canonical conformation in which the ligand main chain adopts an extended β-strand conformation by interacting in an antiparallel fashion with a PDZ β-strand. However, a subset of PDZ domains bind peptides with a bent main-chain conformation and the specificities of these non-canonical domains could not be explained based on canonical structures. Our analysis provides a structural portrait of the PDZ domain family, which serves as a guide in understanding the structural basis for the diverse specificities across the family.     

Comprehensive analysis of all evolutionary paths between two divergent PDZ domain specificities

To understand the molecular evolution of functional diversity in protein families, we comprehensively investigated the consequences of all possible mutation combinations separating two peptide‐binding domains with highly divergent specificities. We analyzed the Erbin PDZ domain (Erbin‐PDZ), which exhibits canonical type I specificity, and a synthetic Erbin‐PDZ variant (E‐14) that differs at six positions and exhibits an atypical specificity that closely resembles that of the natural Pdlim4 PDZ domain (Pdlim4‐PDZ). We constructed a panel of 64 PDZ domains covering all possible transitions between Erbin‐PDZ and E‐14 (i.e., the panel contained variants with all possible combinations of either the Erbin‐PDZ or E‐14 sequence at the six differing positions). We assessed the specificity profiles of the 64 PDZ domains using a C‐terminal phage‐displayed peptide library containing all possible genetically encoded heptapeptides. The specificity profiles clustered into six distinct groups, showing that intermediate domains can be nodes for the evolution of divergent functions. Remarkably, three substitutions were sufficient to convert the specificity of Erbin‐PDZ to that of Pdlim4‐PDZ, whereas Pdlim4‐PDZ contains 71 differences relative to Erbin‐PDZ. X‐ray crystallography revealed the structural basis for specificity transition: a single substitution in the center of the binding site, supported by contributions from auxiliary substitutions, altered the main chain conformation of the peptide ligand to resemble that of ligands bound to Pdlim4‐PDZ. Our results show that a very small set of mutations can dramatically alter protein specificity, and these findings support the hypothesis whereby complex protein functions evolve by gene duplication followed by cumulative mutations.     

Solution structure of the PDZ2 domain from cytosolic human phosphatase hPTP1E complexed with a peptide reveals contribution of the beta2-beta3 loop to PDZ domain-ligand interactions

The solution structure of the second PDZ domain from human phosphatase hPTP1E in complex with a C-terminal peptide from the guanine nucleotide exchange factor RA-GEF-2 has been determined using 2D and 3D heteronuclear NMR experiments. Compared to previously solved structures, the hPTP1E complex shows an enlarged interaction surface with the C terminus of the bound peptide. Novel contacts were found between the long structured beta2/beta3 loop of the PDZ domain and the sixth amino acid residue from the C terminus of the peptide. This work underlines the importance of the beta2/beta3 loop for ligand selection by PDZ domains.     

Peptide selectivity between the PDZ domains of human pregnancy‐related serine proteases (HtrA1, HtrA2, HtrA3, and HtrA4) can be reshaped by different halogen probes

The human HtrA family of serine proteases (HtrA1, HtrA2, HtrA3, and HtrA4) are the key enzymes associated with pregnancy and closely related to the development and progression of many pathological events. Previously, it was found that halogen substitution at the indole moiety of peptide Trp‐1 residue can form a geometrically satisfactory halogen bond with the Drosophila discs large, zona occludens‐1 (PDZ) domain of HtrA proteases. Here, we attempt to systematically investigate the effect of substitution with 4 halogen types and 2 indole positions on the binding affinity and specificity of peptide ligands to the 4 HtrA PDZ domains. The complex structures, interaction energies, halogen‐bonding strength, and binding affinity of domain‐peptide systems were modeled, analyzed, and measured via computational modeling and fluorescence‐based assay. It is revealed that there is a compromise between the local rearrangement of halogen bond involving different halogen atoms and the global optimization of domain‐peptide interaction; the substitution position is fundamentally important for peptide‐binding affinity, while the halogen type can effectively shift peptide selectivity between the 4 domains. The HtrA1‐PDZ and HtrA4‐PDZ as well as HtrA2‐PDZ and HtrA3‐PDZ respond similarly to different halogen substitutions of peptide; –Br substitution at R2‐position and –I substitution at R4‐position are most effective in improving peptide selectivity for HtrA1‐PDZ/HtrA4‐PDZ and HtrA2‐PDZ/HtrA3‐PDZ, respectively; –F substitution would not address substantial effect on peptide selectivity for all the 4 domains. Consequently, the binding affinities of a native peptide ligand DSRIWWV^–COOH as well as its 4 R2‐halogenated counterparts were determined as 1.9, 1.4, 0.5, 0.27, and 0.92 μM, which are basically consistent with computational analysis. This study would help to rationally design selective peptide inhibitors of HtrA family members by using different halogen substitutions.     

Quantification of PDZ domain specificity, prediction of ligand affinity and rational design of super-binding peptides

Transient macromolecular complexes are often formed by protein-protein interaction domains (e.g. PDZ, SH2, SH3, WW) which recognize linear sequence motifs with in vitro affinities typically in the micromolar range. The analysis of the resulting interaction networks requires a quantification of domain specificity and selectivity towards all possible ligands with physiologically relevant affinity. As representative examples, we determined specificity as a function of ligand sequence-dependent affinity contributions by statistical analysis of peptide library screens for the AF6, ERBIN and SNA1 (alpha-1-syntrophin) PDZ domains. For this purpose, the three PDZ domains were first screened for binding with a peptide library comprising 6223 human C termini created by SPOT synthesis. Based on the detected ligand preferences, we designed focused peptide libraries (profile libraries). These libraries were used to quantify the affinity contributions of the four C-terminal ligand residues by means of ANOVA models (analysis of variance) relating the C-terminal ligand sequences to the corresponding dissociation constants. Our models agreed well with experimentally determined dissociation constants and allowed us to design super binding peptides. The latter were shown experimentally to bind to their cognate PDZ domains with the highest affinity. In addition, we determined structure-activity relationships and thereby rationalized the position-specific affinity contributions. Furthermore, we used the statistical models to predict the dissociation constants for the complete ligand sequence space and thus determined the specificity overlap for the three investigated PDZ domains (). Altogether, we present an efficient method for profiling protein-protein interaction domains that provides a biophysical picture of specificity and selectivity. This approach allows the rational design of functional experiments and provides a basis for simulating interaction networks in the field of systems biology.     

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

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

Structure-Based Prediction of the Peptide Sequence Space Recognized by Natural and Synthetic PDZ Domains

Protein-protein recognition, frequently mediated by members of large families of interaction domains, is one of the cornerstones of biological function. Here, we present a computational, structure-based method to predict the sequence space of peptides recognized by PDZ domains, one of the largest families of recognition proteins. As a test set, we use a considerable amount of recent phage display data that describe the peptide recognition preferences for 169 naturally occurring and engineered PDZ domains. For both wild-type PDZ domains and single point mutants, we find that 70-80% of the most frequently observed amino acids by phage display are predicted within the top five ranked amino acids. Phage display frequently identified recognition preferences for amino acids different from those present in the original crystal structure. Notably, in about half of these cases, our algorithm correctly captures these preferences, indicating that it can predict mutations that increase binding affinity relative to the starting structure. We also find that we can computationally recapitulate specificity changes upon mutation, a key test for successful forward design of protein-protein interface specificity. Across all evaluated data sets, we find that incorporation backbone sampling improves accuracy substantially, irrespective of using a crystal or NMR structure as the starting conformation. Finally, we report successful prediction of several amino acid specificity changes from blind tests in the DREAM4 peptide recognition domain specificity prediction challenge. Because the foundational methods developed here are structure based, these results suggest that the approach can be more generally applied to specificity prediction and redesign of other protein-protein interfaces that have structural information but lack phage display data.