Reconstruction of a genome-scale metabolic model and in silico analysis of the polymalic acid producer Aureobasidium pullulans CCTCC M2012223

Aureobasidium pullulans is a yeast-like fungus used for producing biopolymers e.g. polymalic acid (PMA) and pullulan. A high PMA producing strain, A. pullulans CCTCC M2012223, was isolated and sequenced in our previous study. To understand its metabolic performance, a genome-scale metabolic model, iZX637, consisting of 637 genes, 1347 reactions and 1133 metabolites, was reconstructed based on genome annotation and literature mining studies. The iZX637 model was validated by simulating cell growth, utilization of carbon and nitrogen sources, and gene essentiality analysis in A. pullulans. We further validated our model, designed a simulation program for the prediction of PMA production using experimental data, and further analyzed the carbon flux distribution and change with increasing PMA synthesis rates. Through the calculated flux distribution, NADPH- and NADH-dependent methylenetetrahydrofolate dehydrogenase (MTHFD) were found to be associated with the transfer of reducing equivalents from NADPH to NADH for supplementing NADH in the metabolic network. Furthermore, under the high PMA synthesis rate, a large amount of carbon flux was through pyruvate into malic acid via the reductive TCA cycle. Thus, pyruvate carboxylase, which can convert pyruvate to oxaloacetate with CO2 fixation, may also be an important target for PMA synthesis. These results illustrated that the model iZX637 was a powerful tool for optimization of A. pullulans metabolism and identification of targets for guiding metabolic engineering.