No abstract

To investigate the relationship between the yield of 1,3-propanediol (1,3-PD) and the flux variation in metabolic pathways of Klebsiella pneumoniae, an optimized calculation method was constructed on basis of dynamic flux balance analysis by combining genome-scale flux balance analysis with a kinetic model of extracellular metabolites. Through optimizing calculations, a more completely expanded metabolic pathway was obtained, which includes the previously reported metabolic pathway and additional three pathways or site: a pentose phosphate pathway (PPP) elicited at the dihydroxyacetone (DHA) node to provide more reducing equivalents; a branch of synthetic amino acids at the 3-phosphoglycerate (3PG) node; and the α-ketoglutarate site in the tricarboxylic acid (TCA) cycle leading to anabolic pathways for glutamate and other amino acids. On this basis, the relationships between the dynamic flux distribution of the important nodes in the metabolic pathway and the yield of 1,3-propanediol were analyzed. First, dynamic flux change from DHA to the PPP is positively correlated with the yield. Second, variation in flux in the TCA cycle is also positively correlated with the yield of 1,3-propanediol. In addition, the influence of the feedback loop formed by the cofactor tetrahydrofolate on the flux change of TCA in the amino acid anabolic pathway was examined. These results are of important reference value and have guiding significance for the extension of the glycerol metabolism pathway in K. pneumoniae, the rational transformation of genetic engineering in bacteria, and the optimization of metabolic pathways for industrial production.

In this study, a combined optimization method was developed to optimize the N‐terminal site‐specific PEGylation of recombinant hirudin variant‐2 (HV2) with different molecular weight mPEG‐propionaldehyde (mPEG‐ALD), which is a multifactor‐influencing process. The HV2‐PEGylation with 5 kDa mPEG‐ALD was first chosen to screen significant factors and determine the locally optimized conditions for maximizing the yield of mono‐PEGylated product using combined statistical methods, including the Plackett–Burman design, steepest ascent path analysis, and central composition design for the response surface methodology (RSM). Under the locally optimized conditions, PEGylation kinetics of HV2 with 5, 10, and 20 kDa mPEG‐ALD were further investigated. The molar ratio of polyethylene glycol to HV2 and reaction time (the two most significant factors influencing the PEGylation efficiency) were globally optimized in a wide range using kinetic analysis. The data predicted by the combined optimization method using RSM and kinetic analysis were in good agreement with the corresponding experiment data. PEGylation site analysis revealed that almost 100% of the obtained mono‐PEGylated‐HV2 was modified at the N‐terminus of HV2. This study demonstrated that the developed method is a useful tool for the optimization of the N‐terminal site‐specific PEGylation process to obtain a homogeneous mono‐PEGylated protein with desirable yield.