Investigation of two metabolic engineering approaches for (R,R)-2,3-butanediol production from glycerol in Bacillus subtilis

Flux Balance Analysis (FBA) is a well-known bioinformatics tool for metabolic engineering design. Previously, we have successfully used single-level FBA to design metabolic fluxes in Bacillus subtilis to enhance (R,R)-2,3-butanediol (2,3-BD) production from glycerol. OptKnock is another powerful technique for devising gene deletion strategies to maximize microbial growth coupling with improved biochemical production. It has never been used in B. subtilis. In this study, we aimed to compare the use of single-level FBA and OptKnock for designing enhanced 2,3-BD production from glycerol in B. subtilis.Single-level FBA and OptKnock were used to design metabolic engineering approaches for B. subtilis to enhance 2,3-BD production from glycerol. Single-level FBA indicated that deletion of ackA, pta, lctE, and mmgA would improve the production of 2,3-BD from glycerol, while OptKnock simulation suggested the deletion of ackA, pta, mmgA, and zwf. Consequently, strains LM01 (single-level FBA-based) and MZ02 (OptKnock-based) were constructed, and their capacity to produce 2,3-BD from glycerol was investigated. The deletion of multiple genes did not negatively affect strain growth and glycerol utilization. The highest 2,3-BD production was detected in strain LM01. Strain MZ02 produced 2,3-BD at a similar level as the wild type, indicating that the OptKnock prediction was erroneous. Two-step FBA was performed to examine the reason for the erroneous OptKnock prediction. Interestingly, we newly found that zwf gene deletion in strain MZ02 improved lactate production, which has never been reported to date. The predictions of single-level FBA for strain MZ02 were in line with experimental findings.We showed that single-level FBA is an effective approach for metabolic design and manipulation to enhance 2,3-BD production from glycerol in B. subtilis. Further, while this approach predicted the phenotypes of generated strains with high precision, OptKnock prediction was not accurate. We suggest that OptKnock modelling predictions be evaluated by using single-level FBA to ensure the accuracy of metabolic pathway design. Furthermore, the zwf gene knockout resulted in the change of metabolic fluxes to enhance the lactate productivity.