Sequence-specific long range networks in PSD-95/discs large/ZO-1 (PDZ) domains tune their binding selectivity

Protein-protein interactions mediated by modular protein domains are critical for cell scaffolding, differentiation, signaling, and ultimately, evolution. Given the vast number of ligands competing for binding to a limited number of domain families, it is often puzzling how specificity can be achieved. Selectivity may be modulated by intradomain allostery, whereby a remote residue is energetically connected to the functional binding site via side chain or backbone interactions. Whereas several energetic pathways, which could mediate intradomain allostery, have been predicted in modular protein domains, there is a paucity of experimental data to validate their existence and roles. Here, we have identified such functional energetic networks in one of the most common protein-protein interaction modules, the PDZ domain. We used double mutant cycles involving site-directed mutagenesis of both the PDZ domain and the peptide ligand, in conjunction with kinetics to capture the fine energetic details of the networks involved in peptide recognition. We performed the analysis on two homologous PDZ-ligand complexes and found that the energetically coupled residues differ for these two complexes. This result demonstrates that amino acid sequence rather than topology dictates the allosteric pathways. Furthermore, our data support a mechanism whereby the whole domain and not only the binding pocket is optimized for a specific ligand. Such cross-talk between binding sites and remote residues may be used to fine tune target selectivity.     

A conformational switch in the CRIB-PDZ module of Par-6

Here, we report a novel mechanism of PDZ (PSD-95/Dlg/ZO-1) domain regulation that distorts?a conserved element of PDZ ligand recognition. The polarity regulator Par-6 assembles a conserved multiprotein complex and is directly modulated by the?Rho GTPase Cdc42. Cdc42 binds the adjacent Cdc42/Rac interactive binding (CRIB) and PDZ domains of Par-6, increasing C-terminal ligand binding affinity by 10-fold. By solving structures of the isolated PDZ domain and a disulfide-stabilized CRIB-PDZ, we detected a conformational switch that controls affinity by altering the configuration of the conserved "GLGF" loop. As a result, lysine 165 is displaced from the PDZ core by an adjacent hydrophobic residue, disrupting coordination of the PDZ ligand-binding cleft. Stabilization of the CRIB:PDZ interface restores K165 to its canonical location in the binding pocket. We conclude that a unique "dipeptide switch" in the Par-6 PDZ transmits a signal for allosteric activation to the ligand-binding pocket.     

Solution structure of AF-6 PDZ domain and its interaction with the C-terminal peptides from Neurexin and Bcr

AF-6 is a key molecule essential for structure organization of cell-cell junction of polarized epithelia. It belongs to a novel cell-cell adhesion system. The AF-6 PDZ domain mediates interactions by binding to a specific amino acid sequence in target proteins. Here we report the solution structure of the AF-6 PDZ domain determined by NMR. Previously, the AF-6 PDZ domain was considered to be a class II PDZ domain. However we found that a unique hydrophilic amino acid, Gln70, at position alphaB1 makes the alphaB/betaB groove of the AF-6 PDZ domain significantly different from that of the canonical class II PDZ domain. The AF-6 PDZ domain does not have the second hydrophobic binding pocket, and the N-terminal end of alphaB is closer to betaB. Using BIACORE and NMR chemical shift perturbation experiments, we have studied the binding characteristics of the PDZ domain to the C-terminal peptide of Neurexin, KKNKDKEYYV, and that of Bcr, KRQSILFSTEV. The C-terminal peptide of Neurexin is a class II ligand, whereas that of Bcr is a class I ligand. The dissociation constants of these ligands were 4.08 x 10(-7) and 2.23 x 10(-6) m, respectively. Each of the four C-terminal positions in Neurexin and Bcr may contribute to the interaction. The three-dimensional models of the AF-6 PDZ-Neurexin C-terminal peptide complex and the AF-6 PDZ-Bcr C-terminal peptide complex were built up by molecular dynamics simulations. Unlike the canonical class II PDZ domain, Ala74 at alphaB5 rather than the residue at alphaB1 makes direct hydrophobic contact with the side chain of Tyr at the -2 position of the ligand.     

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.     

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.     

An In Silico Analysis of the Binding Modes and Binding Affinities of Small Molecule Modulators of PDZ-Peptide Interactions

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known K_i or K_d values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high K_i or K_d moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.     

Solution structure of the extended neuronal nitric oxide synthase PDZ domain complexed with an associated peptide.

The PDZ domain of neuronal nitric oxide synthase (nNOS) functions as a scaffold for organizing the signal transduction complex of the enzyme. The NMR structure of a complex composed of the nNOS PDZ domain and an associated peptide suggests that a two-stranded beta-sheet C-terminal to the canonical PDZ domain may mediate its interaction with the PDZ domains of postsynaptic density-95 and alpha-syntrophin. The structure also provides the molecular basis of recognition of Asp-X-Val-COOH peptides by the nNOS PDZ domain. The role of the C-terminal extension in Asp-X-Val-COOH peptide binding is investigated. Additionally, NMR studies further show that the Asp-X-Val-COOH peptide and a C-terminal peptide from a novel cytosolic protein named CAPON bind to the same pocket of the nNOS PDZ domain.     

Investigation of the PDZ domain ligand binding site using chemically modified peptides

Several chemically modified analogues to a tightly binding ligand for the second PDZ domain of MAGI-3 were synthesized and evaluated for their ability to compete with native peptide ligands. N-methyl scanning of the ligand backbone amides revealed the energetically important hydrogen bonds between the ligand backbone and the PDZ domain. Analogues to the ligand's conserved threonine/serine(-2) residue, involved in a side chain to side chain hydrogen bond with a conserved histidine in the PDZ domain, revealed that the interaction is highly sensitive to the steric structure around the hydroxyl group of this residue. Analogues of the ligand carboxy terminus revealed that the full hydrogen bond network of the GLGF loop is important in ligand binding.     

Rapid characterization of the binding property of HtrA2/Omi PDZ domain by validation screening of PDZ ligand library

HtrA2/Omi is a mammalian mitochondrial serine protease, and was found to have dual roles in mammalian cells, not only acting as an apoptosis-inducing protein but also maintaining mitochondrial homeostasis. PDZ domain is one of the most important protein-protein interaction modules and is involved in a variety of important cellular functions, such as signal transduction, degradation of proteins, and formation of cytoskeleton. Recently, it was reported that the PDZ domain of HtrA2/Omi might regulate proteolytic activity through its interactions with ligand proteins. In this study, we rapidly characterized the binding properties of HtrA2/Omi PDZ domain by validation screening of the PDZ ligand library with yeast two-hybrid approach. Then, we predicted its novel ligand proteins in human proteome and reconfirmed them in the yeast two-hybrid system. Finally, we analyzed the smallest networks bordered by the shortest path length between the protein pairs of novel interactions to evaluate the confidence of the identified interactions. The results revealed some novel binding properties of HtrA2/Omi PDZ domain. Besides the reported Class II PDZ motif, it also binds to Class I and Class III motifs, and exhibits restricted variability at P⁻³, which means that the P⁻³ residue is selected according to the composition of the last three residues. Seven novel ligand proteins of HtrA2/Omi PDZ domain were discovered, providing significant clues for further clarifying the roles of HtrA2/Omi. Moreover, this study proves the high efficiency and practicability of the newly developed validation screening of candidate ligand library method for binding property characterization of peptide-binding domains.     

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.     

基于弹性网络模型的蛋白质变构路径与关键残基识别研究

目的 变构效应在蛋白质生物学功能执行过程中发挥着重要的调控作用，如何基于蛋白质空间结构，有效识别变构信号的传播路径和关键的残基位点是蛋白质结构-功能关系研究领域的热点科学问题。方法 本研究利用基于弹性网络模型（elastic network model，ENM）的力分布计算方法，通过分析蛋白质对外力的响应过程，来识别体系的变构路径以及变构过程中的关键残基。在该方法中，对蛋白质的关键变构位点施加外力，通过对体系形变以及内力分布情况的分析，有效识别与外力承载区域形变相耦合的关键残基，从而得到力信号在蛋白质结构内的传播路径。结果 利用该方法研究了人类磷酸甘油酸激酶（human phosphoglycerate kinase，hPGK）和蛋白质酪氨酸磷酸酶（protein tyrosine phosphatase，PTP）PDZ2结构域的变构调控路径和关键残基。对于hPGK，识别出从底物结合位点到铰链区的两条变构信号传导路径。对于PTP PDZ2，也成功识别出从配体结合位点传递到蛋白质远端的两条长程变构调控路径。计算结果与实验和分子动力学（molecular dynamics，MD）模拟得到的结果一致。结论 本研究为蛋白质体系关键残基识别及变构路径研究提供了有效的分析方法。     

Molecular determinants for the complex binding specificity of the PDZ domain in PICK1

PICK1 (protein interacting with C kinase 1) contains a single PDZ domain known to mediate interaction with the C termini of several receptors, transporters, ion channels, and kinases. In contrast to most PDZ domains, the PICK1 PDZ domain interacts with binding sequences classifiable as type I (terminating in (S/T)XPhi; X, any residue) as well as type II (PhiXPhi; Phi, any hydrophobic residue). To enable direct assessment of the affinity of the PICK1 PDZ domain for its binding partners we developed a purification scheme for PICK1 and a novel quantitative binding assay based on fluorescence polarization. Our results showed that the PICK1 PDZ domain binds the type II sequence presented by the human dopamine transporter (-WLKV) with an almost 15-fold and >100-fold higher affinity than the type I sequences presented by protein kinase Calpha (-QSAV) and the beta(2)-adrenergic receptor (-DSLL), respectively. Mutational analysis of Lys(83) in the alphaB1 position of the PDZ domain suggested that this residue mimics the function of hydrophobic residues present in this position in regular type II PDZ domains. The PICK1 PDZ domain was moreover found to prefer small hydrophobic residues in the C-terminal P(0) position of the ligand. Molecular modeling predicted a rank order of (Val > Ile > Leu) that was verified experimentally with up to a approximately 16-fold difference in binding affinity between a valine and a leucine in P(0). The results define the structural basis for the unusual binding pattern of the PICK1 PDZ domain by substantiating the critical role of the alphaB1 position (Lys(83)) and of discrete side chain differences in position P(0) of the ligands.     

Structural basis for NHERF1 PDZ domain binding

The Na(+)/H(+) exchange regulatory factor-1 (NHERF1) is a scaffolding protein that possesses two tandem PDZ domains and a carboxy-terminal ezrin-binding domain (EBD). The parathyroid hormone receptor (PTHR), type II sodium-dependent phosphate cotransporter (Npt2a), and β2-adrenergic receptor (β2-AR), through their respective carboxy-terminal PDZ-recognition motifs, individually interact with NHERF1 forming a complex with one of the PDZ domains. In the basal state, NHERF1 adopts a self-inhibited conformation, in which its carboxy-terminal PDZ ligand interacts with PDZ2. We applied molecular dynamics (MD) simulations to uncover the structural and biochemical basis for the binding selectivity of NHERF1 PDZ domains. PDZ1 uniquely forms several contacts not present in PDZ2 that further stabilize PDZ1 interactions with target ligands. The binding free energy (ΔG) of PDZ1 and PDZ2 with the carboxy-terminal, five-amino acid residues that form the PDZ-recognition motif of PTHR, Npt2a, and β2-AR was calculated and compared with the calculated ΔG for the self-association of NHERF1. The results suggest that the interaction of the PTHR, β2-adrenergic, and Npt2a involves competition between NHERF1 PDZ domains and the target proteins. The binding of PDZ2 with PTHR may also compete with the self-inhibited conformation of NHERF1, thereby contributing to the stabilization of an active NHERF1 conformation.     

The kinetics of PDZ domain-ligand interactions and implications for the binding mechanism

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

Unusual binding interactions in PDZ domain crystal structures help explain binding mechanisms

PDZ domains most commonly bind the C‐terminus of their protein targets. Typically the C‐terminal four residues of the protein target are considered as the binding motif, particularly the C‐terminal residue (P0) and third‐last residue (P‐2) that form the major contacts with the PDZ domain's “binding groove”. We solved crystal structures of seven human PDZ domains, including five of the seven PDLIM family members. The structures of GRASP, PDLIM2, PDLIM5, and PDLIM7 show a binding mode with only the C‐terminal P0 residue bound in the binding groove. Importantly, in some cases, the P‐2 residue formed interactions outside of the binding groove, providing insight into the influence of residues remote from the binding groove on selectivity. In the GRASP structure, we observed both canonical and noncanonical binding in the two molecules present in the asymmetric unit making a direct comparison of these binding modes possible. In addition, structures of the PDZ domains from PDLIM1 and PDLIM4 also presented here allow comparison with canonical binding for the PDLIM PDZ domain family. Although influenced by crystal packing arrangements, the structures nevertheless show that changes in the positions of PDZ domain side‐chains and the αB helix allow noncanonical binding interactions. These interactions may be indicative of intermediate states between unbound and fully bound PDZ domain and target protein. The noncanonical “perpendicular” binding observed potentially represents the general form of a kinetic intermediate. Comparison with canonical binding suggests that the rearrangement during binding involves both the PDZ domain and its ligand.     

Convergent and divergent ligand specificity among PDZ domains of the LAP and zonula occludens (ZO) families

We present a detailed comparative analysis of the PDZ domains of the human LAP proteins Erbin, Densin-180, and Scribble and the MAGUK ZO-1. Phage-displayed peptide libraries and in vitro affinity assays were used to define ligand binding profiles for each domain. The analysis reveals the importance of interactions with all four C-terminal residues of the ligand, which constitute a core recognition motif, and also the role of interactions with more upstream ligand residues that support and modulate the core binding interaction. In particular, the results highlight the importance of site(-1), which interacts with the penultimate residue of ligand C termini. Site(-1) was found to be monospecific in the Erbin PDZ domain (accepts tryptophan only), bispecific in the first PDZ domain of ZO-1 (accepts tryptophan or tyrosine), and promiscuous in the Scribble PDZ domains. Furthermore, it appears that the level of promiscuity within site(-1) greatly influences the range of potential biological partners and functions that can be associated with each protein. These findings show that subtle changes in binding specificity can significantly alter the range of biological partners for PDZ domains, and the insights enhance our understanding of this diverse family of peptide-binding modules.     

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.     

Rapid characterization of the binding property of HtrA2/Omi PDZ domain by validation screening of PDZ ligand library

HtrA2/Omi is a mammalian mitochondrial serine protease, and was found to have dual roles in mammalian cells, not only acting as an apoptosis-inducing protein but also maintaining mitochondrial homeostasis. PDZ domain is one of the most important protein-protein interaction modules and is involved in a variety of important cellular functions, such as signal transduction, degradation of proteins,and formation of cytoskeleton. Recently, it was reported that the PDZ domain of HtrA2/Omi might regulate proteolytic activity through its interactions with ligand proteins. In this study, we rapidly characterized the binding properties of HtrA2/Omi PDZ domain by validation screening of the PDZ ligand library with yeast two-hybrid approach. Then, we predicted its novel ligand proteins in human proteome and reconfirmed them in the yeast two-hybrid system. Finally, we analyzed the smallest networks bordered by the shortest path length between the protein pairs of novel interactions to evaluate the confidence of the identified interactions. The results revealed some novel binding properties of HtrA2/Omi PDZ domain. Besides the reported Class Ⅱ PDZ motif, it also binds to Class Ⅰ and Class Ⅲ motifs, and exhibits restricted variability at P-3, which means that the P-3 residue is selected according to the composition of the last three residues. Seven novel ligand proteins of HtrA2/Omi PDZ domain were discovered, providing significant clues for further clarifying the roles of HtrA2/Omi.Moreover, this study proves the high efficiency and practicability of the newly developed validation screening of candidate ligand library method for binding property characterization of peptide-binding domains.     

Solution structure and backbone dynamics of the AF-6 PDZ domain/Bcr peptide complex

The human AF-6, a scaffold protein between cell membrane-associated proteins and the actin cytoskeleton, plays an important role in special cell-cell junctions and signal transduction. It can be phosphorylated by the protein kinase Bcr, which allows efficient binding of the C terminus of Bcr to the PDZ domain of AF-6 and consequently enhances the binding affinity of AF-6 to Ras. Formation of the AF-6, Bcr, and Ras ternary complex results in down-regulation of the Ras-mediated signal transduction pathway. To better understand the molecular basis for the recognition of the AF-6 PDZ domain and Bcr, we solve the solution structure of the AF-6 PDZ domain complexed with the C-terminal peptide of Bcr and explore the interactions between them in detail. Compared with previously reported structures, the complex exhibits a noncanonical binding mode of PDZ/peptide. Owing to the distinct residues involved in the AF-6 PDZ domain and Bcr peptide interaction, the interaction mode does not adapt to the existing classification rules that have been put forward, based on the ligand or the PDZ domain specificity. Furthermore, the PDZ domain of AF-6 can bind to the C terminus of Bcr efficiently after phosphorylation of AF-6 by the Bcr kinase. The phosphorylation may induce a conformational change of AF-6, which makes the binding surface on the PDZ domain accessible to Bcr for efficient binding. This study not only characterizes the structural details of the AF-6 PDZ/Bcr peptide complex, but also provides a potential target for future drug design and disease therapy.