Clustering under the line graph transformation: application to reaction network
 Jose C Nacher^{1}Email author,
 Nobuhisa Ueda^{1},
 Takuji Yamada^{1},
 Minoru Kanehisa^{1} and
 Tatsuya Akutsu^{1}
DOI: 10.1186/147121055207
© Nacher et al; licensee BioMed Central Ltd. 2004
Received: 14 October 2004
Accepted: 24 December 2004
Published: 24 December 2004
Abstract
Background
Many 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.
Results
We have applied the line graph transformation to a hierarchical network and the degreedependent 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 k^{0.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.
Conclusions
The weak scaling found for the transformed network indicates that the reaction network can be identified as a degreeindependent clustering network. By using this result, the hierarchical classification of the reaction network is discussed.
Background
Recent studies on network science demonstrate that cellular networks are described by universal features, which are also present in nonbiological complex systems, as for example social networks or WWW. Most networks encountered in real world have scalefree topology, in particular networks of fundamental elements of cells as proteins and chemical substrates [1–4]. In these networks, the distribution of node degree follows a powerlaw as P(k) ~ k^{γ}(i.e., frequency of the nodes that are connected to k other nodes). The degree of a node is the number of other nodes to which it is connected.
One of the most successful models for explaining that scalefree topology was proposed by BarabásiAlbert [5], which introduced a meanfield method to simulate the growth dynamics of individual nodes in a continuum theory framework. However, although that model was a milestone to understand the behavior of real complex networks, it could not reproduce all the observed features in real networks such as clustering dependence. The observed properties of networks with N nodes are: scalefree of degree distribution P(k) ~ k^{γ}, powerlaw scaling of clustering coefficient C(k) ~ k^{1} and a high value for the average of the clustering coefficient <C > and its independence with network size. In particular, the dependence of C(k) ~ k^{1} was one of the results obtained by [6]. In order to bring under a single framework all these observed properties in real networks Ravasz et al. (the RSMOB model in what follows) suggested successfully a hierarchical and modular topology [7, 8]. In [8], a network with the above mentioned properties was called hierarchical network. We note that this deterministic model is an extension of the original model shown in [9]. It is also worth noticing that this modular topology was also suggested in biological networks by [10, 11]. Interestingly, these properties of networks have been found in many non biological and biological networks. One of them, which is the subject of our study, is the metabolic network.
It is interesting to note that the metabolic network is an example of bipartite networks [12]. In a bipartite network there are two kinds of nodes and edges only connect nodes of different kinds. In the metabolic network these nodes are chemical compounds and reactions. The network generated by the chemical compounds (reactions) is called compound (reaction) projection. A line graph transformation (i.e., each edge between two nodes becomes a node of the transformed network) may relate both projections.
However, although the line graph transformation works fine on bipartite networks, the transformed network (in the particular case of metabolic networks) may not be totally the same as the reaction projection. This issue is discussed in detail later. In addition, we will show by comparing with the experimental data, that this fact does not affect our qualitative results. Furthermore, a detailed analysis of the line graph transformation focused on the degree distribution P(k) and applied to some real networks can be found in [13]. In that work, similarities and differences between the line graph transformation and the metabolic network are also discussed. There it was found that if the initial network follows a powerlaw P(k) ~ k^{γ}, the transformed network preserves the scalefree topology and in most cases the exponent is increased by one unit as P(k) ~ k^{γ+1}.
It is also worth noting that the line graph transformation has recently been applied with success by PereiraLeal et al. [15] on the protein interaction network with the aim to detect functional modules. In that work, the edges (interactions) between two proteins become the nodes of the transformed network (interaction network). By means of the line graph transformation, the interaction network has a higher clustering coefficient than the protein network. By using the TribeMCL algorithm [16] they are able to detect clusters in the more highly clustered interaction network. These clusters are transformed back to the initial proteinprotein network to identify which proteins can form functional clusters. At this point, we note that the aim of our study is not to detect functional modules from the metabolic network. In our work the line graph transformation is used successfully to evoke general topological properties related to the clustering degree of the reaction network.
The observed topological properties related to the clustering degree of the metabolic network (in particular, the chemical compound network) have been properly described by means of the RSMOB model. In the present work, our aim is to study the clustering coefficients C(k) and <C > of the reaction network by using two approaches: Firstly, we derive mathematical equations of those coefficients in the transformed network. Secondly, we apply the line graph transformation to a hierarchical network. The results from both methods are compared with experimental data of reactions from KEGG database [14] showing a good agreement. Though we started this work motivated by theoretical interest in the line graph transformation, the results provide explanation for the difference of C(k) between the compound network and the reaction network.
In our work, the hierarchical network is generated by the RSMOB model, where the nodes correspond to chemical compounds and the edges correspond to reactions. While the RSMOB model reproduces successfully the hierarchical properties of the compound network, here we show that this hierarchical model also stores adequate information to reproduce the experimental data of the reaction network. Our study indicates that it is enough to apply the line graph transformation to the hierarchical network to extract that information. While C(k) follows the powerlaw k^{1.1} for the initial hierarchical network (compound network), C(k) scales weakly as k^{0.08} for the transformed network (reaction network). Consequently, we conclude that the reaction network may not be classified as a hierarchical network, as it is defined in [8].
Remark
In [8], a network with scalefree topology, scaling law of C(k) ~ k^{1}, and high degree of clustering was called hierarchical network. Consequently, the RSMOB model shown in [7, 8] was developed to bring these properties under a single roof. Furthermore, in [7, 8, 17] some networks (in particular metabolic network) were classified as hierarchical network according to the above definition. To be precise, it was argued that the signature of the intrinsic hierarchy (or hierarhical modularity) is the scaling law of C(k). Moreover, in a more recent work [18], it was claimed that traditional random and scalefree models do not have a hierarchical topology because C(k) is independent of k (i.e., flat plot of C(k)). In addition, analyses of C(k) were recently carried out in [19] to uncover the structural organization and hierarchy of nonbiological weighted networks. At this point, we must note that we have followed the research done by Barabási et al, [7, 8, 18] and consequently, we have used its definition of hierarchical network in the present work.
However, it is also worth noticing that another way to quantify the hierarchical topology of a network is recently introduced by [20, 21]. It is based on the concept of a hierarchical path: a path between nodes i and j is called hierarchical if (1) the node degrees grow monotonously ("up path"), and it is followed by a path where the node degrees decrease monotonously ("down path") or (2) the node degrees along this path changes monotonously from one node to the other. The fraction of shortest paths in a network, which are also hierarchical paths is called H. If H is very close to 1, the network shows a hierarchical organization. This definition seems interesting, and consequently, as a future work it would be worth to examine some biological networks (in particular, metabolic networks) by using this approach. One remark about this concept is that it focuses on hierarchy and may not contain enough information about modularity or clustering. For a brief discussion of these issues, we refer to some very useful notes written by Dorogovtsev et al. [21] (see also [22] for further information about related topics in networks).
Results and discussion
Clustering coefficients C(k) and <C>
Recent analyses have demonstrated that the metabolic network has a hierarchical organization, with properties as: scalefree degree distribution P(k) ~ k^{γ}, powerlaw dependence of clustering coefficient C(k) ~ k^{1} and independence with network size of the average clustering coefficient <C >, where N is the total number of nodes in a network [7]. The clustering coefficient can be defined for each node i as:
Geometrically, n_{ i }gives the number of triangles that go through node i. The factor k_{ i }(k_{ i } 1)/2 gives the total number of triangles that could go through node i (i.e., total number of triangles obtained when all the neighbors of node i are connected to each other). In the case of Fig. 1A, there is one triangle that contains node 1 (dashdotted lines), and a total of 6 triangles could be generated as the maximum. Hence, the clustering coefficient of node 1 is C_{1} = 1/6.
On the other hand, the average clustering coefficient <C > characterizes the overall tendency of nodes to form clusters as a function of the total size of the network N. The mathematical expression is:
The structure of the network is given by the function C(k), which is defined as the average clustering coefficient over nodes with the same node degree k. This function is written as:
where N_{ k }is the number of nodes with degree k, and the sum runs over the N_{ k }nodes with degree k. A scaling law k^{1} for this magnitude is an indication of the hierarchical topology of a network.
Once the theoretical definitions have been introduced, our aim is to analyse how the coefficients <C > and C(k) are modified under the line graph transformation.
Line graph transformation to metabolic networks: spurious nodes
Given an undirected graph G, defined by a set of nodes V(G) and a set of edges E(G), we associate another graph L(G), called the line graph of G, in which V(L(G)) = E(G), and two nodes are adjacent if and only if they have a common endpoint in G (i.e., E(L(G)) = {{(u, v), (v, w)}(u, v) ∈ E(G), (v, w) ∈ E(G)}). This construction of graph L(G) from the initial graph G is called line graph transformation [23].
It is worth noting that in a previous work [13] the degree distribution P(k) was studied by applying line graph transformation to synthetic and real networks. There it is assumed an initial graph G with scalefree topology as P(k) ≈ k^{γ}. As the degree of each transformed node (i.e., an edge in G) will be roughly around k, the distribution of the line graph L(G) should be k·k^{γ}= k^{γ+1}with degree around k. Therefore, it is concluded that if we have a graph G with a probability distribution following a powerlaw as k^{γ}, then L(G) will follow a powerlaw as k^{γ+1}. The real networks under study were proteinprotein interaction, WWW, and metabolic networks. In Fig. 1B, we can see an example of the line graph transformation applied to a subgraph of the metabolic network.
Equations of C(k) and <C> under the line graph transformation
The clustering coefficient for the node a in the transformed network can be written by using Eq. (1) as:
where k = k' + k"  2, because the edge a vanishes in the graph L(G). This equation ignores cases where edges in the graph G, b and b' for example, have a common node as endpoint (i.e., existence of triangles or loops in Fig. 3C). However, we can quantify these cases by using a new parameter l. As we can see in Fig. 3C–D, edges with one common node as endpoint in the graph G means one additional edge in the graph L(G). This additional edge in L(G) connects two neighbors of node a. By following definition of Eq. (1), it means that n_{ a }increases its value by one unit. We can consider these cases by increasing one unit the parameter l for each common node as endpoint of two edges in the graph G (for example, l = 1 means one common node). We write Eq. (4) after introducing the parameter l as:
where if l = 0 means that there are not loops and we recover Eq. (4). Though it is more realistic to consider the parameter l as a function of k' and k", we have considered l as an independent parameter. However, this simplification does not affect the qualitative features of our results. It should be noted that l always contributes to increasing the value of C_{ a }(k) and C_{ a }(k) ≤ 1 always holds from the definition. In order to study the limits of Eq. (5) we consider the following two cases:

a) k' = k": We analyse the case where both degrees have the same value. We also consider the cases when l = 0 and l ≠ 0 in order to study the effect of triangles. We show the results in Fig. 4. For large k', Eq. (5) goes asymptotically to 1/2 for l = 0 and l ≠ 0. We also see that for k' ≥ 25, all lines are very close to 1/2. For low k' and l = 0, C_{ a }(k) takes values from 0.33 (k' = 3) to 0.48 (k' = 20). Hence, we see in Fig. 4 that higher values of l (more triangles) increase the values of C_{ a }(k).

b) k" = constant, k' >>k": We plot in Fig. 5 three cases. k" is fixed with constant values as k" = 5 (black), k" = 10 (red), k" = 20 (blue) and k' is a free parameter. We see that C_{ a }(k) approaches to 1 when k' takes large values. For low k', the case k" = 5 shows a minimum with a few values of k' below 1/2. As we can see with dotted and dashdotted lines in Fig. 5, the presence of triangles (l ≠ 0) increases the value of C_{ a }(k). Finally, for k" = 10 and k" = 20, we see that only a few values of C_{ a }(k) are slightly below 1/2 for low k'. This analysis is complemented by calculating the minimum value of C_{ a }(k) analytically as: . The value of k', where the function C_{ a }(k) takes the minimum value, is given by:
where positive solution of the square root is written. By substituting this equation into Eq. (5), it is possible to calculate the minimum value of C_{ a }(k) for each configuration of l and k".
From these two cases, we can conclude that for hubs (i.e., those nodes with high degree (k' and k" >> 1)) and for highly clustered networks (many triangles l >> 1), the values of C_{ a }(k) in the transformed network are between around [ , 1].
To calculate the distribution of C(k) in the transformed space (C^{ T }(k)) we introduce the concept of assortativity. By assortative (disassortative) mixing in networks we understand the preference for nodes with high degree to connect to other high (low) degree nodes [24]. By following Newman [24], we define the probability distribution to choose a randomly edge with two nodes at either end with degrees k' and k" as e_{k'k"}. We also assume that the nodes of the initial network are following a powerlaw distribution k^{γ}and have no assortative mixing. Under these assumptions, the probability distribution e_{k'k"}of edges that link together nodes with degree k' + k" can be written as:
We make a convolution between Eq. (4) and Eq. (7), by summing for all the possible degrees of the two nodes at either end of edges (k', k"), which can generate transformed nodes with degree k = k' + k"  2. Thus, we obtain:
According to the structure of C^{ T }(k) and the behavior of C_{ a }(k) exposed above, C^{ T }(k) will grow smoothly for large k, i.e., scaling weakly with the node degree k. We have calculated numerically this expression and the results are discussed later in Fig. 7.
We have also calculated the analytical expression for <C >, and we have found that <C > has a sizeindependent behavior before and after the line graph transformation is done. We can write the number of nodes with degree k as:
and we assume that C(k) = A·k ^{α}, where A is a constant. This constant changes when we consider hierarchical networks with different number of nodes in the initial cluster [7]. But it seems natural because in that case the degree distribution P(k) ~ k^{γ}of the network also changes. For <C > before the transformation we can write:
Note that the summation in the denominator begins with k = 1 because we renormalize over all the probability distribution.
Furthermore, we can obtain <C > by using the RSMOB model (explained in next section in detail). This model starts by generating a fully connected cluster of m nodes, such that the connectivity of each node is k = m  1. In the following iteration, m  1 replicas of the initial cluster are generated, and linked to the central node of the original cluster in such a way that the central node of the original cluster gains (m  1)·(m  1) edges, and its total connectivity being k = m  1 + (m  1)^{2}. By iterating these procedure, it is easy to see that hubs (i.e., central nodes of each replica) will have connectivities , with j = 1, ..., log_{ m }N being the iteration number. Therefore, assuming that the degree distribution P(k) and the clustering coefficient C(k) are powerlaws with exponents γ' and α respectively, the expression for <C > for the hubs reads as:
where A' is a constant adjusted so that <C >=1 holds for j = 1. The upper limit of the summation log_{ m }N is obtained by means of the expression m^{ j }= N, which gives the total number of nodes in the network and denotes the exponent of the powerlaw distribution of hubs in the RSMOB model. We must note that in a hierarchical network, the number of nodes with different degree k is scarce, therefore the probability distribution of node degree is properly defined as P(k) = (1/N_{ tot })(N_{ k }/Δk), where N_{ k }is the number of nodes with degree k, N_{ tot }is the total number of nodes, and Δk means that nodes with degree k are binned into intervals. In addition, we note that for the hierarchical model, Δk changes linearly with k. Hence, the exponent of the powerlaw is given by γ = 1 + γ', with where m is the number of nodes in the initial module.
Results of <C > evaluated by using Eq. (10) and the needed parameters in that calculation for 3 different setups: γ = 1 + γ', where (P(k) ~ k^{γ}), α (C(k) ~ k^{α}), A (C(k) = A·k^{α}). Eq. (10) is a general expression of <C >.
m initial nodes  γ  α  A  <C > 

3  2.58  1.1  2.34  0.20 
4  2.26  1.1  3.68  0.36 
5  2.16  1.1  5.18  0.54 
Results of <C > evaluated by using Eq. (11) for 3 different setups. The exponent of the powerlaw distribution of hubs is given by . The parameter α has same meaning as in Table 1. We also notice that in Eq. (11), A' is adjusted so that <C >=1 holds for j = 1. Eq. (11) is the particular expression of <C > applied to the RSMOB model.
m initial nodes  γ'  α  <C > (Eq. (11) 

3  1.58  1.1  0.78 
4  1.26  1.1  0.81 
5  1.16  1.1  0.83 
Definitions of functions and their values before and after the line graph transformation is applied to the hierarchical network. N_{ k }: number of nodes of degree k. The † symbol means that these dependences were analyzed in the present work, while the * symbol means that it was studied in our previous work [13].
Func.  Definition  Dependence before  Dependence after (Eq. (11) 

P(k)  N_{ k }/N  k ^{γ}  k^{γ+1}* 
C_{ i }(k)  2n/[k_{ i }(k_{ i } 1)]  k ^{1.1}  k ^{0.08†} 
<C > 
 sizeindependent^{†}  sizeindependent^{†} 
Line graph transformation to a hierarchical network: numerical results
The RSMOB model [8] is able to reproduce the main topological features of the metabolic network. We follow the method described in [8] and generate a hierarchical network. Then, we apply the line graph transformation to that network.
Fig. 6 illustrates the hierarchical network generated by the RSMOB model. The network is made of densely linked 5node modules (it is worth noticing that the number of nodes in the initial module can be different than 5) that are assembled into larger 25node modules (iteration n = 1, 5^{2} = 25 nodes). In the next step four replicas are created and the peripheral nodes are connected again to produce 125node modules (iteration n = 2, 5^{3} = 125 nodes). This process can be repeated indefinitely [8].
To evaluate C(k), we have constructed three hierarchical networks with 3, 4, and 5 initial number of nodes. These networks were generated up to 7 (6561), 5 (4096), and 4 (3125) iterations (nodes), repectively. Once we have constructed these three networks, we apply the line graph transformation to them, and we calculate the C^{ T }(k) clustering coefficient for the transformed networks. In Fig. 7A we show the results of the clustering coefficient of the transformed network. Circles, triangles and squares indicate the values of C^{ T }(k) for the transformed network with 3, 4, and 5 initial nodes, respectively. In Fig. 7A we also plot with continuous lines the values of C^{ T }(k) obtained from Eq. (8). From top to bottom the lines correspond to the networks of 3, 4 and 5 initial nodes, respectively. In Fig. 7A, we see that the lines show an acceptable agreement with the overall tendency of data generated by the transformed network. In Fig. 7B, we see that the results from theoretical calculation of C^{ T }(k) via Eq. (8) (lines) are in good agreement with the experimental data (diamonds) from the KEGG database [14]. Moreover, in order to have enough statistics to compare with the analytical expression for the C^{ T }(k), we have binned into seven intervals the experimental data according to degree k (1 <k ≤ 8 < ,..., 128 <k ≤ 256, 256 <k ≤ 512), and averaged over the C^{ T }(k)'s obtained in that range (red circles). It shows a better agreement between KEGG results and the analytical curves. The only disagreement comes at k = 2. This is easy to understand because in the hierarchical model depicted in Fig. 6, we can only find C(k = 2) = 1 for 3 initial nodes by construction of the network. However, in real networks, we could find nodes which have only two neighbors and, in some cases, these neighbors could be connected. In these cases the clustering coefficient takes value one.
We show in Fig. 9A the values of <C > calculated for networks generated by 3 initial nodes (circles) and 4 initial nodes (triangles) by using the RSMOB model. We see that <C > approaches asymptotically to constant values around ~ 0.66 (circles) and ~ 0.74 (triangles), being independent of the size of the network. Once the line graph transformation is applied, we see that the corresponding values of <C^{ T }> also approach asymptotically to constant values. Hence, <C^{ T }> also is sizeindependent for large N (empty circles and triangles). In addition, we have averaged the experimental value of the clustering coefficient for reactions of 163 organisms found in KEGG database and we have obtained the value of <C^{ T }>= 0.74. We see that the experimental value <C^{ T }> for reactions is in good agreement with the asymptotic values obtained by the transformed network (empty triangles and circles).
Furthermore, we have also calculated <C > by using Eq. (11). This equation should reproduce the results of <C > calculated by using the RSMOB model (dark circles and triangles in Fig. 9A). In Fig. 9B, we see that the results are qualitatively similar to those shown in Fig. 9A (dark circles and triangles).
We remark that the theoretical analysis of <C > and <C^{ T }> done here has also been useful to prove that they are independent of network size.
Conclusions
We have studied here the clustering coefficients C(k) and <C > of the reaction network by applying the line graph transformation to a hierarchical network. This hierarchical network was generated by using the RSMOB model, which reproduces properly the topological features of the metabolic network, in particular the compound network. Our results indicate that by applying the line graph transformation to the hierarchical network, it is possible to extract topological properties of the reaction network, which is embedded in the metabolic network. The RSMOB model stores the adequate information of the reaction network and the line graph transformation is one useful technique to evoke it.
While C(k) scales as k^{1.1} for the initial hierarchical network (compound network), we find C(k) ~ k^{0.08} for the transformed network (reaction network). This theoretical prediction was compared with the experimental data from the KEGG database, finding a good agreement. Our results indicate that the reaction network is a degreeindependent clustering network. Furthermore, the weak scaling of C(k) for the reaction network suggests us that this network may not have hierarchical organization. However, further analyses of this network, and in general of all biological networks, by following the concept of hierarchical path are encouraged [20, 21].
On the other hand, we have also conducted an analytical derivation for the clustering coefficients C(k) and <C >. Expressions for these coefficients were calculated before and after the line graph transformation is applied to the hierarchical network. The agreement obtained by using these expressions was found acceptable, and consequently, they could be useful for further analyses in different networks (biological and nonbiological).
The line graph transformation has recently been applied on metabolic networks [13] to study the scalefree topology of the reaction network, and on the proteinprotein interaction network to detect functional clusters [15]. The work done here is another important application of this interesting technique.
Declarations
Acknowledgements
This work was partially supported by GrantinAid for Scientific Research on Priority Areas (C) "Genome Information Science" from MEXT of Japan.
Authors’ Affiliations
References
 Wuchty S: Scalefree behavior in protein domain networks. Mol Biol Evol 2001, 18: 1694–1702.View ArticlePubMedGoogle Scholar
 Titz B, Schlesner M, Uetz P: What do we learn from highthroughput protein interaction data and networks? Expert Reviews in Proteomics 2004, 1(1):89–99.View ArticleGoogle Scholar
 Wagner A, Fell DA: The small world inside large metabolic networks. Proc Roy Soc London Series B 2001, 268: 1803–1810. 10.1098/rspb.2001.1711View ArticleGoogle Scholar
 Jeong H, Tombor B, Albert R, Oltvai ZN, Barabási AL: The large organization of metabolic networks. Nature 2000, 407: 651–654. 10.1038/35036627View ArticlePubMedGoogle Scholar
 Barabási AL, Albert R: Emergence of scaling in random networks. Science 1999, 286: 509–512. 10.1126/science.286.5439.509View ArticlePubMedGoogle Scholar
 Dorogovtsev SN, Goltsev AV, Mendes JFF: Pseudofractal scalefree web. Physical Review E 2002., 65: 066122–14Google Scholar
 Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabási AL: Hierarchical organization of modularity in metabolic networks. Science 2002, 297: 1551–1555. 10.1126/science.1073374View ArticlePubMedGoogle Scholar
 Ravasz E, Barabási AL: Hierarchical organization in complex networks. Phys Rev E 2003., 67: 026112–17Google Scholar
 Barabási AL, Ravasz E, Vicsek T: Deterministic scalefree networks. Physica A 2001, 299: 559–564.View ArticleGoogle Scholar
 Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature 1999, 402: C47–52. 10.1038/35011540View ArticlePubMedGoogle Scholar
 Milo R, ShenOrr S, Itzkovitz S, Kashtan N, Alon U: Network motifs: simple building blocks of complex networks. Science 2002, 298: 824–827. 10.1126/science.298.5594.824View ArticlePubMedGoogle Scholar
 Holme P, Liljeros F, Edling CR, Kim BJ: Network bipartivity. Phys Rev E 2003., 68: 056107–112Google Scholar
 Nacher JC, Yamada T, Goto S, Kanehisa M, Akutsu T: Two complementary representations of scalefree networks. Physica A 2004, in press. eprint archive, physics/0402072Google Scholar
 Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res 2004, 32: D277D280. The KEGG Database [http://www.genome.ad.jp/KEGG/] 10.1093/nar/gkh063PubMed CentralView ArticlePubMedGoogle Scholar
 PereiraLeal JB, Enright AJ, Ouzounis CA: Detection of functional modules from protein interaction networks. Proteins 2004, 54: 49–57. 10.1002/prot.10505View ArticlePubMedGoogle Scholar
 Enright AJ, Van Dongen S, Ouzounis CA: An efficient algorithm for largescale detection of protein families. Nucleic Acids Res 2002, 30: 1575–1584. 10.1093/nar/30.7.1575PubMed CentralView ArticlePubMedGoogle Scholar
 Barabási AL, Deszo Z, Ravasz E, Yook SH, Oltvai Z: to appear in Sitges Proceedings on Complex Networks. 2004.Google Scholar
 Barabási AL, Oltvai ZN: Network biology: Understanding the cell's functional organization. Nature Genetics Reviews 2004, 5: 101–114. 10.1038/nrg1272View ArticleGoogle Scholar
 Barrat A, Barthelemy M, PastorSatorras R, Vespignani A: The architecture of complex weighted networks. Proceedings of the Nat'l Academy of Sciences 2004, 101: 3747–3752. 10.1073/pnas.0400087101View ArticleGoogle Scholar
 Trusina A, Maslov S, Minnhagen P, Sneppen K: Hierarchy measures in complex networks. 2004.Google Scholar
 Dorogovtsev SN, Mendes JFF: The shortest path to complex networks. 2004.Google Scholar
 Mendes JFF, Dorogovtsev SN, Ioffe AF: Evolution of Networks: From Biological Nets to the Internet and WWW. Oxford Press; 2003.Google Scholar
 Hemminger RL, Beineke LW: Selected Topics in Graph Theory I. Academic Press, London; 1978.Google Scholar
 Newman MEJ: Assortative mixing in networks. Phys Rev Lett 2002., 89: 208701–15Google Scholar
 The algorithm package is available at2003. [http://vlado.fmf.unilj.si/pub/networks/pajek]
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.