Novel topological descriptors for analyzing biological networks
© Dehmer et al; licensee BioMed Central Ltd. 2010
Received: 29 October 2009
Accepted: 17 June 2010
Published: 17 June 2010
Skip to main content
© Dehmer et al; licensee BioMed Central Ltd. 2010
Received: 29 October 2009
Accepted: 17 June 2010
Published: 17 June 2010
Topological descriptors, other graph measures, and in a broader sense, graph-theoretical methods, have been proven as powerful tools to perform biological network analysis. However, the majority of the developed descriptors and graph-theoretical methods does not have the ability to take vertex- and edge-labels into account, e.g., atom- and bond-types when considering molecular graphs. Indeed, this feature is important to characterize biological networks more meaningfully instead of only considering pure topological information.
In this paper, we put the emphasis on analyzing a special type of biological networks, namely bio-chemical structures. First, we derive entropic measures to calculate the information content of vertex- and edge-labeled graphs and investigate some useful properties thereof. Second, we apply the mentioned measures combined with other well-known descriptors to supervised machine learning methods for predicting Ames mutagenicity. Moreover, we investigate the influence of our topological descriptors - measures for only unlabeled vs. measures for labeled graphs - on the prediction performance of the underlying graph classification problem.
Our study demonstrates that the application of entropic measures to molecules representing graphs is useful to characterize such structures meaningfully. For instance, we have found that if one extends the measures for determining the structural information content of unlabeled graphs to labeled graphs, the uniqueness of the resulting indices is higher. Because measures to structurally characterize labeled graphs are clearly underrepresented so far, the further development of such methods might be valuable and fruitful for solving problems within biological network analysis.
Major reasons for the emergence of biological network analysis [1–4] are the extensive use of computer systems during the last decade and the availability of highly demanding and complex biological data sets. For instance, important types of such biological networks are protein-protein interaction networks [5–7], transcriptional regulatory networks [8, 9], and metabolic networks [7, 10, 11]. Note that vertices in such biological networks can represent, e.g., proteins, transcription factors or metabolites which are connected by edges representing interactions, concentrations or reactions, respectively [3, 12]. Thus, vertex-and edge-labeled graphs is an important graph class [13, 14] and useful for modeling biological networks . To name only some well-known examples or methods which have often been applied within biological network analysis, we briefly mention graph classes like scale-free and small-world networks [15, 16], network centralities [12, 17], module and motif detection [18–20], and complexity measures for exploring biological networks structurally [21, 22].
Taking into account that a large number of graph-theoretical methods have been developed so far, approaches to process and meaningfully analyze labeled graphs are clearly underrepresented in the scientific literature. In particular, this holds for chemical graph analysis where various graph-theoretical methods and topological indices have been intensely used, see, e.g., [23–34]. Yet, we state a few examples where such graphs appear in the context of biological network analysis: Structure descriptors to determine the complexity of pathways representing labeled graphs have been used to examine the relationship between metabolic and phylogenetic information, see . Another challenging task relates to determine the similarity between graphs or subgraphs [35–38]. For instance, YANG et al.  recently developed path-and graph matching methods involving vertex-and edge-labeled graphs which turned out to be useful for biological network comparison . Finally, to utilize graph-theoretical concepts for investigating graphs and labeled graphs within molecular biology, HUBER et al.  reviewed several existing software packages and outlined concrete applications .
In this paper, we restrict our analysis to a set of bio-chemical graphs which have already been used for predicting Ames mutagenicity, see . To perform this study, we develop and investigate entropic descriptors for vertex- and edge-labeled graphs. Before sketching the main contributions of our paper, we state some facts about topological descriptors which have been used in mathematical chemistry, drug design, and QSPR/QSAR.
As already mentioned, topological indices have been proven to be powerful tools in drug design, chemometrics, bioinformatics, and mathematical and medicinal chemistry [23, 24, 26, 28, 29, 34, 41–43]. Certainly, one reason for their success can be understood by the fact that there is a strong need to apply empirical models to solve QSPR (Quantitative structure-property relationship)/QSAR (Quantitative structure-activity relationship) problems [24, 28, 29, 44] and related tasks in the just mentioned areas. In this paper, we put the emphasis on developing novel molecular descriptors for tackling a problem in QSAR: We will use structural property descriptors of molecules based on SHANNON' s entropy for predicting Ames mutagenicity, see [40, 45–47]. Generally, we note that the problem of detecting mutagenicity in vitro is based on the bacterial reverse mutation assay (Ames test) and often serves as a crucial tool in drug design and discovery [40, 45–47].
Further, topological descriptors have often been combined with other techniques from statistical data analysis, e.g., clustering methods [26, 48] to infer correlations between the used indices. Besides using topological descriptors for characterizing chemical graphs [27, 32, 49], they have also been applied to quantify the structural similarity of chemicals representing networks [50, 51]. Among the large number of existing topological indices, an important class of such measures relies on SHANNON's entropy to characterize graphs by determining their structural information content [27, 52–54]. Until now, especially these measures have been intensely applied within biology, ecology, and mathematical chemistry [27, 52, 54–60], in particular, to measure the complexity of biological and chemical systems [27, 52, 61]. Recently, we already developed a novel procedure to infer such information-theoretic measures for graphs that results in so-called partition-independent measures [57, 62]. More precisely, we mean that we do not induce partitions using the procedure manifested by Equation (2), (3) in . In this work, partitions using graph invariants and equivalence criteria have been explicitly induced, see, e.g., [27, 52, 53]. Note that we already placed a comment on this problem in the first paragraph of the section 'Partition-Independent Information Measures for Graphs'. In contrast to partition-independent measures, classical partition-based information measures often rely on the problem to group elements manifested by an arbitrary graph invariant according to an equivalence criterion [27, 53, 54, 63].
The contribution of our paper is twofold: First, we develop some novel information-theoretic descriptors having the ability to incorporate vertex- and edge-labels when measuring the information content of a chemical structure. Because we already mentioned that there is a lack of graph measures which can process vertex-and edge-labeled graphs meaningfully, such descriptors need to be further developed. In terms of analyzing chemical structures, that means they can only be adequately represented by graphs if different types of atoms (vertices) and different types of bonds (edges) are considered. Hence, there is a strong need to exploring such labeled networks. Besides developing the novel information-theoretic measures for vertex- and edge-labeled graphs, we will investigate some of their properties thereof (see section 'Properties of the Novel Information-Theoretic Descriptors') [40, 47]. Second, the paper also deals with evaluating the ability of the mentioned descriptors to predict Ames mutagenicity when applying well-known machine learning methods like random forests [64, 65] (RF) and support vector machines [64, 66] (SVM). Starting from chemical structures represented as vectors composed of topological descriptors, we will analyze the prediction performance by focussing on the underlying supervised graph classification problem. We want to emphasize that beside our novel descriptors, we also combine them with other well known information-theoretic and non-information-theoretic measures which turned out to be useful in QSPR/QSAR, see, e.g., . Further, we examine the influence on the prediction performance when taking semantical (labels) and structural information of the graphs into account. Finally, we want to point out that considerable related work has been done so far that deals with investigating multifaceted problems when applying molecular descriptors to machine learning algorithms [67–69]. For example, DESHPANDE et al.  developed an approach to find discriminating substructures of chemical graphs. Then, by using a vector representation model for these graphs, they applied several machine learning methods to chemical databases for classifying these structures meaningfully. Another interesting study was done by XUE et al.  that deals with applying a variety of molecular descriptors to characterize structural and physicochemical properties of molecules . Particularly, they used a feature selection method for automatically selecting molecular descriptors for SVM-prediction of P-glycoprotein substrates and others. As an important result, XUE et al.  determined the reduction of noise and its influence on the prediction accuracy of a statistical learning system . The last contribution we want to sketch in brief is due to MAHÉ et al. . In this work, a graph kernel approach [64, 69] was validated for structure-activity-relationship analysis where special kernels based on random walks were used and optimized. Note that more related work can be found in [40, 71–74].
Definition 1 is a finite, undirected graph. In this paper, we always assume that the considered graphs are connected and do not have loops.
Definition 2 Let G be a finite and undirected graph. δ(v) is called the degree of a vertex v ∈ V and equals the number of edges e ∈ E which are incident with v.
Definition 30 d(u, v) stands for the distance between u ∈ V and v ∈ V expressed as the minimum length of a path between u,v. Further, the quantity σ(v) = maxu∈Vd(u, v) is called the eccentricity of v ∈ V. ρ(G) = maxv∈Vσ(v) is called the diameter of G.
the j-sphere of a vertex v i regarding G.
be unique (finite) vertex and edge alphabets, respectively. and are the corresponding edge and vertex labeling functions. G := (V,E,l V ,l E ) is called a finite, labeled graph.
Clearly, denotes the cardinality of the set of vertices whose distances, starting from v i , are equal to j and possess the vertex label .
To finalize this section, we repeat the definition  of a so-called local information graph of an undirected graph G. In the following, we will use this definition to derive an advanced information functional for incorporating edge- and vertex-labels when measuring the structural information content of a labeled network.
As already outlined, the majority of classical information measures for graphs are based on determining partitions by using an arbitrary graph invariant and an equivalence criterion, see, e.g., [27, 48, 53, 54]. However, DEHMER et al. [57, 62] recently proposed another method for quantifying the structural information content of a graph. The key principle of this approach is to assign a probability value to every vertex in a graph using different information functionals [57, 62]. This results in partition-independent information measures to determine the entropy of the underlying graph topology. We already explained why we call our measures partition-independent (see also the section 'Background'). In a narrow sense, one might argue that to calculate the information functionals and f E (see next section), we also deal with certain graph partitions for quantifying the information content of a vertex- and edge-labeled graph because we have to compute all local information graphs (local subgraphs). But nonetheless, the construction of our information measures basically differs from the ones mentioned in  (see Equation (2), (3)). In fact, we end up with probability values for every vertex of a given graph. Now, in order to start developing the new measures, we briefly recall the most important definitions. A recent review on information-theoretic descriptors to quantify structural information of unlabeled graphs can be found in .
f represents an arbitrary information functional.
Now, we repeat the definition of an information functional for quantifying the structural complexity of unlabeled and unweighted chemical graphs . Generally, this relates to measure the structural information content of a graph that is interpreted as the entropy of the underlying graph topology.
Remark 1 We want to point out that further information functionals have been developed so far. The appropriateness of such a functional that captures structural information of a graph strongly depends on the graph class and on the specific problem under consideration.
Another measure to determine the structural information content is the following one. Until now, it has been used  to perform a statistical analysis when determining structural complexity of real chemical structures and investigating correlations with other molecular descriptors . Mathematical properties thereof were also described in .
λ > 0 is a scaling constant.
In this section, we present novel information measures to quantify structural information of labeled (weighted) chemical structures by adapting the just shown approach. Because the majority of the developed topological indices is only defined for the underlying skeleton of a chemical structure, the further development of descriptors for processing chemical graphs containing heteroatoms and multiple bonds is generally of great importance. Before we start expressing the new definitions, we first point out some related work in this area.
Note that earlier contributions to infer measures for labeled graphs are often based on special distance matrices and polynomial methods [78–80]. Another attempt in this direction was done by IVANCIUC et al.  where this approach is based on defining weighted matrices incorporating special weighting schemes . For example, a definition of a connectivity, adjacency, distance, and reciprocal distance matrix by applying several weighting schemes incorporating chemical information like the atomic bond number, electronegativity, and the covalent radius have been investigated . Then, such matrices have been used to define molecular descriptors for quantifying information of weighted chemical graphs, e.g., organic compounds. To name some concrete examples, we first mention the WIENER index  for vertex-and edge-labeled graphs when applying the known formula for calculating this index with a special weighting scheme as mentioned above . Further, starting from the mentioned weighted matrices, the well-known information indices U, V, X, Y have been extended to determine the structural information content of labeled (weighted) graphs . As a result, IVANCIUC et al. [81, 84] obtained information-theoretic topological descriptors for vertex- and edge-labeled graphs where the underlying (weighted) matrix may contain negative elements and those between zero and one.
We now start by stating the novel partition-independent information-based descriptors to determine the information content of vertex- and edge-labeled graphs. The first definition represents an information functional to account for vertex labels of a chemical structure. For this, we adapt the idea [57, 62] of determining the topological neighborhoods (using j-spheres) for all involved atoms (vertices) of the molecule. By now considering labeled graphs, our first attempt results in an information functional with the property that every vertex in each j-sphere possessing a certain vertex label (atom type) will be weighted differently.
Because it is not always clear how to choose the involved parameter in practice, we further derive an information functional to overcome this problem.
This quantity denotes the number of vertices of ℒ G (v i , j) possessing vertex label .
By determining all local information graphs for the remaining vertices of G, the just shown calculation can be performed analogously.
Next, we are able to derive an information functional that takes the edge labels of a graph G into account. The main idea is to use weighted paths which can be directly determined by calculating the local information graphs.
Now, we present an example how to apply this definition to the local information graphs shown in Fig. 2.
In order to incorporate both edge and vertex labels when determining the topological entropy of a labeled graph, we also derive
Finally, we obtain the following entropy measures for measuring the structural information content of labeled graphs.
Remark 6 We emphasize that according to the above stated definition and the definitions of the underlying information functionals, the resulting information measures are obviously parametric. This property generalizes classical information measures which have often been used in mathematical chemistry, see, e.g., [27, 29, 53, 83]. As already pointed out in, such measures establish a link to machine learning because the parameters could be learned using appropriate datasets. However, we won't study this problem in the present paper.
This section aims to evaluate the just presented (see previous section) information measures for labeled graphs numerically. Also, we will calculate some known information indices to tackle the second part of our study when applying these measures to machine learning algorithms. Our study will be twofold: First, we examine some properties of the measures for labeled graphs when applying them to a large set of real chemical structures. Second, we analyze a QSAR problem by applying supervised machine learning methods [64, 85] using our novel molecular descriptors.
We created the database AG 3982 from the benchmark database called Ames mutagenicity [40, 47] originally used for the evaluation and prediction of the mutagenicity of chemical compounds . The Ames database was created from six different public sources [40, 47] and each chemical structure possesses a class label (0 and 1) that results from the Ames test indicating the genetoxicity of a substance. By starting from the original database Ames mutagenicity [40, 47] containing 6512 chemical compounds, we created AG 3982 by filtering out isomorphic graphs based on the software SubMat . Finally, this procedure resulted in 3982 structurally different skeletons, that is, all atoms and all bonds are considered as equal. Among these 3982 graphs, 1794 possess class label 0 and 2188 possess 1. It holds 2 ≤ |V| ≤ 109; 1 ≤ ρ(G) ≤ 47 ∀ G ∈ AG 3982. To evaluate the novel descriptors for labeled graphs, we then considered these structures as vertex- and edge-labeled graphs. Evidently, for calculating the descriptors of the unlabeled graph versions (skeletons), the corresponding descriptors were used which take only topological information into account.
To generate and process the underlying graph structures, we used the known Molfile format . The graphs from AG 3982 were originally available in Smiles format that we converted to Molfile format (SDF) using a Python procedure. The implementation of all topological descriptors based on the Molfile format was performed by Python using freely available libraries like Networkx, Openbabel and Pybel packages . To perform the graph classification using random forests (RF) [64, 66] and support vector machines (SVM) [64, 66], we used the implementations provided by the Python package Orange . The feature selection was done by Weka .
Before starting to evaluate our novel molecular descriptors, we define some concrete information measures by choosing special weighting schemes for the coefficients.
the measures and follow correspondingly. Further, if the underlying c i will be chosen linearly or exponentially decreasing (in both functionals and f E ); see also that the measures and follow correspondingly (Equation (50), (35), (52), (53)).
Definition 19 Let G = (V, E, l V , l E ) be an undirected finite labeled graph, . If we choose the coefficients of information functional (see Equation (29)) linearly or exponentially decreasing, we call the resulting information measures and .
Sensitivity for AG 3982
Topological Index I
Here, we assume that the distance of a value i in the distance matrix appears 2k i times . μ denotes the cyclomatic number . To evaluate the discrimination power of the novel descriptors for vertex- and edge-labeled graphs, we look at Table 1. When applying the partition-independent measures only to skeletons of AG 3982, we see that the sensitivity values are very high, i.e., the corresponding measures possess a high uniqueness. Further, by incorporating edge- and vertex labels, the underlying measures are able to discriminate all graphs uniquely and, hence, . This corresponds to our anticipation that if we incorporate semantical information like edge- and vertex labels, this leads to an increase of the sensitivity measure expressing the uniqueness of the molecular descriptor. We remark that the partition-based measure I W also discriminates the graphs of AG 3982 quite well. In contrast, the discrimination power of W and I orb is comparably very low.
In the following, we evaluate our novel and other descriptors by applying them to supervised machine learning methods [64, 66]. First, our aim is to determine the classification performance of the underlying graph classification problem, i.e., to predict mutagenicity when applying topological descriptors for unlabeled and labeled graphs using SVM and random forests. Second, we examine the influence on the prediction performance when taking semantical (labels) and structural information of the graphs into account. As expressed in a previous section, AG 3982 can be divided into two classes because every graph possesses a unique label (zero or one). Thus, we here deal with a two-class classification problem. Note that by starting from the same underlying benchmark dataset Ames mutagenicity [40, 47], a related study has already been recently performed . However, HANSEN et al.  used the full database (Ames mutagenicity) containing 6512 compounds, molecular descriptors (Dragon ) based on the 3D structure, and supervised machine learning methods (Gaussian processes, RF, SVM, KKN) to predict mutagenicity. In fact, the main goal of this study was to evaluate the prediction performance based on different implementations of the mentioned machine learning algorithms.
Now, before discussing the classification results, we first state some definitions.
Now, based on the SI-representation (see Equation (62)) of a chemical graph, we tackle the mentioned graph classification problem using RF and SVM. The main steps were as follows:
We performed 10-fold crossvalidation for both classification methods.
When doing cross validation, we did a parameter optimization on the corresponding training sets. By using different kernels like linear polynomials, polynomials of higher degree etc., we found that the RBF kernel give the best results.
The random forest was composed by fifty different trees.
We performed the classification both with all features (information measures) and with only seven features determined by running a feature selection algorithm based on greedy stepwise regression .
The results of classification using RF and SVM.
Comparison of the graph classification using unlabeled and labeled graphs
5U + 2L
5U + 2L
Note that in our previously presented classification, we used eleven indices for unlabeled graphs and five for vertex- and edge-labeled graphs. From this feature set, we generated ten subsets composed of seven randomly selected measures for unlabeled graphs (among the eleven), and ten subsets composed of five randomly selected measures for unlabeled graphs and two measures for vertex- and edge-labeled graphs (among five available).
Based on these sets, we again performed 10-fold cross validation with RF and SVM and averaged the classification results.
As a result, Table 3 reflects that if we apply the information-theoretic descriptors for vertex- and edge-labeled graphs, this leads to very similar results (e.g., by considering F-Measure) as in case of only measuring skeletal (structural) information. The calculated standard deviations support this hypothesis. Based on our intuition, we would normally expect that by additionally incorporating semantical information (labels), the graphs can be distinguished more meaningfully. Therefore, the results from Table 3 are astonishing because incorporating the information-theoretic descriptors for vertex- and edge-labeled graphs did not lead to a significant improvement of the prediction performance.
Classification using the substructure method
By looking at the performance evaluation in Table 4, we see again that the representation model based on the superindex led to prediction results which are similar to the ones by applying the model using the appearance or non-appearance of a substructure (see Table 2). From Table 2 and Table 4, we see that if we apply RF and SVM to perform the graph classification, it seems that the used information indices to create the underlying superindex captures structural information of the graphs (contained in AG 3982) similarly than the model that is based on the substructures. But to give a reason why most of the performance measures (mainly F-Measure) in Table 2 are slightly higher than in Table 4, it is plausible to conjecture that the used topological descriptors might measure more complex structural features like branching and other types of structural complexity than only counting the contained substructures.
This paper dealt with investigating several aspects of information-theoretic measures for vertex- and edge-labeled chemical structures. We now summarize the main results of the paper as follows:
We already mentioned that the majority of the topological indices which have been developed so far are only suitable to characterize unlabeled graphs. By adapting the approach of deriving partition-independent information measures, we developed families of information-theoretic descriptors to incorporate vertex- and edge labels when measuring the structural information content of graphs. First, we did this by calculating spherical neighborhoods and distinguishing atom types for every sphere. For the resulting measures, we presented a weighting scheme for the vertices which takes chemical information of the graphs into account. Second, to reduce the number of parameters, we developed a simplified version based on the so-called local information graphs. Generally, these graphs are induced by shortest paths and provide information about the local information spread in a network. We here assume that information spreads out via shortest paths in the network . By using this principle, we defined an information functional (see Equation (29)) that relies on calculating the occurrences of existing and unique vertex labels within the local information graphs and on determining weighted paths. In this paper, we did not give a formal analysis of the computational complexity of the underlying algorithm to compute the corresponding information measures. However, we point out that it is easy to prove that their computation requires polynomial time.
Using the benchmark database AG 3982, we evaluated the novel information-theoretic descriptors to see how they capture structural information of the chemical graphs. Based on some characteristic properties  of the measures, we found that the higher the value of the final measure is, the more disordered is the label configuration of a graph in question. Another aspect we have studied relates to determine their high uniqueness, that is, their ability to discriminate graphs as unique as possible. As a result, we derived that the measures for calculating the information content of vertex- and edge-labeled graphs possess a very high discrimination power. In particular, the computation of two of those led to sensitivity values equal to one, i.e., the measures distinguished all the graphs uniquely.
Another aim was to predict Ames mutagenicity when using supervised machine learning methods (RF and SVM) and representing the graphs by a vector consisting of topological descriptors (superindex). First, we performed the graph classification based on 10-fold crossvalidation and evaluated the quality of the learned models. Taking into account that we only used (i) 16 and (ii) 7 information measures for classifying the graphs, we obtained feasible results (by using SVM, we reached F-Measures of over seventy percent). However, another goal was to examine the influence of incorporating vertex- and edge-labels when measuring the prediction performance of the underlying graph classification problem. Here, we obtained the result that the prediction performance (by calculating the statistical performance measures) was very similar to the one we obtained by only measuring skeletal (structural) information. From this, interesting future work arises as follows: Because of the obtained results, it would be important to explore the developed measures for determining the structural information content (structural complexity) of the underlying vertex- and edge-labeled graphs in depth. This aims to investigate the measures such that the prediction performance could be significantly improved when applying them to the machine learning methods we have used in this paper. Another reason for the results shown in Table 3 could be certain characteristics of the underlying graphs which need to be analyzed more deeply. As further future work, we will use different datasets to determine the prediction performance of the novel measures. Moreover, we want to perform similar analyses by applying our novel descriptors combined with a large number of other well-known molecular descriptors to the same benchmark database. But this goes beyond the scope of this paper.
As already mentioned (see section 'Introduction'), labeled graphs play an important role when analyzing biological networks. But because the theory of labeled graphs is not well developed so far (compared to the contributions which have been done towards unlabeled graphs), see, e.g., , a thorough investigation of methods for analyzing these graphs is therefore crucial. On the other hand, to gain information about the basic biological understanding when investigating biological networks, the problem of exploring their topology is essential [5–7]. Hence, there is a strong need to further investigate methods to analyze labeled graphs for solving problems in bioinformatics and systems biology [22, 38, 39].
Inspired from this study, we think that especially the development of further measures for labeled graphs can be an interesting and valuable attempt not merely to analyze QSPR/QSAR problems. Besides applying these measures to machine learning methods, we believe that the measures itself might be valuable for those who will investigate biological networks, see, e.g., . In fact, if we incorporate also semantical information of the graphs (instead of only considering structural information), this may lead to more meaningful results when developing methods for characterizing graphs or predictive models to tackle problems in bioinformatics, systems biology, and drug design.
As a conclusive remark, we argue from a mathematical point of view that a further development of the theory of labeled graphs will surely help to develop more sophisticated methods for analyzing biological networks, see, e.g., [2, 22, 38, 39]. The next important step is to prove mathematical properties of such measures and to investigate their relatedness. In addition, there is a need to examine correlations to other existing topological indices numerically.
We thank Stephan Borgert and Abbe Mowshowitz for fruitful discussions. In particular, we thank Frank Emmert-Streib for valuable discussions and for helping to improve the present paper. Also, thanks to Katja Hansen for providing the Ames databases and calling our attention to it. This work was supported by the COMET Center ONCOTYROL and funded by the Federal Ministry for Transport Innovation and Technology (BMVIT) and the Federal Ministry of Economics and Labour/the Federal Ministry of Economy, Family and Youth (BMWA/BMWFJ), the Tiroler Zukunftsstiftung (TZS), and the State of Styria represented by the Styrian Business Promotion Agency (SFG) [and supported by the University for Health Sciences, Medical Informatics and Technology and BIOCRATES Life Sciences AG]. Also, funding from the FIRB ITALBIONET Project is gratefully acknowledged.
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.