Hot spots--a review of the protein-protein interface determinant amino-acid residues

Proteins tendency to bind to one another in a highly specific manner forming stable complexes is fundamental to all biological processes. A better understanding of complex formation has many practical applications, which include the rational design of new therapeutic agents, and the analysis of metabolic and signal transduction networks. Alanine-scanning mutagenesis made possible the detection of the functional epitopes, and demonstrated that most of the protein-protein binding energy is related only to a group of few amino acids at intermolecular protein interfaces: the hot spots. The scope of this review is to summarize all the available information regarding hot spots for a better atomic understanding of their structure and function. The ultimate objective is to improve the rational design of complexes of high affinity and specificity as well as that of small molecules, which can mimic the functional epitopes of the proteic complexes.     

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.     

Detailed microscopic study of the full zipA:FtsZ interface

Protein-protein interaction networks are very important for a wide range of biological processes. Crystallographic structures and mutational studies have generated a large number of information that allowed the discovery of energetically important determinants of specificity at intermolecular protein interfaces and the understanding of the structural and energetic characteristics of the binding hot spots. In this study we have used the improved MMPB/SA (molecular mechanics/Poisson-Boltzmann surface area) approach that combining molecular mechanics and continuum solvent permits to calculate the free energy differences upon alanine mutation. For a better understanding of the binding determinants of the complex formed between the FtsZ fragment and ZipA we extended the alanine scanning mutagenesis study to all interfacial residues of this complex. As a result, we present new mutations that allowed the discovery of residues for which the binding free energy differences upon alanine mutation are higher than 2.0 kcal/mol. We also observed the formation of a hydrophobic pocket with a high warm spot spatial complementarity between FtsZ and ZipA. Small molecules could be designed to bind to these amino acid residues hindering the binding of FtsZ to ZipA. Hence, these mutational data can be used to design new drugs to control more efficiently bacterial infections.     

Wrapping mimicking in drug-like small molecules disruptive of protein-protein interfaces

The discovery of small-molecule drugs aimed at disrupting protein-protein associations is expected to lead to promising therapeutic strategies. The small molecule binds to the target protein thus replacing its natural protein partner. Noteworthy, structural analysis of complexes between successful disruptive small molecules and their target proteins has suggested the possibility that such ligands might somehow mimic the binding behavior of the protein they replace. In these cases, the molecules show a spatial and "chemical" (i.e., hydrophobicity) similarity with the residues of the partner protein involved in the protein-protein complex interface. However, other disruptive small molecules do not seem to show such spatial and chemical correspondence with the replaced protein. In turn, recent progress in the understanding of protein-protein interactions and binding hot spots has revealed the main role of intermolecular wrapping interactions: three-body cooperative correlations in which nonpolar groups in the partner protein promote dehydration of a two-body electrostatic interaction of the other protein. Hence, in the present work, we study some successful complexes between already discovered small disruptive drug-like molecules and their target proteins already reported in the literature and we compare them with the complexes between such proteins and their natural protein partners. Our results show that the small molecules do in fact mimic to a great extent the wrapping behavior of the protein they replace. Thus, by revealing the replacement the small molecule performs of relevant wrapping interactions, we convey precise physical meaning to the mimicking concept, a knowledge that might be exploited in future drug-design endeavors.     

Understanding the molecular basis of MK2-p38α signaling complex assembly: insights into protein-protein interaction by molecular dynamics and free energy studies

The formation of a p38 MAPK and MAPK-activated protein kinase 2 (MK2) signaling complex is physiologically relevant to cellular responses such as the proinflammatory cytokine production. The interaction between p38α isoform and MK2 is of great importance for this signaling. In this study, molecular dynamics simulation and binding free energy calculation were performed on the MK2-p38α signaling complex to investigate the protein-protein interaction between the two proteins. Dynamic domain motion analyses were performed to analyze the conformational changes between the unbound and bound states of proteins during the interaction. The activation loop, αF-I helices, and loops among α helices in the C-lobe of MK2 are found to be highly flexible and exhibit significant changes upon p38α binding. The results also show that after the binding of p38α, the N- and C-terminal domains of MK2 display an opening and twisting motion centered on the activation loop. The molecular mechanics Poisson-Boltzmann and generalized-Born surface area (MM-PB/GBSA) methods were used to calculate binding free energies between MK2 and p38α. The analysis of the components of binding free energy calculation indicates that the van der Waals interaction and the nonpolar solvation energy provide the driving force for the binding process, while the electrostatic interaction contributes critically to the specificity, rather than to MK2-p38α binding affinity. The contribution of each residue at the interaction interface to the binding affinity of MK2 with p38α was also analyzed by free energy decomposition. Several important residues responsible for the protein-protein interaction were also identified.     

Poisson-Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies

The salt dependence of the binding free energy of five protein-protein hetero-dimers and two homo-dimers/tetramers was calculated from numerical solutions to the Poisson-Boltzmann equation. Overall, the agreement with experimental values is very good. In all cases except one involving the highly charged lactoglobulin homo-dimer, increasing the salt concentration is found both experimentally and theoretically to decrease the binding affinity. To clarify the source of salt effects, the salt-dependent free energy of binding is partitioned into screening terms and to self-energy terms that involve the interaction of the charge distribution of a monomer with its own ion atmosphere. In six of the seven complexes studied, screening makes the largest contribution but self-energy effects can also be significant. The calculated salt effects are found to be insensitive to force-field parameters and to the internal dielectric constant assigned to the monomers. Nonlinearities due to high charge densities, which are extremely important in the binding of proteins to negatively charged membrane surfaces and to nucleic acids, make much smaller contributions to the protein-protein complexes studied here, with the exception of highly charged lactoglobulin dimers. Our results indicate that the Poisson-Boltzmann equation captures much of the physical basis of the nonspecific salt dependence of protein-protein complexation.     

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.     

The role of hydrogen bonding via interfacial water molecules in antigen-antibody complexation. The HyHEL-10-HEL interaction

To study the role of hydrogen bonding via interfacial water molecules in protein-protein interactions, we examined the interaction between hen egg white lysozyme (HEL) and its HyHEL-10 variable domain fragment (Fv) antibody. We constructed three antibody mutants (l-Y50F, l-S91A, and l-S93A) and investigated the interactions between the mutant Fvs and HEL. Isothermal titration calorimetry indicated that the mutations significantly decreased the negative enthalpy change (8-25 kJ mol(-1)), despite some offset by a favorable entropy change. X-ray crystallography demonstrated that the complexes had nearly identical structures, including the positions of the interfacial water molecules. Taken together, the isothermal titration calorimetric and x-ray crystallographic results indicate that hydrogen bonding via interfacial water enthalpically contributes to the Fv-HEL interaction despite the partial offset because of entropy loss, suggesting that hydrogen bonding stiffens the antigen-antibody complex.     

Prediction of hot spot residues at protein-protein interfaces by combining machine learning and energy-based methods

Background  

Alanine scanning mutagenesis is a powerful experimental methodology for investigating the structural and energetic characteristics of protein complexes. Individual amino-acids are systematically mutated to alanine and changes in free energy of binding (ΔΔG) measured. Several experiments have shown that protein-protein interactions are critically dependent on just a few residues ("hot spots") at the interface. Hot spots make a dominant contribution to the free energy of binding and if mutated they can disrupt the interaction. As mutagenesis studies require significant experimental efforts, there is a need for accurate and reliable computational methods. Such methods would also add to our understanding of the determinants of affinity and specificity in protein-protein recognition.     

Interaction of monoclonal anti-peptide antibodies with lysozyme

The interaction of monoclonal anti-peptide antibodies with the free peptide and its protein counterpart has been evaluated for hen egg white lysozyme and the peptide constituting residues 38 to 45. Fluorescence methodology has been developed for the measurement of association constants based on resonance energy transfer between the excited tryptophan of antibody and bound peptide ligand conjugated to a fluorescent probe. Five antibodies, four IgM and one IgG, have been assayed by ELISA, and have demonstrated binding to the adsorbed peptide alone, to the adsorbed lysozyme alone, or to both. Multivalent interaction with the adsorbed ligand is a key factor in the efficacy of binding. Measurement of binding constants in homogeneous solution, by equilibrium dialysis and energy transfer, demonstrated that lysozyme was bound to an IgG antipeptide antibody with an association constant (4 X 10(2) M-1) 200-fold less than that for the free peptide (8 X 10(4) M-1). It was also inferred for IgM that an association constant of the order of 10(2) M-1 was sufficient to effect selective interaction in a system providing multivalent interaction. The shared conformations between protein and peptide, implied by the specific reactivity of the anti-peptide antibody with the protein, points to structural fluctuations of the surface regions and residues of globular proteins.     

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.     

On the attribution of binding energy in antigen-antibody complexes McPC 603, D1.3, and HyHEL-5

Using X-ray coordinates of antigen-antibody complexes McPC 603, D1.3, and HyHEL-5, we made semiquantitative estimates of Gibbs free energy changes (delta G) accompanying noncovalent complex formation of the McPC 603 Fv fragment with phosphocholine and the D1.3 or HyHEL-5 Fv fragments with hen egg white lysozyme. Our empirical delta G function, which implicitly incorporates solvent effects, has the following components: hydrophobic force, solvent-modified electrostatics, changes in side-chain conformational entropy, translational/overall rotational entropy changes, and the dilutional (cratic) entropy term. The calculated delta G ranges matched the experimentally determined delta G of McPC 603 and D1.3 complexes and overestimated it (i.e., gave a more negative value) in the case of HyHEL-5. Relative delta G contributions of selected antibody residues, calculated for HyHEL-5 complexes, agreed with those determined independently in site-directed mutagenesis experiments. Analysis of delta G attribution in all three complexes indicated that only a small number of amino acids probably contribute actively to binding energetics. These form a subset of the total antigen-antibody contact surface. In the antibodies, the bottom part of the antigen binding cavity dominated the energetics of binding whereas in lysozyme, the energetically most important residues defined small (2.5-3 nm2) "energetic" epitopes. Thus, a concept of protein antigenicity emerges that involves the active, attractive contributions mediated by the energetic antigenic epitopes and the passive surface complementarity contributed by the surrounding contact area. The D1.3 energetic epitope of lysozyme involved Gly 22, Gly 117, and Gln 121; the HyHEL-5 epitope consisted of Arg 45 and Arg 68. These are also the essential antigenic residues determined experimentally. The above positions belong to the most protruding parts of the lysozyme surface, and their backbones are not exceptionally flexible. Least-squares analysis of six different antibody binding regions indicated that the geometry of the VH-VL interface beta-barrel is well conserved, giving no indication of significant changes in domain-domain contacts upon complex formation.     

X-ray snapshots of the maturation of an antibody response to a protein antigen

The process whereby the immune system generates antibodies of higher affinities during a response to antigen (affinity maturation) is a prototypical example of molecular evolution. Earlier studies have been confined to antibodies specific for small molecules (haptens) rather than for proteins. We compare the structures of four antibodies bound to the same site on hen egg white lysozyme (HEL) at different stages of affinity maturation. These X-ray snapshots reveal that binding is enhanced, not through the formation of additional hydrogen bonds or van der Waals contacts or by an increase in total buried surface, but by burial of increasing amounts of apolar surface at the expense of polar surface, accompanied by improved shape complementarity. The increase in hydrophobic interactions results from highly correlated rearrangements in antibody residues at the interface periphery, adjacent to the central energetic hot spot. This first visualization of the maturation of antibodies to protein provides insights into the evolution of high affinity in other protein-protein interfaces.     

Amino Acid Sequence and Activity of Green Turtle (Chelonia mydas) Lysozyme

Green turtle lysozyme purified from egg white was sequenced and analyzed its activity. Lysozyme was reduced and pyridylethylated or carboxymethylated to digest with trypsin, chymotrypsin and V8 protease. The peptides yielded were purified by RP-HPLC and sequenced. Every trypsin peptide was overlapped by chymotrypsin peptides and V8 protease peptides. This lysozyme is composed of 130 amino acids including an insertion of a Gly residue between 47 and 48 residues when compared with chicken lysozyme. The amino acid substitutions were found at subsites E and F. Namely Phe34, Arg45, Thr47, and Arg114 were replaced by Tyr, Tyr, Pro, and Asn, respectively. The time course using N-acetylglucosamine pentamer as a substrate showed a reduction of the rate constant of glycosidic cleavage and transglycosylation and increase of binding free energy for subsite E, which proved the contribution of amino acids mentioned above for substrate binding at subsites E and F.     

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.     

Role of bound water in protein-ligand association processes

The role of water molecules on the protein-ligand interface during macromolecular association has been determined. The free energy of association of insulin has been calculated by the methods of molecular mechanics and continual electrostatics (Poisson-Boltzmann model). The previously developed scheme of the decomposition of association free energy onto contributions from individual interactions has been used to calculate intermolecular interactions, the solvation free energy, and the entropies of the process of macromolecular association. An analysis of the calculated oscillation spectra indicated that the presence of water molecules on the protein-protein interface promotes an increase in the contribution of vibration entropy to the free energy of association due to the enhancement of the flexibility of the complex. It was shown that water molecules involved in the formation of protein-water-ligand hydrogen bond network change the balance of forces in the system.     

Crystal structure of anti-Hen egg white lysozyme antibody (HyHEL-10) Fv-antigen complex. Local structural changes in the protein antigen and water-mediated interactions of Fv-antigen and light chain-heavy chain interfaces.

In order to address the recognition mechanism of the fragments of antibody variable regions, termed Fv, toward their target antigen, an x-ray crystal structure of an anti-hen egg white lysozyme antibody (HyHEL-10) Fv fragment complexed with its cognate antigen, hen egg white lysozyme (HEL), was solved at 2.3 A. The overall structure of the complex is similar to that reported in a previous article dealing with the Fab fragment-HEL complex (PDB ID code,). However, the areas of Fv covered by HEL upon complex formation increased by about 100 A(2) in comparison with the Fab-HEL complex, and two local structural differences were observed in the heavy chain of the variable region (VH). In addition, small but significant local structural changes were observed in the antigen, HEL. The x-ray data permitted the identification of two water molecules between the VH and HEL and six water molecules retained in the interface between the antigen and the light chain complementarity determining regions (CDRs) 2 and 3 (CDR-L2 and CDR-L3). These water molecules bridge the antigen-antibody interface through hydrogen bond formation in the VL-HEL interface. Eleven water molecules were found to complete the imperfect VH-VL interface, suggesting that solvent molecules mediate the stabilization of interaction between variable regions. These results suggest that the unfavorable effect of deletion of constant regions on the antigen-antibody interaction is compensated by an increase in favorable interactions, including structural changes in the antigen-antibody interface and solvent-mediated hydrogen bond formation upon complex formation, which may lead to a minimum decreased affinity of the antibody Fv fragment toward its antigen.