From: Exploiting the pathway structure of metabolism to reveal high-order epistasis
MCS = NetKO(S, J, k, R, PFL) |
---|
1: // Arguments: |
2: // S, stoichiometry matrix |
3: // J ⊂ N, a set of objective reactions |
4: // k, cardinality limit for MCS |
5: // R, reactions from which MCS should be chosen |
6: // PFL, pathway fragment limit for tableau algorithm |
7: // Output: |
8: // MCS, a collection of MCS for J |
9: E = tableau(S, PFL) |
10: = {} |
11: for all j ∈ J do |
12: Ej= subset of j-containing pathway fragments in E |
13: = MinHit(Ej, k, R); |
14: = ∪ |
15: end for |
16: remove sets H ∈ from for which there exists a subset in |
17: MCS = Reduce(, J, S) |
18: return MCS |