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.     

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.     

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.     

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.     

Protein subunit interfaces: a statistical analysis of hot spots in Sm proteins

The distinguishing property of Sm protein associations is very high stability. In order to understand this property, we analyzed the interfaces and compared the properties of Sm protein interfaces with those of a test set, the Binding Interface Database (BID). The comparison revealed that the main differences between the interfaces of Sm proteins and those of the BID set are the content of charged residues, the coordination numbers of the residues, knowledge-based pair potentials, and the conservation scores of hot spots. In Sm proteins, the interfaces have more hydrophobic and fewer charged residues than the surfaces, which is also the case for the BID test set and other proteins. However, in the interfaces, the content of charged residues in Sm proteins (26%) is substantially larger than that in the BID set (22%). Hot spots are residues that make up a small fraction of the interfaces, but they contribute most of the binding energy. These residues are critical to protein–protein interactions. Analyses of knowledge-based pair potentials of hot spot and non-hot spot residues in Sm proteins show that they are significantly different; their mean values are 31.5 and 11.3, respectively. In the BID set, this difference is smaller; in this case, the mean values for hot spot and non-hot spot residues are 20.7 and 12.4, respectively. Hence, the pair potentials of hot spots differ significantly for the Sm and BID data sets. In the interfaces of Sm proteins, the amino acids are tightly packed, and the coordination numbers are larger in Sm proteins than in the BID set for both hot spots and non-hot spots. At the same time, the coordination numbers are higher for hot spots; the average coordination number of the hot spot residues in Sm proteins is 7.7, while it is 6.1 for the non-hot spot residues. The difference in the calculated average conservation score for hot spots and non-hot spots in Sm proteins is significantly larger than it is in the BID set. In Sm proteins, the average conservation score for the hot spots is 7.4. Hot spots are surrounded by residues that are moderately conserved (5.9). The average conservation score for the other interface residues is 5.6. The conservation scores in the BID set do not show a significant distinction between hot and non-hot spots: the mean values for hot and non-hot spot residues are 5.5 and 5.2, respectively. These data show that structurally conserved residues and hot spots are significantly correlated in Sm proteins.     

Protein-DNA interactions: structural, thermodynamic and clustering patterns of conserved residues in DNA-binding proteins

Amino acid residues, which play important roles in protein function, are often conserved. Here, we analyze thermodynamic and structural data of protein-DNA interactions to explore a relationship between free energy, sequence conservation and structural cooperativity. We observe that the most stabilizing residues or putative hotspots are those which occur as clusters of conserved residues. The higher packing density of the clusters and available experimental thermodynamic data of mutations suggest cooperativity between conserved residues in the clusters. Conserved singlets contribute to the stability of protein-DNA complexes to a lesser extent. We also analyze structural features of conserved residues and their clusters and examine their role in identifying DNA-binding sites. We show that about half of the observed conserved residue clusters are in the interface with the DNA, which could be identified from their amino acid composition; whereas the remaining clusters are at the protein-protein or protein-ligand interface, or embedded in the structural scaffolds. In protein-protein interfaces, conserved residues are highly correlated with experimental residue hotspots, contributing dominantly and often cooperatively to the stability of protein-protein complexes. Overall, the conservation patterns of the stabilizing residues in DNA-binding proteins also highlight the significance of clustering as compared to single residue conservation.     

Charged residues at protein interaction interfaces: unexpected conservation and orchestrated divergence

Protein-protein interactions play an essential role in the functioning of cell. The importance of charged residues and their diverse role in protein-protein interactions have been well studied using experimental and computational methods. Often, charged residues located in protein interaction interfaces are conserved across the families of homologous proteins and protein complexes. However, on a large scale, it has been recently shown that charged residues are significantly less conserved than other residue types in protein interaction interfaces. The goal of this work is to understand the role of charged residues in the protein interaction interfaces through their conservation patterns. Here, we propose a simple approach where the structural conservation of the charged residue pairs is analyzed among the pairs of homologous binary complexes. Specifically, we determine a large set of homologous interactions using an interaction interface similarity measure and catalog the basic types of conservation patterns among the charged residue pairs. We find an unexpected conservation pattern, which we call the correlated reappearance, occurring among the pairs of homologous interfaces more frequently than the fully conserved pairs of charged residues. Furthermore, the analysis of the conservation patterns across different superkingdoms as well as structural classes of proteins has revealed that the correlated reappearance of charged residues is by far the most prevalent conservation pattern, often occurring more frequently than the unconserved charged residues. We discuss a possible role that the new conservation pattern may play in the long-range electrostatic steering effect.     

蛋白质-蛋白质界面热点残基预测及其在线工具

蛋白质-蛋白质结合热点是界面中对结合自由能有着显著贡献的一小簇残基。捕捉和揭示这类热点残基可以加深对蛋白质间相互作用机制的理解,为蛋白质工程和药物设计提供指导。但实验技术费时费力且代价昂贵。计算工具可用于辅助和补充实验上的尝试。该文较详细、系统地介绍了蛋白质界面热点的特性、计算预测的策略与技术,并应用实例进一步说明这些方法学的特征;还介绍了界面热点的数据库和一些主要的在线预测工具,旨在为设计、挑选和应用这类工具解决特定问题的研究人员提供指南。     

Z-score biological significance of binding hot spots of protein interfaces by using crystal packing as the reference state

Characterization of binding hot spots of protein interfaces is a fundamental study in molecular biology. Many computational methods have been proposed to identify binding hot spots. However, there are few studies to assess the biological significance of binding hot spots. We introduce the notion of biological significance of a contact residue for capturing the probability of the residue occurring in or contributing to protein binding interfaces. We take a statistical Z-score approach to the assessment of the biological significance. The method has three main steps. First, the potential score of a residue is defined by using a knowledge-based potential function with relative accessible surface area calculations. A null distribution of this potential score is then generated from artifact crystal packing contacts. Finally, the Z-score significance of a contact residue with a specific potential score is determined according to this null distribution. We hypothesize that residues at binding hot spots have big absolute values of Z-score as they contribute greatly to binding free energy. Thus, we propose to use Z-score to predict whether a contact residue is a hot spot residue. Comparison with previously reported methods on two benchmark datasets shows that this Z-score method is mostly superior to earlier methods. This article is part of a Special Issue entitled: Computational Methods for Protein Interaction and Structural Prediction.     

Identification of hot spot residues at protein-protein interface

It is known that binding free energy of protein-protein interaction is mainly contributed by hot spot (high energy) interface residues. Here, we investigate the characteristics of hot spots by examining inter-atomic sidechain-sidechain interactions using a dataset of 296 alanine-mutated interface residues. Results show that hot spots participate in strong and energetically favorable sidechain-sidechain interactions. Subsequently, we describe a novel, yet simple 'hot spot' prediction model with an accuracy that is similar to many available approaches. The model is also shown to efficiently distinguish specific protein-protein interactions from non-specific interactions.     

Conservation of polar residues as hot spots at protein interfaces

A number of studies have addressed the question of which are the critical residues at protein-binding sites. These studies examined either a single or a few protein-protein interfaces. The most extensive study to date has been an analysis of alanine-scanning mutagenesis. However, although the total number of mutations was large, the number of protein interfaces was small, with some of the interfaces closely related. Here we show that although overall binding sites are hydrophobic, they are studded with specific, conserved polar residues at specific locations, possibly serving as energy "hot spots." Our results confirm and generalize the alanine-scanning data analysis, despite its limited size. Previously Trp, Arg, and Tyr were shown to constitute energetic hot spots. These were rationalized by their polar interactions and by their surrounding rings of hydrophobic residues. However, there was no compelling reason as to why specifically these residues were conserved. Here we show that other polar residues are similarly conserved. These conserved residues have been detected consistently in all interface families that we have examined. Our results are based on an extensive examination of residues which are in contact across protein interfaces. We utilize all clustered interface families with at least five members and with sequence similarity between the members in the range of 20-90%. There are 11 such clustered interface families, comprising a total of 97 crystal structures. Our three-dimensional superpositioning analysis of the occurrences of matched residues in each of the families identifies conserved residues at spatially similar environments. Additionally, in enzyme inhibitors, we observe that residues are more conserved at the interfaces than at other locations. On the other hand, antibody-protein interfaces have similar surface conservation as compared to their corresponding linear sequence alignment, consistent with the suggestion that evolution has optimized protein interfaces for function.     

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.     

Artificial intelligence based methods for hot spot prediction

Proteins interact through their interfaces to fulfill essential functions in the cell. They bind to their partners in a highly specific manner and form complexes that have a profound effect on understanding the biological pathways they are involved in. Any abnormal interactions may cause diseases. Therefore, the identification of small molecules which modulate protein interactions through their interfaces has high therapeutic potential. However, discovering such molecules is challenging. Most protein–protein binding affinity is attributed to a small set of amino acids found in protein interfaces known as hot spots. Recent studies demonstrate that drug-like small molecules specifically may bind to hot spots. Therefore, hot spot prediction is crucial. As experimental data accumulates, artificial intelligence begins to be used for computational hot spot prediction. First, we review machine learning and deep learning for computational hot spot prediction and then explain the significance of hot spots toward drug design.     

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.     

Restricted mobility of conserved residues in protein-protein interfaces in molecular simulations

Conserved residues in protein-protein interfaces correlate with residue hot-spots. To obtain insight into their roles, we have studied their mobility. We have performed 39 explicit solvent simulations of 15 complexes and their monomers, with the interfaces varying in size, shape, and function. The dynamic behavior of conserved residues in unbound monomers illustrates significantly lower flexibility as compared to their environment, suggesting that already before binding they are constrained in a boundlike configuration. To understand this behavior, we have analyzed the inter- and intrachain hydrogen-bond residence-time in the interfaces. We find that conserved residues are not involved significantly in hydrogen bonds across the interface as compared to nonconserved. However, the monomer simulations reveal that conserved residues contribute dominantly to hydrogen-bond formation before binding. Packing of conserved residues across the trajectories is significantly higher before and after the binding, rationalizing their lower mobility. Backbone torsional angle distributions show that conserved residues assume restricted regions of space and the most visited conformations in the bound and unbound trajectories are similar, suggesting that conserved residues are preorganized. Combined with previous studies, we conclude that conserved residues, hot spots, anchor, and interface-buried residues may be similar residues, fulfilling similar roles.     

Computational Mapping of Anchoring Spots on Protein Surfaces

Protein-protein and protein-peptide interactions are often controlled by few strong contacts that involve hot spot residues. Computational detection of such contacts, termed here anchoring spots, is important for understanding recognition processes and for predicting interactions; it is an essential step in designing interaction interfaces and therapeutic agents. We describe ANCHORSMAP, an algorithm for computational mapping of amino acid side chains on protein surfaces. The algorithm consists of two stages: A geometry based stage (LSMdet), in which sub-pockets adequate for binding single side chains are detected and amino acid probes are scattered near them, and an energy based stage in which optimal positions of the probes are determined through repeated energy minimization and clustering of nearby poses and their ΔG are calculated. ANCHORSMAP employs a new function for ΔG calculations, which is specifically designed for the context of protein-protein recognition by introducing a correction in the electrostatic energy term that compensates for the dielectric shielding exerted by a hypothetical protein bound to the probe.The algorithm successfully detects known anchoring sites and accurately positions the probes. The calculated ΔG rank high the correct anchoring spots in maps produced for unbound proteins. We find that Arg, Trp, Glu and Tyr, which are favorite hot spot residues, are also more selective of their binding environment. The usefulness of anchoring spots mapping is demonstrated by detecting the binding surfaces in the protein-protein complex barnase/barstar and the protein-peptide complex kinase/PKI, and by identifying phenylalanine anchoring sites on the surface of the nuclear transporter NTF2, C-terminus anchors on PDZ domains and phenol anchors on thermolysin. Finally, we discuss the role of anchoring spots in molecular recognition processes.     

Evolvability of yeast protein-protein interaction interfaces

The functional importance of protein-protein interactions indicates that there should be strong evolutionary constraint on their interaction interfaces. However, binding interfaces are frequently affected by amino acid replacements. Change due to coevolution within interfaces can contribute to variability but is not ubiquitous. An alternative explanation for the ability of surfaces to accept replacements may be that many residues can be changed without affecting the interaction. Candidates for these types of residues are those that make interchain interaction only through the protein main chain, β-carbon, or associated hydrogen atoms. Since almost all residues have these atoms, we hypothesize that this subset of interface residues may be more easily substituted than those that make interactions through other atoms. We term such interactions "residue type independent." Investigating this hypothesis, we find that nearly a quarter of residues in protein interaction interfaces make exclusively interchain residue-type-independent contacts. These residues are less structurally constrained and less conserved than residues making residue-type-specific interactions. We propose that residue-type-independent interactions allow substitutions in binding interfaces while the specificity of binding is maintained.     

Statistical analysis of predominantly transient protein-protein interfaces

A non-redundant set of 170 protein-protein interfaces of known structure was statistically analyzed for residue and secondary-structure compositions, pairing preferences and side-chain-backbone interaction frequencies. By focussing mainly on transient protein-protein interfaces, the results underline previous findings for protein-protein interfaces but also show some new interesting aspects of transient interfaces. The residue compositions at interfaces found in this study correlate well with the results of other studies. On average, contacts between pairs of hydrophobic and polar residues were unfavorable, and the charged residues tended to pair subject to charge complementarity. Secondary structure composition analysis shows that neither helices nor beta-sheets are dominantly populated at interfaces. Analyzing the pairing preferences of the secondary structure elements revealed a higher affinity within the same elements and alludes to tight packings. In addition, the results for the side-chain and backbone interaction frequencies, which were measured under more stringent conditions, showed a high occurrence of side-chain-backbone interactions. Taking a closer look at the helix and beta-sheet binding frequencies for a given side-chain and backbone interaction underlined the relevance of tight packings. The polarity of interfaces increased with decreasing interface size. These types of information may be useful for scoring complexes in protein-protein docking studies or for prediction of protein-protein interfaces from the sequences alone.     

Protein domain interfaces: characterization and comparison with oligomeric protein interfaces

The physical and chemical properties of domain-domain interactions have been analysed in two-domain proteins selected from the protein classification, CATH. The two-domain structures were divided into those derived from (i) monomeric proteins, or (ii) oligomeric or complexed proteins. The size, polarity, hydrogen bonding and packing of the intra-chain domain interface were calculated for both sets of two-domain structures. The results were compared with inter-chain interface parameters from permanent and non-obligate protein-protein complexes. In general, the intra-chain domain and inter-chain interfaces were remarkably similar. Many of the intra-chain interface properties are intermediate between those calculated for permanent and non-obligate inter-chain complexes. Residue interface propensities were also found to be very similar, with hydrophobic residues playing a major role, together with positively charged arginine residues. In addition, the residue composition of the domain interfaces were found to be more comparable with domain surfaces than domain cores. The implications of these results for domain swapping and protein folding are discussed.     

Studies of protein-protein interfaces: a statistical analysis of the hydrophobic effect.

Data sets of 362 structurally nonredundant protein-protein interfaces and of 57 symmetry-related oligomeric interfaces have been used to explore whether the hydrophobic effect that guides protein folding is also the main driving force for protein-protein associations. The buried nonpolar surface area has been used to measure the hydrophobic effect. Our analysis indicates that, although the hydrophobic effect plays a dominant role in protein-protein binding, it is not as strong as that observed in the interior of protein monomers. Comparison of interiors of the monomers with those of the interfaces reveals that, in general, the hydrophobic amino acids are more frequent in the interior of the monomers than in the interior of the protein-protein interfaces. On the other hand, a higher proportion of charged and polar residues are buried at the interfaces, suggesting that hydrogen bonds and ion pairs contribute more to the stability of protein binding than to that of protein folding. Moreover, comparison of the interior of the interfaces to protein surfaces indicates that the interfaces are poorer in polar/charged than the surfaces and are richer in hydrophobic residues. The interior of the interfaces appears to constitute a compromise between the stabilization contributed by the hydrophobic effect on the one hand and avoiding patches on the protein surfaces that are too hydrophobic on the other. Such patches would be unfavorable for the unassociated monomers in solution. We conclude that, although the types of interactions are similar between protein-protein interfaces and single-chain proteins overall, the contribution of the hydrophobic effect to protein-protein associations is not as strong as to protein folding. This implies that packing patterns and interatom, or interresidue, pairwise potential functions, derived from monomers, are not ideally suited to predicting and assessing ligand associations or design. These would perform adequately only in cases where the hydrophobic effect at the binding site is substantial.