通过筛选随机多肽文库研究ZO-1中PDZ1结构域的配体结合特点

选择ZO-1的PDZ1结构域作为研究对象,以酵母双杂交为筛选系统,筛选随机多肽文库和与其它PDZ结构域的配体进行相互作用,阐明ZO-1 PDZ1的配体结合特性.ZO-1 PDZ1识别配体C末端保守的氨基酸序列通式可以表示为：［S/T］［F/Y/W］［V/I/L/C］-COOH、［S/T］［K/R］V-COOH、V［F/Y/W］［L/C］-COOH、EYV-COOH.研究发现ZO-1 PDZ1的配体同时具有3种传统PDZ结构域配体的特点,不同的是其结合配体-1位对芳香族氨基酸具有强烈的偏好性.并且某些PDZ结构域配体的-1位和-3位对结构域与配体相互作用的特异性和亲和力有重要的作用.随后通过生物信息学的方法在Swiss-Prot数据库找到与此识别规律相符合的天然人类蛋白质.根据蛋白质的功能和细胞定位等性质选择10个配体用酵母双杂交验证相互作用.证实的相互作用配体有4个.本研究希望用这样的研究策略建立一种有效的研究蛋白质相互作用的方法,通过在全蛋白质组规模上对含有结合配体保守氨基酸序列的蛋白质的查询,理论上可以找到现有数据库中所有可能与目的结构域结合的潜在配体蛋白,特别是那些筛选cDNA文库不容易获得的低丰度配体.     

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.     

Binding free energy landscape of domain-peptide interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step.     

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.     

A thermodynamic ligand binding study of the third PDZ domain (PDZ3) from the mammalian neuronal protein PSD-95

The thermodynamic parameters associated with the binding of several series of linear peptides to the third PDZ domain (PDZ3) of the postsynaptic density 95 protein (PSD-95) have been measured using isothermal titration calorimetry (ITC). Two strategies were pursued in developing these binding ligands: (1) systematic N-terminal truncation of sequences derived from the C-terminal regions of identified PDZ3-binding proteins (CRIPT, neuroligin-1, and citron) and (2) selective mutation of specific positions within a consensus hexapeptide (KKETEV) known to bind PDZ3. Each synthetically prepared peptide was used to titrate PDZ3, which yielded the changes in Gibbs free energy (DeltaG), enthalpy (DeltaH), and entropy (TDeltaS) for the binding event. Selected peptides were subjected to additional analysis, which entailed (1) measuring the change in heat capacity (DeltaCp) upon association, to assess the character of the binding interface, and (2) constructing thermodynamic double mutant cycles, to determine the presence of cooperative effects. From the first series, the CRIPT protein proved to be the better source for higher affinity sequences. From the second series, enhanced binding was associated with peptides that closely adhered to the established motif for class I PDZ domain C-termini, X-(T/S)-X-(V/I/L), and more specifically to a narrower motif of X-T-X-V. Further, in both series a length of six residues was necessary and sufficient to capture maximal affinity. In addition, there were significant influences upon binding by modifying the abutting "X" positions. The cumulative results provide greater detail into the specific nature of ligand binding to PDZ3 and will assist in the development of selective molecular probes for the study of this and structurally homologous PDZ domains.     

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. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

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.     

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.     

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.     

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.     

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

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

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.     

Role of electrostatic interactions in PDZ domain ligand recognition

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

A physical model for PDZ-domain/peptide interactions

The PDZ domain is an interaction motif that recognizes and binds the C-terminal peptides of target proteins. PDZ domains are ubiquitous in nature and help assemble multiprotein complexes that control cellular organization and signaling cascades. We present an optimized energy function to predict the binding free energy (ΔΔG) of PDZ domain/peptide interactions computationally. Geometry-optimized models of PDZ domain/peptide interfaces were built using Rosetta, and protein and peptide side chain and backbone degrees of freedom are minimized simultaneously. Using leave-one-out cross-validation, Rosetta’s energy function is adjusted to reproduce experimentally determined ΔΔG values with a correlation coefficient of 0.66 and a standard deviation of 0.79 kcal mol⁻¹. The energy function places an increased weight on hydrogen bonding interactions when compared to a previously developed method to analyze protein/protein interactions. Binding free enthalpies (ΔΔH) and entropies (ΔS) are predicted with reduced accuracies of R = 0.60 and R = 0.17, respectively. The computational method improves prediction of PDZ domain specificity from sequence and allows design of novel PDZ domain/peptide interactions.     

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.     

Origins of PDZ domain ligand specificity. Structure determination and mutagenesis of the Erbin PDZ domain

The LAP (leucine-rich repeat and PDZ-containing) family of proteins play a role in maintaining epithelial and neuronal cell size, and mutation of these proteins can have oncogenic consequences. The LAP protein Erbin has been implicated previously in a number of cellular activities by virtue of its PDZ domain-dependent association with the C termini of both ERB-B2 and the p120-catenins. The present work describes the NMR structure of Erbin PDZ in complex with a high affinity peptide ligand and includes a comprehensive energetic analysis of both the ligand and PDZ domain side chains responsible for binding. C-terminal phage display has been used to identify preferred ligands, whereas binding affinity measurements provide precise details of the energetic importance of each ligand side chain to binding. Alanine and homolog scanning mutagenesis (in a combinatorial phage display format) identifies Erbin side chains that make energetically important contacts with the ligand. The structure of a phage-optimized peptide (Ac-TGW(-4)ETW(-1)V; IC(50) = approximately 0.15 microm) in complex with Erbin PDZ provides a structural context to understand the binding energetics. In particular, the very favorable interactions with Trp(-1) are not Erbin side chain-mediated (and therefore may be generally applicable to many PDZ domains), whereas the beta2-beta3 loop provides a binding site for the Trp(-4) side chain (specific to Erbin because it has an unusually long loop). These results contribute to a growing appreciation for the importance of at least five ligand C-terminal side chains in determining PDZ domain binding energy and highlight the mechanisms of ligand discrimination among the several hundred PDZ domains present in the human genome.     

Solution structure of the RIM1alpha PDZ domain in complex with an ELKS1b C-terminal peptide

PDZ domains are widespread protein modules that commonly recognize C-terminal sequences of target proteins and help to organize macromolecular signaling complexes. These sequences usually bind in an extended conformation to relatively shallow grooves formed between a beta-strand and an alpha-helix in the corresponding PDZ domains. Because of this binding mode, many PDZ domains recognize primarily the C-terminal and the antepenultimate side-chains of the target protein, which commonly conform to motifs that have been categorized into different classes. However, an increasing number of PDZ domains have been found to exhibit unusual specificities. These include the PDZ domain of RIMs, which are large multidomain proteins that regulate neurotransmitter release and help to organize presynaptic active zones. The RIM PDZ domain binds to the C-terminal sequence of ELKS with a unique specificity that involves each of the four ELKS C-terminal residues. To elucidate the structural basis for this specificity, we have determined the 3D structure in solution of an RIM/ELKS C-terminal peptide complex using NMR spectroscopy. The structure shows that the RIM PDZ domain contains an unusually deep and narrow peptide-binding groove with an exquisite shape complementarity to the four ELKS C-terminal residues in their bound conformation. This groove is formed, in part, by a set of side-chains that is conserved selectively in RIM PDZ domains and that hence determines, at least in part, their unique specificity.     

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.     

Crystal structure of the second PDZ domain of SAP97 in complex with a GluR-A C-terminal peptide

Synaptic targeting of GluR-A subunit-containing glutamate receptors involves an interaction with synapse-associated protein 97 (SAP97). The C-terminus of GluR-A, which contains a class I PDZ ligand motif (-x-Ser/Thr-x-phi-COOH where phi is an aliphatic amino acid) associates preferentially with the second PDZ domain of SAP97 (SAP97(PDZ2)). To understand the structural basis of this interaction, we have determined the crystal structures of wild-type and a SAP97(PDZ2) variant in complex with an 18-mer C-terminal peptide (residues 890-907) of GluR-A and of two variant PDZ2 domains in unliganded state at 1.8-2.44 A resolutions. SAP97(PDZ2) folds to a compact globular domain comprising six beta-strands and two alpha-helices, a typical architecture for PDZ domains. In the structure of the peptide complex, only the last four C-terminal residues of the GluR-A are visible, and align as an antiparallel beta-strand in the binding groove of SAP97(PDZ2). The free carboxylate group and the aliphatic side chain of the C-terminal leucine (Leu907), and the hydroxyl group of Thr905 of the GluR-A peptide are engaged in essential class I PDZ interactions. Comparison between the free and complexed structures reveals conformational changes which take place upon peptide binding. The betaAlpha-betaBeta loop moves away from the C-terminal end of alphaB leading to a slight opening of the binding groove, which may better accommodate the peptide ligand. The two conformational states are stabilized by alternative hydrogen bond and coulombic interactions of Lys324 in betaAlpha-betaBeta loop with Asp396 or Thr394 in betaBeta. Results of in vitro binding and immunoprecipitation experiments using a PDZ motif-destroying L907A mutation as well as the insertion of an extra alanine residue between the C-terminal Leu907 and the stop codon are also consistent with a 'classical' type I PDZ interaction between SAP97 and GluR-A C-terminus.     

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.