Resilience of protein-protein interaction networks as determined by their large-scale topological features

The relationship between the structure and function of biological networks constitutes a fundamental issue in systems biology. Particularly, the structure of protein-protein interaction networks is related to important biological functions. In this work, we investigated how such a resilience is determined by the large scale features of the respective networks. Four species are taken into account, namely yeast Saccharomyces cerevisiae, worm Caenorhabditis elegans, fly Drosophila melanogaster and Homo sapiens. We adopted two entropy-related measurements (degree entropy and dynamic entropy) in order to quantify the overall degree of robustness of these networks. We verified that while they exhibit similar structural variations under random node removal, they differ significantly when subjected to intentional attacks (hub removal). As a matter of fact, more complex species tended to exhibit more robust networks. More specifically, we quantified how six important measurements of the networks topology (namely clustering coefficient, average degree of neighbors, average shortest path length, diameter, assortativity coefficient, and slope of the power law degree distribution) correlated with the two entropy measurements. Our results revealed that the fraction of hubs and the average neighbor degree contribute significantly for the resilience of networks. In addition, the topological analysis of the removed hubs indicated that the presence of alternative paths between the proteins connected to hubs tend to reinforce resilience. The performed analysis helps to understand how resilience is underlain in networks and can be applied to the development of protein network models.     

Topographic Hub Maps of the Human Structural Neocortical Network

Hubs within the neocortical structural network determined by graph theoretical analysis play a crucial role in brain function. We mapped neocortical hubs topographically, using a sample population of 63 young adults. Subjects were imaged with high resolution structural and diffusion weighted magnetic resonance imaging techniques. Multiple network configurations were then constructed per subject, using random parcellations to define the nodes and using fibre tractography to determine the connectivity between the nodes. The networks were analysed with graph theoretical measures. Our results give reference maps of hub distribution measured with betweenness centrality and node degree. The loci of the hubs correspond with key areas from known overlapping cognitive networks. Several hubs were asymmetrically organized across hemispheres. Furthermore, females have hubs with higher betweenness centrality and males have hubs with higher node degree. Female networks have higher small-world indices.     

Network analysis of protein dynamics

The network paradigm is increasingly used to describe the topology and dynamics of complex systems. Here, we review the results of the topological analysis of protein structures as molecular networks describing their small-world character, and the role of hubs and central network elements in governing enzyme activity, allosteric regulation, protein motor function, signal transduction and protein stability. We summarize available data how central network elements are enriched in active centers and ligand binding sites directing the dynamics of the entire protein. We assess the feasibility of conformational and energy networks to simplify the vast complexity of rugged energy landscapes and to predict protein folding and dynamics. Finally, we suggest that modular analysis, novel centrality measures, hierarchical representation of networks and the analysis of network dynamics will soon lead to an expansion of this field.     

Why do hubs in the yeast protein interaction network tend to be essential: reexamining the connection between the network topology and essentiality

The centrality-lethality rule, which notes that high-degree nodes in a protein interaction network tend to correspond to proteins that are essential, suggests that the topological prominence of a protein in a protein interaction network may be a good predictor of its biological importance. Even though the correlation between degree and essentiality was confirmed by many independent studies, the reason for this correlation remains illusive. Several hypotheses about putative connections between essentiality of hubs and the topology of protein-protein interaction networks have been proposed, but as we demonstrate, these explanations are not supported by the properties of protein interaction networks. To identify the main topological determinant of essentiality and to provide a biological explanation for the connection between the network topology and essentiality, we performed a rigorous analysis of six variants of the genomewide protein interaction network for Saccharomyces cerevisiae obtained using different techniques. We demonstrated that the majority of hubs are essential due to their involvement in Essential Complex Biological Modules, a group of densely connected proteins with shared biological function that are enriched in essential proteins. Moreover, we rejected two previously proposed explanations for the centrality-lethality rule, one relating the essentiality of hubs to their role in the overall network connectivity and another relying on the recently published essential protein interactions model.     

Structure networks of E. coli glutaminyl-tRNA synthetase: effects of ligand binding

It is well known that proteins undergo backbone as well as side chain conformational changes upon ligand binding, which is not necessarily confined to the active site. Both the local and the global conformational changes brought out by ligand-binding have been extensively studied earlier. However, the global changes have been reported mainly at the protein backbone level. Here we present a method that explicitly takes into account the side chain interactions, yet providing a global view of the ligand-induced conformational changes. This is achieved through the analysis of Protein Structure Networks (PSN), constructed from the noncovalent side chain interactions in the protein. Here, E. coli Glutaminyl-tRNA synthetase (GlnRS) in the ligand-free and different ligand-bound states is used as a case study to assess the effect of binding of tRNA, ATP, and the amino acid Gln to GlnRS. The PSNs are constructed on the basis of the strength of noncovalent interactions existing between the side chains of amino acids. The parameters like the size of the largest cluster, edge to node ratio, and the total number of hubs are used to quantitatively assess the structure network changes. These network parameters have effectively captured the ligand-induced structural changes at a global structure network level. Hubs, the highly connected amino acids, are also identified from these networks. Specifically, we are able to characterize different types of hubs based on the comparison of structure networks of the GlnRS system. The differences in the structure networks in both the presence and the absence of the ligands are reflected in these hubs. For instance, the characterization of hubs that are present in both the ligand-free and all the ligand-bound GlnRS (the invariant hubs) might implicate their role in structural integrity. On the other hand, identification of hubs unique to a particular ligand-bound structure (the exclusive hubs) not only highlights the structural differences mediated by ligand-binding at the structure network level, but also highlights significance of these amino acids hubs in binding to the ligand and catalyzing the biochemical function. Further, the hubs identified from this study could be ideal targets for mutational studies to ascertain the ligand-induced structure-function relationships in E. coli GlnRS. The formalism used in this study is simple and can be applied to other protein-ligands in general to understand the allosteric changes mediated by the binding of ligands.     

Network structure, and vaccination strategy and effort interact to affect the dynamics of influenza epidemics

There is growing interest in understanding and controlling the spread of diseases through realistically structured host populations. We investigate how network structures, ranging from circulant, through small-world networks, to random networks, and vaccination strategy and effort interact to influence the proportion of the population infected, the size and timing of the epidemic peak, and the duration of the epidemic. We found these three factors, and their higher-order interactions, significantly influenced epidemic development and extent. Increasing vaccination effort (from 0% to 90%) decreased the number of hosts infected while increasing network randomness worked to increase disease spread. On average, vaccinating hosts based on degree (hubs) resulted in the smallest epidemics while vaccinating hosts with the highest clustering coefficient resulted in the largest epidemics. In a targeted test of five vaccination strategies on a small-world network (probability of rewiring edges=5%) with 10% vaccination effort we found that vaccinating hosts preferentially with high-clustering coefficients (similar to real-world strategies) resulted in twice the number of hosts infected as random vaccinations and nearly a 30-fold higher number of cases than our strategy targeting hubs (highest degree hosts). Our model suggests how vaccinations might be implemented to minimize the extent of an epidemic (e.g., duration and total number infected) as well as the timing and number of hosts infected at a given time over a wide range of structured host networks.     

Oligomeric protein structure networks: insights into protein-protein interactions

Background  

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

Phylogenetic imputation of plant functional trait databases

Continental‐scale maps of plant functional diversity are a fundamental piece of data of interest to ecosystem modelers and ecologists, yet such maps have been exceedingly hard to generate. The large effort to compile global plant functional trait databases largely for the purpose of mapping and analyzing the spatial distribution of function has resulted in very sparse data matrices thereby limiting progress. Identifying robust methodologies to gap fill or impute trait values in these databases is an important objective. Here I argue that existing statistical tools from phylogenetic comparative methods can be used to rapidly impute values into global plant functional trait databases due to the large amount of phylogenetic signal often in trait data. In particular, statistical models of phylogenetic signal in traits can be generated from existing data and used to predict missing values of closely related species often with a high degree of accuracy thereby facilitating the continental‐scale mapping of plant function. Despite the promise of this approach, I also discuss potential pitfalls and future challenges that will need to be addressed.     

The dichotomy in degree correlation of biological networks

Most complex networks from different areas such as biology, sociology or technology, show a correlation on node degree where the possibility of a link between two nodes depends on their connectivity. It is widely believed that complex networks are either disassortative (links between hubs are systematically suppressed) or assortative (links between hubs are enhanced). In this paper, we analyze a variety of biological networks and find that they generally show a dichotomous degree correlation. We find that many properties of biological networks can be explained by this dichotomy in degree correlation, including the neighborhood connectivity, the sickle-shaped clustering coefficient distribution and the modularity structure. This dichotomy distinguishes biological networks from real disassortative networks or assortative networks such as the Internet and social networks. We suggest that the modular structure of networks accounts for the dichotomy in degree correlation and vice versa, shedding light on the source of modularity in biological networks. We further show that a robust and well connected network necessitates the dichotomy of degree correlation, suggestive of an evolutionary motivation for its existence. Finally, we suggest that a dichotomous degree correlation favors a centrally connected modular network, by which the integrity of network and specificity of modules might be reconciled.     

Light signaling induces anthocyanin biosynthesis via AN3 mediated COP1 expression

Light signaling plays a pivotal role in controlling plant morphogenesis, metabolism, growth and development. The central process of light signaling pathway is to build the link between light signals and the expression of genes involved. Although studies focused on light signaling toward metabolism have been documented well in the past several decades, most regulation networks of light signaling in a specific metabolic production largely remained unknown. Anthocyanin accumulation in plant tissues depends on the availability of light signals, but only little is known about the potential regulation network underlying light signal controls anthocyanin biosynthesis. Here, we briefly review the recent progress on the light-triggered anthocyanin biosynthesis via ANGUSTIFOLIA3 (AN3) and CONSTITUTIVE PHOTOMORPHOGENIC1 (COP1) network in Arabidopsis.     

Discrete domains mediate the light-responsive nuclear and cytoplasmic localization of Arabidopsis COP1

The Arabidopsis CONSTITUTIVE PHOTOMORPHOGENIC1 (COP1) protein plays a critical role in the repression of photomorphogenesis during Arabidopsis seedling development. We investigated the control of COP1 partitioning between nucleus and cytoplasm, which has been implicated in the regulation of COP1 activity, by using fusion proteins between COP1 and beta-glucuronidase or the green fluorescent protein. Transient expression assays using onion epidermal cells and data from hypocotyl cells of stably transformed Arabidopsis demonstrated that COP1 carries a single, bipartite nuclear localization signal that functions independently of light. Nuclear exclusion was mediated by a novel and distinct signal, bordering the zinc-finger and coiled-coil motifs, that was able to redirect a heterologous nuclear protein to the cytoplasm. The cytoplasmic localization signal functioned in a light-independent manner. Light regulation of nuclear localization was reconstituted by combining the individual domains containing the nuclear localization signal and the cytoplasmic localization signal; the WD-40 repeat domain of COP1 was not required. However, phenotypic analysis of transgenic seedlings suggested that the constitutively nuclear-localized WD-40 repeat domain was able to mimic aspects of COP1 function, as indicated by exaggerated hypocotyl elongation under light conditions.     

Dynamical analysis of yeast protein interaction network during the sake brewing process

Proteins interact with each other for performing essential functions of an organism. They change partners to get involved in various processes at different times or locations. Studying variations of protein interactions within a specific process would help better understand the dynamic features of the protein interactions and their functions. We studied the protein interaction network of Saccharomyces cerevisiae (yeast) during the brewing of Japanese sake. In this process, yeast cells are exposed to several stresses. Analysis of protein interaction networks of yeast during this process helps to understand how protein interactions of yeast change during the sake brewing process. We used gene expression profiles of yeast cells for this purpose. Results of our experiments revealed some characteristics and behaviors of yeast hubs and non-hubs and their dynamical changes during the brewing process. We found that just a small portion of the proteins (12.8 to 21.6%) is responsible for the functional changes of the proteins in the sake brewing process. The changes in the number of edges and hubs of the yeast protein interaction networks increase in the first stages of the process and it then decreases at the final stages.     

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.     

Protein chips for high-throughput doping screening in athletes

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.     

Binding similarity network of ligand

The protein and ligand interaction takes an important part in protein function. Both ligand and its binding site are essential components for understanding how the protein-ligand complex functions. Until now, there have been many studies about protein function and evolution, but they usually lacked ligand information. Accordingly, in this study, we tried to answer the following questions: how much ligand and binding site are associated with protein function, and how ligands themselves are related to each other in terms of binding site. To answer the questions, we presented binding similarity network of ligand. Through the network analysis, we attempted to reveal systematic relationship between the ligand and binding site. The results showed that ligand binding site and function were closely related (conservation ratio, 81%). We also showed conservative tendency of function in line with ligand structure similarity with some exceptional cases. In addition, the binding similarity network of ligand revealed scale-free property to some degree like other biological networks. Since most nodes formed highly connected cluster, a clustering coefficient was very high compared with random. All the highly connected ligands (hubs) were involved in various functions forming large cluster and tended to act as a bridge between modular clusters in the network.     

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.