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.

Robust generation of elementary flux modes

Application of dynamic metabolic flux analysis for process modeling: Robust flux estimation with regularization, confidence bounds, and selection of elementary modes

Robustness analysis of the Escherichia coli metabolic network

Fluxomers: a new approach for 13C metabolic flux analysis

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

An Improved Multi-parametric Programming Algorithm for Flux Balance Analysis of Metabolic Networks

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

Validation-based model selection for 13C metabolic flux analysis with uncertain measurement errors

Minimizing Enzyme Mass to Decompose Flux Distribution for Identifying Biologically Relevant Elementary Flux Modes.

Flux cost functions and optimal metabolic states

Thermodynamics-based metabolic flux analysis

Expanded Flux Variability Analysis on Metabolic Network of Escherichia Coli

Robust Mutant Strain Design by Pessimistic Optimization

Constraint-based inverse modeling of metabolic networks: a proof of concept

The solution space of metabolic networks: producibility, robustness and fluctuations

Model Validation and Selection in Metabolic Flux Analysis and Flux Balance Analysis

A principal components method constrained by elementary flux modes: analysis of flux data sets

Estimating the size of the solution space of metabolic networks

Conservation of high-flux backbone in alternate optimal and near-optimal flux distributions of metabolic networks

E-Flux2 and SPOT: Validated Methods for Inferring Intracellular Metabolic Flux Distributions from Transcriptomic Data

A Novel Method For Flux Distribution Computation In Metabolic Networks