What do we learn from high-throughput protein interaction data?

The biological significance of protein interactions, their method of generation and reliability is briefly reviewed. Protein interaction networks adopt a scale-free topology that explains their error tolerance or vulnerability, depending on whether hubs or peripheral proteins are attacked. Networks also allow the prediction of protein function from their interaction partners and therefore, the formulation of analytical hypotheses. Comparative network analysis predicts interactions for distantly related species based on conserved interactions, even if sequences are only weakly conserved. Finally, the medical relevance of protein interaction analysis is discussed and the necessity for data integration is emphasized.     

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.     

Classification of proteins based on the properties of the ligand-binding site: the case of adenine-binding proteins

Comparative analysis of protein binding sites for similar ligands yields information about conserved interactions, relevant for ligand affinity, and variable interactions, which are important for specificity. The pattern of variability can indicate new targets for a pharmacologically validated class of compounds binding to a similar site. A particularly vast group of therapeutically interesting proteins using the same or similar substrates are those that bind adenine-containing ligands. Drug development is focusing on compounds occupying the adenine-binding site and their specificity is an issue of paramount importance. We use a simple scheme to characterize and classify the adenine-binding sites in terms of their intermolecular interactions, and show that this classification does not necessarily correspond to protein classifications based on either sequence or structural similarity. We find that only a limited number of the different hydrogen bond patterns possible for adenine-binding is used, which can be utilized as an effective classification scheme. Closely related protein families usually share similar hydrogen patterns, whereas non-polar interactions are less well conserved. Our classification scheme can be used to select groups of proteins with a similar ligand-binding site, thus facilitating the definition of the properties that can be exploited to design specific inhibitors.     

RIP death domain structural interactions implicated in TNF-mediated proliferation and survival

Death domain (DD)-containing proteins are involved in both apoptosis and survival/proliferation signaling induced by activated death receptors. Here, a phylogenetic and structural analysis was performed to highlight differences in DD domains and their key regulatory interaction sites. The phylogenetic analysis shows that receptor DDs are more conserved than DDs in adaptors. Adaptor DDs can be subdivided into those that activate or inhibit apoptosis. Modeling of six homotypic DD interactions involved in the TNF signaling pathway implicates that the DD of RIP (Receptor interacting protein kinase 1) is capable of interacting with the DD of TRADD (TNFR1-associated death domain protein) in two different, exclusive ways: one that subsequently recruits CRADD (apoptosis/inflammation) and another that recruits NFkappaB (survival/proliferation).     

Finding biologically relevant protein domain interactions: conserved binding mode analysis

Proteins evolved through the shuffling of functional domains, and therefore, the same domain can be found in different proteins and species. Interactions between such conserved domains often involve specific, well-determined binding surfaces reflecting their important biological role in a cell. To find biologically relevant interactions we developed a method of systematically comparing and classifying protein domain interactions from the structural data. As a result, a set of conserved binding modes (CBMs) was created using the atomic detail of structure alignment data and the protein domain classification of the Conserved Domain Database. A conserved binding mode is inferred when different members of interacting domain families dock in the same way, such that their structural complexes superimpose well. Such domain interactions with recurring structural themes have greater significance to be biologically relevant, unlike spurious crystal packing interactions. Consequently, this study gives lower and upper bounds on the number of different types of interacting domain pairs in the structure database on the order of 1000-2000. We use CBMs to create domain interaction networks, which highlight functionally significant connections by avoiding many infrequent links between highly connected nodes. The CBMs also constitute a library of docking templates that may be used in molecular modeling to infer the characteristics of an unknown binding surface, just as conserved domains may be used to infer the structure of an unknown protein. The method's ability to sort through and classify large numbers of putative interacting domain pairs is demonstrated on the oligomeric interactions of globins.     

Steric Block High Affinity Oligonucleotide Analogues: A New Tool for Mapping RNA-Protein Binding Sites

Steric-block ON analogues are efficient inhibitors of RNA-protein interaction and therefore have potential to probe RNA sequences for putative protein binding sites and to investigate mechanisms of protein binding. The packaging process of HIV-1 is highly specific involving an interaction between the Gag protein and a conserved sequence that is only present on genomic viral RNA. Using oligonucleotide probes we have confirmed that the terminal purine loop is the major Gag binding site on SL3 and that a secondary Gag binding site exists at an internal purine bulge. We also demonstrate direct binding of oligonucleotide to their binding sites and confirm this interaction does not alter global RNA conformation, making them highly specific, nondisruptive probes of RNA protein interactions.     

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.     

Network representation of protein interactions: Theory of graph description and analysis

A methodological framework is presented for the graph theoretical interpretation of NMR data of protein interactions. The proposed analysis generalizes the idea of network representations of protein structures by expanding it to protein interactions. This approach is based on regularization of residue‐resolved NMR relaxation times and chemical shift data and subsequent construction of an adjacency matrix that represents the underlying protein interaction as a graph or network. The network nodes represent protein residues. Two nodes are connected if two residues are functionally correlated during the protein interaction event. The analysis of the resulting network enables the quantification of the importance of each amino acid of a protein for its interactions. Furthermore, the determination of the pattern of correlations between residues yields insights into the functional architecture of an interaction. This is of special interest for intrinsically disordered proteins, since the structural (three‐dimensional) architecture of these proteins and their complexes is difficult to determine. The power of the proposed methodology is demonstrated at the example of the interaction between the intrinsically disordered protein osteopontin and its natural ligand heparin.     

Evolutionary conservation and network structure characterize genes of phenotypic relevance for mitosis in human

The impact of gene silencing on cellular phenotypes is difficult to establish due to the complexity of interactions in the associated biological processes and pathways. A recent genome-wide RNA knock-down study both identified and phenotypically characterized a set of important genes for the cell cycle in HeLa cells. Here, we combine a molecular interaction network analysis, based on physical and functional protein interactions, in conjunction with evolutionary information, to elucidate the common biological and topological properties of these key genes. Our results show that these genes tend to be conserved with their corresponding protein interactions across several species and are key constituents of the evolutionary conserved molecular interaction network. Moreover, a group of bistable network motifs is found to be conserved within this network, which are likely to influence the network stability and therefore the robustness of cellular functioning. They form a cluster, which displays functional homogeneity and is significantly enriched in genes phenotypically relevant for mitosis. Additional results reveal a relationship between specific cellular processes and the phenotypic outcomes induced by gene silencing. This study introduces new ideas regarding the relationship between genotype and phenotype in the context of the cell cycle. We show that the analysis of molecular interaction networks can result in the identification of genes relevant to cellular processes, which is a promising avenue for future research.     

Bivalent Motif-Ear Interactions Mediate the Association of the Accessory Protein Tepsin with the AP-4 Adaptor Complex

The heterotetrameric (ϵ-β4-μ4-σ4) complex adaptor protein 4 (AP-4) is a component of a non-clathrin coat involved in protein sorting at the trans-Golgi network (TGN). Considerable interest in this complex has arisen from the recent discovery that mutations in each of its four subunits are the cause of a congenital intellectual disability and movement disorder in humans. Despite its physiological importance, the structure and function of this coat remain poorly understood. To investigate the assembly of the AP-4 coat, we dissected the determinants of interaction of AP-4 with its only known accessory protein, the ENTH/VHS-domain-containing protein tepsin. Using a variety of protein interaction assays, we found that tepsin comprises two phylogenetically conserved peptide motifs, [GS]LFXG[ML]X[LV] and S[AV]F[SA]FLN, within its C-terminal unstructured region, which interact with the C-terminal ear (or appendage) domains of the β4 and ϵ subunits of AP-4, respectively. Structure-based mutational analyses mapped the binding site for the [GS]LFXG[ML]X[LV] motif to a conserved, hydrophobic surface on the β4-ear platform fold. Both peptide-ear interactions are required for efficient association of tepsin with AP-4, and for recruitment of tepsin to the TGN. The bivalency of the interactions increases the avidity of tepsin for AP-4 and may enable cross-linking of multiple AP-4 heterotetramers, thus contributing to the assembly of the AP-4 coat. In addition to revealing critical aspects of this coat, our findings extend the paradigm of peptide-ear interactions, previously established for clathrin-AP-1/AP-2 coats, to a non-clathrin coat.     

Two types of transmembrane homomeric interactions in the integrin receptor family are evolutionarily conserved

Integrins are heterodimers, but recent in vitro and in vivo experiments suggest that they are also able to associate through their transmembrane domains to form homomeric interactions. Two fundamental questions are the biological relevance of these aggregates and their form of interaction in the membrane domain. Although in vitro experiments have shown the involvement of a GxxxG-like motif, several crosslinking in vivo data are consistent with an almost opposite form of interaction between the transmembrane alpha-helices. In the present work, we have explored these two questions using molecular dynamics simulations for all available integrin types. We have tested the hypothesis that homomeric interactions are evolutionary conserved, and essential for the cell, using conservative substitutions to filter out nonnative interactions. Our results show that two models, one involving a GxxxG-like motif (model I) and an almost opposite form of interaction (model II) are conserved across all alpha and beta integrin types, both in homodimers and homotrimers, with different specificities. No conserved interaction was found for homotetramers. Our results are completely independent from experimental data, both during molecular dynamics simulations and in the selection of the correct models. We rationalize previous seemingly conflicting findings regarding the nature of integrin interhelical homomeric interactions.     

The Role of Aromatic–Aromatic Interactions in Strand–Strand Stabilization of β-Sheets

Aromatic–aromatic interactions have long been believed to play key roles in protein structure, folding, and binding functions. However, we still lack full understanding of the contributions of aromatic–aromatic interactions to protein stability and the timing of their formation during folding. Here, using an aromatic ladder in the β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), as a case study, we find that aromatic π stacking plays a greater role in the Phe65–Phe71 cross-strand pair, while in another pair, Phe50–Phe65, hydrophobic interactions are dominant. The Phe65–Phe71 pair spans β-strands 4 and 5 in the β-barrel, which lack interstrand hydrogen bonding, and we speculate that it compensates energetically for the absence of strand–strand backbone interactions. Using perturbation analysis, we find that both aromatic–aromatic pairs form after the transition state for folding of CRABP1, thus playing a role in the final stabilization of the β-sheet rather than in its nucleation as had been earlier proposed. The aromatic interaction between strands 4 and 5 in CRABP1 is highly conserved in the intracellular lipid-binding protein (iLBP) family, and several lines of evidence combine to support a model wherein it acts to maintain barrel structure while allowing the dynamic opening that is necessary for ligand entry. Lastly, we carried out a bioinformatics analysis and found 51 examples of aromatic–aromatic interactions across non-hydrogen-bonded β-strands outside the iLBPs, arguing for the generality of the role played by this structural motif.     

Evolution of Protein Binding Modes in Homooligomers

The evolution of protein interactions cannot be deciphered without a detailed analysis of interaction interfaces and binding modes. We performed a large-scale study of protein homooligomers in terms of their symmetry, interface sizes, and conservation of binding modes. We also focused specifically on the evolution of protein binding modes from nine families of homooligomers and mapped 60 different binding modes and oligomerization states onto the phylogenetic trees of these families. We observed a significant tendency for the same binding modes to be clustered together and conserved within clades on phylogenetic trees; this trend is especially pronounced for close homologs with 70% sequence identity or higher. Some binding modes are conserved among very distant homologs, pointing to their ancient evolutionary origin, while others are very specific for a certain phylogenetic group. Moreover, we found that the most ancient binding modes have a tendency to involve symmetrical (isologous) homodimer binding arrangements with larger interfaces, while recently evolved binding modes more often exhibit asymmetrical arrangements and smaller interfaces.     

Alternate interfaces may mediate homomeric and heteromeric assembly in the transmembrane domains of SNARE proteins

The fusion of a vesicle to a target membrane is mediated by temporally and spatially regulated interactions within a set of evolutionarily conserved proteins. Integral to proper fusion is the interaction between proteins originating on both vesicle and target membranes to form a protein bridge between the two membranes, known as the SNARE complex. This protein complex includes the single-pass transmembrane helix proteins: syntaxin and synaptobrevin. Experimental data and amino acid sequence analysis suggest that an interface of interaction is conserved between the transmembrane regions of the two proteins. However, conflicting reports have been presented on the role of the synaptobrevin transmembrane domain in mediating important protein-protein interactions. To address this question, a thermodynamic study was carried out to determine quantitatively the self-association propensities of the transmembrane domains of synaptobrevin and syntaxin. Our results show that the transmembrane domain of synaptobrevin has only a modest ability to self-associate, whereas the transmembrane domain of syntaxin is able to form stable homodimers. Nevertheless, by a single amino acid substitution, synaptobrevin can be driven to dimerize with the same affinity as syntaxin. Furthermore, crosslinking studies show that dimerization of synaptobrevin is promoted by oxidizing agents. Despite the presence of a conserved cysteine residue in the same location as in synaptobrevin, syntaxin dimerization is not promoted by oxidization. This analysis suggests that subtle yet distinct differences are present between the two transmembrane dimer interfaces. A syntaxin/synaptobrevin heterodimer is able to form under oxidizing conditions, and we propose that the interface of interaction for the heterodimer may resemble the homodimer interface formed by the synaptobrevin transmembrane domain. Computational analysis of the transmembrane sequences of syntaxin and synaptobrevin reveal structural models that correlate with the experimental data. These data may provide insight into the role of transmembrane segments in the mechanism of vesicle fusion.     

Sequence and structural analysis of cellular retinoic acid-binding proteins reveals a network of conserved hydrophobic interactions

Proteins in the intracellular lipid-binding protein (iLBP) family show remarkably high structural conservation despite their low-sequence identity. A multiple-sequence alignment using 52 sequences of iLBP family members revealed 15 fully conserved positions, with a disproportionately high number of these (n=7) located in the relatively small helical region. The conserved positions displayed high structural conservation based on comparisons of known iLBP crystal structures. It is striking that the beta-sheet domain had few conserved positions, despite its high structural conservation. This observation prompted us to analyze pair-wise interactions within the beta-sheet region to ask whether structural information was encoded in interacting amino acid pairs. We conducted this analysis on the iLBP family member, cellular retinoic acid-binding protein I (CRABP I), whose folding mechanism is under study in our laboratory. Indeed, an analysis based on a simple classification of hydrophobic and polar amino acids revealed a network of conserved interactions in CRABP I that cluster spatially, suggesting a possible nucleation site for folding. Significantly, a small number of residues participated in multiple conserved interactions, suggesting a key role for these sites in the structure and folding of CRABP I. The results presented here correlate well with available experimental evidence on folding of CRABPs and their family members and suggest future experiments. The analysis also shows the usefulness of considering pair-wise conservation based on a simple classification of amino acids, in analyzing sequences and structures to find common core regions among homologues.     

富含半胱氨酸的GASA小分子蛋白研究进展

GASA蛋白是植物特有的一类富含半胱氨酸的小分子蛋白, 大多定位于细胞壁, 在植物生长发育和激素信号转导过程中发挥重要作用。该蛋白具有富含12个半胱氨酸残基的GASA结构域, 该结构域被认为是GASA蛋白维持空间结构和发挥功能的关键区域。该文重点综述了植物GASA蛋白的分子结构、亚细胞定位和生物学功能, 并对相关领域的研究进行了 展望。     

The NMR solution structure and function of RPA3313: a putative ribosomal transport protein from Rhodopseudomonas palustris

Protein function elucidation often relies heavily on amino acid sequence analysis and other bioinformatics approaches. The reliance is extended to structure homology modeling for ligand docking and protein–protein interaction mapping. However, sequence analysis of RPA3313 exposes a large, unannotated class of hypothetical proteins mostly from the Rhizobiales order. In the absence of sequence and structure information, further functional elucidation of this class of proteins has been significantly hindered. A high quality NMR structure of RPA3313 reveals that the protein forms a novel split ββαβ fold with a conserved ligand binding pocket between the first β‐strand and the N‐terminus of the α‐helix. Conserved residue analysis and protein–protein interaction prediction analyses reveal multiple protein binding sites and conserved functional residues. Results of a mass spectrometry proteomic analysis strongly point toward interaction with the ribosome and its subunits. The combined structural and proteomic analyses suggest that RPA3313 by itself or in a larger complex may assist in the transportation of substrates to or from the ribosome for further processing. Proteins 2016; 85:93–102. © 2016 Wiley Periodicals, Inc.     

A fluorescent reporter for mapping cellular protein‐protein interactions in time and space

We introduce a fluorescent reporter for monitoring protein–protein interactions in living cells. The method is based on the Split‐Ubiquitin method and uses the ratio of two auto‐fluorescent reporter proteins as signal for interaction (SPLIFF). The mating of two haploid yeast cells initiates the analysis and the interactions are followed online by two‐channel time‐lapse microscopy of the diploid cells during their first cell cycle. Using this approach we could with high spatio‐temporal resolution visualize the differences between the interactions of the microtubule binding protein Stu2p with two of its binding partners, monitor the transient association of a Ran‐GTPase with its receptors at the nuclear pore, and distinguish between protein interactions at the polar cortical domain at different phases of polar growth. These examples further demonstrate that protein–protein interactions identified from large‐scale screens can be effectively followed up by high‐resolution single‐cell analysis.