# Team:TU Delft/Modeling/MFA/explanation

### From 2010.igem.org

Metabolic flux analysis (MFA) is a tool to study metabolic networks quantitatively. It can be used to calculate theoretical maximal yields.

### Metabolic networks - definition

A metabolic network is the collection of all the reactions that occur inside a cell. All these reactions are coupled and form a network together. This network has a couple of transport reactions with substrate (food, such as sugar and alkanes) as an input and products (such as biomass, CO_{2}, ethanol) as output. Other products of the metabolic network are energy, which is stored in molecules such as ATP and NADH, and building blocks for the cell (biomass). The substances between the substrates and the products are called intermediates. An example of this is in figure 1.

The change in the concentration of the substrate and products outside the cell can easily be calculated from measurements. This results in uptake and excretion speeds. The intracellelular rates (inside the cell) can not be figured out as easily, but if enough information is known from the uptake en excretion rates, the intracellular rates can often be calculated.

### Units in MFA

The term *rate* has occured in the previous alinea and needs to be specified further. Rates indicate how fast a reaction goes. Rates are usually defined in moles per amount of biomass per time. In CellNetAnalyzer the units *mmol per gram of dry weight biomass per hour* were used. The reaction that genereates biomass, called μ, is an exception and has the unit per hour, which is derived from biomass produced per biomass per hour.

### Steady-state assumption

In order to determine the intracellular rates, it is assumed that all the intermediates inside the cell are in a steady-state. This means that the amount of the substance is constant, meaning that the total production and the total consumption of a substance are equal to each other. This holds also for co-factors such as ATP and NADH which are used by many reactions. This means that all the rates need to have a ratio that balances out the co-factor production and consumption. This steady-state assumption for the cofactors has a very strong influence on how the fluxes are distributed through a metabolic network.

### Full flux determination

With this analysis all the intracellular fluxes are calculated to fit with a set of measured transport fluxes. For figure 1 this would mean that, if it is known that 1 mmol g_{x}^{-1} h^{-1} of glucose is going in and 1.2 mmol g_{x}^{-1} h^{-1} of CO_{2} is coming out, then the other rates can be calculated by working out all the mass balances. The calculation results in 1.32 mmol g_{x}^{-1} h^{-1} of ethanol coming out and 0.36 mmol g_{x}^{-1} h^{-1} goes through intermediate 2 to the biomass. This is essentially the metabolic flux analysis. For larger networks, however it is required to used lineair algebra or tools such as CellNetAnalyzer.

### Flux optimization

This is the analysis that is used in our iGEM project. This analysis does not try to calculated all the intracellular fluxes with a set of measured transport fluxes, instead it tries to optimize a (set of) fluxes to a set of transport fluxes. To calculate a maximal theoretical yield of a product on a substrate, the production rate needs to be optimized toward a certain substrate uptake rate.

Production by a metabolic network can be optimized with a built in function of CellNetAnalyzer. This function will calculate with a given set of substrate inflows the maximal rate for a product. An advantage of this is that not all the transport rates need to be known, the function simply finds an optimal solution.

For example, in figure 1 there is a inflow of 1 mmol g_{x}^{-1} h^{-1} substrate and the rate for product 2 needs to be maximized. The function will then assign values to all the reaction rates to make the outflow of product 2 as large as possible, while still fulfilling the steady-state condition. For this example the outflow of product 2 will be 0.6 mmol g_{x}^{-1} h^{-1}. The rest goes into biomass. The maximal yield is then 0.6 mol product 2 per mol substrate.

It is also possible to optimize the growth by maximizing μ. If there is no feasible solution for this, then the cell is not viable because it can not use the substrate for growth. This may occur for example when a lot of NADH is produced, but cannot be consumed sufficiently somewhere else in the network. Because of this the network can not maintain the steady-state condition. This situation is called redox stress.

Click here to go back to the MFA overview