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.     

Role of intrinsic disorder in transient interactions of hub proteins

Hubs in the protein-protein interaction network have been classified as "party" hubs, which are highly correlated in their mRNA expression with their partners while "date" hubs show lesser correlation. In this study, we explored the role of intrinsic disorder in date and party hub interactions. The data reveals that intrinsic disorder is significantly enriched in date hub proteins when compared with party hub proteins. Intrinsic disorder has been largely implicated in transient binding interactions. The disorder to order transition, which occurs during binding interactions in disordered regions, renders the interaction highly reversible while maintaining the high specificity. The enrichment of intrinsic disorder in date hubs may facilitate transient interactions, which might be required for date hubs to interact with different partners at different times.     

Evolutionary and physiological importance of hub proteins

It has been claimed that proteins with more interaction partners (hubs) are both physiologically more important (i.e., less dispensable) and, owing to an assumed high density of binding sites, slow evolving. Not all analyses, however, support these results, probably because of biased and less-than reliable global protein interaction data. Here we provide the first examination of these issues using a comprehensive literature-curated dataset of well-substantiated protein interactions in Saccharomyces cerevisiae. Whereas use of less reliable yeast two-hybrid data alone can reject the possibility that local connectivity correlates with measures of dispensability, in higher quality datasets a relatively robust correlation is observed. In contrast, local connectivity does not correlate with the rate of protein evolution even in reliable datasets. This perhaps surprising lack of correlation with evolutionary rate appears in part to arise from the fact that hub proteins do not have a higher density of residues associated with binding. However, hub proteins do have at least one other set of unusual features, namely rapid turnover and regulation, as manifest in high mRNA decay rates and a large number of phosphorylation sites. This, we suggest, is an adaptation to minimize unwanted activation of pathways that might be mediated by adventitious binding to hubs, were they to actively persist longer than required at any given time point. We conclude that hub proteins are more important for cellular growth rate and under tight regulation but are not slow evolving.     

Identification of hub proteins from sequence

Identification of hub proteins from sequence is a challenge in molecular biology. Therefore, it is of interest to predict protein hubs in networks. We describe the prediction of protein "hub" using physiochemical, thermodynamic and conformational properties of amino acid residues in sequence. We have used twenty sequence based features to identify hub behaviour. Linear discriminant analysis and normalised Bayesian approach were utilized for identifying hub proteins solely using these sequence features in E. coli/H. sapiens datasets with accuracies of 99.5/98.6, 87.8/89.6 and 90.1/92.6, respectively.     

Disorder and sequence repeats in hub proteins and their implications for network evolution

Protein interaction networks display approximate scale-free topology, in which hub proteins that interact with a large number of other proteins determine the overall organization of the network. In this study, we aim to determine whether hubs are distinguishable from other networked proteins by specific sequence features. Proteins of different connectednesses were compared in the interaction networks of Saccharomyces cerevisiae, Drosophila melanogaster, Caenorhabditis elegans, and Homo sapienswith respect to the distribution of predicted structural disorder, sequence repeats, low complexity regions, and chain length. Highly connected proteins ("hub proteins") contained significantly more of, and greater proportion of, these sequence features and tended to be longer overall as compared to less connected proteins. These sequence features provide two different functional means for realizing multiple interactions: (1) extended interaction surface and (2) flexibility and adaptability, providing a mechanism for the same region to bind distinct partners. Our view contradicts the prevailing view that scaling in protein interactomes arose from gene duplication and preferential attachment of equivalent proteins. We propose an alternative evolutionary network specialization process, in which certain components of the protein interactome improved their fitness for binding by becoming longer or accruing regions of disorder and/or internal repeats and have therefore become specialized in network organization.     

Identification of transient hub proteins and the possible structural basis for their multiple interactions

Proteins that can interact with multiple partners play central roles in the network of protein-protein interactions. They are called hub proteins, and recently it was suggested that an abundance of intrinsically disordered regions on their surfaces facilitates their binding to multiple partners. However, in those studies, the hub proteins were identified as proteins with multiple partners, regardless of whether the interactions were transient or permanent. As a result, a certain number of hub proteins are subunits of stable multi-subunit proteins, such as supramolecules. It is well known that stable complexes and transient complexes have different structural features, and thus the statistics based on the current definition of hub proteins will hide the true nature of hub proteins. Therefore, in this paper, we first describe a new approach to identify proteins with multiple partners dynamically, using the Protein Data Bank, and then we performed statistical analyses of the structural features of these proteins. We refer to the proteins as transient hub proteins or sociable proteins, to clarify the difference with hub proteins. As a result, we found that the main difference between sociable and nonsociable proteins is not the abundance of disordered regions, in contrast to the previous studies, but rather the structural flexibility of the entire protein. We also found greater predominance of charged and polar residues in sociable proteins than previously reported.     

The relationship between relative solvent accessibility and evolutionary rate in protein evolution

Recent work with Saccharomyces cerevisiae shows a linear relationship between the evolutionary rate of sites and the relative solvent accessibility (RSA) of the corresponding residues in the folded protein. Here, we aim to develop a mathematical model that can reproduce this linear relationship. We first demonstrate that two models that both seem reasonable choices (a simple model in which selection strength correlates with RSA and a more complex model based on RSA-dependent amino acid distributions) fail to reproduce the observed relationship. We then develop a model on the basis of observed site-specific amino acid distributions and show that this model behaves appropriately. We conclude that evolutionary rates are directly linked to the distribution of amino acids at individual sites. Because of this link, any future insight into the biophysical mechanisms that determine amino acid distributions will improve our understanding of evolutionary rates.     

Stratus not altocumulus: a new view of the yeast protein interaction network

Systems biology approaches can reveal intermediary levels of organization between genotype and phenotype that often underlie biological phenomena such as polygenic effects and protein dispensability. An important conceptualization is the module, which is loosely defined as a cohort of proteins that perform a dedicated cellular task. Based on a computational analysis of limited interaction datasets in the budding yeast Saccharomyces cerevisiae, it has been suggested that the global protein interaction network is segregated such that highly connected proteins, called hubs, tend not to link to each other. Moreover, it has been suggested that hubs fall into two distinct classes: "party" hubs are co-expressed and co-localized with their partners, whereas "date" hubs interact with incoherently expressed and diversely localized partners, and thereby cohere disparate parts of the global network. This structure may be compared with altocumulus clouds, i.e., cotton ball-like structures sparsely connected by thin wisps. However, this organization might reflect a small and/or biased sample set of interactions. In a multi-validated high-confidence (HC) interaction network, assembled from all extant S. cerevisiae interaction data, including recently available proteome-wide interaction data and a large set of reliable literature-derived interactions, we find that hub-hub interactions are not suppressed. In fact, the number of interactions a hub has with other hubs is a good predictor of whether a hub protein is essential or not. We find that date hubs are neither required for network tolerance to node deletion, nor do date hubs have distinct biological attributes compared to other hubs. Date and party hubs do not, for example, evolve at different rates. Our analysis suggests that the organization of global protein interaction network is highly interconnected and hence interdependent, more like the continuous dense aggregations of stratus clouds than the segregated configuration of altocumulus clouds. If the network is configured in a stratus format, cross-talk between proteins is potentially a major source of noise. In turn, control of the activity of the most highly connected proteins may be vital. Indeed, we find that a fluctuation in steady-state levels of the most connected proteins is minimized.     

Evolution of the 14-3-3 Protein Family: Does the Large Number of Isoforms in Multicellular Organisms Reflect Functional Specificity?

14-3-3 proteins constitute a family of eukaryotic proteins that are key regulators of a large number of processes ranging from mitosis to apoptosis. 14-3-3s function as dimers and bind to particular motifs in their target proteins. To date, 14-3-3s have been implicated in regulation or stabilization of more than 35 different proteins. This number is probably only a fraction of the number of proteins that 14-3-3s bind to, as reports of new target proteins have become more frequent. An examination of 14-3-3 entries in the public databases reveals 153 isoforms, including alleloforms, reported in 48 different species. The number of isoforms range from 2, in the unicellular organism Saccharomyces cerevisiae, to 12 in the multicellular organism Arabidopsis thaliana. A phylogenetic analysis reveals that there are four major evolutionary lineages: Viridiplantae (plants), Fungi, Alveolata, and Metazoa (animals). A close examination of the aligned amino acid sequences identifies conserved amino acid residues and regions of importance for monomer stabilization, dimer formation, target protein binding, and the nuclear export function. Given the fact that 53% of the protein is conserved, including all amino acid residues in the target binding groove of the 14-3-3 monomer, one might expect little to no isoform specificity for target protein binding. However, using surface plasmon resonance we show that there are large differences in affinity between nine 14-3-3 isoforms of A. thaliana and a target peptide representing a novel binding motif present in the C terminus of the plant plasma membrane H⁺ATPase. Thus, our data suggest that one reason for the large number of isoforms found in multicellular organisms is isoform-specific functions. Received: 19 April 2000 / Accepted: 24 July 2000     

Evolution of specific protein-protein interaction sites following gene duplication

Gene duplication is a common evolutionary process that leads to the expansion and functional diversification of protein subfamilies. The evolutionary events that cause paralogous proteins to bind different protein ligands (functionally diverged interfaces) are investigated and compared to paralogous proteins that bind the same protein ligand (functionally preserved interfaces). We find that functionally diverged interfaces possess more subfamily-specific residues than functionally preserved interfaces. These subfamily-specific residues are usually partially buried at the interface rim and achieve specific binding through optimized hydrogen bond geometries. In addition to optimized hydrogen bond geometries, side-chain modeling experiments suggest that steric effects are also important for binding specificity. Residues that are completely buried at the interface hub are also less conserved in functionally diverged interfaces than in functionally preserved interfaces. Consistent with this finding, hub residues contribute less to free energy of binding in functionally diverged interfaces than in functionally preserved interfaces. Therefore, we propose that protein binding is a delicate balance between binding affinity that primarily occurs at the interface hub and binding specificity that primarily occurs at the interface rim.     

Markov models of amino acid substitution to study proteins with intrinsically disordered regions

Background

Intrinsically disordered proteins (IDPs) or proteins with disordered regions (IDRs) do not have a well-defined tertiary structure, but perform a multitude of functions, often relying on their native disorder to achieve the binding flexibility through changing to alternative conformations. Intrinsic disorder is frequently found in all three kingdoms of life, and may occur in short stretches or span whole proteins. To date most studies contrasting the differences between ordered and disordered proteins focused on simple summary statistics. Here, we propose an evolutionary approach to study IDPs, and contrast patterns specific to ordered protein regions and the corresponding IDRs.

Results

Two empirical Markov models of amino acid substitutions were estimated, based on a large set of multiple sequence alignments with experimentally verified annotations of disordered regions from the DisProt database of IDPs. We applied new methods to detect differences in Markovian evolution and evolutionary rates between IDRs and the corresponding ordered protein regions. Further, we investigated the distribution of IDPs among functional categories, biochemical pathways and their preponderance to contain tandem repeats.

Conclusions

We find significant differences in the evolution between ordered and disordered regions of proteins. Most importantly we find that disorder promoting amino acids are more conserved in IDRs, indicating that in some cases not only amino acid composition but the specific sequence is important for function. This conjecture is also reinforced by the observation that for of our data set IDRs evolve more slowly than the ordered parts of the proteins, while we still support the common view that IDRs in general evolve more quickly. The improvement in model fit indicates a possible improvement for various types of analyses e.g. de novo disorder prediction using a phylogenetic Hidden Markov Model based on our matrices showed a performance similar to other disorder predictors.     

Primary structures of pancreatic ribonucleases from Bovidae: Impala,Thomson's gazelle,nilgai and water buffalo

The amino acid sequences of the pancreatic ribonucleases from impala, Thomson's gazelle, nilgai and two types of water buffalo were studied from several enzymic digests. Peptides were positioned by homology with other ribonucleases of known sequence. Only peptides that differ in amino acid composition from the corresponding peptides of ox and goat ribonuclease were sequenced.The primary structures of pancreatic ribonucleases from 11 species of the Bovidae family are known to date. The evolutionary rate of bovid ribonucleases is 2–3 times lower than the average rate observed in all mammals. A possible explanation is that the presence of a stable symbiotic system in the rumen of grass-eating ruminants has caused a slowing down of the evolutionary rate of pancreatic ribonuclease in this taxon. The subfamily of the Bovinae (five species) exhibits a slightly higher evolutionary rate; most replacements in the Bovinae occur near residues 34–36, the most commonly observed carbohydrate attachment site in ribonucleases.     

Amino Acid Residues in the ω-Minus Region Participate in Cellular Localization of Yeast Glycosylphosphatidylinositol-Attached Proteins

The final destination of glycosylphosphatidylinositol (GPI)-attached proteins in Saccharomyces cerevisiae is the plasma membrane or the cell wall. Two kinds of signals have been proposed for their cellular localization: (i) the specific amino acid residues V, I, or L at the site 4 or 5 amino acids upstream of the GPI attachment site (the omega site) and Y or N at the site 2 amino acids upstream of the omega site for cell wall localization and (ii) dibasic residues in the region upstream of the omega site (the omega-minus region) for plasma membrane localization. The relationships between these amino acid residues and efficiencies of cell wall incorporation were examined by constructing fusion reporter proteins from open reading frames encoding putative GPI-attached proteins. The levels of incorporation were high in the constructs containing the specific amino acid residues and quite low in those containing two basic amino acid residues in the omega-minus region. With constructs that contained neither specific residues nor two basic residues, levels of incorporation were moderate. These correlations clearly suggest that GPI-attached proteins have two different signals which act positively or negatively in cell wall incorporation for their cellular localization.     

On the Role of Aggregation Prone Regions in Protein Evolution,Stability, and Enzymatic Catalysis: Insights from Diverse Analyses

The various roles that aggregation prone regions (APRs) are capable of playing in proteins are investigated here via comprehensive analyses of multiple non-redundant datasets containing randomly generated amino acid sequences, monomeric proteins, intrinsically disordered proteins (IDPs) and catalytic residues. Results from this study indicate that the aggregation propensities of monomeric protein sequences have been minimized compared to random sequences with uniform and natural amino acid compositions, as observed by a lower average aggregation propensity and fewer APRs that are shorter in length and more often punctuated by gate-keeper residues. However, evidence for evolutionary selective pressure to disrupt these sequence regions among homologous proteins is inconsistent. APRs are less conserved than average sequence identity among closely related homologues (≥80% sequence identity with a parent) but APRs are more conserved than average sequence identity among homologues that have at least 50% sequence identity with a parent. Structural analyses of APRs indicate that APRs are three times more likely to contain ordered versus disordered residues and that APRs frequently contribute more towards stabilizing proteins than equal length segments from the same protein. Catalytic residues and APRs were also found to be in structural contact significantly more often than expected by random chance. Our findings suggest that proteins have evolved by optimizing their risk of aggregation for cellular environments by both minimizing aggregation prone regions and by conserving those that are important for folding and function. In many cases, these sequence optimizations are insufficient to develop recombinant proteins into commercial products. Rational design strategies aimed at improving protein solubility for biotechnological purposes should carefully evaluate the contributions made by candidate APRs, targeted for disruption, towards protein structure and activity.     

Hubs with network motifs organize modularity dynamically in the protein-protein interaction network of yeast

Background

It has been recognized that modular organization pervades biological complexity. Based on network analysis, ‘party hubs’ and ‘date hubs’ were proposed to understand the basic principle of module organization of biomolecular networks. However, recent study on hubs has suggested that there is no clear evidence for coexistence of ‘party hubs’ and ‘date hubs’. Thus, an open question has been raised as to whether or not ‘party hubs’ and ‘date hubs’ truly exist in yeast interactome.

Methodology

In contrast to previous studies focusing on the partners of a hub or the individual proteins around the hub, our work aims to study the network motifs of a hub or interactions among individual proteins including the hub and its neighbors. Depending on the relationship between a hub''s network motifs and protein complexes, we define two new types of hubs, ‘motif party hubs’ and ‘motif date hubs’, which have the same characteristics as the original ‘party hubs’ and ‘date hubs’ respectively. The network motifs of these two types of hubs display significantly different features in spatial distribution (or cellular localizations), co-expression in microarray data, controlling topological structure of network, and organizing modularity.

Conclusion

By virtue of network motifs, we basically solved the open question about ‘party hubs’ and ‘date hubs’ which was raised by previous studies. Specifically, at the level of network motifs instead of individual proteins, we found two types of hubs, motif party hubs (mPHs) and motif date hubs (mDHs), whose network motifs display distinct characteristics on biological functions. In addition, in this paper we studied network motifs from a different viewpoint. That is, we show that a network motif should not be merely considered as an interaction pattern but be considered as an essential function unit in organizing modules of networks.     

Nature disfavors sequences of alternating polar and non-polar amino acids: implications for amyloidogenesis

Recent experiments with combinatorial libraries of de novo proteins have demonstrated that sequences designed to contain polar and non-polar amino acid residues arranged in an alternating pattern form fibrillar structures resembling beta-amyloid. This finding prompted us to probe the distribution of alternating patterns in the sequences of natural proteins. Analysis of a database of 250,514 protein sequences (79,708,024 residues) for all possible binary patterns of polar and non-polar amino acid residues revealed that alternating patterns occur significantly less often than other patterns with similar compositions. The under-representation of alternating binary patterns in natural protein sequences, coupled with the observation that such patterns promote amyloid-like structures in de novo proteins, suggests that sequences of alternating polar and non-polar amino acids are inherently amyloidogenic and consequently have been disfavored by evolutionary selection.     

Revisiting Date and Party Hubs: Novel Approaches to Role Assignment in Protein Interaction Networks

The idea of “date” and “party” hubs has been influential in the study of protein–protein interaction networks. Date hubs display low co-expression with their partners, whilst party hubs have high co-expression. It was proposed that party hubs are local coordinators whereas date hubs are global connectors. Here, we show that the reported importance of date hubs to network connectivity can in fact be attributed to a tiny subset of them. Crucially, these few, extremely central, hubs do not display particularly low expression correlation, undermining the idea of a link between this quantity and hub function. The date/party distinction was originally motivated by an approximately bimodal distribution of hub co-expression; we show that this feature is not always robust to methodological changes. Additionally, topological properties of hubs do not in general correlate with co-expression. However, we find significant correlations between interaction centrality and the functional similarity of the interacting proteins. We suggest that thinking in terms of a date/party dichotomy for hubs in protein interaction networks is not meaningful, and it might be more useful to conceive of roles for protein-protein interactions rather than for individual proteins.