Characterization of SH2D1A missense mutations identified in X-linked lymphoproliferative disease patients

X-linked lymphoproliferative disease (XLP) is a primary immunodeficiency characterized by extreme susceptibility to Epstein-Barr virus. The XLP disease gene product SH2D1A (SAP) interacts via its SH2 domain with a motif (TIYXXV) present in the cytoplasmic tail of the cell-surface receptors CD150/SLAM, CD84, CD229/Ly-9, and CD244/2B4. Characteristically, the SH2D1A three-pronged interaction with Tyr(281) of CD150 can occur in absence of phosphorylation. Here we analyze the effect of SH2D1A protein missense mutations identified in 10 XLP families. Two sets of mutants were found: (i) mutants with a marked decreased protein half-life (e.g. Y7C, S28R, Q99P, P101L, V102G, and X129R) and (ii) mutants with structural changes that differently affect the interaction with the four receptors. In the second group, mutations that disrupt the interaction between the SH2D1A hydrophobic cleft and Val +3 of its binding motif (e.g. T68I) and mutations that interfere with the SH2D1A phosphotyrosine-binding pocket (e.g. C42W) abrogated SH2D1A binding to all four receptors. Surprisingly, a mutation in SH2D1A able to interfere with Thr -2 of the CD150 binding motif (mutant T53I) severely impaired non-phosphotyrosine interactions while preserving unaffected the binding of SH2D1A to phosphorylated CD150. Mutant T53I, however, did not bind to CD229 and CD224, suggesting that SH2D1A controls several critical signaling pathways in T and natural killer cells. Because no correlation is present between identified types of mutations and XLP patient clinical presentation, additional unidentified genetic or environmental factors must play a strong role in XLP disease manifestations.     

Novel mode of ligand binding by the SH2 domain of the human XLP disease gene product SAP/SH2D1A

BACKGROUND: The Src homology 2 (SH2) domains of cytoplasmic signaling proteins generally bind phosphotyrosine (pTyr) sites in the context of carboxy-terminal residues. SAP (also known as SH2D1A or DSHP), the product of the gene that is mutated in human X-linked lymphoproliferative (XLP) disease, comprises almost exclusively a single SH2 domain, which may modulate T-cell signaling by engaging T-cell co-activators such as SLAM, thereby blocking binding of other signaling proteins that contain SH2 domains. The SAP-SLAM interaction can occur in a phosphorylation-independent manner. RESULTS: To characterize the interaction between SAP and SLAM, we synthesized peptides corresponding to the SAP-binding site at residue Y281 in SLAM. Both phosphorylated and non-phosphorylated versions of an 11-residue SLAM peptide bound SAP, with dissociation constants of 150 nM and 330 nM, respectively. SLAM phosphopeptides that were truncated either at the amino or carboxyl terminus bound with high affinity to SAP, suggesting that the SAP SH2 domain recognizes both amino-terminal and carboxy-terminal sequences relative to the pTyr residue. These results were confirmed by nuclear magnetic resonance (NMR) studies on (15)N- and (13)C-labeled SAP complexed with three SLAM peptides: an amino-terminally truncated phosphopeptide, a carboxy-terminally truncated phosphopeptide and a non-phosphorylated Tyr-containing full-length peptide. CONCLUSIONS: The SAP SH2 domain has a unique specificity. Not only does it bind peptides in a phosphorylation-independent manner, it also recognizes a pTyr residue either preceded by amino-terminal residues or followed by carboxy-terminal residues. We propose that the three 'prongs' of a peptide ligand (the amino and carboxyl termini and the pTyr) can engage the SAP SH2 domain, accounting for its unusual properties. These data point to the flexibility of modular protein-interaction domains.     

Crystal structures of the XLP protein SAP reveal a class of SH2 domains with extended, phosphotyrosine-independent sequence recognition.

SAP, the product of the gene mutated in X-linked lymphoproliferative syndrome (XLP), consists of a single SH2 domain that has been shown to bind the cytoplasmic tail of the lymphocyte coreceptor SLAM. Here we describe structures that show that SAP binds phosphorylated and nonphosphorylated SLAM peptides in a similar mode, with the tyrosine or phosphotyrosine residue inserted into the phosphotyrosine-binding pocket. We find that specific interactions with residues N-terminal to the tyrosine, in addition to more characteristic C-terminal interactions, stabilize the complexes. A phosphopeptide library screen and analysis of mutations identified in XLP patients confirm that these extended interactions are required for SAP function. Further, we show that SAP and the similar protein EAT-2 recognize the sequence motif TIpYXX(V/I).     

Binding of SAP SH2 domain to FynT SH3 domain reveals a novel mechanism of receptor signalling in immune regulation

SAP (or SH2D1A), an adaptor-like molecule expressed in immune cells, is composed almost exclusively of a Src homology 2 (SH2) domain. In humans, SAP is mutated and either absent or non-functional in X-linked lymphoproliferative (XLP) syndrome, a disease characterized by an inappropriate response to Epstein-Barr virus (EBV) infection. Through its SH2 domain, SAP associates with tyrosines in the cytoplasmic domain of the SLAM family of immune cell receptors, and is absolutely required for the function of these receptors. This property results from the ability of SAP to promote the selective recruitment and activation of FynT, a cytoplasmic Src-related protein tyrosine kinase (PTK). Here, we demonstrate that SAP operates in this pathway by binding to the SH3 domain of FynT, through a second region in the SAP SH2 domain distinct from the phosphotyrosine-binding motif. We demonstrate that this interaction is essential for SAP-mediated signalling in T cells, and for the capacity of SAP to modulate immune cell function. These observations characterize a biologically important signalling mechanism in which an adaptor molecule composed only of an SH2 domain links a receptor devoid of intrinsic catalytic activity to the kinase required for its function.     

CD150 association with either the SH2-containing inositol phosphatase or the SH2-containing protein tyrosine phosphatase is regulated by the adaptor protein SH2D1A

CD150 (SLAM/IPO-3) is a cell surface receptor that, like the B cell receptor, CD40, and CD95, can transmit positive or negative signals. CD150 can associate with the SH2-containing inositol phosphatase (SHIP), the SH2-containing protein tyrosine phosphatase (SHP-2), and the adaptor protein SH2 domain protein 1A (SH2D1A/DSHP/SAP, also called Duncan's disease SH2-protein (DSHP) or SLAM-associated protein (SAP)). Mutations in SH2D1A are found in X-linked lymphoproliferative syndrome and non-Hodgkin's lymphomas. Here we report that SH2D1A is expressed in tonsillar B cells and in some B lymphoblastoid cell lines, where CD150 coprecipitates with SH2D1A and SHIP. However, in SH2D1A-negative B cell lines, including B cell lines from X-linked lymphoproliferative syndrome patients, CD150 associates only with SHP-2. SH2D1A protein levels are up-regulated by CD40 cross-linking and down-regulated by B cell receptor ligation. Using GST-fusion proteins with single replacements of tyrosine at Y269F, Y281F, Y307F, or Y327F in the CD150 cytoplasmic tail, we found that the same phosphorylated Y281 and Y327 are essential for both SHP-2 and SHIP binding. The presence of SH2D1A facilitates binding of SHIP to CD150. Apparently, SH2D1A may function as a regulator of alternative interactions of CD150 with SHP-2 or SHIP via a novel TxYxxV/I motif (immunoreceptor tyrosine-based switch motif (ITSM)). Multiple sequence alignments revealed the presence of this TxYxxV/I motif not only in CD2 subfamily members but also in the cytoplasmic domains of the members of the SHP-2 substrate 1, sialic acid-binding Ig-like lectin, carcinoembryonic Ag, and leukocyte-inhibitory receptor families.     

Dual functional roles for the X-linked lymphoproliferative syndrome gene product SAP/SH2D1A in signaling through the signaling lymphocyte activation molecule (SLAM) family of immune receptors

The X-linked lymphoproliferative (XLP) syndrome gene encodes a protein named SAP or SH2D1A that is composed of a single Src homology 2 (SH2) domain. Two models have been proposed for its function in lymphocyte signaling. One postulates that it acts as an inhibitor of interactions between the phosphatase SHP-2 and the immune receptor SLAM. The other suggests that it functions as an adaptor to promote the recruitment of a kinase, FynT, to SLAM. Here, we provide evidence in support of both roles for SAP. Using an array of peptides derived from the SLAM family of receptors, we demonstrate that SAP binds with comparable affinities to the same sites in these receptors as do the SH2 domains of SHP-2 and SHIP, suggesting that these three proteins may compete against one another in binding to a given SLAM family receptor. Furthermore, in vitro and in vivo binding studies indicate that SAP is capable of binding directly to FynT, an interaction mediated by the FynT SH3 domain. In cells, FynT was shown to be indispensable for SLAM tyrosine phosphorylation, which, in turn, was drastically enhanced by SAP. Because SAP also blocked the recruitment of SHP-2 to SLAM in these cells, we propose a dual functional role for SAP in SLAM signaling by acting both as an adaptor for FynT and an inhibitor to SHP-2 binding. The physiological relevance of the dual functional role for SAP is underscored by the observation that disease-causing SAP mutants exhibited significantly reduced affinities to both FynT and SLAM.     

Association of the X-linked lymphoproliferative disease gene product SAP/SH2D1A with 2B4, a natural killer cell-activating molecule, is dependent on phosphoinositide 3-kinase

Natural killer (NK) cells express an activating receptor, 2B4, that enhances cellular cytotoxicity. Upon NK cell activation by ligation of 2B4, the intracellular domain of 2B4 associates with the X-linked lymphoproliferative disease (XLP) gene product, signaling lymphocytic activation molecule-associated protein/SH2D1A (SAP/SH2D1A). Defective intracellular association of 2B4 with mutated SAP/SH2D1A is likely to underlie the defects in cytotoxicity observed in NK cells from patients with XLP. We report here a role for phosphoinositide 3-kinase (PI3K) in the recruitment and association of SAP/SH2D1A to 2B4 in human NK cells. The activation of normal NK cells by ligation of 2B4 leads to the phosphorylation of 2B4, recruitment of SAP/SH2D1A, and association of the p85 regulatory subunit of PI3K. The inhibition of PI3K enzymatic activity with either wortmannin or LY294002 prior to 2B4 ligation does not alter the association of 2B4 with the p85 subunit but prevents the recruitment of SAP/SH2D1A to 2B4. In addition, PI3K inhibitors significantly diminish the cytotoxic function of primary NK cells. This observed inhibition of cytotoxicity, present in normal NK cells, was less apparent or absent in NK cells derived from a patient with XLP. These data indicate that the cytotoxicity of activated NK cells is mediated by the association of 2B4 and SAP/SH2D1A, and that this association is dependent upon the activity of PI3K.     

Prediction of binding sites of peptide recognition domains: an application on Grb2 and SAP SH2 domains

Determination of the binding motif and identification of interaction partners of the modular domains such as SH2 domains can enhance our understanding of the regulatory mechanism of protein-protein interactions. We propose here a new computational method to achieve this goal by integrating the orthogonal information obtained from binding free energy estimation and peptide sequence analysis. We performed a proof-of-concept study on the SH2 domains of SAP and Grb2 proteins. The method involves the following steps: (1) estimating the binding free energy of a set of randomly selected peptides along with a sample of known binders; (2) clustering all these peptides using sequence and energy characteristics; (3) extracting a sequence motif, which is represented by a hidden Markov model (HMM), from the cluster of peptides containing the sample of known binders; and (4) scanning the human proteome to identify binding sites of the domain. The binding motifs of the SAP and Grb2 SH2 domains derived by the method agree well with those determined through experimental studies. Using the derived binding motifs, we have predicted new possible interaction partners for the Grb2 and SAP SH2 domains as well as possible interaction sites for interaction partners already known. We also suggested novel roles for the proteins by reviewing their top interaction candidates.     

Structural basis for SH2D1A mutations in X-linked lymphoproliferative disease

X-linked lymphoproliferative disease (XLP) is a rare and severe immune deficiency, characterized by abnormal immune responses to the Epstein-Barr virus. Recently, the gene responsible for XLP, SH2D1A, has been identified and shown to code for a small cytoplasmic protein with an SH2 domain that interacts with SLAM and 2B4, two receptorial molecules involved in signal transduction in T and NK cells, respectively. A variety of SH2D1A gene mutations have been reported thus far in XLP males. Here we describe a single-strand conformation polymorphism assay for mutation analysis in XLP. Four novel patients with SH2D1A mutations are described. These mutants, and the others previously reported in the literature, have been included in a Registry (SH2D1Abase) that is fully accessible on the World Wide Web. A three-dimensional model of the SH2 domain of the SH2D1A protein has been developed, based on homology with other SH2 domains. The structural consequences of disease-causing SH2D1A mutations are discussed.     

Disease-causing SAP mutants are defective in ligand binding and protein folding

The X-linked lymphoproliferative (XLP) syndrome is caused by mutations or deletions in the SH2D1A gene that encodes an SH2 domain protein named SH2D1A or SAP. The identification of a number of missense mutations within the protein's SH2 domain, each of which can directly cause disease, provides a unique opportunity to investigate the function of an interaction protein module, SH2, in the pathogenesis of XLP. We show here that SAP mutants found in XLP patients are defective in binding its physiological ligands signaling lymphocyte activating molecule (SLAM), a co-receptor in T cell activation, and Fyn, a Src family protein tyrosine kinase. Consequently, these mutants are deficient in signaling through the SLAM receptor. This is reflected by compromised abilities for the mutants to recruit Fyn to SLAM and to activate Fyn, by reduced phosphorylation of the receptor, and by deficiencies for the mutants in blocking binding of SHP-2 to SLAM. Furthermore, all mutants examined are defective in protein folding as manifested by their significantly reduced melting temperatures upon thermal denaturation, compared to that of SAP. Taken together, these results suggest that defects in ligand binding, receptor signaling, and protein folding collectively contribute to the loss of function for disease-causing SAP mutants.     

Side-chain dynamics of the SAP SH2 domain correlate with a binding hot spot and a region with conformational plasticity

X-linked lymphoproliferative disease is caused by mutations in the protein SAP, which consists almost entirely of a single SH2 domain. SAP interacts with the Tyr281 site of the T<-->B cell signaling protein SLAM via its SH2 domain. Interestingly, binding is not dependent on phosphorylation but does involve interactions with residues N-terminal to the Tyr. We have used 15N and 2H NMR relaxation experiments to investigate the motional properties of the SAP SH2 domain backbone amides and side-chain methyl groups in the free protein and complexes with phosphorylated and non-phosphorylated peptides derived from the Tyr281 site of SLAM. The most mobile methyl groups are in side-chains with large RMSD values between the three crystal structures of SAP, suggesting that fast time-scale dynamics in side-chains is associated with conformational plasticity. The backbone amides of two residues which interact with the C-terminal part of the peptides experience fast time-scale motions in the free SH2 domain that are quenched upon binding of either the phosphorylated or non-phosphorylated peptide. Of most importance, the mobility of methyl groups in and around the binding site for residues in the N-terminus of the peptide is significantly restricted in the complexes, underscoring the dominance of this interaction with SAP and demonstrating a correlation between changes in rapid side-chain motion upon binding with local binding energy.     

Diversity in peptide recognition by the SH2 domain of SH2B1

SH2B1 is a multidomain protein that serves as a key adaptor to regulate numerous cellular events, such as insulin, leptin, and growth hormone signaling pathways. Many of these protein‐protein interactions are mediated by the SH2 domain of SH2B1, which recognizes ligands containing a phosphorylated tyrosine (pY), including peptides derived from janus kinase 2, insulin receptor, and insulin receptor substrate‐1 and ?2. Specificity for the SH2 domain of SH2B1 is conferred in these ligands either by a hydrophobic or an acidic side chain at the +3 position C‐terminal to the pY. This specificity for chemically disparate species suggests that SH2B1 relies on distinct thermodynamic or structural mechanisms to bind to peptides. Using binding and structural strategies, we have identified unique thermodynamic signatures for each peptide binding mode, and several SH2B1 residues, including K575 and R578, that play distinct roles in peptide binding. The high‐resolution structure of the SH2 domain of SH2B1 further reveals conformationally plastic protein loops that may contribute to the ability of the protein to recognize dissimilar ligands. Together, numerous hydrophobic and electrostatic interactions, in addition to backbone conformational flexibility, permit the recognition of diverse peptides by SH2B1. An understanding of this expanded peptide recognition will allow for the identification of novel physiologically relevant SH2B1/peptide interactions, which can contribute to the design of obesity and diabetes pharmaceuticals to target the ligand‐binding interface of SH2B1 with high specificity.     

The XLP syndrome protein SAP interacts with SH3 proteins to regulate T cell signaling and proliferation

The gene sap/shd1a, which encodes a 128-residue SH2 domain protein, is frequently deleted or mutated in the X-linked lymphoproliferative syndrome (XLP). The SAP SH2 domain differs from others in the same class in that it is not only capable of binding to a phosphotyrosine-containing peptide, it can also associate with an SH3 domain using a distinct surface. This novel mode of ligand-binding is initially discovered in the SLAM-SAP-Fyn complex that plays a critical role in T cell and natural killer cell activation. To identify additional binding partners for SAP, we screened a panel of 12 SH3 domains derived from regulatory proteins and identified NCK1 as a novel target of SAP in T cells. NMR analysis demonstrated that the NCK1 and Fyn SH3 domains possessed comparable affinities for SAP and engaged the same set of residues on the surface of the SAP SH2 domain. Depletion of SAP by siRNA caused a significant decrease in NCK1 tyrosine phosphorylation as well as the phosphorylation of other T cell receptor (TCR) downstream proteins such as LAT and SLP-76. Moreover, SAP was shown to regulate T cell proliferation through the MAP-kinase Erk. Taken together, our work identifies NCK1 as a novel physiological partner for SAP and a direct regulator of TCR signaling and T cell proliferation.     

Binding specificity of SH2 domains: insight from free energy simulations

Cellular signal transduction pathways are controlled by specific protein-protein interactions mediated by the binding of short peptides to small modular interaction domains. To gain insights into the specificity of these interactions, the association of phosphotyrosine-containing peptides to Src Homology 2 (SH2) domains is characterized using computations. Molecular dynamics simulations based on high-resolution crystal structures complemented by homology models are used to calculate the absolute binding free energies for 25 SH2-peptides pairs. The calculations are carried out using a potential of mean force free energy simulations method with restraining potentials that was developed previously (Woo and Roux, Proc Natl Acad Sci USA 2005;102:6825-6830). The method is utilized in conjunction with an implicit solvent representation to reduce the computational cost to characterize the association of five SH2 domains and five peptides. Specificity is ascertained by directly comparing the affinities of a given SH2 domain binding for any of the different peptides. For three of the five SH2 domains, the computational results rank the native peptides, as the most preferred binding motif. For the remaining two SH2 domains, high affinity binding motifs other than the native peptides are identified. This study illustrates how free energy computations can complement experiments in trying to elucidate complex protein-protein interactions networks.     

Crystal structure of the C-terminal SH2 domain of the p85alpha regulatory subunit of phosphoinositide 3-kinase: an SH2 domain mimicking its own substrate.

The binding properties of Src homology-2 (SH2) domains to phosphotyrosine (pY)-containing peptides have been studied in recent years with the elucidation of a large number of crystal and solution structures. Taken together, these structures suggest a general mode of binding of pY-containing peptides, explain the specificities of different SH2 domains, and may be used to design inhibitors of pY binding by SH2 domain-containing proteins. We now report the crystal structure to 1.8 A resolution of the C-terminal SH2 domain (C-SH2) of the P85alpha regulatory subunit of phosphoinositide 3-kinase (PI3 K). Surprisingly, the carboxylate group of Asp2 from a neighbouring molecule occupies the phosphotyrosine binding site and interacts with Arg18 (alphaA2) and Arg36 (betaB5), in a similar manner to the phosphotyrosine-protein interactions seen in structures of other SH2 domains complexed with pY peptides. It is the first example of a non-phosphate-containing, non-aromatic mimetic of phosphotyrosine binding to SH2 domains, and this could have implications for the design of substrate analogues and inhibitors. Overall, the crystal structure closely resembles the solution structure, but a number of loops which demonstrate mobility in solution are well defined by the crystal packing. C-SH2 has adopted a binding conformation reminiscent of the ligand bound N-terminal SH2 domain of PI3K, apparently induced by the substrate mimicking of a neighbouring molecule in the crystal.     

Potent inhibitory ligands of the GRB2 SH2 domain from recombinant peptide libraries.

We cloned and expressed the SH2 domain of human GRB2 as glutathione S-transferase and maltose binding protein fusion proteins. We screened three phagemid-based fd pVIII-protein phage display libraries against SH2 domain fusion proteins. Sequence analysis of the peptide extensions yielded a variety of related peptides. By examining the ability of the phage clones to bind other SH2 domains, we demonstrated that the phage were specific for the SH2 domain of GRB2. Based on the sequence motif identified in the "random" library screening experiment, we also built and screened a phage display library based on a Tyr-X-Asn motif (X5-Tyr-X-Asn-X8). To examine the affinity of the phage derived peptides for GRB2, we set up a radioligand competition binding assay based on immobilized GRB2 and radiolabelled autophosphorylated EGFR ICD as the radioligand. Results obtained with peptide competitors derived from the phage sequences demonstrated that nonphosphotyrosine-containing peptides identified with the phage display technology had an affinity for the receptor similar to tyrosine-phosphorylated peptides derived from the EGFR natural substrate. Interestingly, when the phage display peptides were then phosphorylated on tyrosine, their affinity for GRB2 increased dramatically. We also demonstrated the ability of the peptides to block the binding of the GRB2 SH2 domain to EGFR in a mammalian cell-based binding assay.     

Two phosphorylation-independent sites on the p85 SH2 domains bind A-Raf kinase.

Src homology 2 (SH2) domains mediate phosphotyrosine (pY)-dependent protein:protein interactions involved in signal transduction pathways. We have found that the SH2 domains of the 85-kDa alpha subunit (p85) of phosphatidylinositol 3-kinase (PI3 kinase) bind directly to the serine/threonine kinase A-Raf. In this report we show that the p85 SH2:A-Raf interaction is phosphorylation-independent. The affinity of the p85 C-SH2 domain for A-Raf and phosphopeptide pY751 was similar, raising the possibility that a p85:A-Raf complex may play a role in the coordinated regulation of the PI3 kinase and Raf-MAP kinase pathways. We further show that the p85 C-SH2 domain contains two distinct binding sites for A-Raf; one overlapping the phosphotyrosine-dependent binding site and the other a separate phosphorylation-independent site. This is the first evidence for a second binding site on an SH2 domain, distinct from the phosphotyrosine-binding pocket.     

Crystal structure of the SH3 domain of betaPIX in complex with a high affinity peptide from PAK2

The p21-activated kinases (PAKs) are important effector proteins of the small GTPases Cdc42 and Rac and control cytoskeletal rearrangements and cell proliferation. The direct interaction of PAKs with guanine nucleotide exchange factors from the PIX/Cool family, which is responsible for the localization of PAK kinases to focal complexes in the cell, is mediated by a 24-residue peptide segment in PAKs and an N-terminal src homology 3 (SH3) domain in PIX/Cool. The SH3-binding segment of PAK contains the atypical consensus-binding motif PxxxPR, which is required for unusually high affinity binding. In order to understand the structural basis for the high affinity and specificity of the PIX-PAK interaction, we solved crystal structures for the N-terminal SH3 domain of betaPIX and for the complex of the atypical binding segment of PAK2 with the N-terminal SH3 domain of betaPIX at 0.92 A and 1.3A resolution, respectively. The asymmetric unit of the crystal contains two SH3 domains and two peptide ligands. The bound peptide adopts a conformation that allows for intimate contacts with three grooves on the surface of the SH3 domain that lie between the n-Src and RT-loops. Most notably, the arginine residue of the PxxxPR motif forms a salt-bridge and is tightly coordinated by a number of residues in the SH3 domain. This arginine-specific interaction appears to be the key determinant for the high affinity binding of PAK peptides. Furthermore, C-terminal residues of the peptide engage in additional interactions with the surface of the RT-loop, which significantly increases binding specificity. Compared to a recent NMR structure of a similar complex, our crystal structure reveals an alternate binding mode. Finally, we compare our crystal structure with the recently published betaPIX/Cbl-b complex structure, and suggest the existence of a molecular switch.     

Computational analysis and prediction of the binding motif and protein interacting partners of the Abl SH3 domain

Protein-protein interactions, particularly weak and transient ones, are often mediated by peptide recognition domains, such as Src Homology 2 and 3 (SH2 and SH3) domains, which bind to specific sequence and structural motifs. It is important but challenging to determine the binding specificity of these domains accurately and to predict their physiological interacting partners. In this study, the interactions between 35 peptide ligands (15 binders and 20 non-binders) and the Abl SH3 domain were analyzed using molecular dynamics simulation and the Molecular Mechanics/Poisson-Boltzmann Solvent Area method. The calculated binding free energies correlated well with the rank order of the binding peptides and clearly distinguished binders from non-binders. Free energy component analysis revealed that the van der Waals interactions dictate the binding strength of peptides, whereas the binding specificity is determined by the electrostatic interaction and the polar contribution of desolvation. The binding motif of the Abl SH3 domain was then determined by a virtual mutagenesis method, which mutates the residue at each position of the template peptide relative to all other 19 amino acids and calculates the binding free energy difference between the template and the mutated peptides using the Molecular Mechanics/Poisson-Boltzmann Solvent Area method. A single position mutation free energy profile was thus established and used as a scoring matrix to search peptides recognized by the Abl SH3 domain in the human genome. Our approach successfully picked ten out of 13 experimentally determined binding partners of the Abl SH3 domain among the top 600 candidates from the 218,540 decapeptides with the PXXP motif in the SWISS-PROT database. We expect that this physical-principle based method can be applied to other protein domains as well.