 Software
 Open Access
PFA toolbox: a MATLAB tool for Metabolic Flux Analysis
 Yeimy Morales^{1}Email authorView ORCID ID profile,
 Gabriel Bosque^{2},
 Josep Vehí^{1},
 Jesús Picó^{2} and
 Francisco Llaneras^{1}
https://doi.org/10.1186/s1291801602841
© The Author(s). 2016
Received: 24 October 2015
Accepted: 1 June 2016
Published: 11 July 2016
Abstract
Background
Metabolic Flux Analysis (MFA) is a methodology that has been successfully applied to estimate metabolic fluxes in living cells. However, traditional frameworks based on this approach have some limitations, particularly when measurements are scarce and imprecise. This is very common in industrial environments. The PFA Toolbox can be used to face those scenarios.
Results
Here we present the PFA (Possibilistic Flux Analysis) Toolbox for MATLAB, which simplifies the use of Interval and Possibilistic Metabolic Flux Analysis. The main features of the PFA Toolbox are the following: (a) It provides reliable MFA estimations in scenarios where only a few fluxes can be measured or those available are imprecise. (b) It provides tools to easily plot the results as interval estimates or flux distributions. (c) It is composed of simple functions that MATLAB users can apply in flexible ways. (d) It includes a Graphical User Interface (GUI), which provides a visual representation of the measurements and their uncertainty. (e) It can use stoichiometric models in COBRA format. In addition, the PFA Toolbox includes a User’s Guide with a thorough description of its functions and several examples.
Conclusions
The PFA Toolbox for MATLAB is a freely available Toolbox that is able to perform Interval and Possibilistic MFA estimations.
Keywords
Background
The problem of estimating unknown metabolic fluxes in living cells has been tackled using several methodologies. MFA is one of the most extensively and successfully applied approaches to estimating fluxes [1]. Usually MFA refers to 13CMFA which uses stable isotopically labeled substrates (e.g., 13Clabeled glucose) combined with stoichiometric balancing to estimate the metabolic fluxes in steady state systems [2, 3]. However, in this study we refer to non13CMFA methods. These methods mainly rely on measurements of external fluxes (uptake and production rates) to estimate the flux state of cells. Traditional MFA methods present some limitations when accounting for irreversible reactions [4], underdetermined problems [5], and lack of measurements [6]. To reduce these limitations we have developed Interval [7] and Possibilistic [8] MFA methods, which are wellsuited methodologies for scenarios with limited available data. Their main benefits are the following [6–10]: (a) They can consider the irreversibility of the reactions and other inequality constraints. (b) They are able to represent the measured fluxes as intervals and even distributions to describe the uncertainty of the system. (c) They provide interval estimates, which are more reliable and more informative than pointwise solutions, particularly when multiple flux values are possible. (d) They are able to perform estimations in scenarios of high uncertainty or lack of measurements, being those estimates as reliable as possible. In addition, (e) Possibilistic MFA allows the detection and handling of inconsistencies between a model and a set of measurements. The PFA Toolbox provides all these features while preserving computational efficiency.
In the last years, several published works have used these methodologies to perform interval estimations of metabolic fluxes [9, 11–18] and consistency analysis with Possibilistic MFA [9, 17, 18]. Interval MFA was also implemented in FASIMU [16]. However, any intermediate user of MATLAB, Mathematica, R, etc. can easily implement Interval MFA. The easily implementation of Interval MFA has led to be used more often than Possibilistic MFA, which requires more mathematical development and additional linear optimizations. The PFA Toolbox presented here simplifies the use of both methods.
The PFA Toolbox provides a comprehensive set of MATLAB functions to easily and quickly apply Interval and Possibilistic MFA. The PFA Toolbox is completely free and open source; users are welcome to modify and adapt the toolbox code to build their own particular functions to fulfill specific requirements under the mild conditions described in the accompanying license. In the following subsections, we briefly describe the methods implemented in the toolbox: Interval MFA and Possibilistic MFA. A detailed description of both methods can be found in [6].
Interval MFA
where, considering a system with n metabolites and r reactions, N ∈ R ^{{nxr}} and D ∈ R ^{{rxr}} is a diagonal matrix with D _{ ii } = 1 if the flux is reversible (0 otherwise), and v ∈ R ^{{r}} is the vector of metabolic fluxes. The values of v that are solution of (1) define a flux distribution.
where v _{ m } ^{ m } and v _{ m } ^{ M } are vectors with the minimum and maximum possible values that the measured fluxes v _{ m } can take due to measurement’s uncertainty.
This procedure provides an interval estimate for any flux of interest. These interval estimates are particularly useful in the two situations of having imprecise measurements and/or when few measures are available. Extra details about Interval MFA can be found in [6, 7, 10].
Possibilistic MFA
Possibilistic MFA may be seen as a more flexible and powerful extension of Interval MFA. The methodology is based on two ideas: (a) Representing knowledge with constraints satisfied to a certain degree, thus transforming the feasibility of a potential solution into a gradual notion of “possibility” that accounts for uncertainty, and (b) using computationally efficient optimizationbased methods, such as Linear Programming, to query for the “most possible” solutions. This methodology is able to face two different problems: (a) To evaluate the consistency between a model and a set of measurements, and (b) to obtain rich estimates of metabolic fluxes. Instead of pointwise estimates, it computes interval estimations for a desired degree of possibility and for entire possibility distributions.
Possibilistic MFA starts with a set of modelbased constraints (MOC) defined in (1).
where v _{ m } is the vector of the actual values of the measured fluxes, and w _{ m } is the vector of the measured values for them. Both differ due to errors and imprecisions. This uncertainty is represented by the slack variables ε _{ 1 } , μ _{ 1 } , ε _{ 2 } and μ _{ 2 }. The bounds ε _{ 2 } and μ _{ 2 } define a band of fully possible values for v _{ m } around the measured values w _{ m }. The components ε _{ 1 } and μ _{ 1 } are penalized in a cost index (5) to assign a decreasing possibility to larger errors. Each candidate solution of (1) and (4) can be denoted as δ = {v, w _{ m } , ε _{ 1 } , μ _{ 1 } , ε _{ 2 } , μ _{ 2 }}.
Where α and β are row vectors of accuracy coefficients or weights that define each measurement’s a priori accuracy. These weights need to be defined by the user, e.g., if sensor error is «symmetric», α and β should be defined to be equal.
From this point, Possibilistic MFA calculates different estimates by solving LP problems. You can compute the set of flux values with maximum possibility (a pointwise estimation) or a more informative estimation with intervals or flux distributions.
Pointwise estimations
The solution flux vector v, that we call v _{ mp, } contains the most possible values that are consistent with both the model and the measurements.
This pointwise estimation may be unreliable when multiple solutions are reasonably possible. In these instances, distributions and interval estimates can be computed instead.
Interval estimates
The upper bound is defined by replacing minimum for maximum.
Distributions as estimates
The complete possibility distribution of a flux can also be obtained for marginal and conditional possibilities. Marginal possibilities provide the degree of possibility of each value for a given flux. Conditional distributions are equivalent to normalizing the marginal possibility distribution to a maximum equal to one.
Possibilistic MFA was casted as a linear optimization problem, for which widely known and efficient tools exist. This great computational performance makes the methodology suitable —in principle— for largescale metabolic networks.
More information about the methods and a deeper discussion about the strengths and limitations of each approach can be found in our previous works [6–8, 10] and in the toolbox User’s Guide (http://kikollan.github.io/PFAToolbox/).
Implementation
The PFA Toolbox has been developed to run in MATLAB. Its core is a set of MATLAB functions that solve each step in a typical MFA problem. The code for all functions is provided with the toolbox. The PFA Toolbox also includes a Graphical User Interface (GUI) to represent the measurements in possibilistic terms. The GUI runs within MATLAB.
The toolbox requires solving LP problems, and those are solved with a flexible and efficient external optimizer, YALMIP [21]. We provide a copy of YALMIP within the PFA Toolbox, but further information about it can be found at the YALMIP website [22]. YALMIP can use different LP solvers, and so does the PFA Toolbox. Three LP solvers were tested: IBM ILOG CPLEX by IBM [23], GLPK [24], and Linprog, the LP solver included in MATLAB. However, we do not recommend the use of Linprog, which proved unreliable, especially for larger MFA problems. Instead, CPLEX or GLPK showed excellent performance. CPLEX has a 90day free evaluation version, and can be used free for research and academic purposes. GLPK is freely available.
Results and discussion
List of functions in the PFA Toolbox
Initialization  

initPFAtoolbox  It starts the PFA Toolbox 
1: MFA problem formulation  
define_MOC  It defines the modelbased constraints 
define_PossMeasurements  It represents the measured fluxes 
define_MEC  It defines the measuredbased constraints 
2: Computing estimations  
solve_maxPoss  It calculates the most possible set of flux values 
solve_maxPossIntervals  It calculates the interval of most possible flux values 
solve_PossInterval  It calculates the interval of flux values with the desired possibility 
3: Plotting the estimations  
plot_PossMeasurements  It plots measurements in possibilistic terms 
plot_distribution  It plots the distribution of a given flux 
plot_intervals  It plots interval estimates of a given flux 
4: Other  
Solve_possintervalYMP  Advanced function; read its help. 
solve_Interval  It solves an Interval MFA problem 

» It gives reliable MFA estimations even in uncertain or underdetermined scenarios (those where only a few fluxes can be measured).

» It provides MFA estimations accounting for measurement’s imprecision.

» It provides functions to plot interval estimates and distributions.

» It is composed of simple, free and open functions.
In addition, the toolbox is developed to use stoichiometric models with the format of the COBRA Toolbox (ConstraintBased Reconstruction and Analysis). This format is widely used due to the popularity of COBRA. As an alternative, the user can simply define a model by providing a stoichiometric matrix.
The main features of the toolbox are shown in the next three examples. Additional examples and a thorough description of all functionalities of the toolbox are provided in the User’s Guide. The details about the mathematical methods implemented in the toolbox can be found in [7, 8, 10], and in [6].
Example of flux estimation under data scarcity
The MFA problem consists in the estimation of all six fluxes. Notice, however, that traditional MFA cannot be performed because the problem is undetermined: any pointwise estimate will be only a particular solution of a group of possible ones [5]. The methods in the PFA Toolbox tackle this situation and provide reliable and informative estimates.
In this case, we choose to apply Possibilistic MFA to estimate the fluxes. The first step to solve the problem is to define the modelbased constraints (MOC). Stoichiometric model can be directly defined in the code or be provided in COBRA format.
The next step is the addition of measurements and their uncertainties (in this example, we assume that the measurement w_{4} is very accurate, but w_{6} is not. In agreement with the problem formulation, we assign values to the slack variables μ_{2} and ε_{2}, and the weights α and β (details about this process can be found in the User’s Guide).
Once the MOC and MEC constraints have been defined, the third step is to obtain the estimates. Possibilistic MFA methodology calculates three types of estimations. In this case, we compute three interval estimates for each flux, for conditional possibilities of 0.5, 0.8 and 1.
This same procedure can be applied to obtain other types of estimates, such as the complete possibility distribution for a flux. Those computations can be performed using the function solve_PossInterval. The obtained distributions are for conditional possibilities (see [8] for a detailed explanation of the notion of conditional possibility). These possibilistic distributions can be plotted with the fuction plot_distribution. As an example, Fig. 3b shows the distribution estimation for all the six fluxes. The results show, for instance, that the most possible value for v_{1} is 2.75 mmol/h (π = 1), that v_{1} being equal to 6.1 mmol/h is a less possible situation (π = 0.6), and that a v_{1} being larger than 18 mmol/h is very unlikely (π < 0.1).
The model and the code for all the computations are provided as (Additional file 1a).
Note: to apply Interval MFA a similar protocol can be followed. The main difference is that the measures will be represented as intervals instead of being represented in possibilistic terms.
Example of flux estimation: biomass growth of Pichia pastoris
In this example, we estimate the growth of several chemostat cultures of P. pastoris. For each chemostat only a few extracellular fluxes are measured (mainly substrates uptakes and secretion rates) and the aim is to estimate the cellular growth.
The constraintbased model for P. pastoris used is presented in [18] (see Additional file 2). It is a relatively small representation including only the main catabolic pathways considering the uptake of the usual carbon sources: methanol, glucose and glycerol. The stoichiometric model contains 37 metabolites and 48 reactions, with reversibility accounted for. The stoichiometric matrix and all the measurements can be found in the (Additional file 3) [3135].
The estimations show good agreement with the experimental growth rates (as expected, since this model and the data have been tested previously). Notice that the interval estimates not only predict the growth rates but also provide an indication of the estimation reliability. The complete code for all computations can be found in the (Additional file 1b).
Example of flux estimation: growth of Escherichia coli
Here we use a wellknown model of E. coli, taken from [25] and illustrated in the (Additional file 4). It is a relatively compact model containing 72 metabolites and 95 reactions. We consider six chemostat experiments of E. coli growing in glucose [26]. The datasets contain information only for a handful of extracellular measurements (growth rate, substrate uptake, oxygen uptake, CO_{2} production and acetate and pyruvate secretion). The model and the measurements can be found in the (Additional file 5).
Possibilistic MFA is applied again to estimate the growth rate for all six scenarios. The problem is similar to the previous one, and we assume the same uncertainty for each measurement. However, we now consider a larger model for a different and widely used organism. The computation procedure is analogous to the one previously described. The complete code for all computations can be found in the (Additional file 1c).
The flux estimates computed with the toolbox are compatible with the actual growth rate in all scenarios (Fig. 4b). Notice, however, that the estimates are wider than in the first example (nogrowth is possible in all of them, but the maximum possible growth is near the actual one). The model is larger and the available measurements are not enough to determine completely the flux state of cells. This illustrates one limitation of Interval and Possibilistic MFA: the estimates are only as precise as the uncertainty and the available measurements allow.
Example of consistency analysis with P. pastoris
Experimental data for six chemostat experiments with Pichia pastoris and an analysis of its consistency against a model
Reference  μ  Q_{Glu}  Q_{Glyc}  Q_{Met}  Q_{Pyr}  Q_{Cit}  Q_{EtOH}  OUR  CPR  π_{mp} ^{b} 

Cmmol^{a}/g/h  mmol/g/h  mmol/g/h  mmol/g/h  mmol/g/h  mmol/g/h  mmol/g/h  mmol/g/h  mmol/g/h  
[31]  6.17  0.00  2.75  0.00  0.00  0.00  0.00  3.62  2.35  0.16 
[32]  3.27  0.81  0.00  1.09  0.00  0.00  0.00  4.02  2.68  1.00 
[33]  2.38  0.00  1.21  0.00  N.A.  N.A.  N.A.  1.65  1.22  1.00 
[33]  4.89  0.00  2.40  0.00  N.A.  N.A  N.A.  3.12  2.29  1.00 
[34]  1.40  0.00  0.00  2.55  N.A.  N.A.  N.A.  2.16  1.15  1.00 
[34]  0.94  0.00  0.00  1.87  N.A.  N.A.  N.A.  1.67  0.93  1.00 
We start as in previous examples by defining MOC and MEC constraints. The next step is to compute the estimation. In this example, we compute the most possible solution for each experiment with the solve_maxPoss function. This provides the maximum possibility flux vector and the associated degree of possibility (π_{mp}) between [0, 1] of the most possible solution. This value provides an indication of the agreement between the modelbased constraints (MOC) and the measurements constraints (MEC).
A possibility equal to one is interpreted as a complete consistency; a lower value implies that there are errors in one (or more) of the measurements or in the model. The complete MATLAB code for this computation can be found in (Additional file 1b).
The results presented in Table 2 show that all datasets except one are highly consistent with the model. The dataset 1 has a low degree of possibility (lower 0.2). This suggests that one or more of the measured fluxes in that experiment is unreliable and may contain errors.
All the computations of these four examples were performed with the PFA Toolbox. The computations take approximately 13 s in a 64bit Windows PC (Intel Core™ i5 2.5 GHz processor), using MATLAB R2012a with IBM ILOG CPLEX Optimizer as the solver for Linear Programming problems.
Notes on computational efficiency and large networks
The methods used by the PFA Toolbox, Possibilistic MFA and interval MFA, have been cast as linear optimization problems, and thus they can be solved with computational efficiency. This makes these methodologies suitable for largescale metabolic networks. For instance, when tested on a genomescale E. coli model (iJO1366) that contains 2583 reactions [27], the PFA Toolbox is able to get estimates for all 2507 fluxes with three degrees of possibility (i.e., solving 3x2507 LP problems). Computing those estimates required 120 min in an AMD A10–5800 K with Radeon HD graphic (3.80 GHz) PC and 8 GB of RAM with GLPK optimizer. This suggests that the PFA Toolbox may be able to solve MFA flux estimations of large models with good results and reasonable computational cost.
There is, however, a limitation regarding MFAwise methods when estimating fluxes in large networks: there may be too many flux vectors compatible with the (few) available measurements [28]. Unlike traditional methods, those proposed here may still be of use in this situation. Possibilistic MFA and Interval MFA capture all the equally possible flux states (or “similarly” possible) by means of possibilistic distributions or intervals. If there is a wide range of candidates, however, the estimation may be only slightly informative. If this is the case, one could decide to incorporate a rational assumption, as done in FBA methods [29, 30].
Conclusions
We have presented the PFA Toolbox for MATLAB. This toolbox provides a set of MATLAB functions to apply Interval MFA and Possibilistic MFA in a simple and flexible way. The PFA Toolbox is completely free and open source, and can be modified by its users. The toolbox implements MFAwise methods to perform metabolic flux estimations that are particularly well suited to deal with scenarios of high uncertainty and scarce measurements, which are common in industry.
Availability and requirements
Project name: PFA Toolbox version 1.0.0.
Project home page: http://kikollan.github.io/PFAToolbox/
Operating systems: platform independent.
Programming language: MATLAB
Other requirements: −
License: Own license.
Any restriction to use by nonacademics: none.
Abbreviations
CB, constraintbased model; COBRA, ConstraintBased Reconstruction and Analysis; FASIMU, Fluxbalance Analysis based Simulations; GLPK, GNU Linear Programming kit; GUI, Graphical User Interface; IBM ILOG CPLEX, Highperformance mathematical programming solver for linear programming; LP, Linear Programming; MEC, Measurement constraints; MFA, Metabolic Flux Analysis; MOC, modelbased constraints; PFA, Possibilistic Flux Analysis; YALMIP, Modelling language for advanced modeling and solution of optimization problems
Declarations
Acknowledgements
We acknowledge Ignacio Ribelles for contributing to programming the MATLAB functions and for writing the GUI.
Funding
This research has been partially supported by the Spanish Government (FEDERCICYT: DPI 2014–55276C5–1R). Yeimy Morales is grateful for the BR Grants of the University of Girona (BR2012/26). Gabriel Bosque Chacón is recipient of a doctoral fellowship from the Spanish Government (BES2012–053772).
Availability of data and materials
All data are included in the manuscript, and the associated supplementary material and links provided.
Authors’ contributions
FLL and JP developed the idea for the toolbox, and with JV, they designed the research and coordinated the project. FLL designed the toolbox implementation and wrote the first version of the code. YM contributed to the code, documented it and wrote the user’s documentation. YM and GB developed the examples and debugged the toolbox. YM drafted the first manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Not applicable.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
 Sauer U, Hatzimanikatis V, Bailey J, Hochuli M, Szyperski T, Wuethrich K. Metabolic fluxes in riboflavinproducing Bacillus subtilis. Nature biotechnology. 1997;15(5):448–52.View ArticlePubMedGoogle Scholar
 Wittmann C. Metabolic flux analysis using mass spectrometry. In: Tools and Applications of Biochemical Engineering Science. Berlin: Springer; 2002. p. 39–64.View ArticleGoogle Scholar
 Antoniewicz M. Methods and advances in metabolic flux analysis: a minireview. J Ind Microbiol Biot. 2015;42(3):317–25.View ArticleGoogle Scholar
 AraúzoBravo MR, Shimizu JK. An improved method for statistical analysis of metabolic flux analysis using isotopomermapping matrices with analytical expressions. J Biotech. 2003;05:117–33.View ArticleGoogle Scholar
 Klamt S, Schuster S, Gilles D. Calculability analysis in underdetermined metabolic networks illustrated by a model of the central metabolism in purple nonsulfur bacteria. Biotechnol Bioeng. 2002;77(7):734–51.View ArticlePubMedGoogle Scholar
 Llaneras F. Interval and possibilistic methods for constraintbased metabolic models, PhD Thesis. Universidad Politécnica de Valencia: Departamento de Ingeniería de Sistemas y Automática; 2011.View ArticleGoogle Scholar
 Llaneras F, Picó J. An interval approach for dealing with flux distributions and elementary modes activity patterns. J Theor Biol. 2007;246(2):290–308.View ArticlePubMedGoogle Scholar
 Llaneras F, Sala A, Picó J. A possibilistic framework for constraintbased metabolic flux analysis. BMC Syst Biol. 2009;3(1):79.View ArticlePubMedPubMed CentralGoogle Scholar
 Tortajada M, Llaneras F, Picó J. Validation of a constraintbased model of Pichia pastoris metabolism under data scarcity. BMC Syst Biol. 2010;4(1):115.View ArticlePubMedPubMed CentralGoogle Scholar
 Llaneras F, Picó J. A procedure for the estimation over time of metabolic fluxes in scenarios where measurements are uncertain and/or insufficient. BMC Bioinformatics. 2007;8(1):421.View ArticlePubMedPubMed CentralGoogle Scholar
 Iyer VV, Ovacik MA, Androulakis IP, Roth CM, Ierapetritou MG. Transcriptional and metabolic flux profiling of triadimefon effects on cultured hepatocytes. Toxicology and applied pharmacology. 2010;248(3):165–77.View ArticlePubMedGoogle Scholar
 Zamorano F, Wouwer A, Bastin G. Detailed metabolic flux analysis of an underdetermined network of CHO cells. J Biotechnol. 2010;150(4):497–508.View ArticlePubMedGoogle Scholar
 Iyer V, Yang H, Ierapetritou M, Roth C. Effects of glucose and insulin on HepG2‐C3A cell metabolism. Biotechnol Bioeng. 2010;107(2):347–56.View ArticlePubMedGoogle Scholar
 Iyer V, Androulakis I, Roth C, Ierapetritou M. Effects of Triadimefon on the Metabolism of Cultured Hepatocytes. In: BioInformatics and BioEngineering (BIBE), IEEE International Conference on. 2010. p. 118–23.Google Scholar
 Orman MA, Arai K, Yarmush ML, Androulakis IP, Berthiaume F, Ierapetritou MG. Metabolic flux determination in perfused livers by mass balance analysis: effect of fasting. Biotechnology and bioengineering. 2010;107(5):825–35.View ArticlePubMedGoogle Scholar
 Hoppe A, Hoffmann S, Gerasch A, Gille C, Holzhütter H. FASIMU: flexible software for fluxbalance computation series in large metabolic networks. BMC bioinformatics. 2011;12(1):28.View ArticlePubMedPubMed CentralGoogle Scholar
 González J, FolchFortuny A, Llaneras F, Tortajada M, Picó J, Ferrer A. Metabolic flux understanding of Pichia pastoris grown on heterogenous culture media. Chemometr Intell Lab. 2014;134:89–99.View ArticleGoogle Scholar
 Morales Y, Tortajada M, Picó J, Vehí J, Llaneras F. Validation of an FBA model for Pichia pastoris in chemostat cultures. BMC System Biol. 2014;8(1):142.View ArticleGoogle Scholar
 Stephanopoulos GN, Aristidou AA, Nielsen J. Metabolic Engineering: Principles and Methodologies. San Diego, USA: Academic; 1998.Google Scholar
 Heijden R, Romein B, Heijnen J, Hellinga C, Luyben K. Linear constraint relations in biochemical reaction systems: I & II. Biotech Bioeng. 1994;43(1):3–10.View ArticleGoogle Scholar
 Lofberg J. YALMIP: A toolbox for modeling and optimization in MATLAB. In: IEEE International Symposium on Computer Aided Control Systems Design. 2004. p. 284–9.Google Scholar
 YALMIP Home Page [http://users.isy.liu.se/johanl/yalmip/]. Accessed 11 May 2016.
 IBM ILOG CPLEX Highperformance mathematical programming engine. [http://www01.ibm.com/software/commerce/optimization/cplexoptimizer/]. Accessed 11 May 2016.
 GLPK (GNU Linear programming kit) [http://www.gnu.org/software/glpk/]. Accessed 11 May 2016.
 Orth D, Fleming M, Palsson B. Reconstruction and use of microbial metabolic networks: the core Escherichia coli metabolic model as an educational guide. EcoSal Plus. 2010;4:1.View ArticleGoogle Scholar
 Emmerling M, Dauner M, Ponti A, Fiaux J, Hochuli M, Szyperski T, Wüthrich K, Bailey J, Sauer U. Metabolic flux responses to pyruvate kinase knockout in Escherichia coli. Journal of bacteriology. 2002;184(1):152–64.View ArticlePubMedPubMed CentralGoogle Scholar
 Orth J, Conrad T, Na J, Lerman J, Nam H, Feist A, Palsson B. A comprehensive genome‐scale reconstruction of Escherichia coli metabolism—2011. Molecular systems biology. 2011;7(1):535.View ArticlePubMedPubMed CentralGoogle Scholar
 Bonarius H, Schmid G, Tramper J. Flux analysis of underdetermined metabolic networks: the quest for the missing constraints. Trends in Biotechnology. 1997;15(8):308–14.View ArticleGoogle Scholar
 Palsson BØ. Systems biology: properties of reconstructed networks. New York: Cambridge University Press; 2006.View ArticleGoogle Scholar
 Schilling C, Covert M, Famili I, Church G, Edwards J, Palsson B. Genomescale metabolic model of Helicobacter pylori 26695. Journal of Bacteriology. 2002;184(16):4582–93.View ArticlePubMedPubMed CentralGoogle Scholar
 Solà A, Jouhten P, Maaheimo H, SánchezFerrando F, Szyperski T, Ferrer P. Metabolic flux profiling of Pichia pastoris grown on glycerol/methanol mixtures in chemostat cultures at low and high dilution rates. Microbiol. 2007;153:281–90.View ArticleGoogle Scholar
 Solà A. Estudi del metabolisme central del carboni de Pichia pastoris, PhD Thesis. Universitat Autònoma de Barceloana: Escola Tècnica Superior d’Enginyeria; 2004.Google Scholar
 Jungo C, Rerat C, Marison IW, von Stockar U. Quantitative characterization of the regulation of the synthesis of alcohol oxidase and of the expression of recombinant avidin in a Pichia pastoris Mut + strain. Enzyme Microb Technol. 2006;39:936–44.View ArticleGoogle Scholar
 Tortajada M. Process development for the obtention and use of recombinant glycosidases: expression, modelling and immobilization, PhD Thesis. Universidad Politécnica de Valencia: Departamento de Ingeniería de Sistemas y Automática; 2012.Google Scholar
 Jordà J, de Jesus SS, Peltier S, Ferrer P, Albiol J. Metabolic flux analysis of recombinant Pichia pastoris growing on different glycerol/methanol mixtures by iterative fitting of NMRderived 13Clabelling data from proteinogenic amino acids. New Biotecnol. 2014;31(1):120–32.View ArticleGoogle Scholar