Crystal structure and site 1 binding energetics of human placental lactogen

In primates, placental lactogen (PL) is a pituitary hormone with fundamental roles during pregnancy involving fetal growth, metabolism, and stimulating lactation in the mother. Human placental lactogen (hPL) is highly conserved with human growth hormone (hGH) and both hormones bind to the hPRLR extracellular domain (ECD), the first step in receptor homodimerization, in a Zn2+-dependent manner. A modified surface plasmon resonance method was developed to measure the kinetics for hPL and hGH binding to the hPRLR ECD, with and without Zn2+ and showed that hPL has about a tenfold higher affinity for the hPRLR ECD1 than hGH. The crystal structure of the free state of hPL has been determined to 2.0 A resolution showing the molecule possesses an overall structure similar to other long chain four-helix bundle cytokines. Comparison of the free hPL structure with the 1:1 complex structure of hGH bound to the hPRLR ECD1 suggests that two surface loops undergo conformational changes >10 A upon binding. An 18 residue Ala-scan was used to characterize the binding energy epitope for the site 1 interface of hPL. Individual alanine substitutions at five positions reduced binding affinity by a DeltaDeltaG > or = 3 kcal mol(-1). A comparison of the hPL site 1 epitope with that previously determined for hGH indicates contributions of individual residues track reasonably well between hPL and hGH. In particular, residues involved in the zinc-binding site and Lys172 constitute the principal binding determinants for both hormones. However, several residues that are identical between hPL and hGH contribute quite differently to the binding of the hPRLR ECD1. Additionally, the overall magnitudes of the DeltaDeltaG changes observed from the Ala-scan of hPL were markedly larger than those determined in the comparative scan of hGH to the hPRLR ECD1. The structural and biophysical data presented here show that subtle changes in the structural context of an interaction can lead to significantly different effects at the individual residue level.     

Monte Carlo simulations of the peptide recognition at the consensus binding site of the constant fragment of human immunoglobulin G: the energy landscape analysis of a hot spot at the intermolecular interface

Monte Carlo simulations of molecular recognition at the consensus binding site of the constant fragment (Fc) of human immunoglobulin G (Ig) protein have been performed to analyze structural and thermodynamic aspects of binding for the 13-residue cyclic peptide DCAWHLGELVWCT. The energy landscape analysis of a hot spot at the intermolecular interface using alanine scanning and equilibrium-simulated tempering dynamics with the simplified, knowledge-based energy function has enabled the role of the protein hot spot residues in providing the thermodynamic stability of the native structure to be determined. We have found that hydrophobic interactions between the peptide and the Met-252, Ile-253, His-433, and His-435 protein residues are critical to guarantee the thermodynamic stability of the crystallographic binding mode of the complex. Binding free energy calculations, using a molecular mechanics force field and a solvation energy model, combined with alanine scanning have been conducted to determine the energetic contribution of the protein hot spot residues in binding affinity. The conserved Asn-434, Ser-254, and Tyr-436 protein residues contribute significantly to the binding affinity of the peptide-protein complex, serving as an energetic hot spot at the intermolecular interface. The results suggest that evolutionary conserved hot spot protein residues at the intermolecular interface may be partitioned in fulfilling thermodynamic stability of the native binding mode and contributing to the binding affinity of the complex.     

Structure of a phage display-derived variant of human growth hormone complexed to two copies of the extracellular domain of its receptor: evidence for strong structural coupling between receptor binding sites

The structure of the ternary complex between the phage display- optimized, high-affinity Site 1 variant of human growth hormone (hGH) and two copies of the extracellular domain (ECD) of the hGH receptor (hGHR) has been determined at 2.6 A resolution. There are widespread and significant structural differences compared to the wild-type ternary hGH hGHR complex. The hGH variant (hGH(v)) contains 15 Site 1 mutations and binds>10(2) tighter to the hGHR ECD (hGH(R1)) at Site 1. It is biologically active and specific to hGHR. The hGH(v) Site 1 interface is somewhat smaller and 20% more hydrophobic compared to the wild-type (wt) counterpart. Of the ten hormone-receptor H-bonds in the site, only one is the same as in the wt complex. Additionally, several regions of hGH(v) structure move up to 9A in forming the interface. The contacts between the C-terminal domains of two receptor ECDs (hGH(R1)- hGH(R2)) are conserved; however, the large changes in Site 1 appear to cause global changes in the domains of hGH(R1) that affect the hGH(v)-hGH(R2) interface indirectly. This coupling is manifested by large changes in the conformation of groups participating in the Site 2 interaction and results in a structure for the site that is reorganized extensively. The hGH(v)- hGH(R2) interface contains seven H-bonds, only one of which is found in the wt complex. Several groups on hGH(v) and hGH(R2) undergo conformational changes of up to 8 A. Asp116 of hGH(v) plays a central role in the reorganization of Site 2 by forming two new H-bonds to the side-chains of Trp104(R2) and Trp169(R2), which are the key binding determinants of the receptor. The fact that a different binding solution is possible for Site 2, where there were no mutations or binding selection pressures, indicates that the structural elements found in these molecules possess an inherent functional plasticity that enables them to bind to a wide variety of binding surfaces.     

A binding free energy hot spot in the ankyrin repeat protein GABPbeta mediated protein-protein interaction

The frequently observed ankyrin repeat motif represents a structural scaffold evolved for mediating protein-protein interactions. As such, these repeats modulate a diverse range of cellular functions. We thermodynamically characterized the heterodimeric GA-binding protein (GABP) alphabeta complex and focused specifically on the interaction mediated by the ankyrin repeat domain of the GABPbeta. Our isothermal titration calorimetric analysis of the interaction between the GABP subunits determined an association constant (K(A)) of 6.0 x 10(8) M(-1) and that the association is favorably driven by a significant change in enthalpy (DeltaH) and a minor change in entropy (-TDeltaS). A total of 16 GABPbeta interface residues were chosen for alanine scanning mutagenesis. The calorimetrically measured differences in the free energy of binding were compared to computationally calculated values resulting in a correlation coefficient r = 0.71. We identified three spatially contiguous hydrophobic and aromatic residues that form a binding free energy hot spot (DeltaDeltaG > 2.0 kcal/mol). One residue provides structural support to the hot spot residues. Three non-hot spot residues are intermediate contributors (DeltaDeltaG approximately 1.0 kcal/mol) and create a canopy-like structure over the hot spot residues to possibly occlude solvent and orientate the subunits. The remaining interface residues are located peripherally and have weak contributions. Finally, our mutational analysis revealed a significant entropy-enthalpy compensation for this interaction.     

NMR study of the differential contributions of residues of transforming growth factor alpha to association with its receptor

A heteronuclear NMR study of human transforming growth factor alpha (TGFalpha) in complex with the epidermal growth factor receptor extracellular domain (EGFR-ED) provided an effective method for delineating the relative contributions of the residues of the ligand to its affinity for the receptor. In conjunction with previously obtained mutagenesis data, these results indicate that while a large number of residues are involved in complex formation and make up the binding interface, a small subset contribute most of the binding energy. They also show that while the residues which contribute to receptor binding are localized on one face of the molecule, the specific residues that play the major role in the affinity of TGFalpha in the complex are in two distinct regions of TGFalpha. This suggests that two small functional epitopes each composed of two residues exist within a larger structural epitope presented on the binding face. These results give the most detailed picture to date of the receptor binding determinants and yield further insight into the formation of the ligand-receptor complex.     

Dissecting cooperative and additive binding energetics in the affinity maturation pathway of a protein-protein interface

When two proteins associate they form a molecular interface that is a structural and energetic mosaic. Within such interfaces, individual amino acid residues contribute distinct binding energies to the complex. In combination, these energies are not necessarily additive, and significant positive or negative cooperative effects often exist. The basis of reliable algorithms to predict the specificities and energies of protein-protein interactions depends critically on a quantitative understanding of this cooperativity. We have used a model protein-protein system defined by an affinity maturation pathway, comprising variants of a T cell receptor Vbeta domain that exhibit an overall affinity range of approximately 1500-fold for binding to the superantigen staphylococcal enterotoxin C3, in order to dissect the cooperative and additive energetic contributions of residues within an interface. This molecular interaction has been well characterized previously both structurally, by x-ray crystallographic analysis, and energetically, by scanning alanine mutagenesis. Through analysis of group and individual maturation and reversion mutations using surface plasmon resonance spectroscopy, we have identified energetically important interfacial residues, determined their cooperative and additive energetic properties, and elucidated the kinetic and thermodynamic bases for molecular evolution in this system. The summation of the binding free energy changes associated with the individual mutations that define this affinity maturation pathway is greater than that of the fully matured variant, even though the affinity gap between the end point variants is relatively large. Two mutations in particular, both located in the complementarity determining region 2 loop of the Vbeta domain, exhibit negative cooperativity.     

Dissecting the energetics of an antibody-antigen interface by alanine shaving and molecular grafting.

Alanine-scanning mutagenesis on human growth hormone (hGH) identified 5 primary determinants (Arg 8, Asn 12, Arg 16, Asp 112, and Asp 116) for binding to a monoclonal antibody (MAb 3) (Jin L, Fendly BM, Wells JA, 1992, J Mol Biol 226:851-865). To further analyze the energetic importance of residues surrounding these five, we mutated all neighboring residues to alanine in groups of 7-16 (a procedure we call alanine shaving). Even the most extremely mutated variant, with 16 alanine substitutions, caused less than a 10-fold reduction in binding affinity to MAb3. By comparison, mutating any 1 of the 5 primary determinants to alanine caused a 6- to > 500-fold reduction in affinity. Replacing any of the 4 charged residues (Arg 8, Arg 16, Asp 112, and Asp 116) with a homologous residue (i.e., Arg to Lys or Asp to Glu) caused nearly as large a reduction in affinity as the corresponding alanine replacement. It was possible to graft the 5 primary binding determinants onto a nonbinding homologue of hGH, human placental lactogen (hPL), which has 86% sequence identity to hGH. The grafted hPL mutant bound 10-fold less tightly than hGH to MAb3 but bound as well as hGH when 2 additional framework mutations were introduced. Attempts to recover binding affinity by grafting the MAb3 epitope onto more distantly related scaffolds having a similar 4-helix bundle motif, such as human prolactin (23% sequence identity) or granulocyte colony-stimulating factor, were unsuccessful.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

The thrombin epitope recognizing thrombomodulin is a highly cooperative hot spot in exosite I

The functional epitope of thrombin recognizing thrombomodulin was mapped using Ala-scanning mutagenesis of 54 residues located around the active site, the Na(+) binding loop, the 186-loop, the autolysis loop, exosite I, and exosite II. The epitope for thrombomodulin binding is shaped as a hot spot in exosite I, centered around the buried ion quartet formed by Arg(67), Lys(70), Glu(77), and Glu(80), and capped by the hydrophobic residues Tyr(76) and Ile(82). The hot spot is a much smaller subset of the structural epitope for thrombomodulin binding recently documented by x-ray crystallography. Interestingly, the contribution of each residue of the epitope to the binding free energy shows no correlation with the change in its accessible surface area upon formation of the thrombin-thrombomodulin complex. Furthermore, residues of the epitope are strongly coupled in the recognition of thrombomodulin, as seen for the interaction of human growth hormone and insulin with their receptors. Finally, the Ala substitution of two negatively charged residues in exosite II, Asp(100) and Asp(178), is found unexpectedly to significantly increase thrombomodulin binding.     

Increased site 1 affinity improves biopotency of porcine growth hormone. Evidence against diffusion dependent receptor dimerization

Based on phage display optimization studies with human growth hormone (GH), it is thought that the biopotency of GH cannot be increased. This is proposed to be a result of the affinity of the first receptor for hormone far exceeding that which is required to trap the hormone long enough to allow diffusion of the second receptor to form the ternary complex, which initiates signaling. We report here that despite similar site 1 kinetics to the hGH/hGH receptor interaction, the potency of porcine GH for its receptor can be increased up to 5-fold by substituting hGH residues involved in site 1 binding into pGH. Based on extensive mutations and BIAcore studies, we show that the higher potency and site 1 affinity of hGH for the pGHR is primarily a result of a decreased off-rate associated with residues in the extended loop between helices 1 and 2 that interact with the two key tryptophans Trp104 and Trp169 in the receptor binding hot spot. Our mutagenic analysis has also identified a second determinant (Lys165), which in addition to His169, restricts the ability of non-primate hormones to activate hGH receptor. The increased biopotency of GH that we observe can be explained by a model for GH receptor activation where subunit alignment is critical for effective signaling.     

Computational detection of the binding-site hot spot at the remodeled human growth hormone-receptor interface

A hierarchical computational approach is used to identify the engineered binding-site cavity at the remodeled intermolecular interface between the mutants of human growth hormone (hGH) and the extracellular domain of its receptor (hGHbp). Multiple docking simulations are conducted with the remodeled hGH-hGHbp complex for a panel of potent benzimidazole-containing inhibitors that can restore the binding affinity of the wild-type complex, and for a set of known nonactive small molecules that contain different heterocyclic motifs. Structural clustering of ligand-bound conformations and binding free-energy calculations, using the AMBER force field and a continuum solvation model, can rapidly locate and screen numerous ligand-binding modes on the protein surface and detect the binding-site hot spot at the intermolecular interface. Structural orientation of the benzimidazole motif in the binding-site cavity closely mimics the position of the hot spot residue W104 in the crystal structure of the wild-type complex, which is recognized as an important structural requirement for restoring binding affinity. Despite numerous pockets on the protein surface of the mutant hGH-hGHbp complex, the binding-site cavity presents the energetically favorable hot spot for the benzimidazole-containing inhibitors, whereas for a set of nonactive molecules, the lowest energy ligand conformations do not necessarily bind in the engineered cavity. The results reveal a dominant role of the intermolecular van der Waals interactions in providing favorable ligand-protein energetics in the redesigned interface, in agreement with the experimental and computational alanine scanning of the hGH-hGHbp complex.     

Exploring the capacity of minimalist protein interfaces: interface energetics and affinity maturation to picomolar KD of a single-domain antibody with a flat paratope

A major architectural class in engineered binding proteins ("antibody mimics") involves the presentation of recognition loops off a single-domain scaffold. This class of binding proteins, both natural and synthetic, has a strong tendency to bind a preformed cleft using a convex binding interface (paratope). To explore their capacity to produce high-affinity interfaces with diverse shape and topography, we examined the interface energetics and explored the affinity limit achievable with a flat paratope. We chose a minimalist paratope limited to two loops found in a natural camelid heavy-chain antibody (VHH) that binds to ribonuclease A. Ala scanning of the VHH revealed only three "hot spot" side chains and additional four residues important for supporting backbone-mediated interactions. The small number of critical residues suggested that this is not an optimized paratope. Using selection from synthetic combinatorial libraries, we enhanced its affinity by >100-fold, resulting in variants with Kd as low as 180 pM with no detectable loss of binding specificity. High-resolution crystal structures revealed that the mutations induced only subtle structural changes but extended the network of interactions. This resulted in an expanded hot spot region including four additional residues located at the periphery of the paratope with a concomitant loss of the so-called "O-ring" arrangement of energetically inert residues. These results suggest that this class of simple, single-domain scaffolds is capable of generating high-performance binding interfaces with diverse shape. More generally, they suggest that highly functional interfaces can be designed without closely mimicking natural interfaces.     

Mutational analysis of the complex of human RNase inhibitor and human eosinophil-derived neurotoxin (RNase 2)

RNase inhibitor (RI) binds diverse proteins in the pancreatic RNase superfamily with extremely high avidity. Previous studies showed that tight binding of RNase A and angiogenin (Ang) is achieved primarily through interactions of hot spot residues in the 434-460 C-terminal segment of RI with the enzymatic active site; Asp435 of RI forms key hydrogen bonds with the catalytic lysine in both complexes, whereas the other contacts are largely distinctive. Here we have investigated the structural basis for recognition of a third ligand, eosinophil-derived neurotoxin (EDN), by single-site and multisite mutagenesis. Surprisingly, Ala replacement of Asp435 decreases affinity for EDN only by 14-fold, as compared to the several hundred-fold decreases with RNase A and Ang, and individual mutations of three other hot spot residues-Tyr434, Tyr437, and Ser460-have essentially no effect. Ala substitutions of nine additional residues, selected by examining a computational model of the RI.EDN complex, also have no marked impact. Overall, the losses in affinity for the single-residue variants examined account for only approximately 25% of the free energy of binding for the complex. However, multisite mutagenesis of RI reveals strong superadditivity of mutational effects, indicating that part of this shortfall reflects negative cooperativity. Replacement of Tyr434 together with Asp435 or Tyr437 increases K(i) by 540- and 290-fold, respectively. Thus, the C-terminal region of RI again plays an important role in ligand recognition, although probably smaller than for binding RNase A and Ang. Simultaneous substitutions of three neighboring tryptophans (261, 263, and 318) on RI attenuate affinity even more dramatically (by 4900-fold), indicating that the interactions of this RI region also contribute a considerable amount of the binding energy for the EDN complex. These findings highlight the potential importance of cooperativity in protein-protein interactions and the consequent limitations of single-site mutagenesis for assessing interface energetics.     

Alternative views of functional protein binding epitopes obtained by combinatorial shotgun scanning mutagenesis

Combinatorial shotgun scanning mutagenesis was used to analyze two large, related protein binding sites to assess the specificity and importance of individual side chain contributions to binding affinity. The strategy allowed for cost-effective generation of a plethora of functional data. The ease of the technology promoted comprehensive investigations, in which the classic alanine-scanning approach was expanded with two additional strategies, serine- and homolog-scanning. Binding of human growth hormone (hGH) to the hGH receptor served as the model system. The entire high affinity receptor-binding sites (site 1) of wild-type hGH (hGH_wt) and of an affinity-improved variant (hGH_v) were investigated and the results were compared. The contributions that 35 residue positions make to binding were assessed on each hormone molecule by both serine- and homolog-scanning. The hormone molecules were displayed on the surfaces of bacteriophage, and the 35 positions were randomized simultaneously to allow equal starting frequencies of the wild-type residue and either serine or a homologous mutation in separate libraries. Functional selections for binding to the hGH receptor shifted the relative wild-type/mutant frequencies at each position to an extent characteristic of the functional importance of the side chain. Functional epitope maps were created and compared to previous maps obtained by alanine-scanning. Comparisons between the different scans provide insights into the affinity maturation process that produced hGH_v. The serine and homolog-scanning results expand upon and complement the alanine-scanning results and provide additional data on the robustness of the high affinity receptor-binding site of hGH.     

Sequence determinants of E2-E6AP binding affinity and specificity

The conjugation of ubiquitin to substrates requires a series of enzymatic reactions consisting of an activating enzyme (E1), conjugating enzymes (E2) and ligases (E3). Tagging the appropriate substrate with ubiquitin is achieved by specific E2-E3 and E3-substrate interactions. E6AP, a member of the HECT family of E3s, has been previously shown to bind and function with the E2s UbcH7 and UbcH8. To decipher the sequence determinants of this specificity we have developed a quantitative E2-E3 binding assay based on fluorescence polarization and used this assay to measure the affinity of wild-type and mutant E2-E6AP interactions. Alanine scanning of the E6AP-UbcH7 binding interface identified four side-chains on UbcH7 and six side-chains on E6AP that contribute more than 1 kcal/mol to the binding free energy. Two of the hot spot residues from UbcH7 (K96 and K100) are conserved in UbcH8 but vary across other E2s. To determine if these are key specificity determining residues, we attempted to induce a tighter association between the E2 UbcH5b and E6AP by mutating the corresponding positions in UbcH5b to lysine residues. Surprisingly, the mutations had little effect, but rather a mutation at UbcH7 position 4, which is not at a hot spot on the UbcH7-E6AP interface, significantly strengthened UbcH5bs affinity for E6AP. This result indicates that E2-E3 binding specificities are a function of both favorable interactions that promote binding, and unfavorable interactions that prevent binding with unwanted partners.     

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.     

Shotgun alanine scanning shows that growth hormone can bind productively to its receptor through a drastically minimized interface

The high affinity binding site (Site1) of the human growth hormone (hGH) binds to its cognate receptor (hGHR) via a concave surface patch containing about 35 residues. Using 167 sequences from a shotgun alanine scanning analysis of Site1, we have determined that over half of these residues can be simultaneously changed to an alanine or a non-isosteric amino acid while still retaining a high affinity interaction. Among these hGH variants the distribution of the mutation is highly variable throughout the interface, although helix 4 is more conserved than the other binding elements. Kinetic and thermodynamic analyses were performed on 11 representative hGH Site1 variants that contained 14-20 mutations. Generally, the tightest binding variants showed similar associated rate constants (k(on)) as the wild-type (wt) hormone, indicating that their binding proceeds through a similar transition state intermediate. However, calorimetric analyses indicate very different thermodynamic partitioning: wt-hGH binding exhibits favorable enthalpy and entropy contributions, whereas the variants display highly favorable enthalpy and highly unfavorable entropy contributions. The heat capacities (DeltaCp) on binding measured for wt-hGH and its variants are significantly larger than normally seen for typical protein-protein interactions, suggesting large conformational or solvation effects. The multiple Site1 mutations are shown to indirectly affect binding of the second receptor at Site2 through an allosteric mechanism. We show that the stability of the ternary hormone-receptor complex reflects the affinity of the Site2 binding and is surprisingly exempt from changes in Site1 affinity, directly demonstrating that dissociation of the active signaling complex is a stepwise process.     

Analysis of serine proteinase-inhibitor interaction by alanine shaving

We analyzed the energetic importance of residues surrounding the hot spot (the P(1) position) of bovine pancreatic trypsin inhibitor (BPTI) in interaction with two proteinases, trypsin and chymotrypsin, by a procedure called molecular shaving. One to eight residues of the structural epitope, composed of two extended and exposed loops, were mutated to alanine(s). Although truncation of the side chains of residues surrounding the P(1) position to methyl groups caused a decrease in Delta G(den) values up to 6.4 kcal mole(-1), it did not influence the overall conformation of the inhibitor. We found that the replacement of up to six residues with alanines was fully additive at the level of protein stability. To analyze the influence of the structural epitope on the association energy, we determined association constants for BPTI variants and both enzymes and applied the additivity analysis. Shaving of two binding loops led to a progressive drop in the association energy, more pronounced for trypsin (decrease up to 9.6 kcal mole(-1)) than chymotrypsin (decrease up to 3.5 kcal mole(-1)). In the case of extensively mutated variants interacting with chymotrypsin, the association energies agreed very well with the values calculated from single mutational effects. However, when P(1)-neighboring residues were shaved to alanine(s), their contribution to the association energy was not fully removed because of the presence of methyl groups and main chain-main chain intermolecular hydrogen bonds. Moreover, the hot spot had a different contribution to the complex stability in the fully shaved BPTI variant compared with the wild type, which was caused by perturbations of the P(1)-S(1) electrostatic interaction.     

Analysis of the interactions of human ribonuclease inhibitor with angiogenin and ribonuclease A by mutagenesis: importance of inhibitor residues inside versus outside the C-terminal "hot spot"

Ribonuclease inhibitor (RI) binds diverse mammalian RNases with extraordinary avidity. Here, we have investigated the structural basis for this tight binding and broad specificity by mutational analysis of the complexes of RI with angiogenin (Ang) and RNase A (K(D)=0.5 fM and 43 fM, respectively). Both crystal structures are known; the interfaces are large, and the ligands dock similarly, although few of the specific interactions formed are analogous. Our previous mutagenesis studies focused primarily on one contact region, containing RI 434-438 and the enzymatic active site. Many single-residue replacements produced extensive losses of binding energy (2.3-5.9 kcal/mol), suggesting that this region constitutes a "hot spot" in both cases. We have now explored the roles of most of the remaining RI residues that interact with Ang and/or RNase A. One major cluster in each complex lies in a Trp-rich area of RI, containing Trp261, Trp263, Trp318, and Trp375. Although the energy losses from individual replacements in this portion of the Ang complex were small-to-moderate (0-1.5 kcal/mol), the changes from multiple substitutions were much greater than additive, and the binding energy provided by this region is estimated to be approximately 6 kcal/mol (30 % of total). Effects of replacing combinations of hot spot components had also been found to be superadditive, and this negative cooperativity is now shown to extend to the neighboring contact residue RI Ser460. The overall contribution of the hot spot, taking superadditivity into account, is then approximately 14-15 kcal/mol. The hot spot and Trp-rich regions, although spatially well separated, are themselves functionally linked. No other parts of the RI-Ang interface appear to be energetically important. Binding of RNase A is more sensitive to substitutions throughout the interface, with free energy losses>/=1 kcal/mol produced by nearly all replacements examined, so that the sum of losses greatly exceeds the binding energy of the complex. This discrepancy can be explained, in part, by positive cooperativity, as evident from the subadditive effects observed when combinations of residues in either the hot spot or Trp-rich region are replaced. These findings suggest that the binding energy may be more widely distributed in the RNase A complex than in the Ang complex.     

Magnitude of the hydrophobic effect at central versus peripheral sites in protein-protein interfaces

Hydrophobic interactions are essential for stabilizing protein-protein complexes, whose interfaces generally consist of a central cluster of hot spot residues surrounded by less important peripheral residues. According to the O-ring hypothesis, a condition for high affinity binding is solvent exclusion from interacting residues. This hypothesis predicts that the hydrophobicity at the center is significantly greater than at the periphery, which we estimated at 21 cal mol(-1) A(-2). To measure the hydrophobicity at the center, structures of an antigen-antibody complex where a buried phenylalanine was replaced by smaller hydrophobic residues were determined. By correlating structural changes with binding free energies, we estimate the hydrophobicity at this central site to be 46 cal mol(-1) A(-2), twice that at the periphery. This context dependence of the hydrophobic effect explains the clustering of hot spots at interface centers and has implications for hot spot prediction and the design of small molecule inhibitors.