Hot regions in protein--protein interactions: the organization and contribution of structurally conserved hot spot residues

Structurally conserved residues at protein-protein interfaces correlate with the experimental alanine-scanning hot spots. Here, we investigate the organization of these conserved, computational hot spots and their contribution to the stability of protein associations. We find that computational hot spots are not homogeneously distributed along the protein interfaces; rather they are clustered within locally tightly packed regions. Within the dense clusters, they form a network of interactions and consequently their contributions to the stability of the complex are cooperative; however the contributions of independent clusters are additive. This suggests that the binding free energy is not a simple summation of the single hot spot residue contributions. As expected, around the hot spot residues we observe moderately conserved residues, further highlighting the crucial role of the conserved interactions in the local densely packed environment. The conserved occurrence of these organizations suggests that they are advantageous for protein-protein associations. Interestingly, the total number of hydrogen bonds and salt bridges contributed by hot spots is as expected. Thus, H-bond forming residues may use a "hot spot for water exclusion" mechanism. Since conserved residues are located within highly packed regions, water molecules are easily removed upon binding, strengthening electrostatic contributions of charge-charge interactions. Hence, the picture that emerges is that protein-protein associations are optimized locally, with the clustered, networked, highly packed structurally conserved residues contributing dominantly and cooperatively to the stability of the complex. When addressing the crucial question of "what are the preferred ways of proteins to associate", these findings point toward a critical involvement of hot regions in protein-protein interactions.     

Complementarity of stability patches at the interfaces of protein complexes: Implication for the structural organization of energetic hot spots

A rational design of protein complexes with defined functionalities and of drugs aimed at disrupting protein–protein interactions requires fundamental understanding of the mechanisms underlying the formation of specific protein complexes. Efforts to develop efficient small‐molecule or protein‐based binders often exploit energetic hot spots on protein surfaces, namely, the interfacial residues that provide most of the binding free energy in the complex. The molecular basis underlying the unusually high energy contribution of the hot spots remains obscure, and its elucidation would facilitate the design of interface‐targeted drugs. To study the nature of the energetic hot spots, we analyzed the backbone dynamic properties of contact surfaces in several protein complexes. We demonstrate that, in most complexes, the backbone dynamic landscapes of interacting surfaces form complementary “stability patches,” in which static areas from the opposing surfaces superimpose, and that these areas are predominantly located near the geometric center of the interface. We propose that a diminished enthalpy–entropy compensation effect augments the degree to which residues positioned within the complementary stability patches contribute to complex affinity, thereby giving rise to the energetic hot spots. These findings offer new insights into the nature of energetic hot spots and the role that backbone dynamics play in facilitating intermolecular recognition. Mapping the interfacial stability patches may provide guidance for protein engineering approaches aimed at improving the stability of protein complexes and could facilitate the design of ligands that target complex interfaces.     

Free energy contributions to direct readout of a DNA sequence

The energetic contributions of individual DNA-contacting side chains to specific DNA recognition in the human papillomavirus 16 E2C-DNA complex is small (less than 1.0 kcal mol(-1)), independent of the physical and chemical nature of the interaction, and is strictly additive. The sum of the individual contributions differs 1.0 kcal mol(-1) from the binding energy of the wild-type protein. This difference corresponds to the contribution from the deformability of the DNA, known as "indirect readout." Thus, we can dissect the energetic contribution to DNA binding into 90% direct and 10% indirect readout components. The lack of high energy interactions indicates the absence of "hot spots," such as those found in protein-protein interfaces. These results are compatible with a highly dynamic and "wet" protein-DNA interface, yet highly specific and tight, where individual interactions are constantly being formed and broken.     

Protein hot spots: the islands of stability

Understanding the structural basis of protein-protein interactions (PPIs) may shed light on the organization and functioning of signal transduction and metabolic networks and may assist in structure-based design of ligands (drugs) targeting protein-protein interfaces. The residues at the bimolecular interface, designated as the hot spots, contribute most of the free binding energy of PPI. To date, there is no conclusive atomistic explanation for the unique functional properties of the hot spots. We hypothesized that backbone compliance may play a role in protein-protein recognition and in the mechanism of binding of small-molecule compounds to protein surfaces. We used a steered molecular dynamics simulation to explore the compliance properties of the backbone of surface-exposed residues in several model proteins: interleukin-2, mouse double minute protein 2 and proliferating cell nuclear antigen. We demonstrated that protein surfaces exhibit distinct patterns in which highly immobile residues form defined clusters ("stability patches") alternating with areas of moderate to high mobility. These "stability patches" tend to localize in functionally important regions involved in protein-protein recognition. We propose a mechanism by which the distinct structural organization of the hot spots may contribute to their role in mediating PPI and facilitating binding of structurally diverse small-molecule compounds to protein surfaces.     

Small-world network approach to identify key residues in protein-protein interaction

We show that protein complexes can be represented as small-world networks, exhibiting a relatively small number of highly central amino-acid residues occurring frequently at protein-protein interfaces. We further base our analysis on a set of different biological examples of protein-protein interactions with experimentally validated hot spots, and show that 83% of these predicted highly central residues, which are conserved in sequence alignments and nonexposed to the solvent in the protein complex, correspond to or are in direct contact with an experimentally annotated hot spot. The remaining 17% show a general tendency to be close to an annotated hot spot. On the other hand, although there is no available experimental information on their contribution to the binding free energy, detailed analysis of their properties shows that they are good candidates for being hot spots. Thus, highly central residues have a clear tendency to be located in regions that include hot spots. We also show that some of the central residues in the protein complex interfaces are central in the monomeric structures before dimerization and that possible information relating to hot spots of binding free energy could be obtained from the unbound structures.     

Contributions of amino acid side chains to the kinetics and thermodynamics of the bivalent binding of protein L to Ig kappa light chain

Protein L is a bacterial surface protein with 4-5 immunoglobulin (Ig)-binding domains (B1-B5), each of which appears to have two binding sites for Ig, corresponding to the two edges of its beta-sheet. To verify these sites biochemically and to probe their relative contributions to the protein L-Ig kappa light chain (kappa) interaction, we compared the binding of PLW (the Y47W mutant of the B1 domain) to that of mutants designed to disrupt binding to sites 1 and 2, using gel filtration, BIAcore surface plasmon resonance, fluorescence titration, and solid-phase radioimmunoassays. Gel filtration experiments show that PLW binds kappa both in 1:1 complexes and multivalently, consistent with two binding sites. Covalent dimers of the A20C and V51C mutants of PLW were prepared to eliminate site 1 and site 2 binding, respectively; both the A20C and V51C dimers bind kappa in 1:1 complexes and multivalently, indicating that neither site 1 nor site 2 is solely responsible for kappa binding. The A20R mutant was designed computationally to eliminate site 1 binding while preserving site 2 binding; consistent with this design, the A20R mutant binds kappa in 1:1 complexes but not multivalently. To probe the contributions of amino acid side chains to binding, we prepared 75 point mutants spanning nearly every residue of PLW; BIAcore studies of these mutants revealed two binding-energy "hot spots" consistent with sites 1 and 2. These data indicate that PLW binds kappa at both sites with similar affinities (high nanomolar), with the strongest contributions to the binding energy from Tyr34 (site 2) and Tyr36 (site 1). Compared to other protein-protein complexes, the binding is insensitive to amino acid substitutions at these sites, consistent with the large number of main chain interactions relative to side chain interactions. The strong binding of protein L to Ig kappa light chains of various species may result from the ambidextrous binding of the B1-B5 domains and the unimportance of specific side chain interactions.     

How similar are protein folding and protein binding nuclei? Examination of vibrational motions of energy hot spots and conserved residues

The underlying physico-chemical principles of the interactions between domains in protein folding are similar to those between protein molecules in binding. Here we show that conserved residues and experimental hot spots at intermolecular binding interfaces overlap residues that vibrate with high frequencies. Similarly, conserved residues and hot spots are found in protein cores and are also observed to vibrate with high frequencies. In both cases, these residues contribute significantly to the stability. Hence, these observations validate the proposition that binding and folding are similar processes. In both packing plays a critical role, rationalizing the residue conservation and the experimental alanine scanning hot spots. We further show that high-frequency vibrating residues distinguish between protein binding sites and the remainder of the protein surface.     

Evidence of Conformational Selection Driving the Formation of Ligand Binding Sites in Protein-Protein Interfaces

Many protein-protein interactions (PPIs) are compelling targets for drug discovery, and in a number of cases can be disrupted by small molecules. The main goal of this study is to examine the mechanism of binding site formation in the interface region of proteins that are PPI targets by comparing ligand-free and ligand-bound structures. To avoid any potential bias, we focus on ensembles of ligand-free protein conformations obtained by nuclear magnetic resonance (NMR) techniques and deposited in the Protein Data Bank, rather than on ensembles specifically generated for this study. The measures used for structure comparison are based on detecting binding hot spots, i.e., protein regions that are major contributors to the binding free energy. The main tool of the analysis is computational solvent mapping, which explores the surface of proteins by docking a large number of small “probe” molecules. Although we consider conformational ensembles obtained by NMR techniques, the analysis is independent of the method used for generating the structures. Finding the energetically most important regions, mapping can identify binding site residues using ligand-free models based on NMR data. In addition, the method selects conformations that are similar to some peptide-bound or ligand-bound structure in terms of the properties of the binding site. This agrees with the conformational selection model of molecular recognition, which assumes such pre-existing conformations. The analysis also shows the maximum level of similarity between unbound and bound states that is achieved without any influence from a ligand. Further shift toward the bound structure assumes protein-peptide or protein-ligand interactions, either selecting higher energy conformations that are not part of the NMR ensemble, or leading to induced fit. Thus, forming the sites in protein-protein interfaces that bind peptides and can be targeted by small ligands always includes conformational selection, although other recognition mechanisms may also be involved.     

Protein-protein interactions; coupling of structurally conserved residues and of hot spots across interfaces. Implications for docking

Hot spot residues contribute dominantly to protein-protein interactions. Statistically, conserved residues correlate with hot spots, and their occurrence can distinguish between binding sites and the remainder of the protein surface. The hot spot and conservation analyses have been carried out on one side of the interface. Here, we show that both experimental hot spots and conserved residues tend to couple across two-chain interfaces. Intriguingly, the local packing density around both hot spots and conserved residues is higher than expected. We further observe a correlation between local packing density and experimental deltadeltaG. Favorable conserved pairs include Gly coupled with aromatics, charged and polar residues, as well as aromatic residue coupling. Remarkably, charged residue couples are underrepresented. Overall, protein-protein interactions appear to consist of regions of high and low packing density, with the hot spots organized in the former. The high local packing density in binding interfaces is reminiscent of protein cores.     

AutoMatch: target-binding protein design and enzyme design by automatic pinpointing potential active sites in available protein scaffolds

Proteins perform their functions mainly via active sites, whereas other parts of the proteins comprise the scaffolds, which support the active sites. One strategy for protein functional design is transplanting active sites, such as catalytic sites for enzyme or binding hot spots for protein-protein interactions, onto a new scaffold. AutoMatch is a new program designed for efficiently elucidating suitable scaffolds and potential sites on the scaffolds. Backrub motions are used to treat backbone flexibility during the design process. A step-by-step checking strategy and cluster-representation examination strategy were developed to solve the large combinatorial problem for the matching of active-site conformations. In addition, a grid-based binding energy scoring method was used to filter the solutions. An enzyme design benchmark and a protein-protein interaction design benchmark were built to test the algorithm. AutoMatch could identify the hot spots in the nonbinding protein and rank them within the top five results for 8 of 10 target-binding protein design cases. In addition, among the 15 enzymes tested, AutoMatch can identify the catalytic active sites in the apoprotein and rank them within the top five results for 13 cases. AutoMatch was also tested for screening scaffold library in designing binding proteins targeting influenza hemagglutinin, HIV gp120, and epidermal growth factor receptor kinase, respectively. AutoMatch, and the two test sets, ActApo and ActFree, are available for noncommercial applications at http://mdl.ipc.pku.edu.cn/cgi-bin/down.cgi.     

The modular organization of domain structures: insights into protein-protein binding

Domains are the building blocks of proteins and play a crucial role in protein-protein interactions. Here, we propose a new approach for the analysis and prediction of domain-domain interfaces. Our method, which relies on the representation of domains as residue-interacting networks, finds an optimal decomposition of domain structures into modules. The resulting modules comprise highly cooperative residues, which exhibit few connections with other modules. We found that non-overlapping binding sites in a domain, involved in different domain-domain interactions, are generally contained in different modules. This observation indicates that our modular decomposition is able to separate protein domains into regions with specialized functions. Our results show that modules with high modularity values identify binding site regions, demonstrating the predictive character of modularity. Furthermore, the combination of modularity with other characteristics, such as sequence conservation or surface patches, was found to improve our predictions. In an attempt to give a physical interpretation to the modular architecture of domains, we analyzed in detail six examples of protein domains with available experimental binding data. The modular configuration of the TEM1-beta-lactamase binding site illustrates the energetic independence of hotspots located in different modules and the cooperativity of those sited within the same modules. The energetic and structural cooperativity between intramodular residues is also clearly shown in the example of the chymotrypsin inhibitor, where non-binding site residues have a synergistic effect on binding. Interestingly, the binding site of the T cell receptor beta chain variable domain 2.1 is contained in one module, which includes structurally distant hot regions displaying positive cooperativity. These findings support the idea that modules possess certain functional and energetic independence. A modular organization of binding sites confers robustness and flexibility to the performance of the functional activity, and facilitates the evolution of protein interactions.     

Structural stability of binding sites: Consequences for binding affinity and allosteric effects

During the course of biological function, proteins interact with other proteins, ligands, substrates, inhibitors, etc. These interactions occur at precisely defined locations within the protein but their effects are sometimes propagated to distal regions, triggering highly specific responses. These effects can be used as signals directed to activate or inhibit other sites, modulate interactions with other molecules, and/or establish inter‐molecular communication networks. During the past decade, it has become evident that the energy of stabilization of the protein structure is not evenly distributed throughout the molecule and that, under native conditions, proteins lack global cooperativity and are characterized by the occurrence of multiple independent local unfolding events. From a biological point of view, it is important to assess if this uneven distribution reflects specific functional requirements. For example, are binding sites more likely to be found in well structured regions, unstable regions, or mixed regions? In this article, we have addressed these questions by performing a structure‐based thermodynamic stability analysis of non‐structurally homologous proteins for which high resolution structures of their complexes with specific ligands are available. The results of these studies indicate that for all 16 proteins considered, the binding sites have a dual character and are characterized by the presence of regions with very low structural stability and regions with high stability. In many cases the low stability regions are loops that become stable and cover a significant portion of low molecular weight ligands upon binding. For enzymes, catalytic residues are usually, but not always, located in regions with high structural stability. It is shown that this arrangement provides significant advantages for the optimization of binding affinity of small ligands. In allosteric enzymes, low stability regions in the regulatory site are shown to play a crucial role in the transmission of information to the catalytic site. Proteins 2000;41:63–71. © 2000 Wiley‐Liss, Inc.     

FTMove: A Web Server for Detection and Analysis of Cryptic and Allosteric Binding Sites by Mapping Multiple Protein Structures

Protein mapping distributes many copies of different molecular probes on the surface of a target protein in order to determine binding hot spots, regions that are highly preferable for ligand binding. While mapping of X-ray structures by the FTMap server is inherently static, this limitation can be overcome by the simultaneous analysis of multiple structures of the protein. FTMove is an automated web server that implements this approach. From the input of a target protein, by PDB code, the server identifies all structures of the protein available in the PDB, runs mapping on them, and combines the results to form binding hot spots and binding sites. The user may also upload their own protein structures, bypassing the PDB search for similar structures. Output of the server consists of the consensus binding sites and the individual mapping results for each structure - including the number of probes located in each binding site, for each structure. This level of detail allows the users to investigate how the strength of a binding site relates to the protein conformation, other binding sites, and the presence of ligands or mutations. In addition, the structures are clustered on the basis of their binding properties. The use of FTMove is demonstrated by application to 22 proteins with known allosteric binding sites; the orthosteric and allosteric binding sites were identified in all but one case, and the sites were typically ranked among the top five. The FTMove server is publicly available at https://ftmove.bu.edu.     

The Modular Organization of Domain Structures: Insights into Protein–Protein Binding

Domains are the building blocks of proteins and play a crucial role in protein–protein interactions. Here, we propose a new approach for the analysis and prediction of domain–domain interfaces. Our method, which relies on the representation of domains as residue-interacting networks, finds an optimal decomposition of domain structures into modules. The resulting modules comprise highly cooperative residues, which exhibit few connections with other modules. We found that non-overlapping binding sites in a domain, involved in different domain–domain interactions, are generally contained in different modules. This observation indicates that our modular decomposition is able to separate protein domains into regions with specialized functions. Our results show that modules with high modularity values identify binding site regions, demonstrating the predictive character of modularity. Furthermore, the combination of modularity with other characteristics, such as sequence conservation or surface patches, was found to improve our predictions. In an attempt to give a physical interpretation to the modular architecture of domains, we analyzed in detail six examples of protein domains with available experimental binding data. The modular configuration of the TEM1-β-lactamase binding site illustrates the energetic independence of hotspots located in different modules and the cooperativity of those sited within the same modules. The energetic and structural cooperativity between intramodular residues is also clearly shown in the example of the chymotrypsin inhibitor, where non–binding site residues have a synergistic effect on binding. Interestingly, the binding site of the T cell receptor β chain variable domain 2.1 is contained in one module, which includes structurally distant hot regions displaying positive cooperativity. These findings support the idea that modules possess certain functional and energetic independence. A modular organization of binding sites confers robustness and flexibility to the performance of the functional activity, and facilitates the evolution of protein interactions.     

Protein-protein interactions: hot spots and structurally conserved residues often locate in complemented pockets that pre-organized in the unbound states: implications for docking

Energetic hot spots account for a significant portion of the total binding free energy and correlate with structurally conserved interface residues. Here, we map experimentally determined hot spots and structurally conserved residues to investigate their geometrical organization. Unfilled pockets are pockets that remain unfilled after protein-protein complexation, while complemented pockets are pockets that disappear upon binding, representing tightly fit regions. We find that structurally conserved residues and energetic hot spots are strongly favored to be located in complemented pockets, and are disfavored in unfilled pockets. For the three available protein-protein complexes with complemented pockets where both members of the complex were alanine-scanned, 62% of all hot spots (DeltaDeltaG>2kcal/mol) are within these pockets, and 60% of the residues in the complemented pockets are hot spots. 93% of all red-hot residues (DeltaDeltaG>/=4kcal/mol) either protrude into or are located in complemented pockets. The occurrence of hot spots and conserved residues in complemented pockets highlights the role of local tight packing in protein associations, and rationalizes their energetic contribution and conservation. Complemented pockets and their corresponding protruding residues emerge among the most important geometric features in protein-protein interactions. By screening the solvent, this organization shields backbone hydrogen bonds and charge-charge interactions. Complemented pockets often pre-exist binding. For 18 protein-protein complexes with complemented pockets whose unbound structures are available, in 16 the pockets are identified to pre-exist in the unbound structures. The root-mean-squared deviations of the atoms lining the pockets between the bound and unbound states is as small as 0.9A, suggesting that such pockets constitute features of the populated native state that may be used in docking.     

Protein binding hot spots and the residue-residue pairing preference: a water exclusion perspective

Background  

A protein binding hot spot is a small cluster of residues tightly packed at the center of the interface between two interacting proteins. Though a hot spot constitutes a small fraction of the interface, it is vital to the stability of protein complexes. Recently, there are a series of hypotheses proposed to characterize binding hot spots, including the pioneering O-ring theory, the insightful 'coupling' and 'hot region' principle, and our 'double water exclusion' (DWE) hypothesis. As the perspective changes from the O-ring theory to the DWE hypothesis, we examine the physicochemical properties of the binding hot spots under the new hypothesis and compare with those under the O-ring theory.     

Exosite interactions determine the affinity of factor X for the extrinsic Xase complex

The initiation of coagulation results from the activation of factor X by an enzyme complex (Xase) composed of the trypsin-like serine proteinase, factor VIIa, bound to tissue factor (TF) on phospholipid membranes. We have investigated the basis for the protein substrate specificity of Xase using TF reconstituted into vesicles of phosphatidylcholine, phosphatidylserine, or pure phosphatidylcholine. We show that occupation of the active site of VIIa within Xase by a reversible inhibitor or an alternate peptidyl substrate is sufficient to exclude substrate interactions at the active site but does not alter the affinity of Xase for factor X. This is evident as classical competitive inhibition of peptidyl substrate cleavage but as classical noncompetitive inhibition of factor X activation by active site-directed ligands. This implies that the productive recognition of factor X by Xase arises from a multistep reaction requiring an initial interaction at sites on the enzyme complex distinct from the active site (exosites), followed by active site interactions and bond cleavage. Exosite interactions determine protein substrate affinity, whereas the second binding step influences the maximum catalytic rate for the reaction. We also show that competitive inhibition can be achieved by interfering with exosite binding using factor X derivatives that are expected to have limited or abrogated interactions with the active site of VIIa within Xase. Thus, substrate interactions at exosites, sites removed from the active site of VIIa within the enzyme complex, determine affinity and binding specificity in the productive recognition of factor X by the VIIa-TF complex. This may represent a prevalent strategy through which distinctive protein substrate specificities are achieved by the homologous enzymes of coagulation.     

TrkA receptor "hot spots" for binding of NT-3 as a heterologous ligand

Neurotrophins signal via Trk tyrosine kinase receptors. Nerve growth factor (NGF) is the cognate ligand for TrkA, the brain-derived neurotrophic factor for TrkB, and NT-3 for TrkC. NT-3 also binds TrkA as a lower affinity heterologous ligand. Because neurotrophin-3 (NT-3) interactions with TrkA are biologically relevant, we aimed to define the TrkA "hot spot" functional docking sites of NT-3. The Trk extracellular domain consists of two cysteine-rich subdomains (D1 and D3), flanking a leucine-rich subdomain (D2), and two immunoglobulin-like subdomains IgC1(D4) and IgC2(D5). Previously, the D5 subdomain was defined as the primary ligand-binding site of neurotrophins for their cognate receptors (e.g. NGF binds and activates through TRKA-D5 hot spots). Here binding studies with truncated and chimeric extracellular subdomains show that TRKA-D5 also includes an NT-3 docking and activation hot spot (site 1), and competition studies show that the NGF and NT-3 hot spots on TRKA-D5 are distinct but partially overlapping. In addition, ligand binding studies provide evidence for an NT-3-binding/allosteric site on TRKA-D4 (site 2). NT-3 docking on sites 1 and/or 2 partially blocks NGF binding. Functional survival studies showed that sites 1 and 2 regulate TrkA activation. NT-3 docking on both sites 1 and 2 affords full agonism, which can be additive with NGF activation of Trk. However, NT-3 docking solely on site 1 is partially agonistic but noncompetitively antagonizes NGF binding and activation of Trk. This study demonstrates that Trk signaling is more complex than previously thought because it involves several receptor subdomains and hot spots.     

Protein flexibility in psychrophilic and mesophilic trypsins. Evidence of evolutionary conservation of protein dynamics in trypsin-like serine-proteases

To shed light on the molecular features related to cold-adaptation in serine-proteases, we have carried out molecular dynamics simulations of homologous mesophilic and psychrophilic trypsins, with particular attention to evaluation of intramolecular interactions and flexibility. Psychrophilic trypsins present fewer interdomain interactions and enhanced localized flexibility in regions close to the catalytic site. Notably, these regions fit well with the pattern of protein flexibility previously reported for psychrophilic elastases. Our results indicate that specific sites within the serine-protease fold can be considered hot spots of cold-adaptation and that psychrophilic trypsins and elastases have independently discovered similar molecular strategies to optimize flexibility at low temperatures.