Flux balance analysis with objective function defined by proteomics data-metabolism of mycobacterium tuberculosis exposed to mefloquine

Abstract

We present a study of the metabolism of the Mycobacterium tuberculosis after exposure to antibiotics using proteomics data and flux balance analysis (FBA). The use of FBA to study prokaryotic organisms is well-established and allows insights into the metabolic pathways chosen by the organisms under different environmental conditions. To apply FBA a specific objective function must be selected that represents the metabolic goal of the organism. FBA estimates the metabolism of the cell by linear programming constrained by the stoichiometry of the reactions in an in silico metabolic model of the organism. It is assumed that the metabolism of the organism works towards the specified objective function. A common objective is the maximization of biomass. However, this goal is not suitable for situations when the bacterium is exposed to antibiotics, as the goal of organisms in these cases is survival and not necessarily optimal growth. In this paper we propose a new approach for defining the FBA objective function in studies when the bacterium is under stress. The function is defined based on protein expression data. The proposed methodology is applied to the case when the bacterium is exposed to the drug mefloquine, but can be easily extended to other organisms, conditions or drugs. We compare our method with an alternative method that uses experimental data for adjusting flux constraints. We perform comparisons in terms of essential enzymes and agreement using enzyme abundances. Results indicate that using proteomics data to define FBA objective functions yields less essential reactions with zero flux and lower error rates in prediction accuracy. With flux variability analysis we observe that overall variability due to alternate optima is reduced with the incorporation of proteomics data. We believe that incorporating proteomics data in the objective function used in FBA may help obtain metabolic flux representations that better support experimentally observed features.