Abstract:Stable isotope-assisted metabolic flux analysis (MFA) is a powerful method to estimate carbon flow and partitioning in metabolic networks. At its core, MFA is a parameter estimation problem wherein the fluxes and metabolite pool sizes are model parameters that are estimated, via optimization, to account for measurements of steady-state or isotopically-nonstationary isotope labeling patterns. As MFA problems advance in scale, they require efficient computational methods for fast and robust convergence. The structure of the MFA problem enables it to be cast as an equality-constrained nonlinear program (NLP), where the equality constraints are constructed from the MFA model equations, and the objective function is defined as the sum of squared residuals (SSR) between the model predictions and a set of labeling measurements. This NLP can be solved by using an algebraic modeling language (AML) that offers state-of-the-art optimization solvers for robust parameter estimation and superior scalability to large networks. When implemented in this manner, the optimization is performed with no distinction between state variables and model parameters. During each iteration of such an optimization, the system state is updated instead of being calculated explicitly from scratch, and this occurs concurrently with improvement in the model parameter estimates. This optimization approach starkly contrasts with traditional "shooting" methods where the state variables and model parameters are kept distinct and the system state is computed afresh during each iteration of a stepwise optimization. Our NLP formulation uses the MFA modeling framework of Wiechert et al. [1], which is amenable to incorporation of the model equations into an NLP. The NLP constraints consist of balances on either elementary metabolite units (EMUs) or cumomers. In this formulation, both the steady-state and isotopically-nonstationary MFA (inst-MFA) problems may be solved as an NLP. For the inst-MFA case, the ordinary differential equation (ODE) system describing the labeling dynamics is transcribed into a system of algebraic constraints for the NLP using collocation. This large-scale NLP may be solved efficiently using an NLP solver implemented on an AML. In our implementation, we used the reduced gradient solver CONOPT, implemented in the General Algebraic Modeling System (GAMS). The NLP framework is particularly advantageous for inst-MFA, scaling well to large networks with many free parameters, and having more robust convergence properties compared to the shooting methods that compute the system state and sensitivities at each iteration. Additionally, this NLP approach supports the use of tandem-MS data for both steady-state and inst-MFA when the cumomer framework is used. We assembled a software, eiFlux, written in Python and GAMS that uses the NLP approach and supports both steady-state and inst-MFA. We demonstrate the effectiveness of the NLP formulation on several examples, including a genome-scale inst-MFA model, to highlight the scalability and robustness of this approach. In addition to typical inst-MFA applications, we expect that this framework and our associated software, eiFlux, will be particularly useful for applying inst-MFA to complex MFA models, such as those developed for eukaryotes (e.g. algae) and co-cultures with multiple cell types.

MOFA: Multi-Objective Flux Analysis for the COBRA Toolbox

Creation and analysis of biochemical constraint-based models: the COBRA Toolbox v3.0

MOF: A Modular Framework for Rapid Application of Optimization Methodologies to General Engineering Design Problems

PyCoMo: a python package for community metabolic model creation and analysis

A robust and efficient method for estimating enzyme complex abundance and metabolic flux from expression data

MOFA+: a statistical framework for comprehensive integration of multi-modal single-cell data

A multiobjective evolutionary algorithm toolbox for computer-aided multiobjective optimization

Nonlinear Multi-Objective Flux Balance Analysis of the Warburg Effect

MOFA+: a probabilistic framework for comprehensive integration of structured single-cell data

Multi-objective resistance-capacitance optimization algorithm: An effective multi-objective algorithm for engineering design problems

A Structured Evaluation of Genome-Scale Constraint-Based Modeling Tools for Microbial Consortia

An efficient chemical reaction algorithm for multi-objective combinatorial bi-level optimization

Bayesian metabolic flux analysis reveals intracellular flux couplings

Elimination of thermodynamically infeasible loops in steady-state metabolic models

MORF: An online tool for exploring microbial cell responses using multi-omics analysis

COVAIN: a toolbox for uni- and multivariate statistics, time-series and correlation network analysis and inverse estimation of the differential Jacobian from metabolomics covariance data

Isotope-assisted metabolic flux analysis as an equality-constrained nonlinear program for improved scalability and robustness

Multi-objective liver cancer algorithm: A novel algorithm for solving engineering design problems

The Metano Modeling Toolbox MMTB: An Intuitive, Web-Based Toolbox Introduced by Two Use Cases

MOANA: Multi-Objective Ant Nesting Algorithm for Optimization Problems