Structural basis for specificity switching of the Src SH2 domain

The Src SH2 domain binds pYEEI-containing phosphopeptides in an extended conformation with a hydrophobic pocket, which includes ThrEF1, binding Ile(pY +3). Mutating ThrEF1 to tryptophan switches specificity to an Asn(pY +2) requirement, yielding a biological mimic of the Grb2 SH2 domain. Here we show that the Src ThrEF1Trp SH2 domain mutant binds pYVNV phosphopeptides in a beta turn conformation, which, despite differing conformations of the interacting tryptophan, closely resembles the native Grb2/pYVNV cognate peptide binding mode. The ThrEF1Trp substitution therefore switches specificity by physically occluding the pTyr +3 binding pocket and by providing additional interaction surface area for Asn(pY +2). This demonstrates structurally how novel SH2 domain specificities may rapidly evolve through single amino acid substitutions and suggests how new signaling pathways may develop.     

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.     

A single amino acid in the SH3 domain of Hck determines its high affinity and specificity in binding to HIV-1 Nef protein.

We have examined the differential binding of Hck and Fyn to HIV-1 Nef to elucidate the structural basis of SH3 binding affinity and specificity. Full-length Nef bound to Hck SH3 with the highest affinity reported for an SH3-mediated interaction (KD 250 nM). In contrast to Hck, affinity of the highly homologous Fyn SH3 for Nef was too weak (KD > 20 microM) to be accurately determined. We show that this distinct specificity lies in a variable loop, the 'RT loop', positioned close to conserved SH3 residues implicated in the binding of proline-rich (PxxP) motifs. A mutant Fyn SH3 with a single amino acid substitution (R96I) in its RT loop had an affinity (KD 380 nM) for Nef comparable with that of Hck SH3. Based on additional mutagenesis studies we propose that the selective recognition of Nef by Hck SH3 is determined by hydrophobic interactions involving an isoleucine residue in its RT loop. Although Nef contains a PxxP motif which is necessary for the interaction with Hck SH3, high affinity binding was only observed for intact Nef protein. The binding of a peptide containing the Nef PxxP motif showed > 300-fold weaker affinity for Hck SH3 than full-length Nef.     

Boundary sequences of the NADPH oxidase p67(phox) C-terminal SH3 domain play on its specificity.

SH3 domains are found in many signal transduction proteins where they mediate protein-protein binding by recognizing specific peptides rich in proline. Based on the analysis of sequence alignment data, the NADPH oxidase p67(phox) C-terminal SH3 domain possesses a typical compact beta-barrel consisting of five beta-strands arranged in two antiparallel beta-sheets of three and two beta-strands. Multiple amino acid substitutions were made at beta e and its flanking residues to determine the role of the boundary sequences in binding activity and conformational specificity of the domain. Analysis of amino acid P55 indicated that all mutants were completely abolished in their binding activities. The substitution of F58 with Y58 showed no effect of the binding, whereas substitution with stop codon abolished activity. Furthermore, when amino acid V59 was substituted with stop codon, activity was also completely abolished. Substitution of E60 with stop codon showed no effect of binding. Moreover, our data show that V59 particularly could not be replaced by Leu. Taken together, these data suggest that V59 may not only contribute the exact boundary site but also play on the specificity for protein-protein interactions in phagocyte NADPH oxidase.     

In vitro evolution of recognition specificity mediated by SH3 domains reveals target recognition rules

We have designed a repertoire of 10(7) different SH3 domains by grafting the residues that are represented in the binding surfaces of natural SH3 domains onto the scaffold of the human Abl-SH3 domain. This phage-displayed library was screened by affinity selection for SH3 domains that bind to the synthetic peptides, APTYPPPLPP and LSSRPLPTLPSP, which are peptide ligands for the human Abl or Src SH3 domains, respectively. By characterizing the isolates, we have observed that as few as two or three amino acid substitutions lead to dramatic changes in recognition specificity. We propose that the ability to shift recognition specificity with a small number of amino acid replacements is an important evolutionary characteristic of protein binding modules. Furthermore, we have used the information obtained by these in vitro evolution experiments to generate a scoring matrix that evaluates the probability that any SH3 domain binds to the peptide ligands for the Abl and Src SH3 domains. A table of predictions for the 28 SH3 domains of baker's yeast is presented.     

Conservation analysis and structure prediction of the SH2 family of phosphotyrosine binding domains.

Src homology 2 (SH2) regions are short (approximately 100 amino acids), non-catalytic domains conserved among a wide variety of proteins involved in cytoplasmic signaling induced by growth factors. It is thought that SH2 domains play an important role in the intracellular response to growth factor stimulation by binding to phosphotyrosine containing proteins. In this paper we apply the techniques of multiple sequence alignment, secondary structure prediction and conservation analysis to 67 SH2 domain amino acid sequences. This combined approach predicts seven core secondary structure regions with the pattern beta-alpha-beta-beta-beta-beta-alpha, identifies those residues most likely to be buried in the hydrophobic core of the native SH2 domain, and highlights patterns of conservation indicative of secondary structural elements. Residues likely to be involved in phosphotyrosine binding are shown and orientations of the predicted secondary structures suggested which could enable such residues to cooperate in phosphate binding. We propose a consensus pattern that encapsulates the principal conserved features of the SH2 domains. Comparison of the proposed SH2 domain of akt to this pattern shows only 12/40 matches, suggesting that this domain may not exhibit SH2-like properties.     

Pleiotropy Can Be Effectively Estimated Without Counting Phenotypes Through the Rank of a Genotype–Phenotype Map

Although pleiotropy, the capability of a gene to affect multiple phenotypes, has been well known as one of the common gene properties, a quantitative estimation remains a great challenge, simply because of the phenotype complexity. Not surprisingly, it is hard for general readers to understand how, without counting phenotypes, gene pleiotropy can be effectively estimated from the genetics data. In this article we extensively discuss the Gu-2007 method that estimated pleiotropy from the protein sequence analysis. We show that this method is actually to estimate the rank (K) of genotype–phenotype mapping that can be concisely written as K = min(r, P_min), where P_min is the minimum pleiotropy among all legitimate measures including the fitness components, and r is the rank of mutational effects of an amino acid site. Together, the effective gene pleiotropy (K_e) estimated by the Gu-2007 method has the following meanings: (i) K_e is an estimate of K = min(r, P_min), the rank of a genotype–phenotype map; (ii) K_e is an estimate for the minimum pleiotropy P_min only if P_min < r; (iii) the Gu-2007 method attempted to estimate the pleiotropy of amino acid sites, a conserved proxy to the true gene pleiotropy; (iv) with a sufficiently large phylogeny such that the rank of mutational effects at an amino acid site is r → 19, one can estimate P_min between 1 and 19; and (v) K_e is a conserved estimate of K because those slightly affected components in fitness have been effectively removed by the estimation procedure. In addition, we conclude that mutational pleiotropy (number of traits affected by a single mutation) cannot be estimated without knowing the phenotypes.     

Characterization of the specific DNA-binding properties of Tnp26, the transposase of insertion sequence IS26

The bacterial insertion sequence (IS) IS26 mobilizes and disseminates antibiotic resistance genes. It differs from bacterial IS that have been studied to date as it exclusively forms cointegrates via either a copy-in (replicative) or a recently discovered targeted conservative mode. To investigate how the Tnp26 transposase recognizes the 14-bp terminal inverted repeats (TIRs) that bound the IS, amino acids in two domains in the N-terminal (amino acids M1–P56) region were replaced. These changes substantially reduced cointegration in both modes. Tnp26 was purified as a maltose-binding fusion protein and shown to bind specifically to dsDNA fragments that included an IS26 TIR. However, Tnp26 with an R49A or a W50A substitution in helix 3 of a predicted trihelical helix–turn–helix domain (amino acids I13–R53) or an F4A or F9A substitution replacing the conserved amino acids in a unique disordered N-terminal domain (amino acids M1–D12) did not bind. The N-terminal M1–P56 fragment also bound to the TIR but only at substantially higher concentrations, indicating that other parts of Tnp26 enhance the binding affinity. The binding site was confined to the internal part of the TIR, and a G to T nucleotide substitution in the TGT at positions 6 to 8 of the TIR that is conserved in most IS26 family members abolished binding of both Tnp26 (M1–M234) and Tnp26 M1–P56 fragment. These findings indicate that the helix–turn–helix and disordered domains of Tnp26 play a role in Tnp26–TIR complex formation. Both domains are conserved in all members of the IS26 family.     

Hierarchy of simulation models in predicting molecular recognition mechanisms from the binding energy landscapes: structural analysis of the peptide complexes with SH2 domains.

Computer simulations using the simplified energy function and simulated tempering dynamics have accurately determined the native structure of the pYVPML, SVLpYTAVQPNE, and SPGEpYVNIEF peptides in the complexes with SH2 domains. Structural and equilibrium aspects of the peptide binding with SH2 domains have been studied by generating temperature-dependent binding free energy landscapes. Once some native peptide-SH2 domain contacts are constrained, the underlying binding free energy profile has the funnel-like shape that leads to a rapid and consistent acquisition of the native structure. The dominant native topology of the peptide-SH2 domain complexes represents an extended peptide conformation with strong specific interactions in the phosphotyrosine pocket and hydrophobic interactions of the peptide residues C-terminal to the pTyr group. The topological features of the peptide-protein interface are primarily determined by the thermodynamically stable phosphotyrosyl group. A diversity of structurally different binding orientations has been observed for the amino-terminal residues to the phosphotyrosine. The dominant native topology for the peptide residues carboxy-terminal to the phosphotyrosine is tolerant to flexibility in this region of the peptide-SH2 domain interface observed in equilibrium simulations. The energy landscape analysis has revealed a broad, entropically favorable topology of the native binding mode for the bound peptides, which is robust to structural perturbations. This could provide an additional positive mechanism underlying tolerance of the SH2 domains to hydrophobic conservative substitutions in the peptide specificity region.     

The relationship between the pY+3 amino acid in the phosphopeptide ligands and the specificity for various SH2 domains

After consideration of the tertiary structure of the pY+3 binding pocket in SH2 domains, isoleucine of Ac-pYEEIE was replaced with various unnatural aliphatic amino acids and the binding affinities for Lck, Src, and Fyn SH2 domains were measured. We determined the characteristics of the amino acid at the pY+3 position in the phosphopeptide ligands for the specificity for the SH2 domains.     

The active site of hydroxynitrile lyase from Prunus amygdalus: Modeling studies provide new insights into the mechanism of cyanogenesis

The FAD-dependent hydroxynitrile lyase from almond (Prunus amygdalus, PaHNL) catalyzes the cleavage of R-mandelonitrile into benzaldehyde and hydrocyanic acid. Catalysis of the reverse reaction-the enantiospecific formation of alpha-hydroxynitriles--is now widely utilized in organic syntheses as one of the few industrially relevant examples of enzyme-mediated C-C bond formation. Starting from the recently determined X-ray crystal structure, systematic docking calculations with the natural substrate were used to locate the active site of the enzyme and to identify amino acid residues involved in substrate binding and catalysis. Analysis of the modeled substrate complexes supports an enzymatic mechanism that includes the flavin cofactor as a mere "spectator" of the reaction and relies on general acid/base catalysis by the conserved His-497. Stabilization of the negative charge of the cyanide ion is accomplished by a pronounced positive electrostatic potential at the binding site. PaHNL activity requires the FAD cofactor to be bound in its oxidized form, and calculations of the pKa of enzyme-bound HCN showed that the observed inactivation upon cofactor reduction is largely caused by the reversal of the electrostatic potential within the active site. The suggested mechanism closely resembles the one proposed for the FAD-independent, and structurally unrelated HNL from Hevea brasiliensis. Although the actual amino acid residues involved in the catalytic cycle are completely different in the two enzymes, a common motif for the mechanism of cyanogenesis (general acid/base catalysis plus electrostatic stabilization of the cyanide ion) becomes evident.     

Structural basis for ubiquitin recognition by SH3 domains

The SH3 domain is a protein-protein interaction module commonly found in intracellular signaling and adaptor proteins. The SH3 domains of multiple endocytic proteins have been recently implicated in binding ubiquitin, which serves as a signal for diverse cellular processes including gene regulation, endosomal sorting, and protein destruction. Here we describe the solution NMR structure of ubiquitin in complex with an SH3 domain belonging to the yeast endocytic protein Sla1. The ubiquitin binding surface of the Sla1 SH3 domain overlaps substantially with the canonical binding surface for proline-rich ligands. Like many other ubiquitin-binding motifs, the SH3 domain engages the Ile44 hydrophobic patch of ubiquitin. A phenylalanine residue located at the heart of the ubiquitin-binding surface of the SH3 domain serves as a key specificity determinant. The structure of the SH3-ubiquitin complex explains how a subset of SH3 domains has acquired this non-traditional function.     

Sequence and structural aspects of the functional diversification of plant alcohol dehydrogenases

The glycolytic proteins in plants are coded by small multigene families, which provide an interesting contrast to the high copy number of gene families studied to date. The alcohol dehydrogenase (Adh) genes encode glycolytic enzymes that have been characterized in some plant families. Although the amino acid sequences of zinc-containing long-chain ADHs are highly conserved, the metabolic function of this enzyme is variable. They also have different patterns of expression and are submitted to differences in nonsynonymous substitution rates between gene copies. It is possible that the Adh copies have been retained as a consequence of adaptative amino acid replacements which have conferred subtle changes in function. Phylogenetic analysis indicates that there have been a number of separate duplication events within angiosperms, and that genes labeled Adh1, Adh2 and Adh3 in different groups may not be homologous. Nonsynonymous/synonymous ratios yielded no signs of positive selection. However, the coefficients of functional divergence (theta) estimated between the Adh1 and Adh2 gene groups indicate statistically significant site-specific shift of evolutionary rates between them, as well as between those of different botanical families, suggesting that altered functional constraints may have taken place at some amino acid residues after their diversification. The theoretical three-dimensional structure of the alcohol dehydrogenase from Arabis blepharophylla was constructed and verified to be stereochemically valid.     

Predicting Protein Function and Binding Profile via Matching of Local Evolutionary and Geometric Surface Patterns

Inferring protein functions from structures is a challenging task, as a large number of orphan protein structures from structural genomics project are now solved without their biochemical functions characterized. For proteins binding to similar substrates or ligands and carrying out similar functions, their binding surfaces are under similar physicochemical constraints, and hence the sets of allowed and forbidden residue substitutions are similar. However, it is difficult to isolate such selection pressure due to protein function from selection pressure due to protein folding, and evolutionary relationship reflected by global sequence and structure similarities between proteins is often unreliable for inferring protein function. We have developed a method, called pevoSOAR (pocket-based evolutionary search of amino acid residues), for predicting protein functions by solving the problem of uncovering amino acids residue substitution pattern due to protein function and separating it from amino acids substitution pattern due to protein folding. We incorporate evolutionary information specific to an individual binding region and match local surfaces on a large scale with millions of precomputed protein surfaces to identify those with similar functions. Our pevoSOAR method also generates a probablistic model called the computed binding a profile that characterizes protein-binding activities that may involve multiple substrates or ligands. We show that our method can be used to predict enzyme functions with accuracy. Our method can also assess enzyme binding specificity and promiscuity. In an objective large-scale test of 100 enzyme families with thousands of structures, our predictions are found to be sensitive and specific: At the stringent specificity level of 99.98%, we can correctly predict enzyme functions for 80.55% of the proteins. The overall area under the receiver operating characteristic curve measuring the performance of our prediction is 0.955, close to the perfect value of 1.00. The best Matthews coefficient is 86.6%. Our method also works well in predicting the biochemical functions of orphan proteins from structural genomics projects.     

Identification of the β-Lactamase Inhibitor Protein-II (BLIP-II) Interface Residues Essential for Binding Affinity and Specificity for Class A β-Lactamases

The interactions between β-lactamase inhibitory proteins (BLIPs) and β-lactamases have been used as model systems to understand the principles of affinity and specificity in protein-protein interactions. The most extensively studied tight binding inhibitor, BLIP, has been characterized with respect to amino acid determinants of affinity and specificity for binding β-lactamases. BLIP-II, however, shares no sequence or structural homology to BLIP and is a femtomolar to picomolar potency inhibitor, and the amino acid determinants of binding affinity and specificity are unknown. In this study, alanine scanning mutagenesis was used in combination with determinations of on and off rates for each mutant to define the contribution of residues on the BLIP-II binding surface to both affinity and specificity toward four β-lactamases of diverse sequence. The residues making the largest contribution to binding energy are heavily biased toward aromatic amino acids near the center of the binding surface. In addition, substitutions that reduce binding energy do so by increasing off rates without impacting on rates. Also, residues with large contributions to binding energy generally exhibit low temperature factors in the structures of complexes. Finally, with the exception of D206A, BLIP-II alanine substitutions exhibit a similar trend of effect for all β-lactamases, i.e., a substitution that reduces affinity for one β-lactamase usually reduces affinity for all β-lactamases tested.     

Six X-linked agammaglobulinemia-causing missense mutations in the Src homology 2 domain of Bruton's tyrosine kinase: phosphotyrosine-binding and circular dichroism analysis

Src homology 2 (SH2) domains recognize phosphotyrosine (pY)-containing sequences and thereby mediate their association to ligands. Bruton's tyrosine kinase (Btk) is a cytoplasmic protein tyrosine kinase, in which mutations cause a hereditary immunodeficiency disease, X-linked agammaglobulinemia (XLA). Mutations have been found in all Btk domains, including SH2. We have analyzed the structural and functional effects of six disease-related amino acid substitutions in the SH2 domain: G302E, R307G, Y334S, L358F, Y361C, and H362Q. Also, we present a novel Btk SH2 missense mutation, H362R, leading to classical XLA. Based on circular dichroism analysis, the conformation of five of the XLA mutants studied differs from the native Btk SH2 domain, while mutant R307G is structurally identical. The binding of XLA mutation-containing SH2 domains to pY-Sepharose was reduced, varying between 1 and 13% of that for the native SH2 domain. The solubility of all the mutated proteins was remarkably reduced. SH2 domain mutations were divided into three categories: 1) Functional mutations, which affect residues presumably participating directly in pY binding (R307G); 2) structural mutations that, via conformational change, not only impair pY binding, but severely derange the structure of the SH2 domain and possibly interfere with the overall conformation of the Btk molecule (G302E, Y334S, L358F, and H362Q); and 3) structural-functional mutations, which contain features from both categories above (Y361C).     

Reconstitution of a native-like SH2 domain from disordered peptide fragments examined by multidimensional heteronuclear NMR

The N-terminal SH2 domain from the p85alpha subunit of phosphatidylinositol 3' kinase is cleaved specifically into 9- and 5-kD fragments by limited proteolytic digestion with trypsin. The noncovalent SH2 domain complex and its constituent tryptic peptides have been investigated using high-resolution heteronuclear magnetic resonance (NMR). These studies have established the viability of the SH2 domain as a fragment complementation system. The individual peptide fragments are predominantly unstructured in solution. In contrast, the noncovalent 9-kD + 5-kD complex shows a native-like (1)H-(15)N HSQC spectrum, demonstrating that the two fragments fold into a native-like structure on binding. Chemical shift analysis of the noncovalent complex compared to the native SH2 domain reveals that the highest degree of perturbation in the structure occurs at the cleavage site within a flexible loop and along the hydrophobic interface between the two peptide fragments. Mapping of these chemical shift changes on the structure of the domain reveals changes consistent with the reduction in affinity for the target peptide ligand observed in the noncovalent complex relative to the intact protein. The 5-kD fragment of the homologous Src protein is incapable of structurally complementing the p85 9-kD fragment, either in complex formation or in the context of the full-length protein. These high-resolution structural studies of the SH2 domain fragment complementation features establish the suitability of the system for further protein-folding and design studies.     

Two Amino Acid Residues Confer Different Binding Affinities of Abelson Family Kinase Src Homology 2 Domains for Phosphorylated Cortactin

The closely related Abl family kinases, Arg and Abl, play important non-redundant roles in the regulation of cell morphogenesis and motility. Despite similar N-terminal sequences, Arg and Abl interact with different substrates and binding partners with varying affinities. This selectivity may be due to slight differences in amino acid sequence leading to differential interactions with target proteins. We report that the Arg Src homology (SH) 2 domain binds two specific phosphotyrosines on cortactin, a known Abl/Arg substrate, with over 10-fold higher affinity than the Abl SH2 domain. We show that this significant affinity difference is due to the substitution of arginine 161 and serine 187 in Abl to leucine 207 and threonine 233 in Arg, respectively. We constructed Abl SH2 domains with R161L and S187T mutations alone and in combination and find that these substitutions are sufficient to convert the low affinity Abl SH2 domain to a higher affinity “Arg-like” SH2 domain in binding to a phospho-cortactin peptide. We crystallized the Arg SH2 domain for structural comparison to existing crystal structures of the Abl SH2 domain. We show that these two residues are important determinants of Arg and Abl SH2 domain binding specificity. Finally, we expressed Arg containing an “Abl-like” low affinity mutant Arg SH2 domain (L207R/T233S) and find that this mutant, although properly localized to the cell periphery, does not support wild type levels of cell edge protrusion. Together, these observations indicate that these two amino acid positions confer different binding affinities and cellular functions on the distinct Abl family kinases.     

In vitro scanning saturation mutagenesis of all the specificity determining residues in an antibody binding site.

For the first time, each specificity determining residue (SDR) in the binding site of an antibody has been replaced with every other possible single amino acid substitution, and the resulting mutants analyzed for binding affinity and specificity. The studies were conducted on a variant of the 26-10 antidigoxin single chain Fv (scFv) using in vitro scanning saturation mutagenesis, a new process that allows the high throughput production and characterization of antibody mutants [Burks,E.A., Chen,G., Georgiou,G. and Iverson,B.L. (1997) Proc. Natl Acad. Sci. USA, 94, 412-417]. Single amino acid mutants of 26-10 scFv were identified that modulated specificity in dramatic fashion. The overall plasticity of the antibody binding site with respect to amino acid replacement was also evaluated, revealing that 86% of all mutants retained measurable binding activity. Finally, by analyzing the physical properties of amino acid substitutions with respect to their effect on hapten binding, conclusions were drawn regarding the functional role played by the wild-type residue at each SDR position. The reported results highlight the value of in vitro scanning saturation mutagenesis for engineering antibody binding specificity, for evaluating the plasticity of proteins, and for comprehensive structure-function studies and analysis.