Enumeration of minimal stoichiometric precursor sets in metabolic networks
- Ricardo Andrade†^{1, 2},
- Martin Wannagat†^{1, 2},
- Cecilia C. Klein^{1, 2},
- Vicente Acuña^{3},
- Alberto Marchetti-Spaccamela^{1, 4},
- Paulo V. Milreu^{5},
- Leen Stougie^{1, 6, 7} and
- Marie-France Sagot^{1, 2}Email author
https://doi.org/10.1186/s13015-016-0087-3
© The Author(s) 2016
Received: 18 February 2016
Accepted: 7 September 2016
Published: 19 September 2016
Abstract
Background
What an organism needs at least from its environment to produce a set of metabolites, e.g. target(s) of interest and/or biomass, has been called a minimal precursor set. Early approaches to enumerate all minimal precursor sets took into account only the topology of the metabolic network (topological precursor sets). Due to cycles and the stoichiometric values of the reactions, it is often not possible to produce the target(s) from a topological precursor set in the sense that there is no feasible flux. Although considering the stoichiometry makes the problem harder, it enables to obtain biologically reasonable precursor sets that we call stoichiometric. Recently a method to enumerate all minimal stoichiometric precursor sets was proposed in the literature. The relationship between topological and stoichiometric precursor sets had however not yet been studied.
Results
Such relationship between topological and stoichiometric precursor sets is highlighted. We also present two algorithms that enumerate all minimal stoichiometric precursor sets. The first one is of theoretical interest only and is based on the above mentioned relationship. The second approach solves a series of mixed integer linear programming problems. We compared the computed minimal precursor sets to experimentally obtained growth media of several Escherichia coli strains using genome-scale metabolic networks.
Conclusions
The results show that the second approach efficiently enumerates minimal precursor sets taking stoichiometry into account, and allows for broad in silico studies of strains or species interactions that may help to understand e.g. pathotype and niche-specific metabolic capabilities. sasita is written in Java, uses cplex as LP solver and can be downloaded together with all networks and input files used in this paper at http://sasita.gforge.inria.fr/.
Keywords
Metabolic network Minimal precursor sets Mixed integer linear programmingBackground
The question of which metabolites an organism needs from its environment (henceforth called the sources) in order to grow or to produce a given set of metabolites (henceforth called the targets) is crucial for both fundamental and applied reasons. This indeed enables to define the growth conditions of organisms in the laboratory, as well as the minimal media necessary for the production of compounds of biotechnological interest (for instance, ethanol). More recently, great interest in establishing which nutrients are exchanged among different organisms in communities such as present in the human gut has also been raised by the interest to develop new strategies for fighting infection that rely on the use of probiotics instead of antibiotics [1]. However the latter requires that: (1) such exchanges are computed in a very efficient way in genome-scale metabolic networks; (2) all possible minimal sets of sources are identified for a given target set of interest in order to fully understand the interactions that may take place among the organisms in a community, as well as the alternative niches that may with time develop for some such organisms.
Early attempts at enumerating all minimal precursor sets (minimal sets of sources) were based only on topology (henceforth called topological precursor sets). Stoichiometry was thus not taken into account, leading to possibly many unfeasible solutions [2–5]. The algorithm of Romero and Karp was based on a backtrack traversing of the metabolic graph from the target compounds to the seeds while Handorf et al. tested the reachability of the target from a heuristically defined collection of sets of sources. Neither enumerated all minimal precursor sets. Cycles, although omnipresent in metabolic networks (e.g. Krebs cycle), were not included until the method of Cottret et al. [4]. However, the latter algorithm could be applied only to small networks due to a high memory requirement; subsequently, Acuña et al. [5] allowed the enumeration of all minimal precursor sets of networks of about 1000 reactions. The authors also pointed out that the enumeration of precursor sets and of precursor cut sets could be done simultaneously in quasi-polynomial total time. Precursor cut sets are a set of sources such that, if they are eliminated, then the target set of interest can no longer be produced by any combination of the remaining sources.
The approach of Zarecki et al. [6] takes stoichiometry into account and consists of two steps. First, the size of a set of sources of minimal cardinality that allows the production of a target is determined solving a mixed integer linear programming problem. In a second step, the authors identify a single set of sources of the determined size such that the sum of the molecular weight of the compounds is minimal.
To our knowledge, there are two algorithms that attempt to enumerate all minimal precursor sets with stoichiometry (henceforth called stoichiometric precursor sets) [7, 8].
Imieliński et al. [7] propose a method that first enumerates all extreme semipositive conservation relations (ESCR), that is the extreme rays of the cone defined by the transposed stoichiometric matrix. The precursor sets are then obtained by the enumeration of hitting sets of the ESCRs. As the authors state, this approach is impractical for genome-scale metabolic networks since it is impossible to enumerate all ESCRs with the current algorithms [7]. Consequently, a method is proposed that enumerates a subset of the ESCRs (those that do not contain water) to obtain (via hitting sets) minimal precursor sets that contain water. These solutions are physiologically minimal (all media contains water), but not necessarily the theoretically minimal.
An approach widely used in flux-balance analysis (FBA) [9] to model an organism’s growth condition is to include a so-called biomass reaction, that consumes in the right amounts every compound needed for a cell to duplicate. Such a reaction has a single product, an artificial compound (representing the duplicated cell) that can be modelled as a target compound in our enumeration approach: a cell can grow and duplicate if it can produce this target compound.
The objective of this paper is twofold. On the theoretical level, we show the relationship between topological and stoichiometric precursor sets and we discuss some complexity results. On the methodological level, we provide two algorithms that enumerate all minimal stoichiometric precursor sets. The first one is based on the above mentioned relationship between topological and stoichiometric precursor sets. Although interesting in terms of theory, it however is not efficient in practice. The second approach, called sasita, uses a similar approach as in [10–12]. Therein the authors enumerate reaction subsets solving recursively mixed integer linear programming (MILP) problems. The reaction subsets correspond to alternate flux distributions to obtain an identical value of the objective function [10], the k-shortest elementary flux modes (EFM) [11] or the smallest minimal cut sets (MCS) [12]. All approaches enumerate minimal reaction sets. Here we consider minimality of the set of compounds.
sasita enables to enumerate all stoichiometrically feasible minimal precursor sets in both models, steady-state and machinery-duplicating, as in [8]. A natural question is to compare our results with those obtained in [8]. Unfortunately, the algorithm of Eker et al. [8] was not publicly available for testing. We nevertheless tried to reproduce the authors’ results by incorporating their definitions in our method, but we were unable to obtain the same results, even using the same Escherichia coli network they made available. More surprisingly, we found a precursor set that is minimal with respect to the solutions found by Eker et al.
In the enumeration of minimal precursor sets for a given target set, Eker et al. [8] are able to work with genome-scale metabolic networks and are exhaustive, but their method is very time and memory consuming. The authors indeed indicate that it required 3 days of execution on a 24-core (with Hyper threading) 2.67 GHz Intel X5650 Xeon CPU-model processor, using the machinery-duplicating model on an Escherichia coli network composed of 2314 unidirectional reactions (no reversible reactions) of which 388 were transport reactions, to enumerate 787 solutions.
We show that we can apply sasita on big networks like the iJO1366 reconstruction of the Escherichia coli K-12 MG1655 with 3646 reactions and 2258 compounds.
In the "Background" section, we provide basic definitions; in particular, we extend ideas from topological precursor sets as defined in [5] in order to incorporate stoichiometry, and we discuss the relationship between topological and stoichiometric solutions. We then describe the sasita algorithm for enumerating all stoichiometrically feasible minimal precursor sets. In the "Results" section, we discuss in detail the comparison with respect to Eker’s et al. proposal. Here we observe that sasita is the first publicly available software to enumerate minimal stoichiometric precursors sets both with the steady state and the machinery duplicating model. Experiments show that sasita can be applied to large genome-scale metabolic networks and we discuss the obtained results.
Methods
Definitions and properties
Topologically, a reaction \(r \in \mathcal {R}\) is defined by its substrates \({\textit{S}ubs(r)}\) and its products \({\textit{P}rod(r)}\), suggesting the interpretation of a reaction as an hyperarc with \({\textit{S}ubs(r)}\) as the set of tail nodes and \({\textit{P}rod(r)}\) as the set of head nodes of a reaction r. In the above example, \({\textit{S}ubs(r_{1})} = \{c_{1},c_{2}\}\) and \({\textit{P}rod(r_{1})} = \{c_{3}\}\). The network \(\mathcal {N}\) can then be seen as a directed hypergraph with stoichiometric coefficients associated with each hyperarc. Given a subset of reactions \(F\subseteq \mathcal {R}\), we denote by \({\textit{S}ubs(F)}\) and \({\textit{P}rod(F)}\) the union of the substrates and products, respectively, of the reactions in F.
Given a network \(\mathcal {N}\) and the sets of source \(\mathcal {X}\) and target \(\mathcal {T}\) compounds, we loosely define a precursor set as a set \(X \subseteq \mathcal {X}\) that can produce all the targets of \(\mathcal {T}\) using a subset of the reactions in \(\mathcal {R}\). The concept of a factory was introduced in [5]; a topological factory is defined as follows:
Definition 1
A set \(F \subseteq \mathcal {R}\) is a topological factory from \(X \subseteq \mathcal {X}\) to \(T \subseteq \mathcal {T}\) if \(T \cup {\textit{S}ubs(F)} \subseteq {\textit{P}rod(F)} \cup X\); i.e., if T and every substrate of every reaction in F is either a source or is produced by some reaction in F.
A set \(X \subseteq \mathcal {X}\) is a topological precursor set (TPS) for \(\mathcal {T}\) if there exists a topological factory from X to \(\mathcal {T}\).
Extending this definition to include stoichiometry requires that any substrate of a reaction in F, should be either produced at least in a same quantity by one or more reactions also in F, or the substrate should be a source compound. Observe that, if the flux vector \(v \in \mathbb {R}^{|\mathcal {R}|}\) denotes the flux of every reaction in the network per time unit, then \(Sv\in \mathbb {R}^{|\mathcal {C}|}\) is the vector of net production of all compounds in the network for the flux v. Furthermore, \((Sv)_{A}\) specifies the net production of the compounds in a set A.
Definition 2
- 1.
\(v_{i} \left\{ \begin{array}{ll} > 0 & \quad i \in F\\ = 0 & \quad \text {otherwise,} \end{array} \right.\)
- 2.
\((Sv)_{\mathcal {C}\setminus X} \ge 0\),
- 3.
\((Sv)_{T} > 0\).
It is possible to adapt Definition 2 to the steady-state assumption by replacing the greater-than-or-equal sign by the equal sign in the second constraint. Thus, the steady-state constraint is put on the compounds in \(\mathcal {C}\) that are neither in X nor in \(\mathcal {T}\).
Definition 3
A set \(X \subseteq \mathcal {X}\) is a stoichiometric precursor set (SPS) of T if there exists a S-factory from X to T. A SPS of T is minimal if it does not contain any other SPS of T.
- 1.
Every S-factory is a topological factory. Every SPS is also a TPS.
- 2.
Not every topological factory is a S-factory. Not every TPS is a SPS.
- 3.
There exist minimal S-factories which are not the union of minimal topological factories.
- 4.
There exist minimal SPSs which do not consist of a union of minimal TPSs.
Figure 2 shows an example where it is indeed not possible to obtain the minimal S-factory that contains r1, r2, r3, and r4 from minimal topological factories. However, it is true that every minimal S-factory is a union of minimal topological factories of the many-to-one transformed network defined as follows:
Definition 4
Given \(\mathcal {N}=(\mathcal {C},\mathcal {R})\), the many-to-one transformation of \(\mathcal {N}\) is the metabolic network \(\Psi (\mathcal {N})=(\mathcal {C},\Psi (\mathcal {R}))\) such that for each reaction \(r\in \mathcal {R}\) and for each metabolite \(a\in {\textit{P}rod(r)}\), there is a reaction \(r_a\) in \(\Psi (\mathcal {R})\) such that \({\textit{S}ubs(r_a)}={\textit{S}ubs(r)}\) and \({\textit{P}rod(r_a)}=\{a\}\).
Given a reaction \(r \in \mathcal {R}\) with \({\textit{P}rod(r)}=\{a_1,\ldots ,a_k\}\), we denote by \(\Psi (r)=\{r_1,\ldots ,r_k\}\) the set of k reactions in \(\Psi (\mathcal {R})\) that correspond to the many-to-one transformation of r, that is, \({\textit{S}ubs(r_i)}={\textit{S}ubs(r)}\) and \({\textit{P}rod(r_i)}=\{a_i\}\). Furthermore, we extend this definition to sets of reactions, that is, if \(R\subseteq \mathcal {R}\), we denote \(\Psi (R)=\cup _{r\in R} \Psi (r)\).
It is clear that all minimal topological factories in \(\mathcal {N}\) are also among the minimal topological factories in the many-to-one network \(\Psi (\mathcal {N})\) if we would retransform the many-to-one reactions into their original hyper-reactions. However, there are additional minimal topological factories in \(\Psi (\mathcal {N})\) that are not minimal in \(\mathcal {N}\) after the retransformation. From this, we claim that every minimal S-factory in \(\mathcal {N}\) is a union of minimal topological factories of \(\Psi (\mathcal {N})\). The following definition and lemmas will provide the basis for the proof of this statement.
Definition 5
- 1.
\(m_0=s\) and \(m_n=t\);
- 2.
\(m_i \in Subs(r_{i+1}),\ i\in \{0,\ldots ,n-1\}\);
- 3.
\(m_i \in Prod(r_i),\ i\in \{1,\ldots ,n\}\).
Lemma 1
Given a minimal S-factory H from \(X \in \mathcal {X}\) to \(T \in \mathcal {T}\) in the network \(\mathcal {N}\) , for every reaction \(r \in H\) , there is always at least one path in H from one of the products of r to some metabolite in \(\mathcal {T}\).
Proof
Let us suppose without loss of generality that \(|T| = 1\). We are going to prove the lemma by contradiction. Suppose that there is a reaction \(r \in H\) such that \({\textit{P}rod(r)}=\left\{ p_{1},\dots ,p_{k}\right\}\), and that for all \(p_{i} \in {\textit{P}rod(r)}\) there is no path from \(p_{i}\) to T in H. Since there is no path to T in H that includes r, if \(r' \in H\) is a reaction that consumes one product of r, then there cannot be a path to T that includes \(r'\) (in fact if such a path exists then there is also a path that includes r). By repeating the same reasoning, consider the set of reactions \(I_{r}\) corresponding to the reactions in H that consume the products of r and of the reactions that consume the products of those reactions and so on. Let \(\overline{H}= H\setminus I_{r}\); we will argue that \(\overline{H}\) remains a stoichiometric factory.
Lemma 2
- 1.
For every metabolite m in \({\textit{S}ubs(H)}\setminus X\) there is exactly one reaction in H that produces m;
- 2.
For every metabolite m in \({\textit{P}rod(H)}\) there exists a path from m to \(t\in \mathcal {T}\) contained in H.
Proof
We first prove that if H is a minimal topological factory from X to \(\mathcal {T}\), then both statements (1 and 2) above hold. By definition, H is a topological factory from X to \(\mathcal {T}\) if and only if any metabolite in \(\mathcal {T}\) and in \({\textit{S}ubs(H)}\) is a source in X or is produced by some reaction in H. Let m be a metabolite in \({\textit{S}ubs(H)}\setminus X\). Then by definition there is a reaction r in H that produces m. Suppose however that there is another reaction \(r'\) in H also producing m. Then \({\textit{P}rod(H\setminus \{r'\})}= {\textit{P}rod(H)}\). Thus, \(H\setminus \{r'\}\) would still be a topological factory from X to \(\mathcal {T}\) which contradicts the minimality and thus proves Statement 1.
Now let m be a metabolite in \({\textit{P}rod(H)}\). We show that there is a path from m to some target in \(\mathcal {T}\). By contradiction, suppose that there is no such path. Consider then the following iterative process. Starting from \(M=\mathcal {T}\) and \(R=\emptyset\), consider all reactions \(H'\) of H that produce the metabolites in M. Then add to R all reactions in \(H'\) and to M all substrates of \(H'\), and repeat the process until no reaction is added. Clearly, for all \(m \in M\), either \(m \in \mathcal {X}\) or m is produced by some reaction in R, and therefore R is a topological factory from X to \(\mathcal {T}\). Since all reactions in \(\mathcal {R}\) are also in H and H is a minimal topological factory, \(R=H\) and m must have been included in M in some iteration. Clearly from that iteration, we can recover a path from m to some metabolite in \(\mathcal {T}\) by going backwards in the described process and therefore Statement 2 above holds.
In order to prove the opposite implication, we first observe that, if both statements are true, then by Definition 1, the set H corresponds to a topological factory from X to \(\mathcal {T}\). Therefore, we only need to show that it corresponds to a minimal topological factory. Let \(H'\subseteq H\) be a minimal topological factory from X to \(\mathcal {T}\). By contradiction, suppose that \(H'\ne H\), then there is a reaction r in \(H\setminus H'\). Let a be the product of r. By hypothesis, there is a path from a to \(\mathcal {T}\) in H. However each reaction in the path is the only one in H producing the metabolites composing its products. Clearly the last reaction in the path (which produces a target) must also belong to \(H'\). Thus, at some point in the path, there is a metabolite which is the product of a reaction in H which is not in \(H'\), and is the substrate of a reaction in \(H'\). There is thus a substrate of a reaction in \(H'\) that is not produced by any reaction in \(H'\), which is a contradiction with the fact that H is a topological factory from X to \(\mathcal {T}\). Therefore, \(H'=H\) and the minimality is proved. \(\square\)
The following theorem shows that any minimal S-factory is the union of minimal topological factories in the many-to-one network.
Theorem 1
- 1.
\(F_{1}, \ldots , F_{k}\subseteq \Psi (H)\);
- 2.
For each reaction r in H there is \(i\in \{1,\ldots ,n\}\) such that \(\Psi (r)\cap F_i\ne \emptyset\).
Proof
From a given minimal S-factory H, we select a reaction \(r \in H\). By Lemma 1, we know that there is a path from at least one product m of r to one target compound t. Clearly, there is a path \(p=(M_p,R_p)\) from m to t in \(\Psi (\mathcal {N})\). Since \(\Psi (\mathcal {N})\) is a many-to-one network, every metabolite in \(M_p\) is produced by only one reaction in \(R_p\). We show that p can be extended to a topological factory from X to t. Starting from the set \(R_0=R_p\), we consider the set of metabolites \(M_0={\textit{S}ubs(R_0)}\setminus {\textit{P}rod(R_0)}\), that is, the set of substrates that are not produced by any reaction in the set. Let c be any metabolite in \(M_0\). In the S-factory H in \(\mathcal {N}\), there exists a reaction h that produces c. Let \(h_c\) be the many-to-one reaction in \(\Psi (h_c)\) that produces c, that is, \({\textit{P}rod(h_c)}=\{c\}\). We define \(R_1=R_0\cup \{h_c\}\) as the new set of reactions and \(M_1={\textit{S}ubs(R_1)}\backslash {\textit{P}rod(R_1)}\). We repeat this process defining \(R_{i+1}\) and \(M_{i+1}\) by choosing any metabolite in \(M_i\) until \(M_{i+1}\) is empty. By construction, the set of reactions \(F=R_{i+1}\) satisfies the two properties of Lemma 2, and therefore F is a minimal topological factory from X to t contained in \(\Psi (H)\). Repeating this process for every reaction \(r\in H\), we obtain a set \(F_{1}, \ldots , F_{k}\) of topological factories from X to t satisfying the desired properties. \(\square\)
Complexity
We now discuss the main complexity results for finding and enumerating \(SPSs\). The next theorem shows that deciding whether a set is a \(SPS\) can be done efficiently.
Theorem 2
Given a network \(\mathcal {N}\) , a subset \(X\subseteq \mathcal {X}\) of sources and a target set \(\mathcal {T}\) , we can decide in polynomial time whether X is a SPS for \(\mathcal {T}\).
Proof
We are going to show that it suffices to solve a linear optimisation problem to decide whether X is a SPS for \(\mathcal {T}\).
If \(v^{*}\) is a solution to M1 and \(f(v^{*}) > 0\) then the support of \(v^{*}\) is a stoichiometric factory from X to \(\mathcal {T}\) and X is a stoichiometric precursor set for \(\mathcal {T}\). \(\square\)
As concerns the problem of enumerating all solutions, we first observe that the proof that enumerating all minimal TPSs cannot be done in polynomial total time (that is, in the size of the input and the number of solutions) unless P = NP given in [5] can be immediately applied to show that enumerating all minimal SPSs cannot be done in polynomial total time unless P=NP. The same observation holds for enumerating all minimal cut sets (SCSs), which we define as follows:
Definition 6
A set \(X \subseteq \mathcal {X}\) is a stoichiometric cut set (SCS) (topological cut set (TCS)), if \(\mathcal {X}\setminus X\) is not a stoichiometric precursor set (topological precursor set).
We now show that the simultaneous enumeration of minimal SPSs and SCSs can be done in quasi-polynomial time. Notice that in [5], a quasi-polynomial time algorithm to simultaneously enumerate all TPSs and TCSs was presented by formulating the problem with a monotone boolean formula and then using a result of [15]. Such approach is possible even in the case of SPSs and SCSs.
Theorem 3
The set of minimal SPSs and the set of minimal SCSs can be enumerated in total quasi-polynomial time.
Proof
Define the Boolean function \(f{:}\,2^\mathcal {S}\rightarrow \{0,1\}\) as \(f(X)=1\) if X is a SPS and \(f(X)=0\) otherwise. Clearly, this function is monotone: if \(f(X)=1\) then \(f(Y)=1\) for any set \(Y\supseteq X\). The collection \(\mathcal {P}\) of minimal SPSs is the collection of all minimal sets in \(\mathcal {S}\) that evaluate to 1 and the collection \(\mathcal {C}\) of minimal SCSs is the collection of all minimal sets whose complement in \(\mathcal {S}\) evaluates to 0. In the context of monotone Boolean functions, minimal SPSs correspond to the prime implicants and minimal SCSs to the prime implicates of f. In [15], a general algorithm is proposed to jointly enumerate prime implicants and prime implicates of any Boolean function. The algorithm and time analysis are rather technical and we only give a brief description of the incremental algorithm applied to our case. Briefly, given two collections of solutions already found, that is, of collections \((\mathcal {P}',\mathcal {C'})\) of SPSs and SCSs, the algorithm finds a set \(X\subseteq \mathcal {S}\) such that X does not contain any minimal SPS in \(\mathcal {P}'\) and \(\mathcal {S}\setminus X\) does not contain any minimal SCS in \(\mathcal {C'}\) (or proves that such set does not exist). Since either X is a SPS or \(\mathcal {S}\setminus X\) is a SCS, we have found a new solution not in \((\mathcal {P}',\mathcal {C'})\). Such a new solution is found in time \(O(n(\tau +n))+m^{O(\log m)}\) where \(n=|\mathcal {S}|\), m is the number of partial solutions already found (i.e. \(m=|\mathcal {P}'| +|\mathcal {C}'|\)) and \(\tau\) is the time needed to evaluate f. Since \(\tau\) is polynomial, we conclude the proof. \(\square\)
Relation to previous work
The paper in the literature that comes closest to ours is [8]. In fact, one of their definitions coincides completely with our definition of SPS. However, their work concentrates on a more restrictive model, which they call machinery-duplicating. The underlying idea of the latter is that each compound involved in a path from the precursor set to the target set should be produced in strictly positive amount, allowing a cell to therefore duplicate itself. We translate their definition by using the concept of factory (cf. Definition 2).
Definition 7
- 1.
\(v_{i} \left\{ \begin{array}{ll} > 0 &{} \quad i \in F\\ = 0 &{} \quad\text {otherwise,} \end{array} \right.\)
- 2.
\((Sv)_{\mathcal {C}\setminus X} \quad \ge 0\),
- 3.
\((Sv)_{T \cup Y} \quad > 0\).
Comparing this definition to Definition 2, clearly any MD-SPS is a SPS, but not the other way around. Moreover, not every minimal MD-SPS is a minimal SPS. In their work, the authors claim that for the growth of a colony of cells, one must consider that not only the biomass compounds should be produced in positive amount, but also all the reactants of every reaction with nonzero flux that are not sources. However, as we already mentioned, cycles like the one in Fig. 1 are considered unfeasible according to the machinery-duplicating model. Yet cycles with this structure are present in real networks and play an important role in metabolism, such as in the urea or the Krebs cycle.
Enumerating precursor sets via MILP
In the following section, we describe how to enumerate all minimal \(SPS\) and MD-SPS using a MILP approach similar to [10–12]. The authors of these papers describe methods that enumerate reaction subsets by recursively solving MILP problems. Therein, solutions obtained in a previous step are excluded from the solution space.
Enumeration of minimal \(SPS\)
We now present a practical method to enumerate all minimal stoichiometric precursor sets that allow to produce the set \(\mathcal {T}\) in a positive amount. We iteratively solve a series of optimisation problems: at each iteration a mixed integer linear programming (MILP) problem is solved to obtain a minimal precursor set X; then we define a new MILP by adding a constraint that removes the obtained solution X and all the sets that contain it from the feasible set. We keep repeating this process until all solutions are found.
We need some additional definitions. For each source compound \(x_j \in \mathcal {X}\), we add to \(\mathcal {R}\) a reaction, which we call source-pool reaction, that produces \(x_j\) from nothing (with stoichiometric coefficient 1). We denote this new set by \(\overline{\mathcal {R}}\) and the set containing all source-pool reactions by \(\overline{\mathcal {R}}_{\mathcal {X}}\). This set of reactions allows to model the availability of the source compounds since the upper bounds on their fluxes are linked to the amount of each source that is available. In the sequel \(\overline{S}\) denotes the stoichiometric matrix S obtained by adding the columns given by the set of reactions \(\overline{\mathcal {R}}_{\mathcal {X}}\) and v the flux vector, U is an upper bound constant for the values of each flux, \(\epsilon\) is a vector of size \(|\mathcal {T}|\) with an arbitrarily small positive real number in all coordinates, b is the vector of binary variables associated with each compound in \(\mathcal {X}\), and we assume that \(b_j=1\) (\(b_j=0\)) implies that compound \(x_j\) is used (not used) to produce the target.
If the obtained problem has no feasible solution, then we claim that we have found all minimal precursor sets.
MILP constraints for MD-SPS
In the work of Eker et al. [8], the machinery-duplicating model is defined through the use of linear constraints and boolean operators. If a set of sources is a MD-SPS, this implies it is a feasible solution according to their model. The authors also present a method to enumerate all MD-SPSs. Suppose we are given a set \(\{X_1,\ldots , X_k\}\) of precursor sets that were already found. Their method consists in finding a minimal subset of sources Y that verifies two conditions: (1) Y has at least one source in common with each precursor set in \(\{X_1,\ldots ,X_k\}\); and (2) the complement of Y must be able to produce the target according to the machinery-duplicating model. If one can find such a subset Y, a minimal precursor set can be obtained by taking the complement \(\overline{Y}\) of Y, and finding one minimal subset of \(\overline{Y}\). If no Y verifying the above conditions can be found, all minimal precursor sets have been enumerated and the algorithm stops.
Results and discussion
In this section, we present the experiments we realised and discuss the results we obtained. We start by comparing the method we developed with the one of Eker et al. [8]. We then show the performance of sasita versus the approach where minimal \(SPSs\) are obtained from combinations of minimal topological factories in the many-to-one network. Finally, we apply sasita to some genome-scale metabolic networks, obtained from Monk et al. [16]. The objective of this last part is both to illustrate how our method can be used and to validate it by reproducing the findings of the authors.
All the experiments were performed using an Intel QuadCore i7-4770 computer with 16 GB of RAM memory. The algorithm sasita is coded in Java (OpenJDK IcedTea) and uses cplex (IBM ILOG AMPL/CPLEX 12.5.1) for solving the MILP models; the constants are fixed as follows: \(\epsilon = 0.5\), \(\overline{\epsilon }= 0.5\), \(U = 1000.0\). The constraints (2) were coded using indicator constraints to avoid numerical instability. The software and all network and input files can be downloaded at http://sasita.gforge.inria.fr/.
Comparison between sasita and Eker et al.’s approach
We start by calling attention to the fact that the comparison with the method of Eker et al. was difficult due to the fact that it is not publicly available. We also were not able to obtain it upon request. We therefore implemented a version of sasita that enumerates all minimal MD-SPS using the constraints given by Eq. (5). As input we took the metabolic network, the set of sources \(\mathcal {X}\) and the set of targets \(\mathcal {T}\) provided in the supplementary material of Eker et al. [8]. The authors provided also a list of “auxiliary compounds” without which, according to them, their model does not work. No auxiliary compound appears in the minimal precursor sets that are enumerated by Eker et al. It is not clear how these compounds are handled in their approach. If we treat such auxiliary compounds as ordinary ones, we are not able to enumerate a single MD-SPS with sasita. If we add a source-pool reaction for each one of the auxiliary compounds, we obtain the minimal MD-SPS X = {CCO-PERI-BAC@SULFATE}. Eker et al. find 787 solutions and all of them contain Sulfate. So the minimal solution X we found is in fact a subset of all their solutions.
We provide in the Additional file 1 a list of reactions F that form a MD-stoichiometric factory from X to \(\mathcal {T}\), the flux values in F, and the stoichiometric matrix restricted to the reactions in F. Furthermore, we show that all substrates of the reactions in F and the target set \(\mathcal {T}\) are produced in a positive amount using the reactions in F. Hence, the minimal MD-SPS fulfils the properties of a precursor set according to the machinery-duplicating model [8]. Such minimal MD-SPS is not found by Eker et al. probably because they do some preprocessing on the network that is not described in their paper and that we were not able to obtain upon request.
Comparison between sasita and combinatorial approach
We ran both approaches, i.e. sasita and a combinatorial approach (called combi) where minimal \(SPSs\) are obtained from combinations of minimal topological factories in the many-to-one network, on several instances. Our objective was to analyse the differences in the running times between both approaches, so we set cplex into single thread mode for sasita. Table 1 shows clearly that the MILP approach is more efficient than the combinatorial one. The networks of S. muelleri, C. ruddii and B. aphidicola were obtained from metexplore, filtering out ubiquitous metabolites and pairs of co-factors. We obtained the E. coli core model from http://www.systemsbiology.ucsd.edu/InSilicoOrganisms/Ecoli/EcoliSBML. As sources we considered all compounds that are not produced by a reaction or those that are produced by reversible reactions only. For the E. coli strains, we used the same networks from Monk et al. [16] and considered as sources the compounds from Table 3.
Our MILP approach (sasita) versus the combinatorial one
Strain | #Compounds/#reactions | #Sources/#targets | \(t_{sasita}\) (s) | \(t_{Combi}\) (s) |
---|---|---|---|---|
S. muelleri | 76/64 | 9/Pyruvate | <1 | <1 |
C. ruddii | 128/126 | 45/Pyruvate | <1 | 1 |
B. aphidicola | 282/245 | 91/l-Histidine | 2 | 2 |
E. coli core^{a} | 72/126 | 14/Biomass core | <1 | 42 |
E. coli core | 78/126 | 14/Biomass core | <1 | * |
E. coli CFT073 | 1911/2949 | 26/Biomass core | 12 | * |
E. coli EDL993 | 1895/2943 | 25/Biomass core | 13 | * |
E. coli K-12 | 1806/2854 | 25/Biomass core | 12 | * |
E. coli Sakai | 1895/2942 | 25/Biomass core | 12 | * |
Enumerating minimal precursor sets in genome-scale metabolic networks
In this case, we based our experiments on the work of Monk et al. [16] who investigated the pan and core metabolic capabilities of 55 Escherichia coli and Shigella strains based on genome-scale reconstructions of their metabolism. By core is meant the elements shared by all strains and by pan the union of the elements from all strains. As concerns the latter in particular, the authors found the pan to be enriched in alternate carbon metabolic pathways. In order to determine the functional differences among the strains, the authors computed by flux balance analysis (FBA) the growth phenotypes of 385 nutrients (henceforth called the test metabolites/compounds), each considered individually as a source of carbon, nitrogen, phosphorus and sulfur, aerobically and anaerobically. To that purpose, an in silico minimal medium that contains a sole carbon, nitrogen, phosphorus and sulfur source was defined. The authors then replaced the sole carbon source by each of the 385 test metabolites one at a time. Whether or not these new media constituted a growth condition was tested by FBA. The procedure was repeated for each source in the minimal media, namely for nitrogen, phosphorus and sulfur, as well as for each strain. The resulting metabolic phenotypes indicated strain-specific adaptation to nutritional environments.
Our first goal was to validate our method: we compared the results obtained with sasita to the ones in Monk et al. [16]. We enumerated and compared the minimal precursor sets allowing for biomass production of the E. coli strains, which included commensals as well as both intestinal and extraintestinal pathogens. We used for this the genome-scale metabolic models from Monk et al. [16]. The strains were E. coli str. K-12 MG1655 (Commensal), E. coli O157:H7 str. Sakai (Enterohemorrhagic E. coli, EHEC), E. coli O157:H7 EDL933 (EHEC), and E. coli CFT073 (Uropathogenic E. coli, UPEC). The same 385 compounds tested in [16] were given as part of the sources for different runs of sasita. Since Monk et al. [16] were not interested in minimal solutions, we wanted to check whether our solutions were subsets of their solutions.
The second goal was to explore some solutions that were only found by sasita in order to illustrate one application of our method. These solutions contain more than one of the test metabolites, after excluding sources from the minimal media (i.e. carbon, nitrogen, phosphorus and sulfur). Such solutions were explored as concerns strain-specific growth and their relation to niches and to pathotypes.
Differences between solutions found by Monk et al. [16] and by sasita
Network | Matches | Not found |
---|---|---|
E. coli CFT073 | 599 | 0 |
E. coli O157:H7 EDL933 | 597 | 0 |
E. coli str. K-12 MG1655 | 607 | 7 |
E. coli O157:H7 str. Sakai | 597 | 0 |
Two experiments were conducted. In both cases, oxygen was always available and we used as target an artificial compound that is added as an extra product of the core biomass reaction, with stoichiometry of 1.0. Also, the minimal media for E. coli CFT073 contained tryptophan, its auxotrophy. We now give a general description of each experiment. The exact list of compounds for each experiment as well as all networks can be found in the sasita website.
Minimal media compounds for each E. coli strain
Strain | Minimal media |
---|---|
E. coli CFT073 | Calcium, cob(I)alamin, chloride, \(\text{Co}^{2+}\), \(\text{Cu}^{2+}\), \(\text{Fe}^{2+}\), \(\text{Fe}^{3+}\), \(\text{H}^{+}\), potassium, magnesium, \(\text{Mn}^{2+}\), molybdate, calcium, nickel, selenate, selenite, tungstate, zinc, 6-Acetyl-d-Glucose, sulfate, ammonium, diphosphate, nicotinate, l-Tryptophan, \(\text{O}_{2}\) |
E. coli EDL993 , K-12, Sakai | Calcium, cob(I)alamin, chloride, \(\text{Co}^{2+}\), \(\text{Cu}^{2+}\), \(\text{Fe}^{2+}\), \(\text{Fe}^{3+}\), \(\text{H}^{+}\), potassium, magnesium, \(\text{Mn}^{2+}\), molybdate, calcium, nickel, selenate, selenite, tungstate, zinc, 6-Acetyl-d-Glucose, sulfate, ammonium, diphosphate, nicotinate, \(\text{O}_{2}\) |
In order to check the ability of each strain to grow on the 385 test compounds as sources of carbon, nitrogen, phosphorus and sulfur, we ran Experiment 2. For each network we considered as input source set a subset of the minimal media compounds plus a subset of the 385 test metabolites.
For subsets of the minimal media compounds, we considered the set of minimal media compounds minus one of the following: glucose, ammonium, phosphate or sulfate respectively. Since we were removing a compound from the minimal media, we included only the test metabolites that could replace the removed one (if we removed glucose, we considered only the set of test compounds that have carbon in their composition and so on). This was done because considering all the 385 test metabolites together leads to a combinatorial explosion of the number of solutions that are unpractical to enumerate with sasita. In one case, namely when we remove glucose, the set of test metabolites to include was also too big and we needed to split it further in two smaller sets. As a side effect, this split of the input compounds can lead to a loss of some solutions, namely those containing compounds that are in different input sets. We thus may lose some solutions that have more than one minimal media compound replaced by two or more test metabolites. However, our split guarantees that at least all the sets considered in Monk et al. [16] are possible combinations of our input compounds because the authors replace glucose, ammonium, phosphate or sulfate from the minimal media with only one of the test metabolites.
Number of solutions found for each E. coli strain
Strain | Solutions |
---|---|
E. coli CFT073 | 837 |
E. coli EDL993 | 11.164 |
E. coli K-12 | 13.732 |
E. coli Sakai | 11.164 |
This difference is remarkable but not surprising: E. coli CFT073 has a tryptophan auxotrophy, and tryptophan itself can be a source of carbon and nitrogen. Since we search for minimal precursor sets and tryptophan is always a source, there is no need for any extra source of carbon and nitrogen, thus reducing the number of solutions for this strain.
Number of solutions with one test metabolite
Sources | E. coli CFT 073 | E. coli EDL933 | E. coli K-12 | E. coli Sakai | ||||
---|---|---|---|---|---|---|---|---|
\(\text {O}_{2}\) | No \(\text {O}_{2}\) | \(\text {O}_{2}\) | No \(\text {O}_{2}\) | \(\text {O}_{2}\) | No \(\text {O}_{2}\) | \(\text {O}_{2}\) | No \(\text {O}_{2}\) | |
C | 0 | 0 | 104 | 51 | 109 | 61 | 104 | 54 |
C, N | 0 | 0 | 51 | 42 | 52 | 44 | 51 | 42 |
C, N, P | 0 | 0 | 37 | 22 | 38 | 22 | 37 | 22 |
C, N, S | 0 | 0 | 8 | 0 | 8 | 0 | 8 | 0 |
C, P | 0 | 0 | 14 | 22 | 13 | 21 | 14 | 22 |
C, S | 0 | 0 | 0 | 0 | 2 | 0 | 0 | 0 |
N | 0 | 0 | 11 | 14 | 14 | 10 | 11 | 14 |
P | 51 | 51 | 6 | 6 | 7 | 7 | 6 | 6 |
S | 22 | 0 | 14 | 0 | 10 | 0 | 14 | 0 |
Minimal media | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
Total | 74 | 52 | 246 | 158 | 254 | 166 | 246 | 161 |
Total | 126 | 404 | 420 | 407 |
These solutions contain one test metabolite plus a subset of the minimal media in which the test metabolite replaces one, or more of the following compounds: glucose, ammonium, phosphate and sulfate. They therefore correspond to minimal sets of the solutions found by Monk et al. [16].
Other differences found in the comparison of our results with those from [16] are presented and discussed below. Some such differences arise for E. coli CFT073 for which, as mentioned, tryptophan can always be a source of carbon and nitrogen. Since we find a minimal solution without any test metabolite and without glucose and ammonium, but with tryptophan, it is a minimal subset of all solutions from [16] considering the replacement of glucose or ammonium by a test metabolite. Furthermore, a few minimal precursor sets for which Monk et al. [16] found no growth are present among our solutions because of the different conditions we allowed in our test, namely some sources were available at bigger amounts and the compounds were allowed to accumulate. There were as well seven solutions for which Monk et al. [16] found a positive flux and we did not (see Table 2). Those solutions are for the K-12 strain. Among these solutions, 6 are related with the compounds 4-Hydroxy-l-threonine and Oxaloacetate from the exchange subsystem, and in fact there are no reactions in the network that use those compounds to produce anything. The remaining solution is the one using Thiosulfate as source of sulfur, and we confirmed no growth by FBA for this condition. There is one last difference, for the solution with the test metabolite Fe(III) dicitrate which allows growth as a carbon source in aerobic and anaerobic conditions for E. coli K-12. We explicitly found only the anaerobic solution. Since the aerobic one can be seen as a superset of the anaerobic, it is not a minimal precursor set. This does not happen in the other test metabolites owing to different iron oxidation states in each solution. We thus found different aerobic and anaerobic solutions when we enumerated the minimal precursor sets.
In conclusion, we obtained almost all of the solutions found by Monk et al. [16], showing that the nutrient sources of alternate catabolic pathways are part of the minimal precursor sets that allow for biomass production in most of the tested E. coli strains.
Second goal: Original results using sasita
Solutions with more than one test metabolite
Test metabolites | Solutions | E. coli strains |
---|---|---|
2 | 710 | CFT073 |
2 | 3.387 | EDL993;K-12;Sakai |
2 | 298 | EDL993;Sakai |
2 | 656 | K-12 |
3 | 4.992 | EDL993;K-12;Sakai |
3 | 703 | EDL993;Sakai |
3 | 2.106 | K-12 |
4 | 1.027 | EDL993;K-12;Sakai |
4 | 299 | EDL993;Sakai |
4 | 1.113 | K-12 |
5 | 19 | EDL993;K-12;Sakai |
5 | 21 | EDL993;Sakai |
5 | 5 | K-12 |
6 | 5 | EDL993;K-12;Sakai |
6 | 4 | EDL993;Sakai |
Most of these solutions actually have two or three test metabolites (4.961 and 7.801 respectively). In both cases, more than 60% of those solutions are aerobic.
Solutions with two or three test metabolites
Test metabolites | Solutions | E. coli strains |
---|---|---|
2 | 674 | CFT073 |
2 | 556 | EDL993; K-12; Sakai |
2 | 6 | EDL993; Sakai |
3 | 2.641 | EDL993; K-12; Sakai |
3 | 143 | EDL993; Sakai |
3 | 511 | K-12 |
From these results, the pairs of test metabolites in the six solutions specific to E. coli EDL993 and E. coli Sakai are: N-Acetyl-d-galactosamine 1-phosphate with butanesulfonate, ethanesulfonate or taurine, respectively; all aerobic and each one in two solutions with different iron states (line in italics in Table 7). N-Acetyl-d-galactosamine 1-phosphate was shown to give extraintestinal pathogenic strains of E. coli a catabolic advantage when compared to commensals, supporting growth in 100% of the cases compared to 67%, respectively [16]. This compound supported growth as a sole carbon source in aerobic and anaerobic conditions for extra and intracellular pathogens (E. coli CTF073 together with tryptophan, E. coli EDL993 and E. coli Sakai), however not for the commensal strain E. coli K-12. Furthermore, the enterohemorrhagic strains E. coli EDL933 and E. coli Sakai were shown to occupy the same niche in the streptomycin-treated mouse intestine [17] while E. coli EDL933 was shown not to colonise the same niche and does not use the same sugars as carbon source as the commensal E. coli K-12 [18, 19]. The three solutions detailed above and the solutions presented in Table 7 therefore represent metabolic capabilities that are specific to the pathogenic strains analysed here when compared to the commensal strain E. coli K-12, in agreement with the niches occupied by such strains.
These results suggest that sasita can depict pathotype and niche-specific metabolic capabilities which allow broad in silico studies of strains or species interactions. For instance, an extension of the analysis presented in this paper to a larger dataset of E. coli strains including both pathogenic and commensal biotypes could help predict in silico sets of commensal strains that would prevent the colonisation of pathogens due to a consumption by the native microbiota of the nutrients required by the pathogen (see the experimental study of mutant phenotypes in Maltby et al. [20]).
Conclusions
We examined the relationship between topological and stoichiometric precursor sets. We highlighted that stoichiometric precursor sets can be obtained from combinations of minimal topological factories in the many-to-one network. However, this does not lead to an efficient method. We then presented sasita, an efficient algorithm for the exhaustive enumeration of minimal precursor sets for a given target that takes into account stoichiometry. To the best of our knowledge, there exists only one previous approach for this problem due to Eker et al. [8] who proposed two different constraint models, steady-state and machinery-duplicating. However, in their computations, the authors use only the latter, that requires a strictly positive net production of the intermediate compounds on the path from the sources to the target. This model may exclude solutions as we showed (Fig. 1).
In our experiments, we enumerated and compared the minimal precursor sets of nutrient sources of alternate catabolic pathways allowing for biomass production of some Escherichia coli strains, comprising commensal and both intestinal and extraintestinal pathogens, using genome-scale metabolic models. We compared our results to those of Monk et al. [16] in order to have a guideline on part of the solutions we generated, since our approach is different from the one that the authors used, and our results go beyond such comparison. We found metabolic capabilities that distinguish the strains compared in their ability to catabolise nutrients, and such were specific to pathotypes and niches of E. coli strains.
Our method can therefore be used in a wide variety of applications in order to study minimal growth conditions as well as strains and/or species interactions based on their catabolic abilities and their nutritional niches. One valuable application in this context would be to predict patterns of colonisation of commensal and pathogenic E. coli strains in the intestine.
Our method can furthermore be used to refine a metabolic network. If growth of an organism is observed for a defined medium in the laboratory, but no minimal precursor set is a subset of such medium, then either the metabolic network lacks reactions, e.g. export reactions, or the biomass function is not well formulated.
The execution of all experiments of the comparison with Monk et al. [16] took altogether around 5 days. The execution times ranged from 20 s to 12 h. Running the experiments in parallel, one could such retrieve the results after less than a day. We cannot claim to be more efficient than Eker et al. [8] as their software is not available for testing. However, we are guaranteed to enumerate all minimal \(SPSs\) and MD-SPS.
Notes
Declarations
Authors' contributions
RA and MW designed and implemented the algorithms, RA, MW and CCK performed the tests and analyzed the results, RA, MW, VA, AMS, PVM, LS, MFS wrote and revised the theorems and proofs, MFS supervised the work. All authors contributed to the writing of the paper. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Funding
RA acknowledges CNPq/Brasil for the financial support and VA for the support of Fondecyt 1140631, CIRIC-INRIA Chile and FONDAP Center for Genome Regulation. CCK and MW were recipients of a grant from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007-2013)/ERC Grant Agreement No [247073]10 SISYPHE.
Open AccessThis 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.
Authors’ Affiliations
References
- Lin C, Chang C-J, Lu C-C, Martel J, Ojcius D, Ko Y-F, Young J, Lai H-C. Impact of the gut microbiota, prebiotics, and probiotics on human health and disease. Biomed J. 2014:37(5):259–268. doi:https://doi.org/10.4103/2319-4170.138314.View ArticlePubMedGoogle Scholar
- Romero P, Karp P. Nutrition-related analysis of pathway/genome databases. In: Pacific symposium on biocomputing’01. 2001. p. 470–82.Google Scholar
- Handorf T, Christian N, Ebenhöh O, Kahn D. An environmental perspective on metabolism. J Theor Biol. 2008;252:530–7. doi:https://doi.org/10.1016/j.jtbi.2007.10.036.View ArticlePubMedGoogle Scholar
- Cottret L, Milreu PV, Acuna V, Marchetti-Spaccamela A, Martinez FV, Sagot MF, Stougie L. Enumerating precursor sets of target metabolites in a metabolic network. In: Lecture notes in bioinformatics, vol 5251. 2007. p. 233–44.Google Scholar
- Acuña V, Milreu PV, Cottret L, Marchetti-Spaccamela A, Stougie L, Sagot M-F. Algorithms and complexity of enumerating minimal precursor sets in genome-wide metabolic networks. Bioinformatics. 2012;28(19):2474–83.View ArticlePubMedGoogle Scholar
- Zarecki R, Oberhardt MA, Reshef L, Gophna U, Ruppin E. A novel nutritional predictor links microbial fastidiousness with lowered ubiquity, growth rate, and cooperativeness. PLoS Comput Biol. 2014;10(7):1–12. doi:https://doi.org/10.1371/journal.pcbi.1003726.View ArticleGoogle Scholar
- Imieliński M, Belta C, Rubin H, Halász Á. Systematic analysis of conservation relations in Escherichia coli genome-scale metabolic network reveals novel growth media. Biophys J. 2006;90(8):2659–72. doi:https://doi.org/10.1529/biophysj.105.069278.View ArticlePubMedPubMed CentralGoogle Scholar
- Eker S, Krummenacker M, Shearer A, Tiwari A, Keseler I, Talcott C, Karp P. Computing minimal nutrient sets from metabolic networks via linear constraint solving. BMC Bioinform. 2013;14(1):114. doi:https://doi.org/10.1186/1471-2105-14-114.View ArticleGoogle Scholar
- Watson MR. Metabolic maps for the apple 2. Biochem Soc Trans. 1984;12(6):1093–4. doi:https://doi.org/10.1042/bst0121093.View ArticleGoogle Scholar
- Lee S, Phalakornkule C, Domach MM, Grossmann IE. Recursive MILP model for finding all the alternate optima in LP models for metabolic networks. Comput Chem Eng. 2000;24(2–7):711–6. doi:https://doi.org/10.1016/S0098-1354(00)00323-9.View ArticleGoogle Scholar
- de Figueiredo LF, Podhorski A, Rubio A, Kaleta C, Beasley JE, Schuster S, Planes FJ. Computing the shortest elementary flux modes in genome-scale metabolic networks. Bioinformatics. 2009;25(23):3158–65. doi:https://doi.org/10.1093/bioinformatics/btp564.View ArticlePubMedGoogle Scholar
- von Kamp A, Klamt S. Enumeration of smallest intervention strategies in genome-scale metabolic networks. PLoS Comput Biol. 2014;10(1):1003378. doi:https://doi.org/10.1371/journal.pcbi.1003378.View ArticleGoogle Scholar
- Gallo G, Longo G, Nguyen S, Pallottino S. Directed hypergraphs and applications. Discrete Appl Math. 1993;42(2–3):177–201.View ArticleGoogle Scholar
- Ausiello G, Franciosa PG, Frigioni D. Directed hypergraphs: problems, algorithmic results, and a novel decremental approach. In: ICTCS. 2001. p. 312–27.Google Scholar
- Gurvich V, Khachiyan L. On generating the irredundant conjunctive and disjunctive normal forms of monotone boolean functions. Discrete Appl Math. 1999;96:363–73.View ArticleGoogle Scholar
- Monk J, Charusanti P, Aziz RK, Lerman JA, Premyodhin N, Orth JD, Feist AM, Palsson BØ. Genome-scale metabolic reconstructions of multiple Escherichia coli strains highlight strain-specific adaptations to nutritional environments. Proc Natl Acad Sci. 2013;110(50):20338–43. doi:https://doi.org/10.1073/pnas.1307797110.View ArticlePubMedPubMed CentralGoogle Scholar
- Meador J, Caldwell ME, Cohen PS, Conway T. Escherichia coli pathotypes occupy distinct niches in the mouse intestine. Infect Immun. 2014;82(5):1931–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Fabich A, Jones SA, Chowdhury FZ, Cernosek A, Anderson A, Smalley D, McHargue JW, Hightower GA, Smith JT, Autieri SM, Leatham MP, Lins JJ, Allen RL, Laux DC, Cohen PS, Conway T. Comparison of carbon nutrition for pathogenic and commensal Escherichia coli strains in the mouse intestine. Infect Immun. 2008;76:1143–52. doi:https://doi.org/10.1128/iai.01386-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Leatham M, Banerjee S, Autieri SM, Mercado-Lubo R, Conway T, Cohen PS. Precolonized human commensal Escherichia coli strains serve as a barrier to E. coli O157:H7 growth in the streptomycin-treated mouse intestine. Infect Immun. 2009;77(7):2876–86. doi:https://doi.org/10.1128/iai.00059-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Maltby R, Leatham-Jensen MP, Gibson T, Cohen PS, Conway T. Nutritional basis for colonization resistance by human commensal Escherichia coli strains HS and Nissle 1917 against E. coli O157:H7 in the mouse intestine. PLoS One. 2013;8(1):53957. doi:https://doi.org/10.1371/journal.pone.0053957.View ArticleGoogle Scholar