Reconstruction and analysis of genome-scale metabolic model of weak Crabtree positive yeast Lachancea kluyveri

Abstract Background Lachancea kluyveri, a weak Crabtree positive yeast, has been extensively studied for its unique URC pyrimidine catabolism pathway. It produces more biomass than Saccharomyces cerevisiae due to the underlying weak Crabtree effect and resorts to optimal fermentation only in oxygen limiting conditions that render it a suitable host for industrial-scale protein production. Ethyl acetate, an important industrial chemical, has been demonstrated to be a major overflow metabolite during aerobic batch cultivation with a specific rate of 0.12 g per g dry weight per hour. Here, we reconstruct a genome-scale metabolic model of the yeast to better explain the observed phenotypes and aid further hypothesis generation. Results We report the first genome-scale metabolic model, iPN730, using Build Fungal Model in KBase workspace. The inconsistencies in the draft model were semi-automatically corrected using literature and published datasets. The curated model comprises of 1235 reactions, 1179 metabolites, and 730 genes distributed in 8 compartments (organelles). The in silico viability in different media conditions and the growth characteristics in various carbon sources show good agreement with experimental data. Dynamic flux balance analysis describes the growth dynamics, substrate utilization and product formation kinetics in various oxygen-limited conditions. The URC pyrimidine degradation pathway incorporated into the model enables it to grow on uracil or urea as the sole nitrogen source. Conclusion The genome-scale metabolic construction of L. kluyveri will provide a better understanding of metabolism, particularly that of pyrimidine metabolism and ethyl acetate production. Metabolic flux analysis using the model will enable hypotheses generation to gain a deeper understanding of metabolism in weakly Crabtree positive yeast and in fungal biodiversity in general.