Mortierella 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. This 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. The simulation results confirmed the model as a useful tool for future research on the metabolism of PUFAs.

Synthetic biology aims to develop programmable tools to perform complex functions such as redistributing metabolic flux in industrial microorganisms. However, development of protein-level circuits is limited by availability of designable, orthogonal, and composable tools. Here, with the aid of engineered viral proteases and proteolytic signals, we build two sets of controllable protein units, which can be rationally configured to three tools. Using a protease-based dynamic regulation circuit to fine-tune metabolic flow, we achieve 12.63 g L⁻¹ shikimate titer in minimal medium without inducer. In addition, the carbon catabolite repression is alleviated by protease-based inverter-mediated flux redistribution under multiple carbon sources. By coordinating reaction rate using a protease-based oscillator in E. coli, we achieve d-xylonate productivity of 7.12 g L⁻¹ h⁻¹ with a titer of 199.44 g L⁻¹. These results highlight the applicability of programmable protein switches to metabolic engineering for valuable chemicals production.