 Research
 Open Access
 Published:
Characterizing building blocks of resource constrained biological networks
BMC Bioinformatics volume 20, Article number: 318 (2019)
Abstract
Background
Identification of motifs–recurrent and statistically significant patterns–in biological networks is the key to understand the design principles, and to infer governing mechanisms of biological systems. This, however, is a computationally challenging task. This task is further complicated as biological interactions depend on limited resources, i.e., a reaction takes place if the reactant molecule concentrations are above a certain threshold level. This biochemical property implies that network edges can participate in a limited number of motifs simultaneously. Existing motif counting methods ignore this problem. This simplification often leads to inaccurate motif counts (over or underestimates), and thus, wrong biological interpretations.
Results
In this paper, we develop a novel motif counting algorithm, Partially Overlapping MOtif Counting (POMOC), that considers capacity levels for all interactions in counting motifs.
Conclusions
Our experiments on real and synthetic networks demonstrate that motif count using the POMOC method significantly differs from the existing motif counting approaches, and our method extends to largescale biological networks in practical time. Our results also show that our method makes it possible to characterize the impact of different stress factors on cell’s organization of network. In this regard, analysis of a S. cerevisiae transcriptional regulatory network using our method shows that oxidative stress is more disruptive to organization and abundance of motifs in this network than mutations of individual genes. Our analysis also suggests that by focusing on the edges that lead to variation in motif counts, our method can be used to find important genes, and to reveal subtle topological and functional differences of the biological networks under different cell states.
Background
Biological networks describe interactions between molecules such as genes and proteins [1, 2]. These networks are often modeled as graphs where nodes and edges represent molecules and their interactions, respectively. Biological networks are involved in many key biological processes including transcriptional regulation, interactions between a cell and its environment, and controlling a cell’s specificity [3–5]. Understanding biological networks is essential for understanding how cells function. Efforts on computational analysis of biological networks have been growing rapidly in recent years as largescale data collection at low cost is now possible.
One of the most fundamental challenges in computational network studies is the motif counting problem. A network motif is a pattern of local interconnections (i.e., a small subnetwork) observed significantly more frequently in a given network than in a random network of the same size [6, 7]. Existing studies have already uncovered existence of network motifs such as feed forward loop and bifan (see Fig. 1) [8, 9]. Motifs utilize the basic control mechanisms to govern biologically important dynamic behaviors, such as oscillations, generation of molecular pulses, and rapid or delayed responses [7, 10]. Thus, the presence or relative abundance of motifs in biological networks is often used to characterize their topology, function, and robustness [11, 12]. Network motifs have been effectively used to study the biological processes that regulate transcription [9], to find the genetic factors that impact various diseases [13, 14] and to discover new drugs [15].
Identifying motifs and counting them in biological networks is a computationally challenging task as it requires solving the subgraph isomorphism problem, which is NPcomplete [16]. Several methods have been developed to count instances of a motif in a given network [17–20]. These methods could be categorized into two classes [17]. The first class counts all instances of a given motif ignoring the fact that some motifs may share edges (F_{1} measure) [19]. The second class counts all nonoverlapping instances of a given motif, i.e., those which do not share any edge (F_{2} measure). Figure 2 shows the difference between these two frequency measures on a hypothetical network G. Consider the motif pattern M in Fig. 2b. Our input network G in Fig. 2a yields six possible embeddings of M shown in Fig. 2di. Thus, F_{1} measure of M in G is six. However, out of these six embeddings at most two can be chosen without picking the same edge multiple times (e.g., Fig. 2d and i). Thus, F_{2} measure of M in G is two.
Notice that both F_{1} and F_{2} measures make opposite, yet very strong assumptions regarding how the cells realize interactions. The former one assumes that the same interaction can participate in an arbitrary number of motif instances at the same time. The latter one limits the participation of an interaction to a single motif instance. Although, these two assumptions simplify the motif counting problem, they rarely reflect how cells operate interactions. Each interaction utilizes the molecules participating in that interaction. Thus the abundance of the interacting molecules makes it possible to include the corresponding edges appear in multiple (yet limited) number of motif instances. For example, depending on the concentration level of a metabolite, an enzymatic reaction may take place in one (low reaction/edge capacity) or many metabolic pathways (high edge capacity) [21–23]. We use the term capacity of an edge to denote the number of motif instances that edge can participate simultaneously. Thus, even if two cells have the same underlying biological network topology, they may yield different number of motifs of the same topology. For example, if we allow partial overlap of the motif M in the network G (see Fig. 2b and c), we find four possible embeddings of motif M in G (Fig. 2d, g, h and i). A motif counting approach that ignores the edge capacities will lead to unrealistic motif counts, and wrong biological interpretations. Thus, motif counting approaches that take edge capacities into account are needed.
In the literature, several approaches have been developed to count motifs while also incorporating edge information [24–27]. These methods are designed for weighted networks, and find motifs with specific weights. As we explain in the subsequent sections in detail, the problem considered in this paper is fundamentally different (See “Discussion and Conclusions” section for details on motifs in weighted networks).
Our contributions. In this study, we build a new motif counting algorithm that allows partial overlap between different embeddings of a given motif on the target network. Briefly, given a target network, a motif topology, and a positive capacity value for each interaction in the target network, we count the maximum number of ways to place the motif on the target network, so that no edge appears in more motif embeddings than its capacity. Notice that the classical counting measures F_{1} and F_{2} are special instances of our measure [28, 29]. If we set the capacity of all edges to one, our motif counting problem reduces to nonoverlapping motif counting with F_{2} measure. Similarly, if we set the edge capacities to infinity our problem reduces to motif counting using F_{1} measure.
We develop a novel motif counting method called Partially Overlapping MOtif Counting (POMOC), that computes the number of partially overlapping instances of a given motif in a given network. POMOC algorithm first finds all instances of a given motif M in the network G, then it chooses the motif instances that are guaranteed to exist without using any edge more than its capacity. For a given motif embedding, if the capacities of its edges are more than the number of embeddings those edges are part of, this embedding exists in our solution. Next, for the remaining motif instances, our algorithm randomly adds some embeddings whose edges are not more than the capacity constraints, into the resulting embedding set. It gradually expands the resulting set by taking one embedding out of the resulting set and inserting another two embeddings to the set.
Our experimental results on synthetic and real datasets demonstrate that our algorithm finds vastly different motif counts than F_{1} and F_{2} measures. Since biological interactions are resource limited, leading to varying edge capacities, we hypothesize that our method provides a more accurate approach to counting network motifs in biological networks than existing methods. Although, the POMOC method is slightly slower than motif counting with F_{1} and F_{2} measures, our method remains to be practical for all network sizes and motifs we test here. Our results on a S. cerevisiae transcriptional regulatory network suggest that oxidative stress is more disruptive to abundance and organization of network motifs than genetic mutations. Our analysis on the yeast network also suggests that our method can be used to find the key genes, which lead to topological and functional differences in biological networks under varying genetic backgrounds and growth conditions.
The rest of the paper is organized as follows. We present our algorithm in “Methods” section. We experimentally evaluate our method in “Results” section and provide a brief conclusion in “Discussion and Conclusions” section.
Methods
Here, we describe our method for counting partially overlapping motifs in networks. Preliminaries and problem definition section provides the preliminaries needed to describe our method. Counting partial overlapping motifs section discusses our algorithm.
Preliminaries and problem definition
We denote a given biological network with graph G=(V,E,c), where V={v_{1},v_{2},…,v_{n}} and E={e_{1},e_{2},…,e_{m}} represent the set of nodes (molecules) and the set of edges (interactions) among those nodes, respectively. The function \(\mathbf {c}:E \rightarrow \mathbb {Z}^{+}\) shows the capacity of the edges. To simplify our notation, in the rest of this paper, ∀e_{i}∈E we use c_{i} to denote c(e_{i}) (i.e., the capacity of the edge e_{i}) and the vector C=(c_{1},c_{2},…,c_{m}) to represent the capacity of all edges in E in sorted order of edge indices. For example, in Fig. 2c, the value in the form x(y) denotes that edge x has capacity y.
Given a motif pattern M, we represent the ith embedding of M in G with H_{i}⊆E (i.e., H_{i} constitutes a subgraph of G, which is topologically isomorphic to M). We denote the set of all possible embeddings of M in G with \(\mathcal {H}(M)\). Given an edge e_{i}∈E and a subset \(\mathcal {H}^{\prime }\) of \(\mathcal {H}(M)\), we denote the set of all embeddings in \(\mathcal {H}^{\prime }\) containing edge e_{i} with \(\mathrm {f_{i}}(M,\mathcal {H}^{\prime }) = \{ H_{j} \in \mathcal {H}^{\prime } e_{i} \in H_{j} \}\). We denote the number of embeddings in the set \(\mathcal {H}^{\prime } \subseteq \mathcal {H}(M)\) containing each interaction in E with the vector \(C_{\mathcal {H}^{\prime }} = (\mathrm {f}_{1}(M,\mathcal {H}^{\prime }), \mathrm {f}_{2}(M,\mathcal {H}^{\prime }), \dots, \mathrm {f}_{m}(M,\mathcal {H}^{\prime }))\). We say that \(\mathcal {H}^{\prime }\) is feasible if no interaction appears in more embeddings in \(\mathcal {H}^{\prime }\) than its capacity that is \(\forall i, \mathrm {f}_{i}(M,\mathcal {H}^{\prime }) \le c_{i}\). Figure 2b to i explain this on an example. Here, the capacity of the edges is C=(1,2,2,2,1,2,2,1). This network yields six embeddings (see Fig. 2d to i). Consider the subset \(\mathcal {H}^{\prime }=\{H_{1},H_{6}\}\). The capacity \(\mathcal {H}^{\prime }\) uses is \(C_{\mathcal {H}^{\prime }} = (1,1,1,0,0,1,1,1)\), which is less than the imposed capacity constraints in C. Thus, \(\mathcal {H}^{\prime }\) is feasible. The subset of embeddings {H_{1},H_{2}}, however, is not feasible as this set contains edge e_{1} twice, which is more than its capacity (c_{1}=1).
Next, we formally define the partially overlapping motif counting problem.
Definition 1.
(PARTIALLY OVERLAPPING MOTIF COUNTING). Consider a graph G=(V,E,c) conditioned with edge capacity constraints. Given a motif pattern M, partially overlapping motif counting problem seeks to find a largest feasible subset of the set of motif embeddings \(\mathcal {H}(M)\).
Notice that Definition 1 provides a general formulation of motif counting problem. When the capacity constraints of each edge is set to infinity it reduces to a motif counting problem using F_{1} measure. When the capacity of all edges are set to one, it counts nonoverlapping motifs (i.e., F_{2} measure). Next, we present the partially overlapping motif counting problem on an example.
Example 1.
Different ways to select embeddings leads to different number of possible embeddings. Consider the network in Fig. 2c with capacity constraints C=(1,2,2,2,1,2,2,1). Embedding set 1 (see Fig. 2j) includes embeddings H_{2},H_{5} and H_{6}, and is feasible as its usage of the capacity is (1,1,0,2,1,2,1,1). Embedding set 2 (see Fig. 2k) includes embeddings H_{1},H_{4},H_{5} and H_{6}. This set is also feasible. The number of embeddings in this set is four. The partially overlapping motif count is thus four as this is the largest set of feasible embeddings.
Counting partially overlapping motifs is an NP complete problem for several reasons. First, it requires solving the subgraph isomorphism problem, which is NPcomplete [30]. Furthermore, as we explain above the F_{2} count is a special instance of the partially overlapping motif counting problem as we can reduce the nonoverlapping motif counting problem to partially overlapping motif counting by setting the capacity of all edges to one. The Maximum Independent Set (MIS) problem, which is NPcomplete [16], reduces to the nonoverlapping motif counting problem [19, 29]. Thus the partially overlapping motif counting problem requires solving at least two NPcomplete problems. In this paper, we develop a scalable method to tackle this problem using the local search strategy.
Counting partial overlapping motifs
In this section, we discuss our POMOC algorithm. Algorithm 1 presents the pseudocode of our method. Our algorithm takes a network G=(V,E,c) and a motif pattern M as input. Briefly, our algorithm has four main steps: (1) We locate all possible embeddings of M in G (line 1). At this step, we ignore the number of embeddings of M sharing each edge. (2) We determine the embeddings, which are guaranteed to exist in the final solution (lines 25). (3) We construct an initial, random yet feasible, solution by including a subset of the remaining embeddings in the set found in Step 2 (lines 56). (4) We iteratively improve this solution by replacing an embedding in the current solution with two or one new embeddings without violating feasibility of the solution (lines 711).
The first step of our algorithm is identical to computing the F_{1} count. This is a well studied problem in the literature. We use the method developed by Elhesha et al [31] for this step as it is one of the most recent and efficient methods. One can however replace this step with another method for F_{1} count without affecting the rest of our algorithm. Below, we explain Steps 2, 3, and 4 in detail.
Step 2. An embedding H_{r} is guaranteed to exist in the solution set if each edge of H_{r} has large enough capacity to realize all embeddings that have this edge. Formally, H_{r} exists in result set if \(\forall e_{i} \in H_{r}, \mathrm {f}_{i}(M,\mathcal {H}(M))\le c_{i}\). For example, in the network in Fig. 2c, embedding H_{6} (Fig. 2i) satisfies this criteria. We prove the correctness of this step at the end of this section.
Step 3. Once we identify the set of all embeddings, which are guaranteed to be in the final solution, we move them from the set \(\mathcal {H}\) into solution set \(\mathcal {A}\). We then update the capacity constraints in the input graph G as follows. For each embedding \(H_{r} \in \mathcal {A}\), we reduce the capacities of all the edges e_{i}∈H_{r} by one as H_{r} is in the solution set. We then build a new graph, called the overlap graph for the remaining embeddings in \(\mathcal {H}\). Each node in this graph corresponds to an embedding in \(\mathcal {H}\). We include an edge between two nodes if their corresponding embeddings share at least one edge. Figure 3 depicts the overlap graph of the remaining five embeddings after moving embedding H_{6} from \(\mathcal {H}\) into the result set. Next, we generate a random feasible solution iteratively using the overlap graph as follows. At each iteration, we randomly pick a node v_{r} from the overlap graph, and include the corresponding embedding (say H_{r}) in the solution. We then reduce the capacity of all the edges of H_{r} in G by one. If the capacity of an edge drops to zero, it means that that edge cannot participate in any other embedding (say H_{s}) without violating the feasibility of solution. If such an embedding H_{s} exists, its corresponding node in the overlap graph is a neighbor of v_{r} as they share that edge in G. Thus, we remove these neighbors of v_{r} in the overlap graph, which denote an embedding with an edge of zero capacity. We repeat these iterations to grow the random feasible solution set until the overlap graph becomes empty.
Step 4. So far, we have generated an initial feasible solution, which consists of those embeddings, which are guaranteed to be in the final result (Step 2) and those chosen randomly (Step 3). Next, we try to increase the size of this solution set iteratively by applying two strategies, namely local search and perturbation. We elaborate on these strategies next.
Local search iteratively increases the size of the solution set by replacing an embedding in the solution set with two other embeddings while maintaining the feasibility of the solution. For each embedding H_{i} in the solution, we study its corresponding node v_{i} in the overlap graph along with the set of neighboring nodes v_{i}. We find all pairs of nodes (v_{j},v_{k}) in the overlap graph, which satisfy all of the following two conditions.
(i) v_{j} and v_{k} are neighbors of v_{i}.
(ii) Replacing the embedding H_{i} corresponding to v_{i} in the solution set with those of H_{j} and H_{k} corresponding to nodes v_{j} and v_{k} respectively does not violate the feasibility of the solution.
Once we identify all such node pairs, we randomly pick one and replace the embedding for node v_{i} with those of v_{j} and v_{k}. After each swap, we update the capacity of the edges contained in the embeddings H_{i},H_{j} and H_{k}. Notice that each swap operation increases the solution set size by one.
Similar to many local search algorithms, POMOC has the potential to get trapped in the local optimum. To escape the local optimum, before each local search, POMOC perturbs the current best solution to explore a slightly different search space. More specifically, POMOC first selects a node v_{i} such that H_{i} is in the solution set and replaces it with an embedding H_{j} such that v_{j} is a neighbor of v_{i} in the overlap graph if the solution set remains feasible. As a result, the size of the solution set remains unchanged after the perturbation. For each node in the solution, POMOC does a Bernoulli trial with a probability value p to determine if a replacement is needed for this node. Given a solution S, POMOC does S×p replacements approximately. If a replacement is decided, POMOC finds all valid nonsolution neighbors and randomly pick one to replace. Here, valid neighbors means the replacement of these nodes is consistent with the capacity constraints.
We repeat updating the solution by applying local search and perturbation until the size of the solution set does not improve for user supplied number of iterations. Finally, we prove the correctness of Step 2 of our algorithm in the following theorem.
Theorem 1.
Given a graph G=(V,E,c) conditioned with edge capacity constraints and a motif pattern M, consider the set of all possible embeddings of M, \(\mathcal {H}(M)\). An embedding \(H_{r} \in \mathcal {H}(M)\) must be included by partially overlapping motif counting problem, if \(\forall e_{i} \in H_{r}, \mathrm {f}_{i}(M,\mathcal {H}(M))\le c_{i}\).
Proof. We prove the theorem by contradiction. Assume that \(\mathcal {H}^{\prime }\) is the largest feasible subset of the set of motif embeddings \(\mathcal {H}(M)\), which does not have H_{r}. That is, \(\mathcal {H}^{\prime } \subseteq \mathcal {H}(M) \backslash H_{r}\). Now we construct a larger subset of embeddings \(\mathcal {H}^{\prime \prime } = \mathcal {H}^{\prime } \cup \{H_{r}\}\). ∀e_{i}∉H_{r}, we have
In addition, ∀e_{i}∈H_{r}, we have
Thus, \(\mathcal {H}^{\prime \prime }\) is also a feasible subset of embeddings, which yields a contradiction to the assumption that \(\mathcal {H}^{\prime }\) which does not include H_{r} is the largest feasible subset of \(\mathcal {H}(M)\). Thus, an embedding H_{r} must be included by partially overlapping motif counting problem if \(\forall e_{i} \in H_{r}, \mathrm {f}_{i}(M,\mathcal {H}(M))\le c_{i}\). □
Results
In this section, we experimentally evaluate performance of our method on synthetic and real datasets. We consider four motif topologies, which are commonly studied in the literature; namely bifan, biparallel, cascade and feed forward loop (see Fig. 1). These motifs have been shown to be overrepresented in biological networks under the F_{1} count [8, 9]. We compare our method’s motif count and running time to two existing approaches: counting with F_{1} and F_{2} measures. First, we describe in detail the synthetic and real datasets used in our experiments. We then present the results.
Datasets
SYNTHETIC DATASET. We generate directed random networks to test robustness and scalability of our algorithm under three network parameters: network size (i.e., number of nodes), edge capacity, and topology model. We generate synthetic networks of varying sizes and edge capacities using three network topology models: ErdosRenyi random graph (ER) [32], WattsStrogatz smallworld (WS) [33] and BarabasiAlbert preferential attachment (BA) [34]. We randomly assign the direction of each edge. We set a capacity value to each edge using two alternative approaches. In the first approach, we set the capacity of all edges to the same value; 1, 2 or 3. In the second approach, we randomly assign each edge capacity to a value between 1 and 3, thus different edges can have different capacities.
REAL DATASET. We use S. cerevisiae transcriptional regulatory network [8, 35]. This network contains 690 nodes and 1081 edges. We use the S. cerevisiae gene expression dataset, GSE26169, obtained from the GEO database to set capacities of interactions [36]. This dataset contains expression data under control and oxidative stress conditions in seven genetic backgrounds; wildtype, and Glr1, Gpx1, Gpx2, Grx1, Grx2 and Yap1 mutants leading to 14 different conditions (i.e., 2×7). We assign the capacity of each network edge using the capacity of the reactant gene. For each condition, we calculate the capacity of each gene as log(e_{g})/κ, where e_{g} and κ represent the expression level of gene g and capacity constant, respectively. In our experiments, we use κ=2. There are several main reasons behind our choice of capacity function. First, one could replace the gene expression levels with the protein abundance values. However, when the protein abundance data is not available, transcription values are often used as an indicator of the protein abundance. Second, logarithmic transformation is commonly used in studying microarray gene expression data [37] as gene expression values are highly skewed. Logarithmic transformation stabilizes the variance, compresses the range of data and makes the data more normally distributed so it allows statistics to be applied to the data.
In order to observe whether there is any correlation between the topology of the genes in the interaction networks and their transcription (thus the capacity), we create a scatter plot for the distribution of capacity levels and degrees of reactant genes for each mutant and wild type under normal and oxidative stress conditions. Figure 4 presents the results for Glr1 mutant under oxidative stress condition. We observe that there is no correlation between node degree and capacity levels of genes. We observe the same pattern for all the remaining 13 conditions (results not shown).
IMPLEMENTATION AND SYSTEM DETAILS. We implement the POMOC algorithm in C++. We perform all the computational experiments on a Linux machine equipped with Intel core i7 processor 3.6 GHz CPU and 12GBs RAM.
Evaluation on synthetic networks
In this section, we compare POMOC to motif counting with F_{1} and F_{2} measures on synthetically generated networks for varying network size, network topology and edge capacity.
Effects of network size
We generate random networks of varying sizes (200, 400, 800 and 1600 nodes) using ER, WS and BA topology models. We set the average degree of nodes to six, and the capacity for each edge to two. For each network size and model we generate 10 networks, and report the average motif count and running time. Figure 5 reports the results.
Effects of network size on motif count (Fig. 5a to c). We observe that motif count using partial overlaps significantly differs from that of F_{1} and F_{2} measures. Motif counting with F_{2} measure finds substantially lower motif counts than our method for each motif type in all network topologies for almost all motif types. In ER networks, only biparallel motif count shows observable difference between partial overlap and F_{1} measure. In WS networks, feed forward loop and biparallel motifs show the highest differences. In BA networks, all motif topologies show difference; the most significant variation is observed for the biparallel motif. The POMOC method and motif counting with F_{1} measure shows the most significant difference in BA networks. This difference can be explained as follows. ER and WS models generate synthetic networks whose nodes have similar degrees leading to uniform motif distributions in these network topologies. However, BA networks contain hub nodes with high degrees, which leads to nonuniform distribution of motifs to edges. Thus all motif embeddings can be realized in ER and WS networks but not in BA without violating our edge capacity constraints.
All three motif counting methods show the least motif count in synthetic networks generated under the ER model. Motif counts in WS and BA networks are comparable, but substantially higher than those in ER networks. This difference can be explained by the fact that ER model generates networks that are more disconnected compared to WS and BA models. Our results also show that regardless of the network model and counting measure the motif count varies dramatically across different motif topologies. More importantly, the motif count distribution exhibits significant variance across different random network models (see Fig. 5a to c). While biparallel motif is observed most in ER and BA networks, feed forward loop is the most abundant motif in WS networks. Cascade motif is found least in ER and BA networks, but bifan is the least abundant motif in WS networks.
Effects of network size on running time (Fig. 5d to f). Our results demonstrate that counting motifs takes the most time with partial overlap constraint, and least time with F_{1} measure for all network sizes, topologies and motif types. This is not surprising since F_{2} count and partial overlap require solving the F_{1} count as the first step. Furthermore, since F_{2} count enforces identifying nonoverlapping set of motifs, it eliminates the motif embeddings identified in F_{1} count more aggressively than the partial overlap. That said, we observe that the running time of our method is either very close to those of F_{1} and F_{2} measures or remains to be practical for all network sizes and motifs we test.
We observe that the running times for all three measures are greatly affected by the underlying network topology. Although all running times are less than 1 s for ER and WS networks, they go up to 10 s for BA networks. The second factor that influences the running time of our method is the motif topology. While our method’s running time for bifan motif is largest on ER and WS networks, that for biparallel motif is largest on BA networks. Cascade motif takes the shortest time in all three network topologies due to its smaller abundance and simple three node topology.
Three motif counting approaches show network topology dependent differences. In ER networks, only running times for feed forward loop and biparallel motifs show observable difference for the three measures. Similarly, in WS networks, running times for feed forward loop, biparallel and cascade motifs differ. The running time for our method differs from that of the F_{2} measure by a larger margin on BA networks for all motif types and network sizes. The most significant gap is observed for the biparallel motif. Finally, we would like to note that running time is independent of the motif count in ER and WS networks. However, there is a positive correlation between motif count and running time for BA networks; running time increases with increasing motif count.
In summary, our experiments suggest that motif distributions are heavily impacted by the network topology. In particular, BA networks show the most difference between POMOC method, and motif counting with F_{1} and F_{2} measures. Since biological networks often have similar topological characteristics as BA networks [2], we conjecture that our method will be crucial in determining the motif counts for real networks. On the other hand, our experiments also suggest that POMOC method is slightly slower than motif counting with F_{1} and F_{2} measures. While all three network topologies show running time differences for three motif counting approaches, the most significant difference is observed in BA networks. In all three network models, we observe that as network size increases running time increases linearly. Our method, however, scales to large networks (it runs in less than 10 s even for a network with 1600 nodes). Thus, our method is scalable to genome scale biological networks.
Effects of capacity
Here, we compare the performance of our method to that of motif counting with F_{1} and F_{2} measures under varying edge capacity levels. We compare these three approaches using two metrics: motif count and running time. Similar to the previous section, we generate random networks with size 800 and average node degree 6 using ER, WS and BA network topology models. The capacity for each network edge is set using two distinct approaches. In the first approach we set capacity of all edges to 2 (C_{2}) or 3 (C_{3}). In the second approach, for each network edge, we randomly set capacity to an integer between 1 and 3 (C_{R}). Note that motif counting with F_{1} and F_{2} measures correspond to setting the edge capacity to infinity (C_{INF}) and one (C_{1}), respectively. For each edge capacity and network topology model we generate 10 random networks, and report the average and two standard errors of motif count and running time. Figure 6 shows the results.
Effects of capacity on motif count (Fig. 6a to c). Our results suggest that the effects of edge capacity depends greatly on the network topology model. The motif counts for the WS and BA networks are substantially higher than that for the ER network. For the ER network, at all capacity levels, biparallel and cascade motifs are observed the most and the least, respectively. In WS networks, while feed forward loop is observed the most, bifan motif is observed the least in this network topology. In BA networks, biparallel and cascade motifs are observed the most and the least, respectively.
Next, we analyze how the motif counts change as edge capacity increases. We observe that as the capacity level increases the number of motifs found by the POMOC method tends to increase. This is not surprising as higher capacity values make it possible to have more motifs share the same edge without violating the capacity constraints. Motif counting with F_{1}(C_{INF}) and F_{2}(C_{1}) measures show the most and least number of motifs, respectively. Motif counting with varying edge capacities (C_{R}) results in motif counts that are comparable to C_{1} and C_{2}. Increasing edge capacity affects the motif counts in different ways for the three network topologies. In ER and WS networks, number of motifs are comparable between motif counting with F_{1} measure and partial motif counting with an edge capacity of 2 and 3, respectively. This suggests that for these two network topologies, a small edge capacity is enough to find all possible embeddings for each motif type. However, in BA networks, the motif counts are vastly different between motif counting with F_{1} measure and POMOC method. For this network topology, even an edge capacity of 3 is not enough to find all network motif embeddings. This observation can be explained by the fact that ER and WS models generate synthetic networks whose nodes have similar degrees. This leads to uniform motif distribution in these two network topologies with small number of motif overlaps. Thus, all motif embeddings could be found by using small edge capacities. However, the scalefree networks generated by BA model contain hub nodes with high degrees. The distribution of motifs to edges are nonuniform in these networks leading to large number of motif overlaps. Thus, all motif embeddings cannot be realized in BA networks without violating edge capacity constraints even with an edge capacity of 3–resulting in much lower motif counts as compared to the F_{1} measure.
Effects of capacity on running time (Fig. 6d to f). Our results demonstrate that our method is very fast for both fixed and variable capacity values. We observe that POMOC method’s running time for four motifs is affected by the network topology. The running time for all network motifs is less than one second in ER and WS networks. Although finding motifs in BA networks takes slightly longer, running time is still less than 5 s. Counting bifan and cascade motifs take the longest and the shortest time in ER and WS networks, respectively. Our algorithm shows different performance for BA networks. While bifan motif’s running time is highest in motif counting with F_{1} and F_{2} measures, biparallel motif count’s running time is highest for motif counting with edge capacity levels; C_{2},C_{3} and C_{R}. For all capacity levels, cascade motif takes the shortest time in this network topology.
Running time for each motif changes similarly for ER and WS networks; in general motif counting with C_{R} and C_{INF} takes the longest and the shortest time, respectively. Interestingly, our method’s running time does not increase or decrease monotonically as edge capacity increases in these two network topologies. While running time increases from edge capacity of 1 to 2, the opposite behavior is observed as edge capacity increases from 2 to 3. Again, our algorithm shows slightly different behavior in BA networks. The running time increases as edge capacity increases from 1 to 2, but it is comparable for the capacity levels 2 and 3. In ER and WS networks, the running time does not correlate the motif counts. However, BA networks show positive association between motif count and running time.
In summary, our experiments suggest that BA networks will benefit from our partial motif counting method the most. Since, BA model generates scalefree networks that resemble the real biological networks, our method will potentially lead to more accurate motif counts and biological interpretations. On the other hand, our experiments also demonstrate that POMOC method is fast (less than 5 s) for a network of 800 nodes. Thus, we conjecture that that our method can scale to large real biological networks. Next, we test this conjecture by running our method on real networks.
Evaluation on real networks
In this section, we evaluate the performance of the POMOC method on a S. cerevisiae (budding yeast) transcriptional regulatory network. In our analysis, we focus on the effects of edge capacity to motif counting in seven genetic backgrounds (wild type, and Gpx1, Gpx2, Grx1, Grx2, Glr1 and Yap1 mutants) under two experimental conditions (normal and oxidative stress). As explained above, the POMOC method uses gene expression levels to calculate the network edge capacities. Thus, we use our method to count motifs of the yeast transcription network for 14 edge capacity vectors. Tables 1, 2 and 3 and Fig. 7 present our results. The running times for motif counting with F_{1} and F_{2} measures are less than 1 s in all real network experiments. POMOC method’s running time–in average 22 s–is comparable to those of F_{1} and F_{2} measures, and remains to be practical for real biological networks.
Average motif counts
Here, we compare the average motif counts for four motif topologies over the 14 input networks using the POMOC method and, motif counting with F_{1} and F_{2} measures. Table 1 reports the mean and standard deviation of motif count for each counting approach. Note that F_{1} and F_{2} measures do not utilize edge capacity, and thus, their motif counts do not change under varying edge capacity vectors. Our results show that POMOC finds substantially different number of motifs than that with F_{1} and F_{2} measures. We observe that there is a massive variation in the abundance of different motif topologies from only one up to over two thousand. There is only one instance of cascade. Thus, this motif instance occurs in all three measures as it has no other instance to overlap with. More importantly, partially overlapping motif counting varies greatly in the interval defined by F_{1} and F_{2} measures for different motif topologies. For instance, for biparallel, it is closer to F_{1} measure, while it is closer to F_{2} measure for bifan, and equally distant to F_{1} and F_{2} measures for feed forward loop. This demonstrates that the amount of overlap among different embeddings of a motif depends not only on the network topology but also the motif topology. Similar to Milo et al. (2002) [8], we find that bifan and feed forward motifs are the most abundant motif types in this network. These two motif topologies also show small variation across the 14 networks for the POMOC method. Finally, our results show that the motif count distribution in yeast network differs substantially from that of all three synthetic network models we test (see Fig. 5a to c). Biparallel motif is the most abundant in ER and BA networks, and feed forward is the most abundant in WS network. On the other hand, bifan is the most abundant motif in yeast network with a huge margin.
In summary, our partial motif counting approach finds vastly different motif counts than F_{1} and F_{2} measures. Since real biological interactions are resource limited, leading to varying interaction (edge) capacities, we hypothesize that our method provides a more accurate approach to count network motifs, and to decipher dynamics of biological systems under varying conditions. Our method shows little variation in motif counts on the yeast transcription network (see the standard deviation values in Table 1), which suggests that variation in gene expression does not substantially change the edge capacity levels for the yeast network. That said, by focusing on the edges, which lead to such minor variation can reveal subtle differences between the impact of different cell states. Next, we further investigate the difference among three measures, F_{1},F_{2} and POMOC.
Comparison with F _{1} and F _{2} measures
Here, we focus on the genes that exhibit differences in motif count among F_{1} measure, F_{2} measure and POMOC. To do this, for each frequency measure, we obtain the set of genes which are included in the resulting embedding set for each genetic background under oxidative stress condition. Notice that the gene sets identified by F_{1} and F_{2} measures are same across all conditions. They however may differ for POMOC. We represent the gene set for F_{1} and F_{2} measures with \(A_{F_{1}}, A_{F_{2}}\) respectively. For each genetic background, say the ith background, we represent the gene set identified by POMOC with \(A_{P_{i}}\). Notice that, for all i, the gene set \(A_{P_{i}}\) is always a subset of \(A_{F_{1}}\). We first investigate the difference among genes identified by F_{1} measure and POMOC respectively. To do this, we do the literature analysis on genes included by motif instances. For each gene \(\phantom {\dot {i}\!}g_{j} \in A_{P_{i}}\), we count the number of publications containing gene g_{j} and each of the three keywords ‘oxidative’, ‘respiration’ and ‘fermentation’ in PubMed. For each gene \(\phantom {\dot {i}\!}g_{j} \in A_{F_{1}}\), we also count the motif embeddings containing g_{j} identified by F_{1} measure and POMOC respectively. Here, we calculate the latter as the mean value of the number of embeddings containing g_{i} for all seven genetic backgrounds under oxidative stress conditions. Then we calculate the Spearman correlation between publication count and gene specific motif count generated by F_{1} and POMOC respectively. We observe that the gene set identified by POMOC exhibits much higher correlation with publication count (feed forward loop: correlation = 0.265, bifan: correlation = 0.248) than that by F_{1} measure (feed forward loop: correlation = 0.158, bifan: correlation = 0.007), which implies that our method has the potential to filter genes that do not exhibit significant publication evidence. To test our hypothesis, we collect all genes that belong to \(A_{F_{1}}\) but not POMOC, that is \(\bigcup _{i}{A_{F_{1}} \backslash A_{P_{i}}}\). We observe that these genes contribute to limited publication counts compared to the number of these genes especially for feed forward loop (i.e., feed forward loop: number of genes percentage = 25.26%, publication count percentage = 4.86%, bifan: number of genes percentage = 1.88%, publication count percentage = 0.24%). Thus, the genes found by F_{1} but not by POMOC have very low publication evidence to the oxidative stress.
Next, we investigate the difference between F_{2} measure and POMOC. Similar with the comparison between F_{1} measure and POMOC, we collect all genes that belong to POMOC but not in \(A_{F_{2}}\), that is \(\bigcup _{i}{A_{P_{i}} \backslash A_{F_{2}}}\). Notice that the gene set identified by POMOC is often larger than that by F_{2} measure (but not necessarily a superset). Table 2 presents the top five genes with the largest publication count for feed forward loop and bifan. We observe that there is substantial publication evidence for these genes. For example, Mitochondrial NAD^{+}specific isocitrate dehydrogenases (IDHs) are key regulators of flux through biosynthetic and oxidative pathways in response to cellular energy levels [38]. One lysosomal protease which when mutated can cause juvenile onset neuronal ceroid lipofuscinoses (NCL) is encoded in CLN2, a serine tripeptidyl protease with no obvious yeast orthologue [39]. Thus, POMOC has the potential to discover significant genes which are missed using F_{2}. Next, we analyze the association between motif count and edge capacity focusing on individual genes.
Gene specific motif counts
Here, we take a closer look into the network contents and analyze the distribution of motif embeddings to individual genes. To do that, we first find the set of partially overlapping motif embeddings for each of the 14 genetic background and experimental condition combination. For each gene in each network, then, we count the number of motifs containing that gene. Figure 7 plots these counts for all reactant genes that we use to assign the edge capacities in the S. cerevisiae transcriptional regulatory network. In what follows, we discuss the motif count differences and effects of the edge capacity levels to motif counting in seven genetic backgrounds under control and oxidative stress conditions.
Effects of experimental condition to motif count. Our results demonstrate that most of the genes exhibit similar motif counts under varying genetic backgrounds and growth conditions. We observe more genes away from the y=x line in control vs. oxidative stress conditions (Fig. 5a) than wild type vs. mutant genetic backgrounds (Fig. 5b and c). This suggests that oxidative stress is more disruptive to the organization and abundance of motifs in yeast network than mutations of individual genes.
To understand our results more clearly, from now on, we focus on the analysis of the six reactant genes; Abf1, Dal80, Msn2, Msn4, Skn7 and Yap1. Among these, Msn2, Msn4, Skn7 and Yap1 play key roles in controlling the transcription of oxidative stress response genes [40–44]. Dal80 is implicated in nitrogen depletion and amino acid starvation stress responses [45, 46]. Similarly, Abf1 is involved in nutritional stress response [47]. Table 3 reports how changes in the edge capacity levels affect the number of motifs these six genes participate. For this analysis, first we divide the edge capacity and motif count variations to two groups: significant or nonsignificant. We designate any variation that is more than 10% as significant. We then divide all reactant genes to four groups (00, 01, 10 and 11). Here, the first digit represents whether observed variations in edge capacity is significant (1) or not (0), respectively. Second digit denotes whether the deviation in motif count is significant.
First, we focus on motif count variation in seven genetic backgrounds under control and oxidative stress conditions (Fig. 5a and second digits in Table 3a). Our analysis shows that while most of the reactant genes including Dal80 and Skn7 have similar motif counts, some genes such as Msn2, Msn4 and Yap1 show substantially different motif counts. We observe that motif count differences are highly dependent on genetic background: Msn2 and Msn4 show different motif counts in all genetic backgrounds except Grx2 mutant; Yap1 shows different motif counts only in Gpx1, Gpx2, Grx1, Grx2 and Glr1 mutants; Abf1 shows different motif counts only in wild type and, Grx2 and Yap1 mutants.
Motif counts show small variation between wild type and mutant conditions under control and oxidative stress conditions (Fig. 5bc and and second digits in Table 3bc). Our results suggest that genetic mutants show similar gene expression levels to wild type, which leads to comparable edge capacity levels and motif counts. Under control conditions: Msn2 and Msn4 show motif count variation in Grx2 mutant; Yap1, Skn7 and Abf1 show motif count variation in Yap1 mutant. We observe slightly different behavior under oxidative stress conditions: Yap1 and Abf1 show motif count variation in all genetic backgrounds; Skn7 shows motif count variation in Yap1 mutant. Msn2, Msn4 and Dal80 show similar motif counts under oxidative stress in all genetic backgrounds.
Effects of edge capacity to motif count. Here, we analyze how an individual gene’s edge capacities affect the number of motifs it participates. Figure 7 shows that while genes with highly different edge capacity (darker circles) might have no motif count variation (circles on the y=x line); genes with similar edge capacity (lighter circles) might have motif count variation (circles off the y=x line). This implies that variation in edge capacity levels is not a sufficient condition for alterations in motif count. Also small variations in edge capacities may lead to substantial changes in motif counts. We will further analyze these observations on six stress response genes. The two digits in Table 3 reports the dependency between edge capacity and motif counts.
First, we report how oxidative stress affects motif counts. We observe that although oxidative stress leads to edge capacity variation for many genes, most of the genes show similar motif counts under control and oxidative stress conditions in seven genetic backgrounds. Msn2 and Msn4 show small edge capacity variation, but high motif variation in all genetic backgrounds except Grx2 mutant (Table 3a). Abf1 shows similar edge capacity in all genetic backgrounds, however, it shows different motif counts in wild type and, Grx2 and Yap1 mutants. In Grx2 mutant, Dal80 shows edge capacity variation, but similar motif counts. For Yap1 we see a more complex story. It shows edge capacity variation in Gpx1, Gpx2, Grx1 and Glr1 mutants, which leads to different motif counts. However, while Yap1 shows similar edge capacity levels in wild type, and Grx2 and Yap1 mutants, it has different motif counts in Grx2 mutant.
In almost all genetic backgrounds, majority of the genes do not show edge capacity variation between wild type and mutant genetic backgrounds under control condition (Table 3b). Yap1 shows edge capacity and motif count variation in Yap1 mutant. Msn2 and Msn4 show similar edge capacity but different motif counts in Grx2 mutant compared with wild type. Similarly, Abf1 and Skn7 show motif count variation in Yap1 mutant. Oxidative stress leads to slightly more edge capacity variation compared to control (Table 3c). Yap1 and Abf1 show edge capacity and motif count variation in Yap1 and Grx1 mutants compared with wild type, respectively. Dal80 shows edge capacity variation in Gpx1, Gpx2 and Grx1 mutants, however, it has similar motif counts. Despite lack of edge capacity variation, Yap1 shows different motif counts between wild type and, Gpx1 Gpx2, Grx1, Grx2 and Glr1 mutants. Similarly, Abf1 does not show edge capacity variation in Gpx1, Gpx2, Grx2, Glr1 and Yap1 mutants, but, it has different motif counts. Finally, Skn7 shows motif count variation in Yap1 mutant although having similar edge capacity.
In summary, partial motif counting approach finds substantially different motif counts than that of F_{1} and F_{2} measures. However, the changes in edge capacity is not enough to explain the observed motif count differences for individual genes. While some genes show different motif counts due to edge capacity changes, most genes have the same motif count. There are also genes that show similar edge capacity levels but different motif counts. For these genes, the changes in motif counts could be possibly explained by the edge capacity changes of its neighboring genes.
Discussion and Conclusions
Motif counting in biological networks is an important tool to decipher the topology of biological networks and its function. Existing motif counting approaches either count all or nonoverlapping instances of a given motif. This results in either over or underestimation of the motif counts, since biological reactions are constrained by many factors such as concentration levels of reactant molecules. In this paper, we introduced a novel frequency measure which considers the abundance of the interacting molecules. It allows partially overlap between different embeddings based on edge capacities. This problem is fundamentally different from the problem of finding weighted network motifs which also incorporates the edge information. This is because the edge weights in weighted network is used to identify if an embedding has strong strength. And the resulting embeddings can be counted using both F_{1} and F_{2} measure. In our problem, the edge information however is used to specify at most how many embeddings can share this edge, leading to motif count different from both F_{1} and F_{2} measure. Moreover, the resulting embeddings in our problem do not contain weight.
To address this problem, we presented a novel motif counting method, POMOC, based on edge capacities of a given network. Our experiments on both synthetic and real networks demonstrate that motif count using our method significantly differs from the existing motif counting approaches, and our approach extends to largescale biological networks in practical time. Application of our method to the S. cerevisiae transcriptional regulatory networks demonstrates that the POMOC method reveals topological differences of biological networks under different genetic backgrounds and experimental conditions.
Abbreviations
 BA:

Barabasialbert
 ER:

Erdosrenyi
 FFL:

Feed forward loop
 IDHs:

Isocitrate dehydrogenases
 MIS:

Maximum independent set
 NCL:

Neuronal ceroid lipofuscinoses
 POMOC:

Partially overlapping motif counting
 WS:

Wattsstrogatz
 WT:

Wild type
References
 1
Alon U. Biological networks: the tinkerer as an engineer. Science. 2003; 301(5641):1866–7.
 2
Barabasi AL, Oltvai ZN. Network biology: understanding the cell’s functional organization. Nat Rev Genet. 2004; 5(2):101.
 3
Neph S, Stergachis AB, Reynolds A, Sandstrom R, Borenstein E, Stamatoyannopoulos JA. Circuitry and dynamics of human transcription factor regulatory networks. Cell. 2012; 150(6):1274–86.
 4
Bray D. Protein molecules as computational elements in living cells. Nature. 1995; 376(6538):307.
 5
Neves SR, Tsokas P, Sarkar A, Grace EA, Rangamani P, Taubenfeld SM, Alberini CM, Schaff JC, Blitzer RD, Moraru II, et al. Cell shape and negative links in regulatory motifs together control spatial information flow in signaling networks. Cell. 2008; 133(4):666–80.
 6
Shoval O, Alon U. Snapshot: network motifs. Cell. 2010; 143(2):326.
 7
Alon U. Network motifs: theory and experimental approaches. Nat Rev Genet. 2007; 8(6):450.
 8
Milo R, ShenOrr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U. Network motifs: simple building blocks of complex networks. Science. 2002; 298(5594):824–7.
 9
ShenOrr SS, Milo R, Mangan S, Alon U. Network motifs in the transcriptional regulation network of escherichia coli. Nat Genet. 2002; 31(1):64.
 10
Tyson JJ, Chen KC, Novak B. Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. Curr Opin Cell Biol. 2003; 15(2):221–31.
 11
Hayes W, Sun K, Pržulj N. Graphletbased measures are suitable for biological network comparison. Bioinformatics. 2013; 29(4):483–91.
 12
Prill RJ, Iglesias PA, Levchenko A. Dynamic properties of network motifs contribute to biological network organization. PLoS Biol. 2005; 3(11):343.
 13
Lee DS, Park J, Kay K, Christakis N, Oltvai Z, Barabasi AL. The implications of human metabolic network topology for disease comorbidity. Proc Natl Acad Sci. 2008; 105(29):9880–5.
 14
Hu Y, Flockhart I, Vinayagam A, Bergwitz C, Berger B, Perrimon N, Mohr SE. An integrative approach to ortholog prediction for diseasefocused and other functional studies. BMC Bioinformatics. 2011; 12(1):357.
 15
Sridhar P, Kahveci T, Ranka S. An iterative algorithm for metabolic networkbased drug target identification. In: Biocomputing 2007. Singapore: World Scientific: 2007. p. 88–99.
 16
Garey MR, Johnson DS. Computers and intractability: A guide to the theory of npcompleteness (series of books in the mathematical sciences), ed.Comput Intractability. 1979;340.
 17
Inokuchi A, Washio T, Motoda H. Complete mining of frequent patterns from graphs: Mining graph data. Mach Learn. 2003; 50(3):321–54.
 18
Kuramochi M, Karypis G. Frequent subgraph discovery. In: Proceedings of IEEE International Conference on Data Mining. IEEE: 2001. p. 313–20.
 19
Schreiber F, Schwöbbermeyer H. Frequency concepts and pattern detection for the analysis of motifs in networks. In: Proceedings of Transactions on Computational Systems Biology III. Berlin: Springer: 2005. p. 89–104.
 20
Alon N, Dao P, Hajirasouliha I, Hormozdiari F, Sahinalp SC. Biomolecular network motif counting and discovery by color coding. Bioinformatics. 2008; 24(13):241–9.
 21
Bader JS, Chaudhuri A, Rothberg JM, Chant J. Gaining confidence in highthroughput protein interaction networks. Nat Biotechnol. 2004; 22(1):78.
 22
YungChi C, Prusoff WH. Relationship between the inhibition constant (ki) and the concentration of inhibitor which causes 50 per cent inhibition (i50) of an enzymatic reaction. Biochem Pharmacol. 1973; 22(23):3099–108.
 23
Castellani GC, Quinlan EM, Cooper LN, Shouval HZ. A biophysical model of bidirectional synaptic plasticity: dependence on ampa and nmda receptors. Proc Natl Acad Sci. 2001; 98(22):12772–7.
 24
Saramäki J, Onnela JP, Kertész J, Kaski K. Characterizing motifs in weighted complex networks. In: AIP Conference Proceedings. AIP: 2005. p. 108–17.
 25
Onnela JP, Saramäki J, Kertész J, Kaski K. Intensity and coherence of motifs in weighted complex networks. Phys Rev E. 2005; 71(6):065103.
 26
Vassilevska V, Williams R. Finding, minimizing, and counting weighted subgraphs. In: Proceedings of the Fortyfirst Annual ACM Symposium on Theory of Computing. ACM: 2009. p. 455–64.
 27
Choobdar S, Ribeiro P, Silva F. Motif mining in weighted networks. In: 2012 IEEE 12th International Conference on Data Mining Workshops. IEEE: 2012. p. 210–7.
 28
Tran NH, Choi KP, Zhang L. Counting motifs in the human interactome. Nat Commun. 2013; 4:2241.
 29
Kuramochi M, Karypis G. Finding frequent patterns in a large sparse graph. Data Min Knowl Disc. 2005; 11(3):243–71.
 30
Cook SA. The complexity of theoremproving procedures. In: Proceedings of ACM Symposium on Theory of Computing. ACM: 1971. p. 151–8.
 31
Elhesha R, Kahveci T. Identification of large disjoint motifs in biological networks. BMC Bioinformatics. 2016; 17(1):408.
 32
Erdös P, Rényi A. On random graphs, i. Publ Math (Debr). 1959; 6:290–7.
 33
Watts DJ, Strogatz SH. Collective dynamics of smallworld networks. Nature. 1998; 393(6684):440.
 34
Barabasi A. L., Albert R. Emergence of scaling in random networks. Science. 1999; 286(5439):509–12.
 35
Costanzo MC, Crawford ME, Hirschman JE, Kranz JE, Olsen P, Robertson LS, Skrzypek MS, Braun BR, Hopkins KL, Kondu P, et al. Ypd, pombepd and wormpd: model organism volumes of the bioknowledge library, an integrated resource for protein information. Nucleic Acids Res. 2001; 29(1):75–9.
 36
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, et al. Ncbi geo: Archive for functional genomics data setsupdate. Nucleic Acids Res. 2012; 41(D1):991–5.
 37
Speed R, Thompson P. Determinants of sports sponsorship response. J Acad Mark Sci. 2000; 28(2):226–38.
 38
Taylor AB, Hu G, Hart PJ, McAlisterHenn L. Allosteric motions in structures of yeast nad+specific isocitrate dehydrogenase. J Biol Chem. 2008; 283(16):10872–80.
 39
Breitenbach M, Ralser M, Perrone GG, Iglseder B, Rinnerthaler M, Dawes IW. Oxidative stress and neurodegeneration: the yeast model system. Front Biosci (Landmark Ed). 2013; 18:1174–93.
 40
MoyeRowley WS. Regulation of the transcriptional response to oxidative stress in fungi: similarities and differences. Eukaryotic Cell. 2003; 2(3):381–9.
 41
Thorpe GW, Fong CS, Alic N, Higgins VJ, Dawes IW. Cells have distinct mechanisms to maintain protection against different reactive oxygen species: oxidativestressresponse genes. Proc Natl Acad Sci. 2004; 101(17):6564–9.
 42
Temple MD, Perrone GG, Dawes IW. Complex cellular responses to reactive oxygen species. Trends Cell Biol. 2005; 15(6):319–26.
 43
Lee J, Godon C, Lagniel G, Spector D, Garin J, Labarre J, Toledano MB. Yap1 and skn7 control two specialized oxidative stress response regulons in yeast. J Biol Chem. 1999; 274(23):16040–6.
 44
Morano KA, Grant CM, MoyeRowley WS. The response to heat shock and oxidative stress in saccharomyces cerevisiae. Genetics. 2012; 190(4):1157–95.
 45
Neklesa TK, Davis RW. A genomewide screen for regulators of torc1 in response to amino acid starvation reveals a conserved npr2/3 complex. PLoS Genet. 2009; 5(6):1000515.
 46
Wu WS, Chen BS. Identifying stress transcription factors using gene expression and tfgene association data. Bioinforma and Biol Insights. 2007; 1:292.
 47
Fermi B, Bosio MC, Dieci G. Promoter architecture and transcriptional regulation of abf1dependent ribosomal protein genes in saccharomyces cerevisiae. Nucleic Acids Res. 2016; 44(13):6113–26.
Acknowledgements
We would like to thank the reviewers for their insightful comments and suggestions.
Funding
Publication costs are funded by US National Science Foundation (NSF) grant DBI1262451 to TK.
Availability of data and materials
The datasets generated and/or analyzed during this study along with the code for our models and instructions for their use are available under open licenses https://github.com/KahveciLab/partialOverlappingMotif.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 20 Supplement 12, 2019: Selected original research articles from the Fifth International Workshop on Computational Network Biology: Modeling, Analysis and Control (CNBMAC 2018): Bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume20supplement12.
Author information
Affiliations
Contributions
Developed the method: YR, AD and TK. Conceived and Designed the experiments: YR, AA, AD and TK. Performed the data analysis: YR, AA and TK. Performed the experiments and interpreted the results: YR, AA, AD and TK. Contributed to the writing of the manuscript: YR, AA and TK. All authors read, provided comment and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Ren, Y., Ay, A., Dobra, A. et al. Characterizing building blocks of resource constrained biological networks. BMC Bioinformatics 20, 318 (2019). https://doi.org/10.1186/s128590192838x
Published:
Keywords
 Biological networks
 Network motif
 Partial overlap