The proposed methodology utilizes transformation of the bounds of the embedded linear programming problem of flux balance analysis via a logarithmic barrier (interior point) approach. By exploiting the first-order optimality conditions of the interior-point problem, and with further transformations, the approach results in a system of implicit ordinary differential equations. Results from four case studies, show that the CPU and wall-times obtained using the proposed method are competitive with existing state-of-the art approaches for solving dFBA simulations, for problem sizes up to genome-scale. The differentiability of the proposed approach allows, using existing commercial packages, its application to the optimal control of dFBA problems at a genome-scale size, thus outperforming existing formulations as shown by two dynamic optimization case studies.
- Dynamic flux balance analysis
- Genome-scale metabolic network
- Linear programming
- Ordinary differential equations with embedded optimization