Intrinsic pKas of ionizable residues in proteins: An explicit solvent calculation for lysozyme

Molecular dynamics simulations of triclinic hen egg white lysozyme in aqueous solution were performed to calculate the intrinsic pK_as of 14 ionizable residues. An all-atom model was used for both solvent and solute, and a single 180 ps simulation in conjunction with a Gaussian fluctuation analysis method was used. An advantage of the Gaussian fluctuation method is that it only requires a single simulation of the system in a reference state to calculate all the pK_as in the protein, in contrast to multiple simulations for the free energy perturbation method. pK_int shifts with respect to reference titratable residues were evaluated and compared to results obtained using a finite difference Poisson-Boltzmann (FDPB) method with a continuum solvent model; overall agreement with the direction of the shifts was generally observed, though the magnitude of the shifts was typically larger with the explicit solvent model. The contribution of the first solvation shell to the total charging free energies of the titratable groups was explicitly evaluated and found to be significant. Dielectric shielding between pairs of titratable groups was examined and found to be smaller than expected. The effect of the approximations used to treat the long-range interactions on the pK_int shifts is discussed. © 1994 Wiley-Liss, Inc.     

Functional tuning of the catalytic residue pKa in a de novo designed esterase

AlleyCatE is a de novo designed esterase that can be allosterically regulated by calcium ions. This artificial enzyme has been shown to hydrolyze p‐nitrophenyl acetate (pNPA) and 4‐nitrophenyl‐(2‐phenyl)‐propanoate (pNPP) with high catalytic efficiency. AlleyCatE was created by introducing a single‐histidine residue (His¹⁴⁴) into a hydrophobic pocket of calmodulin. In this work, we explore the determinants of catalytic properties of AlleyCatE. We obtained the pK_a value of the catalytic histidine using experimental measurements by NMR and pH rate profile and compared these values to those predicted from electrostatics pK_a calculations (from both empirical and continuum electrostatics calculations). Surprisingly, the pK_a value of the catalytic histidine inside the hydrophobic pocket of calmodulin is elevated as compared to the model compound pK_a value of this residue in water. We determined that a short‐range favorable interaction with Glu¹²⁷ contributes to the elevated pK_a of His¹⁴⁴. We have rationally modulated local electrostatic potential in AlleyCatE to decrease the pK_a of its active nucleophile, His¹⁴⁴, by 0.7 units. As a direct result of the decrease in the His¹⁴⁴ pK_a value, catalytic efficiency of the enzyme increased by 45% at pH 6. This work shows that a series of simple NMR experiments that can be performed using low field spectrometers, combined with straightforward computational analysis, provide rapid and accurate guidance to rationally improve catalytic efficiency of histidine‐promoted catalysis. Proteins 2017; 85:1656–1665. © 2017 Wiley Periodicals, Inc.     

Correlation of exchangeable (NH) and nonexchangeable (C2H) histidine resonances in the proton NMR spectrum of ribonuclease A in aqueous solution.

The ionization characteristics of the hydrogen-bonded His 12 N₁ proton observed to titrate between 11 to 13 ppm in the nmr spectrum of ribonuclease A in H₂O solution are compared with the ionization characteristics of the four histidine C₂ protons in the enzyme. Comparison of the pK_a's of the enzyme in H₂O and D₂O in the absence and presence of cytidine monophosphate (?5′, ?3′, and ?2′) inhibitors, line widths in the presence of Cu II at pH 3.6 and 5.6, and chemical shifts in the presence of AgNO₃ permit a correlation of the exchangeable His 12 N₁ proton with the active site histidine C₂ proton exhibiting the lower ionization pK_a. The histidines with pK_a of 5.1 and 5.6 in ribonuclease A in the absence of salt are assigned in this study to His 12 and His 119, respectively.     

Reassignment of the active site histidines in ribonuclease A by selective deuteration studies.

The C₂H resonance of the active site histidine residue designated AS-2, which has the lower pK_a of the two active site histidines, has been correlated in both RNase A and RNase S by comparing the pH 3 to 5.5 regions of the chemical shift titration curves, the effect of the inhibitor CMP-3′ on the chemical shifts at pH 4.0, and the effect of Cu II on the line widths at pH 3.6. It has been demonstrated that resonance AS-2 is absent in the spectrum of RNase S′ reconstituted using S-peptide deuterated at the C₂ of His 12, and in that of the RNase S′-CMP-3′ complex. We thus demonstrate that histidine AS-2 is in fact His 12 in both enzymes. This finding is in agreement with out previous assignment of the exchangeable NH proton in RNase A to His 12, but reverses the assignments of the active site histidine C₂H resonances made earlier by other authors.     

Comparative studies of the interaction of human and bovine platelet factor 4 with heparin using histidine NMR resonances as spectroscopic probes

ThepK _a values of His-38 and His-50 of the heparin-binding protein, bovine platelet factor 4, are 5.6 and 6.5, respectively, as determined by¹H NMR spectroscopy. The¹H NMR resonance of His-38 of bovine platelet factor 4 which exhibits the lowerpK _a value is perturbed upon heparin binding to a greater degree than the resonance of His-50. Human platelet factor 4 contains the homologous residues His-23 and His-35. ThepK _a values of the two histidine residues of human platelet factor 4 are 5.3 and 6.4. The¹H NMR resonance of the histidine of human platelet factor 4 exhibiting the lowerpK _a value also is perturbed upon heparin binding to a greater degree than the histidine resonance exhibiting the higherpK _a , thereby suggesting comparable heparin-protein interactions in bovine and human platelet factor 4.     

The pKa values of two histidine residues in human haemoglobin, the Bohr effect, and the dipole moments of alpha-helices

Studies of abnormal and chemically modified haemoglobins indicate that in 0.1 m-NaCl about 40% of the alkaline Bohr effect of human haemoglobin is contributed by the C-terminal histidine HC3(146)β. In deoxyhaemoglobin, the imidazole of this histidine forms a salt bridge with aspartate FG1(94)β, in oxyhaemoglobin or carbonmonoxyhaemoglobin it accepts a hydrogen bond from its own NH group instead. Kilmartin et al. (1973) showed that in 0.2 m-NaCl + 0.2 m-phosphate this change of ligation lowered the pK_a of the histidine from 8.0 in Hb³ to 7.1 in HbCO, but Russu et al. (1980) claimed that in bis-Tris buffer without added NaCl its pK_a in HbCO dropped no lower than 7.85, and that in this medium the C-terminal histidine made only a negligible contribution to the alkaline Bohr effect.We have compared the histidine resonances of HbCO A with those of three abnormal haemoglobins: HbCO Cowtown (His HC3(146)β → Leu), HbCO Wood (His FG4(97)β → Leu) and HbCO Malmø (His FG4(97)β → Gln). Our results show that the resonance assigned by Russu et al. to His HC3(146)β in fact belongs to His FG4(97)β. Although in Hb the pK_a of His HC3(146)β is 8.05 ± 0.05 independent of ionic strength, in HbCO its pK_a drops sharply with diminishing ionic strength, so that in the buffer employed by Russu et al. it has a pK_a of 6.2 and makes a contribution to the alkaline Bohr effect that is 57% larger than in the phosphate buffer employed by Kilmartin et al. (1973).In HbCO A, His FG4(97)β does not contribute to the Bohr effect, but in HbCO from which His HC3(146)β has been cleaved (HbCO des-His), His FG4(97)β is in equilibrium between two conformations with different pK_a values. This equilibrium varies with ionic strength and pH, and presumably also with degree of ligation of the haem moiety.In HbCO A, His FG4(97)β has a pK_a of 7.8 compared to the pK_a value of about 6.6 characteristic of free histidines at the surface of proteins. This high pK_a is accounted for by its interaction with the negative pole at the C terminus of helices F and FG. It corresponds to a free energy change of the same order as that observed in the interaction of histidines with carboxylate ions and confirms the strongly dipolar character of α-helices, which manifests itself even when they lie on the surface of the protein.     

Assignment of the histidine 12-threonine 45 hydrogen-bonded proton in the NMR spectrum of ribonuclease A in H2O.

Described herein are proton nmr experiments on chemically modified derivatives of ribonuclease A designed to elucidate the origin of an exchangeable resonance, assigned previously to a histidine ring N proton that titrates between 11 to 13 ppm with a pK_a of 6.1 in H₂O solution. Histidines 48 and 105, which are distant from the active site, are eliminated as candidates for this resonance from inhibitor binding studies on the enzyme in acetate–water solutions. This exchangeable resonance titrates with modified pK_a's and constant area over the above pH range in His-119-N₁-carboxymethylated-RNase A and des-(121–124)-RNase A, thus eliminating the imidazole N₃ proton in the His 119-Asp 121 hydrogen bond. In His-12-N₁-carboxymethylated-RNase A, this resonance is also observable, but broadens on raising the pH above 7 and at elevated temperatures above neutrality. It exhibits a pH-independent chemical shift characteristic of the protonated state of histidine. On the basis of these findings, this exchangeable resonance, designated a, is assigned to the imidazole N₁ proton of His 12, which is hydrogen-bonded to the carbonyl oxygen of Thr 45 in the crystal.     

Origin of the pKa shift of the catalytic lysine in acetoacetate decarboxylase

The pK_a value of Lys115, the catalytic residue in acetoacetate decarboxylate, was calculated using atomic coordinates of the X-ray crystal structure with consideration of the protonation states of all titratable sites in the protein. The calculated pK_a value of Lys115 (pK_a(Lys115)) was unusually low (≈6) in agreement with the experimentally measured value. Although charged residues impact pK_a(Lys115) considerably in the native protein, the significant pK_a(Lys115) downshift in the protein with respect to aqueous solution was mainly due to loss of the solvation energy in the catalytic active site relative to bulk water.     

1H nuclear magnetic resonance studies of histidine-containing di- and tripeptides. Estimation of the effects of charged groups on the pKa value of the imidazole ring.

The nmr titration curves of chemical shifts versus pH were observed for the protons of various histidine-containing di- and tripeptides. With these results, the macroscopic pK_a values and the chemical shifts intrinsic to each ionic species were determined by a computer curve-fitting based on a simple acid dissociation sequence. The pK_a value of the imidazole ring in N-acetyl-L -histidine methylamide was assumed to represent the intrinsic (or unperturbed) pK_a of the imidazole rings of histidine having peptide linkages at both the CO and NH sides. The pK_a values of the imidazole rings observed for most di- and tripeptides were reasonably reproduced by simple calculations using the intrinsic value and the perturbations due to the CO₂^? and NH₃⁺ groups located at various positions. Some other factors affecting the pK_a value of the imidazole ring are also discussed.     

An essential role of active site arginine residue in iodide binding and histidine residue in electron transfer for iodide oxidation by horseradish peroxidase

The objective of the present study is to delineate the role of active site arginine and histidine residues of horseradish peroxidase (HRP) in controlling iodide oxidation using chemical modification technique. The arginine specific reagent, phenylglyoxal (PGO) irreversibly blocks iodide oxidation following pseudofirst order kinetics with second order rate constant of 25.12 min^-1 M^-1. Radiolabelled PGO incorporation studies indicate an essential role of a single arginine residue in enzyme inactivation. The enzyme can be protected both by iodide and an aromatic donor such as guaiacol. Moreover, guaiacol-protected enzyme can oxidise iodide and iodide-protected enzyme can oxidise guaiacol suggesting the regulatory role of the same active site arginine residue in both iodide and guaiacol binding. The protection constant (K_p) for iodide and guaiacol are 500 and 10 M respectively indicating higher affinity of guaiacol than iodide at this site. Donor binding studies indicate that guaiacol competitively inhibits iodide binding suggesting their interaction at the same binding site. Arginine-modified enzyme shows significant loss of iodide binding as shown by increased K_d value to 571 mM from the native enzyme (K_d = 150 mM). Although arginine-modified enzyme reacts with H₂O₂ to form compound II presumably at a slow rate, the latter is not reduced by iodide presumably due to low affinity binding.The role of the active site histidine residue in iodide oxidation was also studied after disubstitution reaction of the histidine imidazole nitrogens with diethylpyrocarbonate (DEPC), a histidine specific reagent. DEPC blocks iodide oxidation following pseudofirst order kinetics with second order rate constant of 0.66 min^-1 M^-1. Both the nitrogens (, ) of histidine imidazole were modified as evidenced by the characteristic peak at 222 nm. The enzyme is not protected by iodide suggesting that imidazolium ion is not involved in iodide binding. Moreover, DEPC-modified enzyme binds iodide similar to the native enzyme. However, the modified enzyme does not form compound II but forms compound I only with higher concentration of H₂O₂ suggesting the catalytic role of this histidine in the formation and autoreduction of compound I. Interestingly, compound I thus formed is not reduced by iodide indicating block of electron transport from the donor to the compound I. We suggest that an active site arginine residue regulates iodide binding while the histidine residue controls the electron transfer to the heme ferryl group during oxidation.     

Myoglobin discriminates between O2, NO, and CO by electrostatic interactions with the bound ligand

 Most biological substrates have distinctive sizes, shapes, and charge distributions which can be recognized specifically by proteins. In contrast, myoglobin must discriminate between the diatomic gases O₂, CO, and NO which are apolar and virtually the same size. Selectivity occurs at the level of the covalent Fe-ligand complexes, which exhibit markedly different bond strengths and electrostatic properties. By pulling a water molecule into the distal pocket, His64(E7)¹ inhibits the binding of all three ligands by a factor of ∼10 compared to that observed for protoheme-imidazole complexes in organic solvents. In the case of O₂ binding, this unfavorable effect is overcome by the formation of a strong hydrogen bond between His64(E7) and the highly polar FeO₂ complex. This favorable electrostatic interaction stabilizes the bound O₂ by a factor of ∼1000, and the net result is a 100-fold increase in overall affinity compared to model hemes or mutants with an apolar residue at position 64. Electrostatic interaction between FeCO and His64 is very weak, resulting in only a two- to three-fold stabilization of the bound state. In this case, the inhibitory effect of distal pocket water dominates, and a net fivefold reduction in K _CO is observed for the wild-type protein compared to mutants with an apolar residue at position 64. Bound NO is stabilized ∼tenfold by hydrogen bonding to His64. This favorable interaction with FeNO exactly compensates for the tenfold inhibition due to the presence of distal pocket water, and the net result is little change in K _NO when the distal histidine is replaced with apolar residues. Thus, it is the polarity of His64 which allows discrimination between the diatomic gases. Direct steric hindrance by this residue plays a minor role as judged by: (1) the independence of K _O2, K _CO, and K _NO on the size of apolar residues inserted at position 64, and (2) the observation of small decreases, not increases, in CO affinity when the mobility of the His64 side chain is increased. Val68(E11) does appear to hinder selectively the binding of CO. However, the extent is no more than a factor of 2–5, and much smaller than electrostatic stabilization of bound O₂ by the distal histidine. Received, accepted: 23 May 1997     

Theoretical studies on the redox-Bohr effect in cytochrome c 3 from Desulfovibrio vulgaris Hildenborough

 The pH dependence of the redox potentials in the tetrahemic cytochrome c ₃ from Desulfovibrio vulgaris Hildenborough (redox-Bohr effect) is here investigated using continuum electrostatics methods. The redox-Bohr effect seems to be associated with changes in the protonation state of charged residues in the protein, but the exact residues had not been identified. The global pK _a of this phenomenon is dependent on the redox state of the molecule, and the influence of the pH on the microscopic potential of each heme has been experimentally quantified. The availability of detailed experimental data provides us with important and unique guides to the performance of ab initio pK _a calculations aiming at the identification of the groups involved. These calculations were performed in several redox states along the reduction pathway, with the double objective of finding groups with redox-linked pK _a shifts, and absolute pK _as compatible with the redox-Bohr effect. The group with the largest pK _a shift along the reduction pathway is propionate D from heme I. Its effect on the redox potential of individual hemes, as calculated by electrostatic calculations, correlates very well with the experimental order of influence, making it a likely candidate. Abnormal titration of the same propionate has been experimentally observed on a homologous cytochrome c ₃ from a different strain, thus strengthening the theoretical result. However, its absolute calculated pK _a in the fully oxidised cytochrome is outside the zone where the phenomenon is known to occur, but the calculation shows a strong dependence on small conformational changes, suggesting large uncertainties in the calculated value. A group with a pK _a value within the experimentally observed range is propionate D from heme IV. Its influence on the redox potential of the hemes does not correlate with the experimental order, indicating that, although it may be one of the possible players on the phenomenon, it cannot be solely responsible for it. Mutation of the Lys45 residue is suggested as an indirect way of probing the importance of the propionate D from heme I in the mechanism. Non-heme groups may also be involved in this process; our calculations indicate His67 and the N-terminal as groups that may play a role. Accuracy and applicability of current continuum electrostatic methods are discussed in the context of this system. Received: 27 March 1997 / Accepted: 19 August 1997     

Modulation of an Active-Site Cysteine pKa Allows PDI to Act as a Catalyst of both Disulfide Bond Formation and Isomerization

Protein disulfide isomerase (PDI) plays a central role in disulfide bond formation in the endoplasmic reticulum. It is implicated both in disulfide bond formation and in disulfide bond reduction and isomerization. To be an efficient catalyst of all three reactions requires complex mechanisms. These include mechanisms to modulate the pK_a values of the active-site cysteines of PDI. Here, we examined the role of arginine 120 in modulating the pK_a values of these cysteines. We find that arginine 120 plays a significant role in modulating the pK_a of the C-terminal active-site cysteine in the a domain of PDI and plays a role in determining the reactivity of the N-terminal active-site cysteine but not via direct modulation of its pK_a. Mutation of arginine 120 and the corresponding residue, arginine 461, in the a′ domain severely reduces the ability of PDI to catalyze disulfide bond formation and reduction but enhances the ability to catalyze disulfide bond isomerization due to the formation of more stable PDI-substrate mixed disulfides. These results suggest that the modulation of pK_a of the C-terminal active cysteine by the movement of the side chain of these arginine residues into the active-site locales has evolved to allow PDI to efficiently catalyze both oxidation and isomerization reactions.     

Direct Determination of Protonation States of Histidine Residues in a 2 Å Neutron Structure of Deoxy-Human Normal Adult Hemoglobin and Implications for the Bohr Effect

We have investigated the protonation states of histidine residues (potential Bohr groups) in the deoxy form (T state) of human hemoglobin by direct determination of hydrogen (deuterium) positions with the neutron protein crystallography technique. The reversible binding of protons is key to the allosteric regulation of human hemoglobin. The protonation states of 35 of the 38 His residues were directly determined from neutron scattering omit maps, with 3 of the remaining residues being disordered. Protonation states of 5 equivalent His residues—αHis20, αHis50, αHis89, βHis143, and βHis146—differ between the symmetry-related globin subunits. The distal His residues, αHis58 and βHis63, are protonated in the α1β1 heterodimer and are neutral in α2β2. Buried residue αHis103 is found to be protonated in both subunits. These distal and buried residues have the potential to act as Bohr groups. The observed protonation states of His residues are compared to changes in their pK_a values during the transition from the T to the R state and the results provide some new insights into our understanding of the molecular mechanism of the Bohr effect.     

On the amino acids involved in the ATPase site of mitochondrial F1 and the implication of the subunit C.

The data on the pH dependence of the K_m for Mg-ATP and the V_m of the ATPase of pig heart mitochondrial F₁ indicate the presence of two groups of different pK's which modify the enzyme activity. The first pK at pH 9.6 ± 0.2 may be related to the possible presence of arginine and/or tyrosine residues in the ATPase site; the second pK at pH 7.2 ± 0.2 could be due to the presence of a histidine residue in the ATPase site or to the involvement of amino groups in the ATPase site. The inhibition induced by photooxidation in the presence of Rose Bengal is not pH dependent in the pH range corresponding to the pK of histidine. The inhibition induced by diethylpyrocarbonate cannot be reversed by hydroxylamine and the characteristics of this inhibition rather correspond to the reaction of the inhibitor with amino groups. Pyridoxal phosphate also inhibits the ATPase activity of F₁ by reaction with amino groups. The presence of ATP or phosphate partially protects against the inhibition induced by diethylpyrocarbonate or pyridoxal phosphate, which indicates that amino groups may be directly or indirectly involved in the binding of nucleotide and phosphate to F₁. Glutaraldehyde also inhibits the enzyme by reacting with amino groups and inducing a crosslinking of the subunits. The disappearance of subunit C is well correlated with the decrease of ATPase activity, indicating that subunit C is essential in the ATPase activity.     

Effect of methylation on the side‐chain pKa value of arginine

Arginine methylation is important in biological systems. Recent studies link the deregulation of protein arginine methyltransferases with certain cancers. To assess the impact of methylation on interaction with other biomolecules, the pK_a values of methylated arginine variants were determined using NMR data. The pK_a values of monomethylated, symmetrically dimethylated, and asymmetrically dimethylated arginine are similar to the unmodified arginine (14.2 ± 0.4). Although the pK_a value has not been significantly affected by methylation, consequences of methylation include changes in charge distribution and steric effects, suggesting alternative mechanisms for recognition.     

A Combinatorial Histidine Scanning Library Approach to Engineer Highly pH-Dependent Protein Switches

There is growing interest in the development of protein switches, which are proteins whose function, such as binding a target molecule, can be modulated through environmental triggers. Efforts to engineer highly pH sensitive protein–protein interactions typically rely on the rational introduction of ionizable groups in the protein interface. Such experiments are typically time intensive and often sacrifice the protein's affinity at the permissive pH. The underlying thermodynamics of proton‐linkage dictate that the presence of multiple ionizable groups, which undergo a pK_a change on protein binding, are necessary to result in highly pH‐dependent binding. To test this hypothesis, a novel combinatorial histidine library was developed where every possible combination of histidine and wild‐type residue is sampled throughout the interface of a model anti‐RNase A single domain VHH antibody. Antibodies were coselected for high‐affinity binding and pH‐sensitivity using an in vitro, dual‐function selection strategy. The resulting antibodies retained near wild‐type affinity yet became highly sensitive to small decreases in pH, drastically decreasing their binding affinity, due to the incorporation of multiple histidine groups. Several trends were observed, such as histidine “hot‐spots,” which will help enhance the development of pH switch proteins as well as increase our understanding of the role of ionizable residues in protein interfaces. Overall, the combinatorial approach is rapid, general, and robust and should be capable of producing highly pH‐sensitive protein affinity reagents for a number of different applications.     

Molecular Dynamics Simulations Elucidate the Mechanism of Proton Transport in the Glutamate Transporter EAAT3

The uptake of glutamate in nerve synapses is carried out by the excitatory amino acid transporters (EAATs), involving the cotransport of a proton and three Na⁺ ions and the countertransport of a K⁺ ion. In this study, we use an EAAT3 homology model to calculate the pK_a of several titratable residues around the glutamate binding site to locate the proton carrier site involved in the translocation of the substrate. After identifying E374 as the main candidate for carrying the proton, we calculate the protonation state of this residue in different conformations of EAAT3 and with different ligands bound. We find that E374 is protonated in the fully bound state, but removing the Na2 ion and the substrate reduces the pK_a of this residue and favors the release of the proton to solution. Removing the remaining Na⁺ ions again favors the protonation of E374 in both the outward- and inward-facing states, hence the proton is not released in the empty transporter. By calculating the pK_a of E374 with a K⁺ ion bound in three possible sites, we show that binding of the K⁺ ion is necessary for the release of the proton in the inward-facing state. This suggests a mechanism in which a K⁺ ion replaces one of the ligands bound to the transporter, which may explain the faster transport rates of the EAATs compared to its archaeal homologs.     

Investigation of phosphotyrosine recognition by the SH2 domain of the Src kinase.

The binding of tyrosine phosphorylated targets by SH2 domains is required for propagation of many cellular signals in higher eukaryotes; however, the determinants of phosphotyrosine (pTyr) recognition by SH2 domains are not well understood. In order to identify the attributes of pTyr required for high affinity interaction with SH2 domains, the binding of the SH2 domain of the Src kinase (Src SH2 domain) to a dephosphorylated peptide, a phosphoserine-containing peptide, and the amino acid pTyr was studied using titration calorimetry and compared with the binding of a high affinity tyrosyl phosphopeptide. The dephosphorylated peptide and the phosphoserine containing peptide both bind extremely weakly to the Src SH2 domain (DeltaGo (dephosphorylated)=-3.6 kcal/mol, DeltaGo (phosphoserine) >-3.7 kcal/mol); however, the DeltaGo value of pTyr binding is more favorable (-4.7 kcal/mol, or 50 % of the entire binding free energy of a high affinity tyrosyl phosphopeptide). These results indicate that both the phosphate and the tyrosine ring of the pTyr are critical determinants of high affinity binding. Alanine mutagenesis was also used to evaluate the energetic contribution to binding of ten residues located in the pTyr-binding site. Mutation of the strictly conserved Arg betaB5 resulted in a large increase in DeltaGo (DeltaDeltaGo=3.2 kcal/mol) while elimination of the other examined residues each resulted in a significantly smaller (DeltaDeltaGo<1.4 kcal/mol) reduction in affinity, indicating that Arg betaB5 is the single most important determinant of pTyr recognition. However, mutation of Cys betaC3, a residue unique to the Src SH2 domain, surprisingly increased affinity by eightfold (DeltaDeltaGo=-1.1 kcal/mol). Using a double mutant cycle analysis, it was revealed that residues of the pTyr-binding pocket are not coupled to the peptide residues C-terminal to the pTyr. In addition, comparison of each residue's DeltaDeltaGo value upon mutation with that residue's sequence conservation among SH2 domains revealed only a modest correlation between a residue's energetic contribution to pTyr recognition and its conservation throughout evolution. The results of this investigation highlight the importance of a single critical interaction, the buried ionic bond between the phosphate of the pTyr and Arg betaB5 of the SH2 domain, driving the binding of SH2 domains to tyrosine phosphorylated targets.     

Solution structure of the SH2 domain of Grb2 complexed with the Shc-derived phosphotyrosine-containing peptide.

The solution structure of growth factor receptor-bound protein 2 (Grb2) SH2 complexed with a Shc-derived phosphotyrosine (pTyr)-containing peptide was determined by nuclear magnetic resonance (NMR) spectroscopy. The pTyr binding site of Grb2 SH2 was similar to those of other SH2 domains. In contrast, the amino acid residues C-terminal to pTyr did not form an extended structure because of steric hindrance caused by a bulky side-chain of Trp121 (EF1). As a result, the peptide formed a turn-structure on the surface of Grb2 SH2. The asparagine residue at the pTyr+2 position of the Shc-peptide interacted with the main-chain carbonyl groups of Lys109 and Leu120. The present solution structure was similar to the crystal structure reported for Grb2 SH2 complexed with a BCR-Abl-derived phosphotyrosine-containing peptide. Finally, the structure of Grb2 SH2 domain was compared with those of the complexes of Src and phospholipase C-gamma1 with their cognate peptides, showing that the specific conformation of the peptide was required for binding to the SH2 domains.