Discovering functions and revealing mechanisms at molecular level from biological networks

With the increasingly accumulated data from high-throughput technologies, study on biomolecular networks has become one of key focuses in systems biology and bioinformatics. In particular, various types of molecular networks (e.g., protein-protein interaction (PPI) network; gene regulatory network (GRN); metabolic network (MN); gene coexpression network (GCEN)) have been extensively investigated, and those studies demonstrate great potentials to discover basic functions and to reveal essential mechanisms for various biological phenomena, by understanding biological systems not at individual component level but at a system-wide level. Recent studies on networks have created very prolific researches on many aspects of living organisms. In this paper, we aim to review the recent developments on topics related to molecular networks in a comprehensive manner, with the special emphasis on the computational aspect. The contents of the survey cover global topological properties and local structural characteristics, network motifs, network comparison and query, detection of functional modules and network motifs, function prediction from network analysis, inferring molecular networks from biological data as well as representative databases and software tools.     

Influence of metabolic network structure and function on enzyme evolution

Background  

Most studies of molecular evolution are focused on individual genes and proteins. However, understanding the design principles and evolutionary properties of molecular networks requires a system-wide perspective. In the present work we connect molecular evolution on the gene level with system properties of a cellular metabolic network. In contrast to protein interaction networks, where several previous studies investigated the molecular evolution of proteins, metabolic networks have a relatively well-defined global function. The ability to consider fluxes in a metabolic network allows us to relate the functional role of each enzyme in a network to its rate of evolution.     

Analyzing protein-protein interaction networks

The advent of the "omics" era in biology research has brought new challenges and requires the development of novel strategies to answer previously intractable questions. Molecular interaction networks provide a framework to visualize cellular processes, but their complexity often makes their interpretation an overwhelming task. The inherently artificial nature of interaction detection methods and the incompleteness of currently available interaction maps call for a careful and well-informed utilization of this valuable data. In this tutorial, we aim to give an overview of the key aspects that any researcher needs to consider when working with molecular interaction data sets and we outline an example for interactome analysis. Using the molecular interaction database IntAct, the software platform Cytoscape, and its plugins BiNGO and clusterMaker, and taking as a starting point a list of proteins identified in a mass spectrometry-based proteomics experiment, we show how to build, visualize, and analyze a protein-protein interaction network.     

Interplay between gene expression noise and regulatory network architecture

Complex regulatory networks orchestrate most cellular processes in biological systems. Genes in such networks are subject to expression noise, resulting in isogenic cell populations exhibiting cell-to-cell variation in protein levels. Increasing evidence suggests that cells have evolved regulatory strategies to limit, tolerate or amplify expression noise. In this context, fundamental questions arise: how can the architecture of gene regulatory networks generate, make use of or be constrained by expression noise? Here, we discuss the interplay between expression noise and gene regulatory network at different levels of organization, ranging from a single regulatory interaction to entire regulatory networks. We then consider how this interplay impacts a variety of phenomena, such as pathogenicity, disease, adaptation to changing environments, differential cell-fate outcome and incomplete or partial penetrance effects. Finally, we highlight recent technological developments that permit measurements at the single-cell level, and discuss directions for future research.     

Coevolution by different functional mechanisms modulates the structure and dynamics of antagonistic and mutualistic networks

A central problem in the study of species interactions is to understand the underlying ecological and evolutionary mechanisms that shape and are shaped by trait evolution in interacting assemblages. The patterns of interaction among species (i.e. network structure) provide the pathways for evolution and coevolution, which are modulated by how traits affect individual fitness (i.e. functional mechanisms). Functional mechanisms, in turn, also affect the likelihood of an ecological interaction, shaping the structure of interaction networks. Here, we build adaptive network models to explore the potential role of coevolution by two functional mechanisms, trait matching and exploitation barrier, in driving trait evolution and the structure of interaction networks. We use these models to explore how different scenarios of coevolution and functional mechanisms reproduce the empirical network patterns observed in antagonistic and mutualistic interactions and affect trait evolution. Scenarios assuming coevolutionary feedback with a strong effect of functional mechanism better reproduce the empirical structure of networks. Antagonistic and mutualistic networks, however, are better explained by different functional mechanisms and the structure of antagonisms is better reproduced than that of mutualisms. Scenarios assuming coevolution by strong trait matching between interacting partners better explain the structure of antagonistic networks, whereas those assuming strong barrier effects better reproduce the structure of mutualistic networks. The dynamics resulting from the feedback between strong functional mechanisms and coevolution favor the stability of antagonisms and mutualisms. Selection favoring trait matching reduces temporal trait fluctuation and the magnitude of arms races in antagonisms, whereas selection due to exploitation barriers reduces temporal trait fluctuations in mutualisms. Our results indicate that coevolutionary models better reproduce the network structure of antagonisms than those of mutualisms and that different functional mechanisms may favor the persistence of antagonistic and mutualistic interacting assemblages.     

Generative probabilistic models for protein-protein interaction networks--the biclique perspective

MOTIVATION: Much of the large-scale molecular data from living cells can be represented in terms of networks. Such networks occupy a central position in cellular systems biology. In the protein-protein interaction (PPI) network, nodes represent proteins and edges represent connections between them, based on experimental evidence. As PPI networks are rich and complex, a mathematical model is sought to capture their properties and shed light on PPI evolution. The mathematical literature contains various generative models of random graphs. It is a major, still largely open question, which of these models (if any) can properly reproduce various biologically interesting networks. Here, we consider this problem where the graph at hand is the PPI network of Saccharomyces cerevisiae. We are trying to distinguishing between a model family which performs a process of copying neighbors, represented by the duplication-divergence (DD) model, and models which do not copy neighbors, with the Barabási-Albert (BA) preferential attachment model as a leading example. RESULTS: The observed property of the network is the distribution of maximal bicliques in the graph. This is a novel criterion to distinguish between models in this area. It is particularly appropriate for this purpose, since it reflects the graph's growth pattern under either model. This test clearly favors the DD model. In particular, for the BA model, the vast majority (92.9%) of the bicliques with both sides ≥4 must be already embedded in the model's seed graph, whereas the corresponding figure for the DD model is only 5.1%. Our results, based on the biclique perspective, conclusively show that a na?ve unmodified DD model can capture a key aspect of PPI networks.     

Computational models of signalling networks for non-linear control

Artificial signalling networks (ASNs) are a computational approach inspired by the signalling processes inside cells that decode outside environmental information. Using evolutionary algorithms to induce complex behaviours, we show how chaotic dynamics in a conservative dynamical system can be controlled. Such dynamics are of particular interest as they mimic the inherent complexity of non-linear physical systems in the real world. Considering the main biological interpretations of cellular signalling, in which complex behaviours and robust cellular responses emerge from the interaction of multiple pathways, we introduce two ASN representations: a stand-alone ASN and a coupled ASN. In particular we note how sophisticated cellular communication mechanisms can lead to effective controllers, where complicated problems can be divided into smaller and independent tasks.     

Wood‐inhabiting fungi with tight associations with other species have declined as a response to forest management

Research on mutualistic and antagonistic networks, such as plant–pollinator and host–parasite networks, has shown that species interactions can influence and be influenced by the responses of species to environmental perturbations. Here we examine whether results obtained for directly observable networks generalize to more complex networks in which species interactions cannot be observed directly. As a case study, we consider data on the occurrences of 98 wood‐inhabiting fungal species in managed and natural forests. We specifically ask if and how much the positions of wood‐inhabiting fungal species within the interaction networks influence their responses to forest management. For this, we utilize a joint species distribution model that partitions variation in species occurrences among environmental (i.e. resource availability) and biotic (i.e. species‐to‐species associations) predictors. Our results indicate that in addition to the direct loss of resource‐specialised species, forest management has indirect effects mediated through interactive associations. In particular, species with strong associative links to other species are especially sensitive to forest management.     

An agent based model for real-time signaling induced in osteocytic networks by mechanical stimuli

We recently observed that insertion of unloaded rest between each load cycle substantially enhanced bone formation induced by mild loading regimens. To begin to explore this result, we have developed an agent based model for real-time signaling induced when osteocytic networks are challenged by mechanical stimuli. In the model, activity induced in individual osteocytes were governed by the following cellular functions: (1) threshold levels of tissue strain magnitudes were required to initiate and maximally activate cells, (2) cell activity beyond thresholds were propagated within localized neighborhoods and influenced recipient cell activity, (3) cellular activity was modulated by 'molecular' stores and the rates at which stores were replenished when cells were quiescent. Using this model, the real-time response of osteocyte networks was determined as the average of individual cell activity. While not explicitly embedded within the model, interactions between cellular functions served as positive, negative, and end-point feedback mechanisms and resulted in unique real-time network responses to distinct mechanical stimuli. Specifically, the real-time network response to cyclic stimuli consisted of a large magnitude transient followed by low-level steady state fluctuations, while rest-inserted stimuli induced multiple secondary transients. Analysis of interaction patterns suggested that rest-inserted stimuli induced this enhanced and sustained signaling within osteocytic networks by enabling cell recovery of expended molecular stores and by efficiently utilizing properties inherent to cell-cell communication in bone. Importantly, this emergence based approach suggested mechanisms potentially underlying the benefit of rest-inserted stimuli and provides a unique framework for a broader exploration of mechanotransduction function within bone.     

A comparative social network analysis of wasp colonies and classrooms: Linking network structure to functioning

A major question in current network science is how to understand the relationship between structure and functioning of real networks. Here we present a comparative network analysis of 48 wasp and 36 human social networks. We have compared the centralisation and small world character of these interaction networks and have studied how these properties change over time. We compared the interaction networks of (1) two congeneric wasp species (Ropalidia marginata and Ropalidia cyathiformis), (2) the queen-right (with the queen) and queen-less (without the queen) networks of wasps, (3) the four network types obtained by combining (1) and (2) above, and (4) wasp networks with the social networks of children in 36 classrooms. We have found perfect (100%) centralisation in a queen-less wasp colony and nearly perfect centralisation in several other queen-less wasp colonies. Note that the perfectly centralised interaction network is quite unique in the literature of real-world networks. Differences between the interaction networks of the two wasp species are smaller than differences between the networks describing their different colony conditions. Also, the differences between different colony conditions are larger than the differences between wasp and children networks. For example, the structure of queen-right R. marginata colonies is more similar to children social networks than to that of their queen-less colonies. We conclude that network architecture depends more on the functioning of the particular community than on taxonomic differences (either between two wasp species or between wasps and humans).     

Strategies for dealing with incomplete information in the modeling of molecular interaction networks

Modelers of molecular interaction networks encounter the paradoxical situation that while large amounts of data are available, these are often insufficient for the formulation and analysis of mathematical models describing the network dynamics. In particular, information on the reaction mechanisms and numerical values of kinetic parameters are usually not available for all but a few well-studied model systems. In this article we review two strategies that have been proposed for dealing with incomplete information in the study of molecular interaction networks: parameter sensitivity analysis and model simplification. These strategies are based on the biologically justified intuition that essential properties of the system dynamics are robust against moderate changes in the value of kinetic parameters or even in the rate laws describing the interactions. Although advanced measurement techniques can be expected to relieve the problem of incomplete information to some extent, the strategies discussed in this article will retain their interest as tools providing an initial characterization of essential properties of the network dynamics.