Hydration of protein-protein interfaces

We present an analysis of the water molecules immobilized at the protein-protein interfaces of 115 homodimeric proteins and 46 protein-protein complexes, and compare them with 173 large crystal packing interfaces representing nonspecific interactions. With an average of 15 waters per 1000 A2 of interface area, the crystal packing interfaces are more hydrated than the specific interfaces of homodimers and complexes, which have 10-11 waters per 1000 A2, reflecting the more hydrophilic composition of crystal packing interfaces. Very different patterns of hydration are observed: Water molecules may form a ring around interfaces that remain "dry," or they may permeate "wet" interfaces. A majority of the specific interfaces are dry and most of the crystal packing interfaces are wet, but counterexamples exist in both categories. Water molecules at interfaces form hydrogen bonds with protein groups, with a preference for the main-chain carbonyl and the charged side-chains of Glu, Asp, and Arg. These interactions are essentially the same in specific and nonspecific interfaces, and very similar to those observed elsewhere on the protein surface. Water-mediated polar interactions are as abundant at the interfaces as direct protein-protein hydrogen bonds, and they may contribute to the stability of the assembly.     

A "solvated rotamer" approach to modeling water-mediated hydrogen bonds at protein-protein interfaces

Water-mediated hydrogen bonds play critical roles at protein-protein and protein-nucleic acid interfaces, and the interactions formed by discrete water molecules cannot be captured using continuum solvent models. We describe a simple model for the energetics of water-mediated hydrogen bonds, and show that, together with knowledge of the positions of buried water molecules observed in X-ray crystal structures, the model improves the prediction of free-energy changes upon mutation at protein-protein interfaces, and the recovery of native amino acid sequences in protein interface design calculations. We then describe a "solvated rotamer" approach to efficiently predict the positions of water molecules, at protein-protein interfaces and in monomeric proteins, that is compatible with widely used rotamer-based side-chain packing and protein design algorithms. Finally, we examine the extent to which the predicted water molecules can be used to improve prediction of amino acid identities and protein-protein interface stability, and discuss avenues for overcoming current limitations of the approach.     

Predicting where Small Molecules Bind at Protein-Protein Interfaces

Small molecules that bind at protein-protein interfaces may either block or stabilize protein-protein interactions in cells. Thus, some of these binding interfaces may turn into prospective targets for drug design. Here, we collected 175 pairs of protein-protein (PP) complexes and protein-ligand (PL) complexes with known three-dimensional structures for which (1) one protein from the PP complex shares at least 40% sequence identity with the protein from the PL complex, and (2) the interface regions of these proteins overlap at least partially with each other. We found that those residues of the interfaces that may bind the other protein as well as the small molecule are evolutionary more conserved on average, have a higher tendency of being located in pockets and expose a smaller fraction of their surface area to the solvent than the remaining protein-protein interface region. Based on these findings we derived a statistical classifier that predicts patches at binding interfaces that have a higher tendency to bind small molecules. We applied this new prediction method to more than 10 000 interfaces from the protein data bank. For several complexes related to apoptosis the predicted binding patches were in direct contact to co-crystallized small molecules.     

Revisiting the Voronoi description of protein-protein interfaces

We developed a model of macromolecular interfaces based on the Voronoi diagram and the related alpha-complex, and we tested its properties on a set of 96 protein-protein complexes taken from the Protein Data Bank. The Voronoi model provides a natural definition of the interfaces, and it yields values of the number of interface atoms and of the interface area that have excellent correlation coefficients with those of the classical model based on solvent accessibility. Nevertheless, some atoms that do not lose solvent accessibility are part of the interface defined by the Voronoi model. The Voronoi model provides robust definitions of the curvature and of the connectivity of the interfaces, and leads to estimates of these features that generally agree with other approaches. Our implementation of the model allows an analysis of protein-water contacts that highlights the role of structural water molecules at protein-protein interfaces.     

Strategies to search and design stabilizers of protein-protein interactions: a feasibility study

Since protein-protein interactions play a pivotal role in the communication on the molecular level in virtually every biological system and process, the search and design for modulators of such interactions is of utmost importance. In recent years many inhibitors for specific protein-protein interactions have been developed, however, in only a few cases, small and druglike molecules are able to interfere in the complex formation of proteins. On the other hand, there are several small molecules known to modulate protein-protein interactions by means of stabilizing an already assembled complex. To achieve this goal, a ligand is binding to a pocket, which is located rim-exposed at the interface of the interacting proteins, for example as the phytotoxin Fusicoccin, which stabilizes the interaction of plant H+-ATPase and 14-3-3 protein by nearly a factor of 100. To suggest alternative leads, we performed a virtual screening campaign to discover new molecules putatively stabilizing this complex. Furthermore, we screen a dataset of 198 transient recognition protein-protein complexes for cavities, which are located rim-exposed at their interfaces. We provide evidence for high similarity between such rim-exposed cavities and usual ligands accommodating active sites of enzymes. This analysis suggests that rim-exposed cavities at protein-protein interfaces are druggable binding sites. Therefore, the principle of stabilizing protein-protein interactions seems to be a promising alternative to the approach of the competitive inhibition of such interactions by small molecules.     

Polyamines: Naturally occurring small molecule modulators of electrostatic protein-protein interactions

Modulations of protein-protein interactions are a key step in regulating protein function, especially in networks. Modulators of these interactions are supposed to be candidates for the development of novel drugs. Here, we describe the role of the small, polycationic and highly abundant natural polyamines that could efficiently bind to charged spots at protein interfaces as modulators of such protein-protein interactions. Using the mitochondrial cytochrome P45011A1 (CYP11A1) electron transfer system as a model, we have analyzed the capability of putrescine, spermidine, and spermine at physiologically relevant concentrations to affect the protein-protein interactions between adrenodoxin reductase (AdR), adrenodoxin (Adx), and CYP11A1. The actions of polyamines on the individual components, on their association/dissociation, on electron transfer, and on substrate conversion were examined. These studies revealed modulating effects of polyamines on distinct interactions and on the entire system in a complex way. Modulation via changed protein-protein interactions appeared plausible from docking experiments that suggested favourable high-affinity binding sites of polyamines (spermine > spermidine > putrescine) at the AdR-Adx interface. Our findings imply for the first time that small endogenous compounds are capable of interfering with distinct components of transient protein complexes and might control protein functions by modulating electrostatic protein-protein interactions.     

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.     

Modeling the Alzheimer Abeta17-42 fibril architecture: tight intermolecular sheet-sheet association and intramolecular hydrated cavities

We investigate Abeta(17-42) protofibril structures in solution using molecular dynamics simulations. Recently, NMR and computations modeled the Abeta protofibril as a longitudinal stack of U-shaped molecules, creating an in-parallel beta-sheet and loop spine. Here we study the molecular architecture of the fibril formed by spine-spine association. We model in-register intermolecular beta-sheet-beta-sheet associations and study the consequences of Alzheimer's mutations (E22G, E22Q, E22K, and M35A) on the organization. We assess the structural stability and association force of Abeta oligomers with different sheet-sheet interfaces. Double-layered oligomers associating through the C-terminal-C-terminal interface are energetically more favorable than those with the N-terminal-N-terminal interface, although both interfaces exhibit high structural stability. The C-terminal-C-terminal interface is essentially stabilized by hydrophobic and van der Waals (shape complementarity via M35-M35 contacts) intermolecular interactions, whereas the N-terminal-N-terminal interface is stabilized by hydrophobic and electrostatic interactions. Hence, shape complementarity, or the "steric zipper" motif plays an important role in amyloid formation. On the other hand, the intramolecular Abeta beta-strand-loop-beta-strand U-shaped motif creates a hydrophobic cavity with a diameter of 6-7 A, allowing water molecules and ions to conduct through. The hydrated hydrophobic cavities may allow optimization of the sheet association and constitute a typical feature of fibrils, in addition to the tight sheet-sheet association. Thus, we propose that Abeta fiber architecture consists of alternating layers of tight packing and hydrated cavities running along the fibrillar axis, which might be possibly detected by high-resolution imaging.     

SANS/SAXS study of the BSA solvation properties in aqueous urea solutions via a global fit approach

We report on the solvation properties and intermolecular interactions of a model protein (bovine serum albumine, BSA) in urea aqueous solutions, as obtained by combining small-angle neutron and X-ray scattering experiments. According to a global fit strategy, all the whole set of scattering curves are analysed by considering a unique model which includes the BSA structure, the protein-protein interactions and the thermodynamic exchange process of water/urea molecules at the protein solvent interface. As a main result, the equilibrium constant that accounts for the difference in composition between the bulk solvent and the protein solvation layer is derived. Results confirm that urea preferentially sticks to the protein surface, inducing a noticeable change in both the repulsive and the attractive interaction potentials.     

ISIS: interaction sites identified from sequence

MOTIVATION: Large-scale experiments reveal pairs of interacting proteins but leave the residues involved in the interactions unknown. These interface residues are essential for understanding the mechanism of interaction and are often desired drug targets. Reliable identification of residues that reside in protein-protein interface typically requires analysis of protein structure. Therefore, for the vast majority of proteins, for which there is no high-resolution structure, there is no effective way of identifying interface residues. RESULTS: Here we present a machine learning-based method that identifies interacting residues from sequence alone. Although the method is developed using transient protein-protein interfaces from complexes of experimentally known 3D structures, it never explicitly uses 3D information. Instead, we combine predicted structural features with evolutionary information. The strongest predictions of the method reached over 90% accuracy in a cross-validation experiment. Our results suggest that despite the significant diversity in the nature of protein-protein interactions, they all share common basic principles and that these principles are identifiable from sequence alone.     

Interaction between bound water molecules and local protein structures: A statistical analysis of the hydrogen bond structures around bound water molecules

Water molecules play an important role in protein folding and protein interactions through their structural association with proteins. Examples of such structural association can be found in protein crystal structures, and can often explain protein functionality in the context of structure. We herein report the systematic analysis of the local structures of proteins interacting with water molecules, and the characterization of their geometric features. We first examined the interaction of water molecules with a large local interaction environment by comparing the preference of water molecules in three regions, namely, the protein–protein interaction (PPI) interfaces, the crystal contact (CC) interfaces, and the non‐interfacial regions. High preference of water molecules to the PPI and CC interfaces was found. In addition, the bound water on the PPI interface was more favorably associated with the complex interaction structure, implying that such water‐mediated structures may participate in the shaping of the PPI interface. The pairwise water‐mediated interaction was then investigated, and the water‐mediated residue–residue interaction potential was derived. Subsequently, the types of polar atoms surrounding the water molecules were analyzed, and the preference of the hydrogen bond acceptor was observed. Furthermore, the geometries of the structures interacting with water were analyzed, and it was found that the major structure on the protein surface exhibited planar geometry rather than tetrahedral geometry. Several previously undiscovered characteristics of water–protein interactions were unfolded in this study, and are expected to lead to a better understanding of protein structure and function. Proteins 2016; 84:43–51. © 2015 Wiley Periodicals, Inc.     

Structural similarity and classification of protein interaction interfaces

Interactions between proteins play a key role in many cellular processes. Studying protein-protein interactions that share similar interaction interfaces may shed light on their evolution and could be helpful in elucidating the mechanisms behind stability and dynamics of the protein complexes. When two complexes share structurally similar subunits, the similarity of the interaction interfaces can be found through a structural superposition of the subunits. However, an accurate detection of similarity between the protein complexes containing subunits of unrelated structure remains an open problem. Here, we present an alignment-free machine learning approach to measure interface similarity. The approach relies on the feature-based representation of protein interfaces and does not depend on the superposition of the interacting subunit pairs. Specifically, we develop an SVM classifier of similar and dissimilar interfaces and derive a feature-based interface similarity measure. Next, the similarity measure is applied to a set of 2,806×2,806 binary complex pairs to build a hierarchical classification of protein-protein interactions. Finally, we explore case studies of similar interfaces from each level of the hierarchy, considering cases when the subunits forming interactions are either homologous or structurally unrelated. The analysis has suggested that the positions of charged residues in the homologous interfaces are not necessarily conserved and may exhibit more complex conservation patterns.     

Crystal structural analysis of protein-protein interactions drastically destabilized by a single mutation

The complex of barnase (bn) and barstar (bs), which has been widely studied as a model for quantitative analysis of protein-protein interactions, is significantly destabilized by a single mutation, namely, bs Asp39 --> Ala, which corresponds to a change of 7.7 kcal x mol(-1) in the free energy of binding. However, there has been no structural information available to explain such a drastic destabilization. In the present study, we determined the structure of the mutant complex at 1.58 A resolution by X-ray crystallography. The complex was similar to the wild-type complex in terms of overall and interface structures; however, the hydrogen bond network mediated by water molecules at the interface was significantly different. Several water molecules filled the cavity created by the mutation and consequently caused rearrangement of the hydrated water molecules at the interface. The water molecules were redistributed into a channel-like structure that penetrated into the complex. Furthermore, molecular dynamics simulations showed that the mutation increased the mobility of water molecules at the interface. Since such a drastic change in hydration was not observed in other mutant complexes of bn and bs, the significant destabilization of the interaction may be due to this channel-like structure of hydrated water molecules.     

Oligomeric protein structure networks: insights into protein-protein interactions

Background  

Protein-protein association is essential for a variety of cellular processes and hence a large number of investigations are being carried out to understand the principles of protein-protein interactions. In this study, oligomeric protein structures are viewed from a network perspective to obtain new insights into protein association. Structure graphs of proteins have been constructed from a non-redundant set of protein oligomer crystal structures by considering amino acid residues as nodes and the edges are based on the strength of the non-covalent interactions between the residues. The analysis of such networks has been carried out in terms of amino acid clusters and hubs (highly connected residues) with special emphasis to protein interfaces.     

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.     

Hot-spot residues at the E9/Im9 interface help binding via different mechanisms

Protein-protein association involves many interface interactions, but they do not contribute equally. Ala scanning experiments reveal that only a few mutations significantly lower binding affinity. These key residues, which appear to drive protein-protein association, are called hot-spot residues. Molecular dynamics simulations of the Colicin E9/Im9 complex show Im9 Glu41 and Im9 Ser50, both hot-spots, bind via different mechanisms. The results suggest that Im9 Ser50 restricts Glu41 in a conformation auspicious for salt-bridge formation across the interface. This type of model may be helpful in engineering hot-spot clusters at protein-protein interfaces and, consequently, the design of specificity. (c) 2008 Wiley Periodicals, Inc. Biopolymers 89: 916-920, 2008.This article was originally published online as an accepted preprint. The "Published Online" date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com.     

Modulating protein-protein interactions with small molecules: the importance of binding hotspots

The modulation of protein-protein interactions (PPIs) by small drug-like molecules is a relatively new area of research and has opened up new opportunities in drug discovery. However, the progress made in this area is limited to a handful of known cases of small molecules that target specific diseases. With the increasing availability of protein structure complexes, it is highly important to devise strategies exploiting homologous structure space on a large scale for discovering putative PPIs that could be attractive drug targets. Here, we propose a scheme that allows performing large-scale screening of all protein complexes and finding putative small-molecule and/or peptide binding sites overlapping with protein-protein binding sites (so-called "multibinding sites"). We find more than 600 nonredundant proteins from 60 protein families with multibinding sites. Moreover, we show that the multibinding sites are mostly observed in transient complexes, largely overlap with the binding hotspots and are more evolutionarily conserved than other interface sites. We investigate possible mechanisms of how small molecules may modulate protein-protein binding and discuss examples of new candidates for drug design.