Evolutionary diversification of protein-protein interactions by interface add-ons.
Ontology highlight
ABSTRACT:
SUBMITTER: Plach MG
PROVIDER: S-EPMC5635890 | biostudies-literature | 2017 Oct
REPOSITORIES: biostudies-literature
Publications
Evolutionary diversification of protein-protein interactions by interface add-ons.
Plach Maximilian G MG Semmelmann Florian F Busch Florian F Busch Markus M Heizinger Leonhard L Wysocki Vicki H VH Merkl Rainer R Sterner Reinhard R
Proceedings of the National Academy of Sciences of the United States of America 20170918 40
Cells contain a multitude of protein complexes whose subunits interact with high specificity. However, the number of different protein folds and interface geometries found in nature is limited. This raises the question of how protein-protein interaction specificity is achieved on the structural level and how the formation of nonphysiological complexes is avoided. Here, we describe structural elements called interface add-ons that fulfill this function and elucidate their role for the diversifica ...[more]
Similar Datasets
Project description:BackgroundDistinguishing biologically relevant interfaces from lattice contacts in protein crystals is a fundamental problem in structural biology. Despite efforts towards the computational prediction of interface character, many issues are still unresolved.ResultsWe present here a protein-protein interface classifier that relies on evolutionary data to detect the biological character of interfaces. The classifier uses a simple geometric measure, number of core residues, and two evolutionary indicators based on the sequence entropy of homolog sequences. Both aim at detecting differential selection pressure between interface core and rim or rest of surface. The core residues, defined as fully buried residues (>95% burial), appear to be fundamental determinants of biological interfaces: their number is in itself a powerful discriminator of interface character and together with the evolutionary measures it is able to clearly distinguish evolved biological contacts from crystal ones. We demonstrate that this definition of core residues leads to distinctively better results than earlier definitions from the literature. The stringent selection and quality filtering of structural and sequence data was key to the success of the method. Most importantly we demonstrate that a more conservative selection of homolog sequences - with relatively high sequence identities to the query - is able to produce a clearer signal than previous attempts.ConclusionsAn evolutionary approach like the one presented here is key to the advancement of the field, which so far was missing an effective method exploiting the evolutionary character of protein interfaces. Its coverage and performance will only improve over time thanks to the incessant growth of sequence databases. Currently our method reaches an accuracy of 89% in classifying interfaces of the Ponstingl 2003 datasets and it lends itself to a variety of useful applications in structural biology and bioinformatics. We made the corresponding software implementation available to the community as an easy-to-use graphical web interface at http://www.eppic-web.org.
| S-EPMC3556496 | biostudies-literature
Project description:EvoDesign (https://zhanglab.ccmb.med.umich.edu/EvoDesign) is an online server system for protein design. The method uses evolutionary profiles to guide the sequence search simulation and demonstrated significant advantages over physics-based approaches in terms of more accurately designing proteins that adopt desired target folds. Despite the success, the previous EvoDesign program focused only on monomer protein design, which limited its ability and usefulness in terms of designing functional proteins. In this work, we propose a new EvoDesign server, which extends the principles of evolution-based design to design protein-protein interactions. Starting from a two-chain complex structure, structurally similar interfaces are identified from known protein-protein interaction databases. An interface evolutionary profile is then constructed from a multiple sequence alignment of the interface analogies, which is combined with a newly developed, atomic-level physical energy function to guide the replica-exchange Monte Carlo simulation search. The purpose of the server is to redesign the specified complex chain to increase its stability and binding affinity for the other chain in the complex. With the improved scope and accuracy of the methodology, the new EvoDesign pipeline should become a useful online tool for functional protein design and drug discovery studies.
| S-EPMC6589126 | biostudies-literature
Project description:We investigate the influence of asymmetric interactions on coevolutionary dynamics of a predator-prey system by using the theory of adaptive dynamics. We assume that the defense ability of prey and the attack ability of predators all can adaptively evolve, either caused by phenotypic plasticity or by behavioral choice, but there are certain costs in terms of their growth rate or death rate. The coevolutionary model is constructed from a deterministic approximation of random mutation-selection process. To sum up, if prey's trade-off curve is globally weakly concave, then five outcomes of coevolution are demonstrated, which depend on the intensity and shape of asymmetric predator-prey interactions and predator's trade-off shape. Firstly, we find that if there is a weakly decelerating cost and a weakly accelerating benefit for predator species, then evolutionary branching in the predator species may occur, but after branching further coevolution may lead to extinction of the predator species with a larger trait value. However, if there is a weakly accelerating cost and a weakly accelerating benefit for predator species, then evolutionary branching in the predator species is also possible and after branching the dimorphic predator can evolutionarily stably coexist with a monomorphic prey species. Secondly, if the asymmetric interactions become a little strong, then prey and predators will evolve to an evolutionarily stable equilibrium, at which they can stably coexist on a long-term timescale of evolution. Thirdly, if there is a weakly accelerating cost and a relatively strongly accelerating benefit for prey species, then evolutionary branching in the prey species is possible and the finally coevolutionary outcome contains a dimorphic prey and a monomorphic predator species. Fourthly, if the asymmetric interactions become more stronger, then predator-prey coevolution may lead to cycles in both traits and equilibrium population densities. The Red Queen dynamic is a possible outcome under asymmetric predator-prey interactions.
| S-EPMC5042472 | biostudies-literature
Project description:We define an interface-interaction network (IIN) to capture the specificity and competition between protein-protein interactions (PPI). This new type of network represents interactions between individual interfaces used in functional protein binding and thereby contains the detail necessary to describe the competition and cooperation between any pair of binding partners. Here we establish a general framework for the construction of IINs that merges computational structure-based interface assignment with careful curation of available literature. To complement limited structural data, the inclusion of biochemical data is critical for achieving the accuracy and completeness necessary to analyze the specificity and competition between the protein interactions. Firstly, this procedure provides a means to clarify the information content of existing data on purported protein interactions and to remove indirect and spurious interactions. Secondly, the IIN we have constructed here for proteins involved in clathrin-mediated endocytosis (CME) exhibits distinctive topological properties. In contrast to PPI networks with their global and relatively dense connectivity, the fragmentation of the IIN into distinctive network modules suggests that different functional pressures act on the evolution of its topology. Large modules in the IIN are formed by interfaces sharing specificity for certain domain types, such as SH3 domains distributed across different proteins. The shared and distinct specificity of an interface is necessary for effective negative and positive design of highly selective binding targets. Lastly, the organization of detailed structural data in a network format allows one to identify pathways of specific binding interactions and thereby predict effects of mutations at specific surfaces on a protein and of specific binding inhibitors, as we explore in several examples. Overall, the endocytosis IIN is remarkably complex and rich in features masked in the coarser PPI, and collects relevant detail of protein association in a readily interpretable format.
| S-EPMC3656101 | biostudies-literature
Project description:Proteins involved in interactions throughout the course of evolution tend to co-evolve and compensatory changes may occur in interacting proteins to maintain or refine such interactions. However, certain residue pair alterations may prove to be detrimental for functional interactions. Hence, determining co-evolutionary pairings that could be structurally or functionally relevant for maintaining the conservation of an inter-protein interaction is important. Inter-protein co-evolution analysis in several complexes utilizing multiple existing methodologies suggested that co-evolutionary pairings can occur in spatially proximal and distant regions in inter-protein interactions. Subsequently, the Co-Var (Correlated Variation) method based on mutual information and Bhattacharyya coefficient was developed, validated, and found to perform relatively better than CAPS and EV-complex. Interestingly, while applying the Co-Var measure and EV-complex program on a set of protein-protein interaction complexes, co-evolutionary pairings were obtained in interface and non-interface regions in protein complexes. The Co-Var approach involves determining high degree co-evolutionary pairings that include multiple co-evolutionary connections between particular co-evolved residue positions in one protein with multiple residue positions in the binding partner. Detailed analyses of high degree co-evolutionary pairings in protein-protein complexes involved in cancer metastasis suggested that most of the residue positions forming such co-evolutionary connections mainly occurred within functional domains of constituent proteins and substitution mutations were also common among these positions. The physiological relevance of these predictions suggested that Co-Var can predict residues that could be crucial for preserving functional protein-protein interactions. Finally, Co-Var web server (http://www.hpppi.iicb.res.in/ishi/covar/index.html) that implements this methodology identifies co-evolutionary pairings in intra and inter-protein interactions.
| S-EPMC8271121 | biostudies-literature
Project description:This study measured the rates at which duplicated genes in the Daphnia pulex genome diverge in function, based on significant differences in their gene expression when paralogs of various ages were tested across 8 to 12 contrasting conditions. The conditions were designated A to L. The number and diversity of experiments were designed to uncover condition specific expression. The age of duplicated genes was measured as the level of silent-site nucleotide substitution (Ks) between paralogs. Two methods were used to evaluate the degree to which expression pattern of paralogs differ. (1) The first method investigated variation in the differential expression patterns among duplicates of individual gene families by calculating the M values (log2 treatment – log2 reference) from eight experiments, then calculating Pearson correlations. (2) The second method is a global investigation which defined a distinguishable expression pattern by a significance criterion (p < 0.05) using ANOVA for the simple statistical model of aov( Yab ~ Xab ), for matrices Yab differential expression M values and Xab gene factors, with replicates. We used pr(M) < 0.05 as criterion that expression differs between paralog genes, for zero to twelve treatments. The hypothesis that is tested is one of how many paralog pairs in each Ks category reach the criterion of a distinguishable expression pattern, which is tested for significance with Fisher's exact test for count data. This method is reliable for as few as two probes for one gene and one probe for the other, although a greater number of replicate probes produced more significant results. We compared the gene expression patterns for condition 1 versus condition 2 using four to six biological replicates, including balanced dye-swaps. All probes including random probes were quantile-normalized across chips, subarrays, samples and replicates.
2011-02-01 | E-GEOD-25850 | biostudies-arrayexpress
Project description:This study measured the rates at which duplicated genes in the Daphnia pulex genome diverge in function, based on significant differences in their gene expression when paralogs of various ages were tested across 8 to 12 contrasting conditions. The conditions were designated A to L. The number and diversity of experiments were designed to uncover condition-specific expression. The age of duplicated genes was measured as the level of silent-site nucleotide substitution (Ks) between paralogs. Two methods were used to evaluate the degree to which expression pattern of paralogs differ. (1) The first method investigated variation in the differential expression patterns among duplicates of individual gene families by calculating the M values (log2 treatment – log2 reference) from eight experiments, then calculating Pearson correlations. (2) The second method is a global investigation which defined a distinguishable expression pattern by a significance criterion (p < 0.05) using ANOVA for the simple statistical model of aov( Yab ~ Xab ), for matrices Yab differential expression M values and Xab gene factors, with replicates. We used pr(M) < 0.05 as criterion that expression differs between paralog genes, for zero to twelve treatments. The hypothesis that is tested is one of how many paralog pairs in each Ks category reach the criterion of a distinguishable expression pattern, which is tested for significance with Fisher's exact test for count data. This method is reliable for as few as two probes for one gene and one probe for the other, although a greater number of replicate probes produced more significant results. We compared the gene expression patterns for condition 1 versus condition 2 using four to six biological replicates, including balanced dye-swaps. All probes including random probes were quantile-normalized across chips, subarrays, samples and replicates.
2011-02-01 | E-GEOD-25841 | biostudies-arrayexpress
Project description:The sperm protein associated with nucleus in the X chromosome (SPANX) genes cluster at Xq27 in two subfamilies, SPANX-A/D and SPANX-N. SPANX-A/D is specific for hominoids and is fairly well characterized. The SPANX-N gave rise to SPANX-A/D in the hominoid lineage approximately 7 MYA. Given the proposed role of SPANX genes in spermatogenesis, we have extended studies to SPANX-N gene evolution, variation, regulation of expression, and intra-sperm localization. By immunofluorescence analysis, SPANX-N proteins are localized in post-meiotic spermatids exclusively, like SPANX-A/D. But in contrast to SPANX-A/D, SPANX-N are found in all ejaculated spermatozoa rather than only in a subpopulation, are localized in the acrosome rather than in the nuclear envelope, and are expressed at a low level in several nongametogenic adult tissues as well as many cancers. Presence of a binding site for CTCF and its testis-specific paralogue BORIS in the SPANX promoters suggests, by analogy to MAGE-A1 and NY-ESO-1, that their activation in spermatogenesis is mediated by the programmed replacement of CTCF by BORIS. Based on the relative density of CpG, the more extended expression of SPANX-N compared to SPANX-A/D in nongametogenic tissues is likely attributed to differences in promoter methylation. Our findings suggest that the recent duplication of SPANX genes in hominoids was accompanied by different localization of SPANX-N proteins in post-meiotic sperm and additional expression in several nongonadal tissues. This suggests a corresponding functional diversification of SPANX gene families in hominoids. SPANX proteins thus provide unique targets to investigate their roles in the function of spermatozoa, selected malignancies, and for SPANX-N, in other tissues as well.
| S-EPMC1831492 | biostudies-literature
Project description:This study measured the rates at which duplicated genes in the Daphnia pulex genome diverge in function, based on significant differences in their gene expression when paralogs of various ages were tested across 8 to 12 contrasting conditions. The conditions were designated A to L. The number and diversity of experiments were designed to uncover condition specific expression. The age of duplicated genes was measured as the level of silent-site nucleotide substitution (Ks) between paralogs. Two methods were used to evaluate the degree to which expression pattern of paralogs differ. (1) The first method investigated variation in the differential expression patterns among duplicates of individual gene families by calculating the M values (log2 treatment – log2 reference) from eight experiments, then calculating Pearson correlations. (2) The second method is a global investigation which defined a distinguishable expression pattern by a significance criterion (p < 0.05) using ANOVA for the simple statistical model of "aov( Yab ~ Xab )", for matrices Yab differential expression M values and Xab gene factors, with replicates. We used pr(M) < 0.05 as criterion that expression differs between paralog genes, for zero to twelve treatments. The hypothesis that is tested is one of how many paralog pairs in each Ks category reach the criterion of a distinguishable expression pattern, which is tested for significance with Fisher's exact test for count data. This method is reliable for as few as two probes for one gene and one probe for the other, although a greater number of replicate probes produced more significant results. We compared the gene expression patterns for condition 1 versus condition 2 using four to six biological replicates, including balanced dye-swaps. All probes including random probes were quantile-normalized across chips, subarrays, samples and replicates.
2011-02-01 | E-GEOD-25845 | biostudies-arrayexpress
Project description:This study measured the rates at which duplicated genes in the Daphnia pulex genome diverge in function, based on significant differences in their gene expression when paralogs of various ages were tested across 8 to 12 contrasting conditions. The conditions were designated A to L. The number and diversity of experiments were designed to uncover condition specific expression. The age of duplicated genes was measured as the level of silent-site nucleotide substitution (Ks) between paralogs. Two methods were used to evaluate the degree to which expression pattern of paralogs differ. (1) The first method investigated variation in the differential expression patterns among duplicates of individual gene families by calculating the M values (log2 treatment – log2 reference) from eight experiments, then calculating Pearson correlations. (2) The second method is a global investigation which defined a distinguishable expression pattern by a significance criterion (p < 0.05) using ANOVA for the simple statistical model of aov( Yab ~ Xab ), for matrices Yab differential expression M values and Xab gene factors, with replicates. We used pr(M) < 0.05 as criterion that expression differs between paralog genes, for zero to twelve treatments. The hypothesis that is tested is one of how many paralog pairs in each Ks category reach the criterion of a distinguishable expression pattern, which is tested for significance with Fisher's exact test for count data. This method is reliable for as few as two probes for one gene and one probe for the other, although a greater number of replicate probes produced more significant results. We compared the gene expression patterns for condition 1 versus condition 2 using four to six biological replicates, including balanced dye-swaps. All probes including random probes were quantile-normalized across chips, subarrays, samples and replicates.
2011-02-01 | E-GEOD-25851 | biostudies-arrayexpress