Fluxomers: a new approach for 13C metabolic flux analysis.

BACKGROUND: The ability to perform quantitative studies using isotope tracers and metabolic flux analysis (MFA) is critical for detecting pathway bottlenecks and elucidating network regulation in biological systems, especially those that have been engineered to alter their native metabolic capacities. Mathematically, MFA models are traditionally formulated using separate state variables for reaction fluxes and isotopomer abundances. Analysis of isotope labeling experiments using this set of variables results in a non-convex optimization problem that suffers from both implementation complexity and convergence problems. RESULTS: This article addresses the mathematical and computational formulation of (13)C MFA models using a new set of variables referred to as fluxomers. These composite variables combine both fluxes and isotopomer abundances, which results in a simply-posed formulation and an improved error model that is insensitive to isotopomer measurement normalization. A powerful fluxomer iterative algorithm (FIA) is developed and applied to solve the MFA optimization problem. For moderate-sized networks, the algorithm is shown to outperform the commonly used 13CFLUX cumomer-based algorithm and the more recently introduced OpenFLUX software that relies upon an elementary metabolite unit (EMU) network decomposition, both in terms of convergence time and output variability. CONCLUSIONS: Substantial improvements in convergence time and statistical quality of results can be achieved by applying fluxomer variables and the FIA algorithm to compute best-fit solutions to MFA models. We expect that the fluxomer formulation will provide a more suitable basis for future algorithms that analyze very large scale networks and design optimal isotope labeling experiments.

Hartmann-Shack analysis errors.

A theoretical model for reconstruction errors in the Hartmann-Shack wave-front sensor was constructed. The measured pattern of a regulargrid of spots can be analyzed by their individual centroiding or by a globalFourier demodulation. We investigate pixelization errors in the cameraand Poisson errors in the camera pixels. We show that by the Fourier demodulation technique it is possible to overcome pixelization errors which occur in the traditional centroid technique. By supporting simulations we show that the two methods coincide for an infinite number of pixels per spot.