Multi-parameters of physiological metabolism, environmental state, cell morphology, and manipulated variable, have mutual influences and constraints on coenzyme Q 10 (CoQ 10) biosynthesis of Rhodobater sphaeroides. In this work, a multi-parameter fusional correlation-analysis were applied to illuminate and optimize the crucial metabolism on CoQ 10 biosynthesis. The quantitative response model for accurate living cells concentration, morphological deformation, specific oxygen consumption rate (Q O2), and specific product biosynthesis rate were established, the results revealed that maintain the optimal electrical conductivity(Cond.) and Q O2 at 9.0±0.5 mS/cm and 0.06-0.08±0.01×10 -7mmol/cell/h, respectively, by adjusting the nutrient feeding and oxygen transfer rate, the cell’s morphology and oxygen uptake rate were strengthened. CoQ 10 production in 50L bioreactors was successfully improved by 35% to 3384mg/L than that of control, which would be more powerful and effective for scale up in the large-scale production of CoQ 10.

A well known \emph{in silico} analysis in metabolic flux analysis is the structural identifiability analysis. It comes from the fact that some enrichment measurement sets cannot uniquely elucidate all intracellular fluxes. To the best of our knowledge, the structural identifiability analysis of dynamic isotope experiments is not available in the literature. In this work, it is shown that if one measurement plan makes the dynamic isotopic fractions balance equations structurally identifiable then for any arbitrary small time interval the plan also makes the equations structurally identifiable. Based on the fact, in order to resolve the local structural identifiablity problem of the dynamic isotopic fractions balance equations approximated with piecewise affine intracellular fluxes, one should check the local structural identifiablity for the corresponding cascaded linear time invariant system at each sampling point with the approach proposed in our earlier work (Lin \emph{et al.}, Math Biosci. 2018; 300:122-12).

A well known issue in metabolic flux analysis with labelling experiments is the structural identifiability analysis. It comes from the fact that some enrichment measurement sets cannot uniquely elucidate all intracellular fluxes. To the best of our knowledge, the structural identifiability analysis of dynamic isotope experiments is not available in the literature. In this work, it is shown that if one measurement plan makes the dynamic isotopic fractions balance equations structurally identifiable then for any arbitrary small time interval the plan also makes the equations structurally identifiable. Based on the fact, in order to resolve the local structural identifiablity problem of the dynamic isotopic fractions balance equations approximated with piecewise affine intracellular fluxes, one should check the local structural identifiablity for the corresponding cascaded linear time invariant system at each sampling point with the approach proposed in our earlier work (Lin \emph{et al.}, Math Biosci. 2018; 300:122-129).