Structural basis of affinity maturation and intramolecular cooperativity in a protein-protein interaction

Although protein-protein interactions are involved in nearly all cellular processes, general rules for describing affinity and selectivity in protein-protein complexes are lacking, primarily because correlations between changes in protein structure and binding energetics have not been well determined. Here, we establish the structural basis of affinity maturation for a protein-protein interaction system that we had previously characterized energetically. This model system exhibits a 1500-fold affinity increase. Also, its affinity maturation is restricted by negative intramolecular cooperativity. With three complex and six unliganded variant X-ray crystal structures, we provide molecular snapshots of protein interface remodeling events that span the breadth of the affinity maturation process and present a comprehensive structural view of affinity maturation. Correlating crystallographically observed structural changes with measured energetic changes reveals molecular bases for affinity maturation, intramolecular cooperativity, and context-dependent binding.     

Structure-based rational design of a Toll-like receptor 4 (TLR4) decoy receptor with high binding affinity for a target protein

Repeat proteins are increasingly attracting much attention as alternative scaffolds to immunoglobulin antibodies due to their unique structural features. Nonetheless, engineering interaction interface and understanding molecular basis for affinity maturation of repeat proteins still remain a challenge. Here, we present a structure-based rational design of a repeat protein with high binding affinity for a target protein. As a model repeat protein, a Toll-like receptor4 (TLR4) decoy receptor composed of leucine-rich repeat (LRR) modules was used, and its interaction interface was rationally engineered to increase the binding affinity for myeloid differentiation protein 2 (MD2). Based on the complex crystal structure of the decoy receptor with MD2, we first designed single amino acid substitutions in the decoy receptor, and obtained three variants showing a binding affinity (K(D)) one-order of magnitude higher than the wild-type decoy receptor. The interacting modes and contributions of individual residues were elucidated by analyzing the crystal structures of the single variants. To further increase the binding affinity, single positive mutations were combined, and two double mutants were shown to have about 3000- and 565-fold higher binding affinities than the wild-type decoy receptor. Molecular dynamics simulations and energetic analysis indicate that an additive effect by two mutations occurring at nearby modules was the major contributor to the remarkable increase in the binding affinities.     

Structural,energetic, and functional analysis of a protein-protein interface at distinct stages of affinity maturation

Due to a paucity of studies that synthesize structural, energetic, and functional analyses of a series of protein complexes representing distinct stages in an affinity maturation pathway, the biophysical basis for the molecular evolution of protein-protein interactions is poorly understood. Here, we combine crystal structures and binding-free energies of a series of variant superantigen (SAG)-major histocompatibility complex (MHC) class II complexes exhibiting increasingly higher affinity to reveal that this affinity maturation pathway is controlled largely by two biophysical factors: shape complementarity and buried hydrophobic surface. These factors, however, do not contribute equivalently to the affinity maturation of the interface, as the former dominates the early steps of the maturation process while the latter is responsible for improved binding in later steps. Functional assays reveal how affinity maturation of the SAG-MHC interface corresponds to T cell activation by SAGs.     

Dissecting the binding energy epitope of a high-affinity variant of human growth hormone: cooperative and additive effects from combining mutations from independently selected phage display mutagenesis libraries

Phage display mutagenesis is a widely used approach to engineering novel protein properties and is especially powerful in probing structure-function relationships in molecular recognition processes. The relative contributions of additive and cooperative binding forces and the influence of conformational diversity in producing a novel protein-protein interface is investigated using as a model an ultra-high-affinity receptor binding variant of human growth hormone (hGHv) that has been previously affinity matured. The modular aspect of how the mutations were grouped in the phage display libraries and combined allowed for a systematic probing of the inherent functional cross-talk between the different secondary structure elements that make up the remodeled hGHv binding surface. We performed an alanine scanning analyses of 35 hGHv residues and determined the kinetics of each variant by surface plasmon resonance (SPR). This analysis showed that there is a significant difference between the additive and cooperative binding forces existing among the selected residues in each library module, and the binding advantage of these residues is maximized over the original wild-type residue when in the context of the other mutations in the library. The degree to which residues in a particular mutagenesis library display binding cooperativity characteristics is generally correlated with the conformational plasticity of the polypeptide chain. Additionally, these cooperativity effects change when the mutations from one library are combined with the mutations from one or several of the other separate libraries. This supports the idea that significant functional cross-talk exists between the combined library modules that can affect the binding energetics of individual residues over a large distance.     

Characterization of T cell receptors engineered for high affinity against toxic shock syndrome toxin-1

Superantigens, including bacterial enterotoxins, are a family of proteins that bind simultaneously to MHC class II molecules and the Vbeta regions of T cell receptors. This cross-linking results in the activation of a large population of T cells that release massive amounts of inflammatory cytokines, ultimately causing a condition known as toxic shock syndrome. The staphylococcal superantigen toxic shock syndrome toxin-1 (TSST-1) is a causative agent of this disease, but its structure in complex with the cognate T cell receptor (human Vbeta2.1) has not been determined. To understand the molecular details of the interaction and to develop high affinity antagonists to TSST-1, we used directed evolution to generate a panel of high affinity receptors for TSST-1. Yeast display libraries of random and site-directed hVbeta2.1 mutants were selected for improved domain stability and for higher affinity binding to TSST-1. Stability mutations allowed the individual Vbeta domains to be expressed in a bacterial expression system. Affinity mutations were generated in CDR2 and FR3 residues, yielding improvements in affinity of greater than 10,000-fold (a K(D) value of 180 pmol). Alanine scanning mutagenesis of hVbeta2.1 wild-type and mutated residues allowed us to generate a map of the binding site for TSST-1 and to construct a docking model for the hVbeta2.1-TSST-1 complex. Our experiments suggest that the energetic importance of a single hVbeta2.1 wild-type residue likely accounts for the restriction of TSST-1 specificity to only this human Vbeta region. The high affinity mutants described here thus provide critical insight into the molecular basis of TSST-1 specificity and serve as potential leads toward the development of therapeutic agents for superantigen-mediated disease.     

X-ray snapshots of the maturation of an antibody response to a protein antigen

The process whereby the immune system generates antibodies of higher affinities during a response to antigen (affinity maturation) is a prototypical example of molecular evolution. Earlier studies have been confined to antibodies specific for small molecules (haptens) rather than for proteins. We compare the structures of four antibodies bound to the same site on hen egg white lysozyme (HEL) at different stages of affinity maturation. These X-ray snapshots reveal that binding is enhanced, not through the formation of additional hydrogen bonds or van der Waals contacts or by an increase in total buried surface, but by burial of increasing amounts of apolar surface at the expense of polar surface, accompanied by improved shape complementarity. The increase in hydrophobic interactions results from highly correlated rearrangements in antibody residues at the interface periphery, adjacent to the central energetic hot spot. This first visualization of the maturation of antibodies to protein provides insights into the evolution of high affinity in other protein-protein interfaces.     

High affinity T cell receptors from yeast display libraries block T cell activation by superantigens

The alphabeta T cell receptor (TCR) can be triggered by a class of ligands called superantigens. Enterotoxins secreted by bacteria act as superantigens by simultaneously binding to an MHC class II molecule on an antigen- presenting cell and to a TCR beta-chain, thereby causing activation of the T cell. The cross-reactivity of enterotoxins with different Vbeta regions can lead to stimulation of a large fraction of T cells. To understand the molecular details of TCR-enterotoxin interactions and to generate potential antagonists of these serious hyperimmune reactions, we engineered soluble TCR mutants with improved affinity for staphylococcal enterotoxin C3 (SEC3). A library of randomly mutated, single-chain TCRs (Vbeta-linker-Valpha) were expressed as fusions to the Aga2p protein on the surface of yeast cells. Mutants were selected by flow cytometric cell sorting with a fluorescent-labeled SEC3. Various mutations were identified, primarily in Vbeta residues that are located at the TCR:SEC3 interface. The combined mutations created a remodeled SEC3-binding surface and yielded a Vbeta domain with an affinity that was increased by 1000-fold (K(D)=7 nM). A soluble form of this Vbeta mutant was a potent inhibitor of SEC3-mediated T cell activity, suggesting that these engineered proteins may be useful as antagonists.     

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.     

Molecular basis of TCR selectivity, cross-reactivity, and allelic discrimination by a bacterial superantigen: integrative functional and energetic mapping of the SpeC-Vbeta2.1 molecular interface

Superantigens activate large fractions of T cells through unconventional interactions with both TCR beta-chain V domains (Vbetas) and MHC class II molecules. The bacterial superantigen streptococcal pyrogenic exotoxin C (SpeC) primarily stimulates human Vbeta2(+) T cells. Herein, we have analyzed the SpeC-Vbeta2.1 interaction by mutating all SpeC residues that make contact with Vbeta2.1 and have determined the energetic and functional consequences of these mutations. Our comprehensive approach, including mutagenesis, functional readouts from both bulk T cell populations, and an engineered Vbeta2.1(+) Jurkat T cell, as well as surface plasmon resonance binding analysis, has defined the SpeC "functional epitope" for TCR engagement. Although only two SpeC residues (Tyr(15) and Arg(181)) are critical for activation of virtually all human CD3(+) T cells, a larger cluster of four hot spot residues are required for interaction with Vbeta2.1. Three of these residues (Tyr(15), Phe(75), and Arg(181)) concentrate their binding energy on the CDR2 loop residue Ser(52a), a noncanonical residue insertion found only in Vbeta2 and Vbeta4 chains. Plasticity of this loop is important for recognition by SpeC. Although SpeC interacts with the Vbeta2.1 hypervariable CDR3 loop, our data indicate these contacts have little to no influence on the functional interaction with Vbeta2.1. These studies also provide a molecular basis for selectivity and cross-reactivity of SpeC-TCR recognition and reveal a degree of fine specificity in these interactions, whereby certain SpeC mutants are capable of distinguishing between different alleles of the same Vbeta domain subfamily.     

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.     

Understanding the molecular basis of MK2-p38α signaling complex assembly: insights into protein-protein interaction by molecular dynamics and free energy studies

The formation of a p38 MAPK and MAPK-activated protein kinase 2 (MK2) signaling complex is physiologically relevant to cellular responses such as the proinflammatory cytokine production. The interaction between p38α isoform and MK2 is of great importance for this signaling. In this study, molecular dynamics simulation and binding free energy calculation were performed on the MK2-p38α signaling complex to investigate the protein-protein interaction between the two proteins. Dynamic domain motion analyses were performed to analyze the conformational changes between the unbound and bound states of proteins during the interaction. The activation loop, αF-I helices, and loops among α helices in the C-lobe of MK2 are found to be highly flexible and exhibit significant changes upon p38α binding. The results also show that after the binding of p38α, the N- and C-terminal domains of MK2 display an opening and twisting motion centered on the activation loop. The molecular mechanics Poisson-Boltzmann and generalized-Born surface area (MM-PB/GBSA) methods were used to calculate binding free energies between MK2 and p38α. The analysis of the components of binding free energy calculation indicates that the van der Waals interaction and the nonpolar solvation energy provide the driving force for the binding process, while the electrostatic interaction contributes critically to the specificity, rather than to MK2-p38α binding affinity. The contribution of each residue at the interaction interface to the binding affinity of MK2 with p38α was also analyzed by free energy decomposition. Several important residues responsible for the protein-protein interaction were also identified.     

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.     

Analysis of antibody A6 binding to the extracellular interferon gamma receptor alpha-chain by alanine-scanning mutagenesis and random mutagenesis with phage display

The monoclonal antibody A6 binds a conformational epitope comprising mainly the CC' surface loop on the N-terminal fibronectin type-III domain of the extracellular interferon gamma receptor (IFNgammaR). The crystal structure of an A6 Fab-IFNgammaR complex revealed an interface rich in the aromatic side chains of Trp, Tyr, and His residues. These aromatic side chains appear to interact with both polar and hydrophobic groups at the interface, a property which, in general, may be advantageous for ligand binding. To analyze these interactions in more detail, the affinities of 19 A6 alanine-scanning mutants for the IFNgammaR have been measured, using engineered A6 single chain variable region fragments, and a surface plasmon resonance biosensor. Energetically important side chains (DeltaG(mutant) - DeltaG(wt) > 2.4 kcal/mol), that form distinct hot spots in the binding interface, have been identified on both proteins. These include V(L)W92 in A6, whose benzenoid ring appears well situated for a pi-cation (or pi-amine) interaction with the side chain of receptor residue K47 and simultaneously for T-stacking onto the indole ring of W82 in the receptor. At another site, energetically important residues V(H)W52 and V(H)W53, as well as V(H)D54 and V(H)D56, surround the aliphatic side chain of the hot receptor residue K52. Taken together, the results show that side chains distributed across the interface, including many aromatic ones, make key energetic contributions to binding. In addition, the receptor CC' loop has been subjected to random mutagenesis, and receptor mutants with high affinity for A6 have been selected by phage display. Residues previously identified as important for receptor binding to A6 were conserved in the clones isolated. Some mutants, however, showed a much improved affinity for A6, due to changes at Glu55, a residue that appeared to be energetically unimportant for binding the antibody by alanine-scanning mutagenesis. An E55P receptor mutant bound A6 with a 600-fold increase in affinity (K(D) approximately 20 pM), which is one of the largest improvements in affinity from a single point mutation reported so far at any protein-protein interface.     

Investigating the binding specificity of U1A-RNA by computational mutagenesis

The mammalian spliceosomal protein U1A binds a hairpin RNA with picomolar affinity. To examine the origin of this binding specificity, we carried out computational mutagenesis on protein and RNA residues in the U1A-RNA binding interface. Our computational mutagenesis methods calculate the relative binding affinity between mutant and wild-type as the sum of molecular mechanical energies and solvation free energies estimated with a continuum solvent model. We obtained good agreement with experimental studies and we verified mutations that abolish and improve binding. Therefore, we offer these methods as computationally inexpensive tools for investigating and predicting the effects of site-specific mutagenesis.     

The Structural Pathway of Interleukin 1 (IL-1) Initiated Signaling Reveals Mechanisms of Oncogenic Mutations and SNPs in Inflammation and Cancer

Interleukin-1 (IL-1) is a large cytokine family closely related to innate immunity and inflammation. IL-1 proteins are key players in signaling pathways such as apoptosis, TLR, MAPK, NLR and NF-κB. The IL-1 pathway is also associated with cancer, and chronic inflammation increases the risk of tumor development via oncogenic mutations. Here we illustrate that the structures of interfaces between proteins in this pathway bearing the mutations may reveal how. Proteins are frequently regulated via their interactions, which can turn them ON or OFF. We show that oncogenic mutations are significantly at or adjoining interface regions, and can abolish (or enhance) the protein-protein interaction, making the protein constitutively active (or inactive, if it is a repressor). We combine known structures of protein-protein complexes and those that we have predicted for the IL-1 pathway, and integrate them with literature information. In the reconstructed pathway there are 104 interactions between proteins whose three dimensional structures are experimentally identified; only 15 have experimentally-determined structures of the interacting complexes. By predicting the protein-protein complexes throughout the pathway via the PRISM algorithm, the structural coverage increases from 15% to 71%. In silico mutagenesis and comparison of the predicted binding energies reveal the mechanisms of how oncogenic and single nucleotide polymorphism (SNP) mutations can abrogate the interactions or increase the binding affinity of the mutant to the native partner. Computational mapping of mutations on the interface of the predicted complexes may constitute a powerful strategy to explain the mechanisms of activation/inhibition. It can also help explain how an oncogenic mutation or SNP works.     

Selection and analysis of an optimized anti-VEGF antibody: crystal structure of an affinity-matured Fab in complex with antigen.

The Fab portion of a humanized antibody (Fab-12; IgG form known as rhuMAb VEGF) to vascular endothelial growth factor (VEGF) has been affinity-matured through complementarity-determining region (CDR) mutation, followed by affinity selection using monovalent phage display. After stringent binding selections at 37 degrees C, with dissociation (off-rate) selection periods of several days, high affinity variants were isolated from CDR-H1, H2, and H3 libraries. Mutations were combined to obtain cumulatively tighter-binding variants. The final variant identified here, Y0317, contained six mutations from the parental antibody. In vitro cell-based assays show that four mutations yielded an improvement of about 100-fold in potency for inhibition of VEGF-dependent cell proliferation by this variant, consistent with the equilibrium binding constant determined from kinetics experiments at 37 degrees C. Using X-ray crystallography, we determined a high-resolution structure of the complex between VEGF and the affinity-matured Fab fragment. The overall features of the binding interface seen previously with wild-type are preserved, and many contact residues are maintained in precise alignment in the superimposed structures. However, locally, we see evidence for improved contacts between antibody and antigen, and two mutations result in increased van der Waals contact and improved hydrogen bonding. Site-directed mutants confirm that the most favorable improvements as judged by examination of the complex structure, in fact, have the greatest impact on free energy of binding. In general, the final antibody has improved affinity for several VEGF variants as compared with the parental antibody; however, some contact residues on VEGF differ in their contribution to the energetics of Fab binding. The results show that small changes even in a large protein-protein binding interface can have significant effects on the energetics of interaction.     

An ultra‐high affinity protein–protein interface displaying sequence‐robustness

Protein–protein interactions are crucial in biology and play roles in for example, the immune system, signaling pathways, and enzyme regulation. Ultra‐high affinity interactions (K _d <0.1 nM) occur in these systems, however, structures and energetics behind stability of ultra‐high affinity protein–protein complexes are not well understood. Regulation of the starch debranching barley limit dextrinase (LD) and its endogenous cereal type inhibitor (LDI) exemplifies an ultra‐high affinity complex (K _d of 42 pM). In this study the LD–LDI complex is investigated to unveil how robust the ultra‐high affinity is to LDI sequence variation at the protein–protein interface and whether alternative sequences can retain the ultra‐high binding affinity. The interface of LD–LDI was engineered using computational protein redesign aiming at identifying LDI variants predicted to retain ultra‐high binding affinity. These variants present a very diverse set of mutations going beyond conservative and alanine substitutions typically used to probe interfaces. Surface plasmon resonance analysis of the LDI variants revealed that high affinity of LD–LDI requires interactions of several residues at the rim of the protein interface, unlike the classical hotspot arrangement where key residues are found at the center of the interface. Notably, substitution of interface residues in LDI, including amino acids with functional groups different from the wild‐type, could occur without loss of affinity. This demonstrates that ultra‐high binding affinity can be conferred without hotspot residues, thus making complexes more robust to mutational drift in evolution. The present mutational analysis also demonstrates how energetic coupling can emerge between residues at large distances at the interface.     

Guiding a docking mode by phage display: selection of correlated mutations at the staphylokinase-plasmin interface.

During co-evolution of interacting proteins, functionally disruptive mutations on one side of the interface may be compensated by local amino acid changes on the other to restore binding affinity. This information can be useful for geometry-based docking approaches by reducing the translational and rotational space available to the proteins. Here, we demonstrate that correlated mutations at a protein-protein interface can be rapidly identified by selecting a phage-displayed library of a randomly mutated component of the complex for complementation of mutations that decreased binding in the interacting partner. This approach was used to deduce the binding mode of staphylokinase (Sak), a 15.5 kDa "indirect" plasminogen activator on microplasmin (microPli), the 28 kDa serine protease domain of plasmin. Biopanning indicated that residues Arg94 and Gly174 in microPli are located in close proximity to Glu75 and the Glu88:Ile128 pair in Sak, respectively. The coupled mutations Glu94<-->Lys75 reversed and Gly174<-->Lys88:Val128 introduced a salt bridge, whereby the binding affinities (with coupling energies of 1.8 to 2.3 kcal mol-1, respectively) and the plasminogen activation ability of the mutated complexes were partially restored. These findings suggested a unique docking mode of Sak at the western rim of the active-site cleft of microPli, that is in agreement with the structure of the Sak-microPli complex as recently derived by other methods.     

Dissection of the energetic coupling across the Src SH2 domain-tyrosyl phosphopeptide interface

Src Homology (SH2) domains play critical roles in signaling pathways by binding to phosphotyrosine (pTyr)-containing sequences, thereby recruiting SH2 domain-containing proteins to tyrosine-phosphorylated sites on receptor molecules. Investigations of the peptide binding specificity of the SH2 domain of the Src kinase (Src SH2 domain) have defined the EEI motif C-terminal to the phosphotyrosine as the preferential binding sequence. A subsequent study that probed the importance of eight specificity-determining residues of the Src SH2 domain found two residues which when mutated to Ala had significant effects on binding: Tyr beta D5 and Lys beta D3. The mutation of Lys beta D3 to Ala was particularly intriguing, since a Glu to Ala mutation at the first (+1) position of the EEI motif (the residue interacting with Lys beta D3) did not significantly affect binding. Hence, the interaction between Lys beta D3 and +1 Glu is energetically coupled. This study is focused on the dissection of the energetic coupling observed across the SH2 domain-phosphopeptide interface at and around the +1 position of the peptide. It was found that three residues of the SH2 domain, Lys beta D3, Asp beta C8 and AspCD2 (altogether forming the so-called +1 binding region) contribute to the selection of Glu at the +1 position of the ligand. A double (Asp beta C8Ala, AspCD2Ala) mutant does not exhibit energetic coupling between Lys beta D3 and +1 Glu, and binds to the pYEEI sequence 0.3 kcal/mol tighter than the wild-type Src SH2 domain. These results suggest that Lys beta D3 in the double mutant is now free to interact with the +1 Glu and that the role of Lys beta D3 in the wild-type is to neutralize the acidic patch formed by Asp beta C8 and AspCD2 rather than specifically select for a Glu at the +1 position as it had been hypothesized previously. A triple mutant (Lys beta D3Ala, Asp beta C8Ala, AspCD2Ala) has reduced binding affinity compared to the double (Asp beta C8Ala, AspCD2Ala) mutant, yet binds the pYEEI peptide as well as the wild-type Src SH2 domain. The structural basis for such high affinity interaction was investigated crystallographically by determining the structure of the triple (Lys beta D3Ala, Asp beta C8Ala, AspCD2Ala) mutant bound to the octapeptide PQpYEEIPI (where pY indicates a phosphotyrosine). This structure reveals for the first time contacts between the SH2 domain and the -1 and -2 positions of the peptide (i.e. the two residues N-terminal to pY). Thus, unexpectedly, mutations in the +1 binding region affect binding of other regions of the peptide. Such additional contacts may account for the high affinity interaction of the triple mutant for the pYEEI-containing peptide.