BackgroundMortierella alpina is an oleaginous fungus used in the industrial scale production of arachidonic acid (ARA). In order to investigate the metabolic characteristics at a systems level and to explore potential strategies for enhanced lipid production, a genome-scale metabolic model of M. alpina was reconstructed.ResultsThis model included 1106 genes, 1854 reactions and 1732 metabolites. On minimal growth medium, 86 genes were identified as essential, whereas 49 essential genes were identified on yeast extract medium. A series of sequential desaturase and elongase catalysed steps are involved in the synthesis of polyunsaturated fatty acids (PUFAs) from acetyl-CoA precursors, with concomitant NADPH consumption, and these steps were investigated in this study. Oxygen is known to affect the degree of unsaturation of PUFAs, and robustness analysis determined that an oxygen uptake rate of 2.0 mmol gDW−1 h−1 was optimal for ARA accumulation. The flux of 53 reactions involving NADPH was significantly altered at different ARA levels. Of these, malic enzyme (ME) was confirmed as a key component in ARA production and NADPH generation. When using minimization of metabolic adjustment, a knock-out of ME led to a 38.28% decrease in ARA production.ConclusionsThe simulation results confirmed the model as a useful tool for future research on the metabolism of PUFAs.Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-014-0137-8) contains supplementary material, which is available to authorized users.