Robustness analysis of elementary flux modes generated by column generation

•Considers a robust identification of elementary flux modes and a flux over each mode.•Measurement errors are taken into consideration.•Show that a robust optimization problem can be formulated and solved.•Show how the problem can be solved, even when the elementary flux modes are unknown.

Elementary flux modes (EFMs) are vectors defined from a metabolic reaction network, giving the connections between substrates and products. EFMs-based metabolic flux analysis (MFA) estimates the flux over each EFM from external flux measurements through least-squares data fitting. The measurements used in the data fitting are subject to errors. A robust optimization problem includes information on errors and gives a way to examine the sensitivity of the solution of the EFMs-based MFA to these errors. In general, formulating a robust optimization problem may make the problem significantly harder. We show that in the case of the EFMs-based MFA, when the errors are only in measurements and bounded by an interval, the robust problem can be stated as a convex quadratic programming (QP) problem. We have previously shown how the data fitting problem may be solved in a column-generation framework. In this paper, we show how column generation may be applied also to the robust problem, thereby avoiding explicit enumeration of EFMs. Furthermore, the option to indicate intervals on metabolites that are not measured is introduced in this column generation framework. The robustness of the data is evaluated in a case-study, which indicates that the solutions of our non-robust problems are in fact near-optimal also when robustness is considered, implying that the errors in measurement do not have a large impact on the optimal solution. Furthermore, we showed that the addition of intervals on unmeasured metabolites resulted in a change in the optimal solution.