Survey of the geometric association of domain-domain interfaces

Considering the limited success of the most sophisticated docking methods available and the amount of computation required for systematic docking, cataloging all the known interfaces may be an alternative basis for the prediction of protein tertiary and quaternary structures. We classify domain interfaces according to the geometry of domain-domain association. By applying a simple and efficient method called "interface tag clustering," more than 4,000 distinct types of domain interfaces are collected from Protein Quaternary Structure Server and Protein Data Bank. Given a pair of interacting domains, we define "face" as the set of interacting residues in each single domain and the pair of interacting faces as an "interface." We investigate how the geometry of interfaces relates to a network of interacting protein families, such as how many different binding orientations are possible between two families or whether a family uses distinct surfaces or the same surface when the family has diverse interaction partners from various families. We show there are, on average, 1.2-1.9 different types of interfaces between interacting domains and a significant number of family pairs associate in multiple orientations. In general, a family tends to use distinct faces for each partner when the family has diverse interaction partners. Each face is highly specific to its interaction partner and the binding orientation. The relative positions of interface residues are generally well conserved within the same type of interface even between remote homologs. The classification result is available at http://www.biotec.tu-dresden.de/~wkim/supplement.     

Distinct roles of overlapping and non-overlapping regions of hub protein interfaces in recognition of multiple partners

Cellular functions of an organism are maintained by protein-protein interactions. Those proteins that bind multiple partners asynchronously (date hub proteins) are important to make the interaction network coordinated. It is known that many date hub proteins bind different partners at overlapping (OV) interfaces. To understand how OV interfaces of date hub proteins can recognize multiple partners, we analyzed the difference between OV and non-overlapping (Non-OV) regions of interfaces involved in the binding of different partners. By using the structures of 16 date hub proteins with various interaction partners (ranging from 5 to 33), we compared buried surface area, compositions of amino acid residues and secondary structures, and side-chain orientations. It was found that buried interface residues are important for recognizing multiple partners, while exposed interface residues are important for determining specificity to a particular ligand. In addition, our analyses reveal that residue compositions in OV and Non-OV regions are different and that residues in OV region show diverse side-chain torsion angles to accommodate binding to multiple targets.     

Beauty Is in the Eye of the Beholder: Proteins Can Recognize Binding Sites of Homologous Proteins in More than One Way

Understanding the mechanisms of protein–protein interaction is a fundamental problem with many practical applications. The fact that different proteins can bind similar partners suggests that convergently evolved binding interfaces are reused in different complexes. A set of protein complexes composed of non-homologous domains interacting with homologous partners at equivalent binding sites was collected in 2006, offering an opportunity to investigate this point. We considered 433 pairs of protein–protein complexes from the ABAC database (AB and AC binary protein complexes sharing a homologous partner A) and analyzed the extent of physico-chemical similarity at the atomic and residue level at the protein–protein interface. Homologous partners of the complexes were superimposed using Multiprot, and similar atoms at the interface were quantified using a five class grouping scheme and a distance cut-off. We found that the number of interfacial atoms with similar properties is systematically lower in the non-homologous proteins than in the homologous ones. We assessed the significance of the similarity by bootstrapping the atomic properties at the interfaces. We found that the similarity of binding sites is very significant between homologous proteins, as expected, but generally insignificant between the non-homologous proteins that bind to homologous partners. Furthermore, evolutionarily conserved residues are not colocalized within the binding sites of non-homologous proteins. We could only identify a limited number of cases of structural mimicry at the interface, suggesting that this property is less generic than previously thought. Our results support the hypothesis that different proteins can interact with similar partners using alternate strategies, but do not support convergent evolution.     

Thermodynamic dissection of large‐scale domain motions coupled with ligand binding of enzyme I

Domain motions are central to the biological functions of many proteins. The energetics of the motions, however, is often difficult to characterize when motions are coupled with the ligand binding. Here, we determined the thermodynamic parameters of individual domain motions and ligand binding of enzyme I (EI) using strategic domain-deletion mutants that selectively removed particular motions. Upon ligand binding, EI employs two large-scale domain motions, the hinge motion and the swivel motion, to switch between conformational states of distinct domain−domain orientations. Calorimetric analysis of the EI mutants separated the free energy changes of the binding and motions, demonstrating that the unfavorable hinge motion (ΔG = 1.5 kcal mol⁻¹) was driven by the favorable swivel motion (ΔG = −5.2 kcal mol⁻¹). The large free energy differences could be explained by the physicochemical nature of the domain interfaces associated with the motions; the hinge motion employed much narrower interface than the swivel motion without any hydrogen bonds or salt bridges. The small heat capacity further suggested that the packing of the domain interfaces associated with the hinge motion was less compact than that commonly observed in proteins. Lastly, thermodynamic analysis of phosphorylated EI suggests that the domain motions are regulated by the ligand binding and the phosphorylation states. Taken together, the thermodynamic dissection approach illustrates how multiple motions and ligand binding are energetically connected during the functional cycle of EI.     

Dual roles of Munc18-1 rely on distinct binding modes of the central cavity with Stx1A and SNARE complex

Sec1/Munc18 proteins play a fundamental role in multiple steps of intracellular membrane trafficking. Dual functions have been attributed to Munc18-1: it can act as a chaperone when it interacts with monomeric syntaxin 1A, and it can activate soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) for membrane fusion when it binds to SNARE complexes. Although both modes of binding involve the central cavity of Munc18-1, their precise molecular mechanisms of action are not fully understood. In this paper, we describe a novel Munc18-1 mutant in the central cavity that showed a reduced interaction with syntaxin 1A and impaired chaperone function, but still bound to assembled SNARE complexes and promoted liposome fusion and secretion in neuroendocrine cells. Soluble syntaxin 1A H3 domain partially blocks Munc18-1 activation of liposome fusion by occupying the Munc18-1 central cavity. Our findings lead us to propose a transition model between the two distinct binding modes by which Munc18 can control and assist in SNARE-complex assembly during neurotransmitter release.     

Structural and Functional Analysis of Multi-Interface Domains

A multi-interface domain is a domain that can shape multiple and distinctive binding sites to contact with many other domains, forming a hub in domain-domain interaction networks. The functions played by the multiple interfaces are usually different, but there is no strict bijection between the functions and interfaces as some subsets of the interfaces play the same function. This work applies graph theory and algorithms to discover fingerprints for the multiple interfaces of a domain and to establish associations between the interfaces and functions, based on a huge set of multi-interface proteins from PDB. We found that about 40% of proteins have the multi-interface property, however the involved multi-interface domains account for only a tiny fraction (1.8%) of the total number of domains. The interfaces of these domains are distinguishable in terms of their fingerprints, indicating the functional specificity of the multiple interfaces in a domain. Furthermore, we observed that both cooperative and distinctive structural patterns, which will be useful for protein engineering, exist in the multiple interfaces of a domain.     

Bayesian Modeling of the Yeast SH3 Domain Interactome Predicts Spatiotemporal Dynamics of Endocytosis Proteins

SH3 domains are peptide recognition modules that mediate the assembly of diverse biological complexes. We scanned billions of phage-displayed peptides to map the binding specificities of the SH3 domain family in the budding yeast, Saccharomyces cerevisiae. Although most of the SH3 domains fall into the canonical classes I and II, each domain utilizes distinct features of its cognate ligands to achieve binding selectivity. Furthermore, we uncovered several SH3 domains with specificity profiles that clearly deviate from the two canonical classes. In conjunction with phage display, we used yeast two-hybrid and peptide array screening to independently identify SH3 domain binding partners. The results from the three complementary techniques were integrated using a Bayesian algorithm to generate a high-confidence yeast SH3 domain interaction map. The interaction map was enriched for proteins involved in endocytosis, revealing a set of SH3-mediated interactions that underlie formation of protein complexes essential to this biological pathway. We used the SH3 domain interaction network to predict the dynamic localization of several previously uncharacterized endocytic proteins, and our analysis suggests a novel role for the SH3 domains of Lsb3p and Lsb4p as hubs that recruit and assemble several endocytic complexes.     

Design of multi-specificity in protein interfaces

Interactions in protein networks may place constraints on protein interface sequences to maintain correct and avoid unwanted interactions. Here we describe a “multi-constraint” protein design protocol to predict sequences optimized for multiple criteria, such as maintaining sets of interactions, and apply it to characterize the mechanism and extent to which 20 multi-specific proteins are constrained by binding to multiple partners. We find that multi-specific binding is accommodated by at least two distinct patterns. In the simplest case, all partners share key interactions, and sequences optimized for binding to either single or multiple partners recover only a subset of native amino acid residues as optimal. More interestingly, for signaling interfaces functioning as network “hubs,” we identify a different, “multi-faceted” mode, where each binding partner prefers its own subset of wild-type residues within the promiscuous binding site. Here, integration of preferences across all partners results in sequences much more “native-like” than seen in optimization for any single binding partner alone, suggesting these interfaces are substantially optimized for multi-specificity. The two strategies make distinct predictions for interface evolution and design. Shared interfaces may be better small molecule targets, whereas multi-faceted interactions may be more “designable” for altered specificity patterns. The computational methodology presented here is generalizable for examining how naturally occurring protein sequences have been selected to satisfy a variety of positive and negative constraints, as well as for rationally designing proteins to have desired patterns of altered specificity.     

The trihelical bundle subdomain of the GGA proteins interacts with multiple partners through overlapping but distinct sites

The Golgi-localized, gamma-adaptin ear-containing, ARF-binding (GGA) proteins are monomeric clathrin adaptors that mediate the sorting of cargo at the trans-Golgi network and endosomes. The GGAs contain four different domains named Vps27, Hrs, Stam (VHS); GGAs and TOM1 (GAT); hinge; and gamma-adaptin ear (GAE). The VHS domain recognizes transmembrane cargo, whereas the hinge and GAE regions bind clathrin and accessory proteins, respectively. The GAT domain is a polyfunctional module that interacts with various partners including the small GTPase ARF, the endosomal fusion regulator Rabaptin-5, ubiquitin, and the product of the tumor susceptibility gene 101 (TSG101). Previous x-ray crystallographic analyses showed that the GAT region is composed of two subdomains, an N-terminal helix-loop-helix containing the ARF binding site, and a C-terminal triple alpha-helical (trihelical) bundle. In this study, we define the Rabaptin-5 binding site on the GGA1-GAT domain and its relationship to the binding sites for ubiquitin and TSG101. Our observations show that Rabaptin-5, ubiquitin, and TSG101 bind to overlapping but distinct binding sites on the trihelical bundle. The different GAT binding partners engage in both competitive and cooperative interactions that may be important for the function of the GGAs in protein sorting.     

Molecular principles of the interactions of disordered proteins

Thorough knowledge of the molecular principles of protein-protein recognition is essential to our understanding of protein function at the cellular level. Whereas interactions of ordered proteins have been analyzed in great detail, complexes of intrinsically unstructured/disordered proteins (IUPs) have hardly been addressed so far. Here, we have collected a database of 39 complexes of experimentally verified IUPs, and compared their interfaces with those of 72 complexes of ordered, globular proteins. The characteristic differences found between the two types of complexes suggest that IUPs represent a distinct molecular implementation of the principles of protein-protein recognition. The interfaces do not differ in size, but those of IUPs cover a much larger part of the surface of the protein than for their ordered counterparts. Moreover, IUP interfaces are significantly more hydrophobic relative to their overall amino acid composition, but also in absolute terms. They rely more on hydrophobic-hydrophobic than on polar-polar interactions. Their amino acids in the interface realize more intermolecular contacts, which suggests a better fit with the partner due to induced folding upon binding that results in a better adaptation to the partner. The two modes of interaction also differ in that IUPs usually use only a single continuous segment for partner binding, whereas the binding sites of ordered proteins are more segmented. Probably, all these features contribute to the increased evolutionary conservation of IUP interface residues. These noted molecular differences are also manifested in the interaction energies of IUPs. Our approximation of these by low-resolution force-fields shows that IUPs gain much more stabilization energy from intermolecular contacts, than from folding, i.e. they use their binding energy for folding. Overall, our findings provide a structural rationale to the prior suggestions that many IUPs are specialized for functions realized by protein-protein interactions.     

Use of Staphylococcus aureus 6-P-beta-galactosidase and GFP as fusion partners for lactose-specific IIC domain from Staphylococcus aureus

The hydrophilic part of membrane proteins plays an important role in the formation of 3D crystals. The construction of fusion proteins using well crystallizing proteins as fusion partners is a possibility to increase the hydrophilic part of membrane proteins lacking large hydrophilic domains. These fusion proteins might be easier to crystallize. Two bifunctional fusion proteins containing the membrane-bound, lactose-specific enzyme IIC domain of the lactose transporter (IICB(lac)) from S. aureus as N-terminal fusion partner were constructed by gene fusion. The C-terminal fusion partners were S. aureus 6-P-beta-Galactosidase and GFP, respectively. Both proteins were overexpressed in E. coli, purified to homogeneity and kinetically characterized: In the presence of the components of the lactose phosphotransferase system of S. aureus, the hybrid proteins phosphorylated their substrates, indicating that the fusion partners are sufficiently flexibly linked to allow the interaction of the IIC(lac) domain with the IIB(lac) domain of the lactose transporter. The activity of the 6-P-beta-Galactosidase as well as the fluorescence of GFP were preserved in the fusion proteins. The Vmax values determined for the IIC domain in the fusion proteins were dramatically reduced compared with the values determined for the separate IIC(lac) domain and the complete lactose transporter (IICB(lac)). The Km values were only slightly increased indicating that the Vmax values are much more influenced by the fusion than the substrate affinities. The substrate affinity and the Vmax value determined for the GFP-fused IIC(lac) domain are higher than for the 6-P-beta-Galactosidase-fused IIC(lac). The results suggest that the fusion with GFP enables a better interaction with the IIB(lac) domain than the fusion with 6-P-beta-Galactosidase. Moreover, the GFP-fused IIC(lac) domain proved to be more stable than the 6-P-beta-Galactosidase fusion protein.     

Plasticity of BRCA2 function in homologous recombination: genetic interactions of the PALB2 and DNA binding domains

The breast cancer suppressor BRCA2 is essential for the maintenance of genomic integrity in mammalian cells through its role in DNA repair by homologous recombination (HR). Human BRCA2 is 3,418 amino acids and is comprised of multiple domains that interact with the RAD51 recombinase and other proteins as well as with DNA. To gain insight into the cellular function of BRCA2 in HR, we created fusions consisting of various BRCA2 domains and also introduced mutations into these domains to disrupt specific protein and DNA interactions. We find that a BRCA2 fusion peptide deleted for the DNA binding domain and active in HR is completely dependent on interaction with the PALB2 tumor suppressor for activity. Conversely, a BRCA2 fusion peptide deleted for the PALB2 binding domain is dependent on an intact DNA binding domain, providing a role for this conserved domain in vivo; mutagenesis suggests that both single-stranded and double-stranded DNA binding activities in the DNA binding domain are required for its activity. Given that PALB2 itself binds DNA, these results suggest alternative mechanisms to deliver RAD51 to DNA. In addition, the BRCA2 C terminus contains both RAD51-dependent and -independent activities which are essential to HR in some contexts. Finally, binding the small peptide DSS1 is essential for activity when its binding domain is present, but not when it is absent. Our results reveal functional redundancy within the BRCA2 protein and emphasize the plasticity of this large protein built for optimal HR function in mammalian cells. The occurrence of disease-causing mutations throughout BRCA2 suggests sub-optimal HR from a variety of domain modulations.     

Multiple roles of the vesicular-SNARE TI-VAMP in post-Golgi and endosomal trafficking

SNARE (Soluble N-ethylmaleimide-sensitive factor attachment protein receptor) proteins are the core machinery of membrane fusion. Vesicular SNAREs (v-SNAREs) interact with their target SNAREs (t-SNAREs) to form SNARE complexes which mediate membrane fusion. Here we review the basic properties and functions of the v-SNARE TI-VAMP/VAMP7 (Tetanus neurotoxin insensitive-vesicle associated membrane protein). TI-VAMP interacts with its t-SNARE partners, particularly plasmalemmal syntaxins, to mediate membrane fusion and with several regulatory proteins especially via its amino-terminal regulatory Longin domain. Partners include AP-3, Hrb/(Human immunodeficiency virus Rev binding) protein, and Varp (Vps9 domain and ankyrin repeats containing protein) and regulate TI-VAMP’s function and targeting. TI-VAMP is involved both in secretory and endocytic pathways which mediate neurite outgrowth and synaptic transmission, plasma membrane remodeling and lysosomal secretion.     

The interactome of LIM domain proteins: the contributions of LIM domain proteins to heart failure and heart development

LIM domain proteins all contain at least one double zinc-finger motif. They belong to a large family and here we review those expressed mainly in mammalian hearts, but particularly in cardiomyocytes. These proteins contain between one and five LIM domains and usually these proteins contain other domains that have specific functions such as actin-binding, kinases and nuclear translocation motifs. While several recent reviews have summarised the importance of individual LIM domain proteins, this is the first review of its kind to cover all LIMs associated with the heart. Here we examine 33 LIM proteins (including three that bind to, but do not themselves contain, LIM domains) that are implicated in either the development of the heart, heart disorders and failure, or both. Our analysis is consistent with the view that cardiac LIM domain proteins form multiple extensive networks of multi-protein complexes within the myocardium. This multiplicity of binding partners probably protects the heart as it is challenged to maintain cardiac output, until the imbalance reaches a turning point that results in failure. We believe that the complexity of LIM interactions is properly described by the term LIM interactome.     

Specificity of molecular interactions in transient protein-protein interaction interfaces

In this study, we investigate what types of interactions are specific to their biological function, and what types of interactions are persistent regardless of their functional category in transient protein-protein heterocomplexes. This is the first approach to analyze protein-protein interfaces systematically at the molecular interaction level in the context of protein functions. We perform systematic analysis at the molecular interaction level using classification and feature subset selection technique prevalent in the field of pattern recognition. To represent the physicochemical properties of protein-protein interfaces, we design 18 molecular interaction types using canonical and noncanonical interactions. Then, we construct input vector using the frequency of each interaction type in protein-protein interface. We analyze the 131 interfaces of transient protein-protein heterocomplexes in PDB: 33 protease-inhibitors, 52 antibody-antigens, 46 signaling proteins including 4 cyclin dependent kinase and 26 G-protein. Using kNN classification and feature subset selection technique, we show that there are specific interaction types based on their functional category, and such interaction types are conserved through the common binding mechanism, rather than through the sequence or structure conservation. The extracted interaction types are C(alpha)-- H...O==C interaction, cation...anion interaction, amine...amine interaction, and amine...cation interaction. With these four interaction types, we achieve the classification success rate up to 83.2% with leave-one-out cross-validation at k = 15. Of these four interaction types, C(alpha)--H...O==C shows binding specificity for protease-inhibitor complexes, while cation-anion interaction is predominant in signaling complexes. The amine ... amine and amine...cation interaction give a minor contribution to the classification accuracy. When combined with these two interactions, they increase the accuracy by 3.8%. In the case of antibody-antigen complexes, the sign is somewhat ambiguous. From the evolutionary perspective, while protease-inhibitors and sig-naling proteins have optimized their interfaces to suit their biological functions, antibody-antigen interactions are the happenstance, implying that antibody-antigen complexes do not show distinctive interaction types. Persistent interaction types such as pi...pi, amide-carbonyl, and hydroxyl-carbonyl interaction, are also investigated. Analyzing the structural orientations of the pi...pi stacking interactions, we find that herringbone shape is a major configuration in transient protein-protein interfaces. This result is different from that of protein core, where parallel-displaced configurations are the major configuration. We also analyze overall trend of amide-carbonyl and hydroxyl-carbonyl interactions. It is noticeable that nearly 82% of the interfaces have at least one hydroxyl-carbonyl interactions.     

The many faces of Ras: recognition of small GTP-binding proteins

The structures of over 30 complexes of Ras superfamily small GTP-binding proteins bound to diverse protein partners have been reported. Comparison of these complexes using the sequences of the small GTP-binding proteins to align the contact sites shows that virtually all surface positions make contacts with at least one partner protein. Rather than highlighting a single consensus binding site, these comparisons illustrate the remarkable diversity of contacts of Ras superfamily members. Here, a new analysis technique, the interface array, is introduced to quantify patterns of surface contacts. The interface array shows that small GTP-binding proteins are recognized in at least nine distinct ways. Remarkably, binding partners with similar functions, including those with distinct folds, recognize small GTP-binding proteins in similar ways. These classes of shared surface contacts support the occurrence of both divergent and convergent evolutionary processes and suggest that specific effector functions require particular protein–protein contacts.     

Selective formation of Sed5p-containing SNARE complexes is mediated by combinatorial binding interactions

Sed5p is the only syntaxin family member required for protein transport through the yeast Golgi and it is known to bind up to nine other soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) proteins in vivo. We describe in vitro binding experiments in which we identify ternary and quaternary Sed5p-containing SNARE complexes. The formation of SNARE complexes among these endoplasmic reticulum- and Golgi-localized proteins requires Sed5p and is syntaxin-selective. In addition, Sed5p-containing SNARE complexes form selectively and this selectivity is mediated by Sed5p-containing intermediates that discriminate among subsequent binding partners. Although many of these SNAREs have overlapping distributions in vivo, the SNAREs that form complexes with Sed5p in vitro reflect their functionally distinct locales. Although SNARE-SNARE interactions are promiscuous and a single SNARE protein is often found in more than one complex, both the biochemical as well as genetic analyses reported here suggest that this is not a result of nonselective direct substitution of one SNARE for another. Rather our data are consistent with the existence of multiple (perhaps parallel) trafficking pathways where Sed5p-containing SNARE complexes play overlapping and/or distinct functional roles.     

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.     

Exploring functional roles of multibinding protein interfaces

Cellular processes are highly interconnected and many proteins are shared in different pathways. Some of these shared proteins or protein families may interact with diverse partners using the same interface regions; such multibinding proteins are the subject of our study. The main goal of our study is to attempt to decipher the mechanisms of specific molecular recognition of multiple diverse partners by promiscuous protein regions. To address this, we attempt to analyze the physicochemical properties of multibinding interfaces and highlight the major mechanisms of functional switches realized through multibinding. We find that only 5% of protein families in the structure database have multibinding interfaces, and multibinding interfaces do not show any higher sequence conservation compared with the background interface sites. We highlight several important functional mechanisms utilized by multibinding families. (a) Overlap between different functional pathways can be prevented by the switches involving nearby residues of the same interfacial region. (b) Interfaces can be reused in pathways where the substrate should be passed from one protein to another sequentially. (c) The same protein family can develop different specificities toward different binding partners reusing the same interface; and finally, (d) inhibitors can attach to substrate binding sites as substrate mimicry and thereby prevent substrate binding.