J.C. Nacher

Observing metabolic functions at the genome scale

BackgroundHigh-throughput techniques have multiplied the amount and the types of available biological data, and for the first time achieving a global comprehension of the physiology of biological cells has become an achievable goal. This aim requires the integration of large amounts of heterogeneous data at different scales. It is notably necessary to extend the traditional focus on genomic data towards a truly functional focus, where the activity of cells is described in terms of actual metabolic processes performing the functions necessary for cells to live.ResultsIn this work, we present a new approach for metabolic analysis that allows us to observe the transcriptional activity of metabolic functions at the genome scale. These functions are described in terms of elementary modes, which can be computed in a genome-scale model thanks to a modular approach. We exemplify this new perspective by presenting a detailed analysis of the transcriptional metabolic response of yeast cells to stress. The integration of elementary mode analysis with gene expression data allows us to identify a number of functionally induced or repressed metabolic processes in different stress conditions. The assembly of these elementary modes leads to the identification of specific metabolic backbones.ConclusionThis study opens a new framework for the cell-scale analysis of metabolism, where transcriptional activity can be analyzed in terms of whole processes instead of individual genes. We furthermore show that the set of active elementary modes exhibits a highly uneven organization, where most of them conduct specialized tasks while a smaller proportion performs multi-task functions and dominates the general stress response.

Evolutionary history and functional implications of protein domains and their combinations in eukaryotes.

BackgroundIn higher multicellular eukaryotes, complex protein domain combinations contribute to various cellular functions such as regulation of intercellular or intracellular signaling and interactions. To elucidate the characteristics and evolutionary mechanisms that underlie such domain combinations, it is essential to examine the different types of domains and their combinations among different groups of eukaryotes.ResultsWe observed a large number of group-specific domain combinations in animals, especially in vertebrates. Examples include animal-specific combinations in tyrosine phosphorylation systems and vertebrate-specific combinations in complement and coagulation cascades. These systems apparently underwent extensive evolution in the ancestors of these groups. In extant animals, especially in vertebrates, animal-specific domains have greater connectivity than do other domains on average, and contribute to the varying number of combinations in each animal subgroup. In other groups, the connectivities of older domains were greater on average. To observe the global behavior of domain combinations during evolution, we traced the changes in domain combinations among animals and fungi in a network analysis. Our results indicate that there is a correlation between the differences in domain combinations among different phylogenetic groups and different global behaviors.ConclusionRapid emergence of animal-specific domains was observed in animals, contributing to specific domain combinations and functional diversification, but no such trends were observed in other clades of eukaryotes. We therefore suggest that the strategy for achieving complex multicellular systems in animals differs from that of other eukaryotes.

Correlation between structure and temperature in prokaryotic metabolic networks

BackgroundIn recent years, an extensive characterization of network structures has been made in an effort to elucidate design principles of metabolic networks, providing valuable insights into the functional organization and the evolutionary history of organisms. However, previous analyses have not discussed the effects of environmental factors (i.e., exogenous forces) in shaping network structures. In this work, we investigate the effect of temperature, which is one of the environmental factors that may have contributed to shaping structures of metabolic networks.ResultsFor this, we investigate the correlations between several structural properties characterized by graph metrics like the edge density, the degree exponent, the clustering coefficient, and the subgraph concentration in the metabolic networks of 113 prokaryotes and optimal growth temperature. As a result, we find that these structural properties are correlated with the optimal growth temperature. With increasing temperature, the edge density, the clustering coefficient and the subgraph concentration decrease and the degree exponent becomes large.ConclusionThis result implies that the metabolic networks transit with temperature as follows. The density of chemical reactions becomes low, the connectivity of the networks becomes homogeneous such as random networks and both the network modularity, based on the graph-theoretic clustering coefficient, and the frequency of recurring subgraphs decay. In short, metabolic networks undergo a change from heterogeneous and high-modular structures to homogeneous and low-modular structures, such as random networks, with temperature. This finding may suggest that the temperature plays an important role in the design principles of metabolic networks.

Recent Progress on the Analysis of Power-law Features in Complex Cellular Networks

Complex interactions between different kinds of bio-molecules and essential nutrients are responsible for cellular functions. Rapid advances in theoretical modeling and experimental analyses have shown that drastically different biological and non-biological networks share a common architecture. That is, the probability that a selected node in the network has exactly k edges decays as a power-law. This finding has definitely opened an intense research and debate on the origin and implications of this ubiquitous pattern. In this review, we describe the recent progress on the emergence of power-law distributions in cellular networks. We first show the internal characteristics of the observed complex networks uncovered using graph theory. We then briefly review some works that have significantly contributed to the theoretical analysis of cellular networks and systems, from metabolic and protein networks to gene expression profiles. This prevalent topology observed in so many diverse biological systems suggests the existence of generic laws and organizing principles behind the cellular networks.

Two complementary representations of a scale-free network

Several studies on real complex networks from different fields as biology, economy, or sociology have shown that the degree of nodes (number of edges connected to each node) follows a scale-free power-law distribution like P(k)≈k-γ, where P(k) denotes the frequency of the nodes that are connected to k other nodes. Here we have carried out a study on scale-free networks, where a line graph transformation (i.e., edges in an initial network are transformed into nodes) is applied on a power-law distribution. Our results indicate that a power-law distribution as P(k)≈k-γ+1 is found for the transformed network together with a peak for low-degree nodes. In the present work, we show a parametrization of this behaviour and discuss its application on real networks as metabolic networks, protein–protein interaction network and World Wide Web.

A Constructive Approach to Gene Expression Dynamics

Recently, experiments on mRNA abundance (gene expression) have revealed that gene expression shows a stationary organization described by a scale-free distribution. Here we propose a constructive approach to gene expression dynamics which restores the scale-free exponent and describes the intermediate state dynamics. This approach requires only one assumption: Markov property.

Clustering under the line graph transformation: application to reaction network

BackgroundMany real networks can be understood as two complementary networks with two kind of nodes. This is the case of metabolic networks where the first network has chemical compounds as nodes and the second one has nodes as reactions. In general, the second network may be related to the first one by a technique called line graph transformation (i.e., edges in an initial network are transformed into nodes). Recently, the main topological properties of the metabolic networks have been properly described by means of a hierarchical model. While the chemical compound network has been classified as hierarchical network, a detailed study of the chemical reaction network had not been carried out.ResultsWe have applied the line graph transformation to a hierarchical network and the degree-dependent clustering coefficient C(k) is calculated for the transformed network. C(k) indicates the probability that two nearest neighbours of a vertex of degree k are connected to each other. While C(k) follows the scaling law C(k) ~ k-1.1 for the initial hierarchical network, C(k) scales weakly as k0.08 for the transformed network. This theoretical prediction was compared with the experimental data of chemical reactions from the KEGG database finding a good agreement.ConclusionsThe weak scaling found for the transformed network indicates that the reaction network can be identified as a degree-independent clustering network. By using this result, the hierarchical classification of the reaction network is discussed.

Flexible construction of hierarchical scale-free networks with general exponent.

Extensive studies have been done to understand the principles behind architectures of real networks. Recently, evidence for hierarchical organization in many real networks has also been reported. Here, we present a hierarchical model that reproduces the main experimental properties observed in real networks: scale-free of degree distribution P (k) [frequency of the nodes that are connected to k other nodes decays as a power law P (k) approximately k(-gamma) ] and power-law scaling of the clustering coefficient C (k) approximately k(-1) . The major points of our model can be summarized as follows. (a) The model generates networks with scale-free distribution for the degree of nodes with general exponent gamma>2 , and arbitrarily close to any specified value, being able to reproduce most of the observed hierarchical scale-free topologies. In contrast, previous models cannot obtain values of gamma>2.58 . (b) Our model has structural flexibility because (i) it can incorporate various types of basic building blocks (e.g., triangles, tetrahedrons, and, in general, fully connected clusters of n nodes) and (ii) it allows a large variety of configurations (i.e., the model can use more than n-1 copies of basic blocks of n nodes). The structural features of our proposed model might lead to a better understanding of architectures of biological and nonbiological networks.

A stochastic approach to multi-gene expression dynamics

Abstract In the last years, tens of thousands gene expression profiles for cells of several organisms have been monitored. Gene expression is a complex transcriptional process where mRNA molecules are translated into proteins, which control most of the cell functions. In this process, the correlation among genes is crucial to determine the specific functions of genes. Here, we propose a novel multi-dimensional stochastic approach to deal with the gene correlation phenomena. Interestingly, our stochastic framework suggests that the study of the gene correlation requires only one theoretical assumption— Markov property —and the experimental transition probability, which characterizes the gene correlation system. Finally, a gene expression experiment is proposed for future applications of the model.

Identification of metabolic units induced by environmental signals

MOTIVATION Biological cells continually need to adapt the activity levels of metabolic functions to changes in their living environment. Although genome-wide transcriptional data have been gathered in a large variety of environmental conditions, the connections between the expression response to external changes and the induction or repression of specific metabolic functions have not been investigated at the genome scale. RESULTS We present here a correlation-based analysis for identifying the expression response of genes involved in metabolism to specific external signals, and apply it to analyze the transcriptional response of Saccharomyces cerevisiae to different stress conditions. We show that this approach leads to new insights about the specificity of the genomic response to given environmental changes, and allows us to identify genes that are particularly sensitive to a unique condition. We then integrate these signal-induced expression data with structural data of the yeast metabolic network and analyze the topological properties of the induced or repressed subnetworks. They reveal significant discrepancies from random networks, and in particular exhibit a high connectivity, allowing them to be mapped back to complete metabolic routes.