Skip to main content

Metabolic flux analysis of coenzyme Q10 synthesized by Rhodobacter sphaeroides under the influence of different pH regulators


Coenzyme Q10 (CoQ10) is crucial for human beings, especially in the fields of biology and medicine. The aim of this experiment was to investigate the conditions for increasing CoQ10 production. At present, microbial fermentation is the main production method of CoQ10, and the production process of microbial CoQ10 metabolism control fermentation is very critical. Metabolic flux is one of the most important determinants of cell physiology in metabolic engineering. Metabolic flux analysis (MFA) is used to estimate the intracellular flux in metabolic networks. In this experiment, Rhodobacter sphaeroides was used as the research object to analyze the effects of aqueous ammonia (NH3·H2O) and calcium carbonate (CaCO3) on the metabolic flux of CoQ10. When CaCO3 was used to adjust the pH, the yield of CoQ10 was 274.43 mg·L−1 (8.71 mg·g−1 DCW), which was higher than that of NH3·H2O adjustment. The results indicated that when CaCO3 was used to adjust pH, more glucose-6-phosphate (G6P) entered the pentose phosphate (HMP) pathway and produced more NADPH, which enhanced the synthesis of CoQ10. At the chorismic acid node, more metabolic fluxes were involved in the synthesis of p-hydroxybenzoic acid (pHBA; the synthetic precursor of CoQ10), enhancing the anabolic flow of CoQ10. In addition, Ca2+ produced by the reaction of CaCO3 with organic acids promotes the synthesis of CoQ10. In summary, the use of CaCO3 adjustment is more favorable for the synthesis of CoQ10 by R. sphaeroides than NH3·H2O adjustment. The migration of metabolic flux caused by the perturbation of culture conditions was analyzed to compare the changes in the distribution of intracellular metabolic fluxes for the synthesis of CoQ10. Thus, the main nodes of the metabolic network were identified as G6P and chorismic acid. This provides a theoretical basis for the modification of genes related to the CoQ10 synthesis pathway.


Coenzyme Qn (CoQn) is a lipid-soluble quinone compound formed by conjugation of the benzoquinone group at the head and a polyisoprene hydrophobic side chain at the tail [1]. It exists in the cells of various organisms. The subscript n of CoQn represents the number of isoprene repeating units in the molecule. The number of isoprene units on the side chains of CoQn varies from organism to organism. The CoQn molecule in the human body contains 10 isoprene units, namely CoQ10. CoQ10 was discovered by Frederick Crane and colleagues in 1957 [2]. CoQ10 is the only lipid-soluble antioxidant synthesized endogenously in the human body [3]. Reduced CoQ10 (Fig. 1) is the main form present in the body, accounting for more than 90% of the total CoQ10 in the human circulation [4, 5]. The reduced CoQ10 releases electrons during the electron transfer process, and the electrons are transferred through the respiratory chain to eventually produce ATP to provide energy for the body. CoQ10 deficiency in the human body will cause some diseases. In addition, CoQ10 is also beneficial to the skin. It can improve the skin and enhance the expression of collagen and elastin [6, 7]. It has been shown that CoQ10 levels in skin and skin surface lipids decline with age. CoQ10 is very important to human beings and has great value. It is used in many fields such as medicine, food, and cosmetics.

There are three main methods to produce CoQ10. However, there is a preference for the microbial fermentation method, which is considered the most feasible, over the biological tissue extraction method and chemical synthesis method [8,9,10]. CoQ10 synthesized by microbial fermentation is an all-trans conformation with biological activity, and the method does not produce other optical isomers [11, 12]. In addition, microorganisms grow and reproduce quickly, can utilize inexpensive raw materials, and produce fewer contaminants during the production process. One of the microorganisms that can synthesize CoQ10 efficiently and has been used in industrial production is Rhodobacter sphaeroides [13, 14]. It is a gram-negative purple non sulfur bacterium (PNSB) that can grow aerobically, anaerobically, photoautotrophically, and heterotrophically [15,16,17]. It has many other functions, such as biological hydrogen production, synthesis of carotenoids, fixation of CO2 and N2, remediation of heavy metal pollution and power generation [18,19,20,21].

Extracellular organic acids are continuously produced by R. sphaeroides during the fermentation of CoQ10. This leads to a continuous decrease in the pH of the environment and affects the physiological activity of the bacteria. Organic acids are generally neutralized by adding CaCO3 or NH3·H2O.

Fig. 1
figure 1

Chemical structure of reduced CoQ10 (CoQ10H2)

The process of CoQ10 biosynthesis in R. sphaeroides is complex (Fig. 2). First, the head quinone ring and the tail polyisoprene chain of CoQ10 need to be synthesized separately [22]. The synthesis of the quinone ring is the shikimate pathway. This pathway uses erythrose-4-phosphate (E4P) and phosphoenolpyruvate (PEP) as the starting substrates, which are converted to chorismic acid in various steps. The chorismic acid is then converted to p-hydroxybenzoic acid (pHBA), a precursor of the quinone ring, by the action of the chorismic acid lyase encoded by the UbiC gene. The polyisoprene chain is synthesized via the methylerythritol phosphate (MEP) pathway. The starting substances of this pathway are glyceraldehyde-3-phosphate and pyruvate, which are converted to isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) by various reactions, and then formed into decaprenyl diphosphate (DPP) by various reactions. DPP and pHBA are converted to 3-decaprenyl-4-hydroxybenzoate via the ubiquinone pathway by the action of p-hydroxybenzoic acid polyisoprene transferase encoded by the UbiA gene and enter the quinone ring modification pathway. After decarboxylation, hydroxylation, and methylation processes, CoQ10 is finally formed.

Fig. 2
figure 2

CoQ10 synthesis pathway in R. sphaeroides

Methods and materials

Microorganism and culture media

The R. sphaeroides used in this study was the F3-40 strain (from the College of Life Sciences, Fujian Normal University, Fuzhou, China). The culture media included agar slant culture medium (containing 18 g·L− 1 glucose, 4.5 g·L− 1 yeast powder, 20 g·L− 1 agar powder, 2 g·L− 1 monosodium glutamate, 5.25 g·L− 1 (NH4)2SO4, 4.2 g·L− 1 NaCl, 8.75 g·L− 1 MgSO4, 1.5 g·L− 1 KH2PO4, 1.2 g·L− 1 A), seed culture medium (containing 18 g·L− 1 glucose, 4.5 g·L− 1 yeast powder, 2 g·L− 1 monosodium glutamate, 5.25 g·L− 1 (NH4)2SO4, 4.2 g·L− 1 NaCl, 8.75 g·L− 1 MgSO4, 1.5 g·L− 1 KH2PO4, 1.2 g·L− 1 A), and fermentation culture medium (containing 35.6 g·L− 1 glucose, 6.39 g·L− 1 corn steep powder, 7.39 g·L− 1 monosodium glutamate, 6.28 g·L− 1 (NH4)2SO4, 2.5 g·L− 1 NaCl, 12 g·L− 1 MgSO4, 0.8 g·L− 1 KH2PO4, 0.1 g·L− 1 MnSO4, 1.6 g·L− 1 A, 0.1 g·L− 1 B).

Cultivation method

Two colonies of R. sphaeroides cultured on agar plates for 6 days were picked out and transferred into the primary seed medium (50 mL/250 mL), and incubated at 32 °C with shaking at 220 rpm for 24 h. Then, they were transferred into the secondary seed medium (350 mL/2500 mL) at 10% inoculum and incubated at 32 °C with shaking at 220 rpm for 22 h. Finally, the secondary seed medium was inoculated at 15% inoculum into a 10 L fermenter. The inoculated fermenter was filled with 7 L. The temperature was 34 °C, the stirring speed was 300 rpm, the aeration rate was 3 L·min− 1, and the tank pressure was 0.04 MPa. The pH was adjusted by adding CaCO3 and NH3·H2O during the fermentation process. Samples were taken at 3 h intervals to determine various parameters.

Determination of glucose

A bioanalytical sensor (model SBA-40D, from Biology Institute of Shandong Academy of Sciences) was used for analysis. The fermentation broth was diluted by a certain multiple and then filtered. Twenty-five microliters of filtrate was taken for injection determination. The determination was performed three times in parallel and the average value was taken.

Determination of cell dry weight (DCW)

Five milliliters of the fermentation broth was centrifuged (1 × 104 rpm, 10 min), and the supernatant was removed. After washing with distilled water, the samples were centrifuged again, and the supernatant was removed. The organism in the centrifuge tube was dried in an oven at 100 °C to a constant weight. For fermentation broth containing CaCO3, a drop of 6 mol·L− 1 HCl solution was added to dissolve the residual CaCO3 in the fermentation broth before centrifugation.

Extraction of CoQ10 and determination of its content

CoQ10 was extracted by ultrasonic extraction [23, 24]. Five milliters of the fermentation broth was added to a 50 mL brown volumetric flask. A drop of 6 mol·L− 1 HCl solution was added. Shake well and add 10 mL of acetone. After mixing, 0.5 mL of 30% H2O2 solution was added. Finally, 20–30 mL ethanol was added. The volumetric flask was placed into the ultrasonic cleaner for 1 min, and the volumetric flask was not capped. Subsequently, the volumetric flask was fixed with ethanol and capped. A strip of cardboard was used to separate the cap from the neck of the bottle to retain the gap. The volumetric flask was again placed in an ultrasonic cleaner for 45 min, with the temperature controlled below 35 °C. At the end of the treatment, the flask was well shaken and allowed to stand for 30 min. Finally, the supernatant was filtered through a 0.22 μm organic filter.

The filtrate was collected for the determination of CoQ10 by high-performance liquid chromatography (HPLC). The chromatographic column was Hypersil ODS-SP (4.6 mm × 100 mm, 5 μm). The mobile phase consisted of anhydrous methanol and anhydrous ethanol in a ratio of 65: 35 (v/v). The detection wavelength was set at 275 nm. The column temperature was 30 °C, the injection volume was 20 µL, the flow rate was 1.1 mL·min− 1, and the elution time was 15 min. The calculation formula of CoQ10 concentration in the sample is:

$$c = \frac{{S \cdot {V_o} \cdot K \cdot {c_o}}}{{{S_o} \cdot V}}$$

where c is the concentration of CoQ10 in the test sample (mg·L− 1); co is the concentration of CoQ10 in the standard sample (mg·L− 1); K is the dilution multiple of the test sample; S is the peak area of the test sample (mAU·s); So is the peak area of the standard sample (mAU·s); V is the injection volume of the test sample (µL); and Vo is the injection volume of the standard sample (µL).

Determination of various free amino acid concentrations

The concentration of free amino acids in the fermentation broth was determined by HPLC. A Venusil-AA amino acid analysis column was used. Mobile phases A and B were 0.1 mol·L− 1 CH3COONa solution (pH 6.5) and 80% acetonitrile solution, respectively. The detection wavelength was set at 254 nm. The column temperature was 40 °C, the injection volume was 20 µL, and the flow rate was 1 ml·min− 1. Gradient elution was performed.

Determination of organic acid concentration

The concentration of organic acids in the fermentation broth was determined by HPLC [25]. The chromatographic column was an Inertsil ODS-SP C18. Mobile phases A and B were 0.01 mol·L− 1 KH2PO4 buffer (pH adjusted to 2.8 by H3PO4) and acetonitrile solution, respectively, and were injected into the column at a ratio of 95: 5 (v/v). The detection wavelength was set at 215 nm. The column temperature was 25 °C, the injection volume was 20 µL, and the flow rate was 1 ml·min− 1.


Construction of a metabolic network for CoQ10 synthesis in R. sphaeroides

The synthesis of CoQ10 by R. sphaeroides belongs to growth-coupled fermentation. Therefore, when analyzing the metabolic flux of CoQ10 synthesis, the carbon source consumed by the growth of the bacterium cannot be ignored, and the synthesis pathway of the basic cellular components should be considered. The main cellular components and their contents refer to the experimental data measured by Saheed lmam [26]. The fractions of R. sphaeroides under aerobic conditions are shown in Table 1. It contained 17.6% pHBA. Among fatty acids, oleic acid (C18:1) has the highest content, accounting for 85%. Stearic acid (C18:0), soft fatty acid (C16:0), and palmitoleic acid (C16:1) accounted for 9%, 5%, and 1%, respectively.

Table 1 Distribution of cellular fraction components of R. sphaeroides under aerobic contributions

Since R. sphaeroides uses glucose as the sole carbon source, the glyoxylate cycle pathway is largely inoperative and is not considered in this metabolic network. The intermediate reactions without branching points are simplified into a reaction equation. For intermediate metabolites with branching points, the rate of production was equal to the rate of consumption, and the net accumulation was 0, assuming that these metabolites were in the proposed steady state.

It is assumed that NADPH produced by the HMP pathway and the TCA pathway were not oxidatively phosphorylated but were all used for biosynthesis. The fermentation supernatant of 39–42 h in the late logarithmic growth phase of the organism was analyzed for composition. The results showed that the supernatant basically did not contain free amino acids. The content of organic acid metabolic byproducts such as formic acid, acetic acid, propionic acid, butyric acid, lactic acid, pyruvic acid, and citric acid was also very low. Therefore, the influence of these organic acids on metabolic flux can be ignored in the metabolic network studied.

Modeling of metabolic flux equilibrium for CoQ10 synthesis by R. sphaeroides

Metabolic flux analysis is a metabolic network chemometric used to describe the conversion of substrates into metabolites and cellular composition. The distribution of carbon flow under the condition of carbon source limitation is generally considered. The equilibrium equation of each metabolite can be expressed as:


where Xmet is the metabolite concentration; rmet is the net rate of metabolite formation; and u is the dilution factor of the culture medium.

Since most of the intermediate metabolites have low intracellular concentrations, the dilution effect can be negligible. According to the proposed steady-state assumption that the concentration of intermediate metabolites is in equilibrium (dXmet/dt = 0), it follows that: rmet = 0.

The net rate of intermediate metabolite formation used is expressed in the form of a matrix equation as:

$$S \cdot v = 0$$

where S is the n × m stoichiometric coefficient matrix, and the rows of the matrix represent the stoichiometric coefficients of the rate reactions of the intermediate metabolites; v is an m-dimensional column vector; n is the number of intermediate metabolites; and m is the total number of metabolic reaction rates.

The degree of freedom to be solved is F = m - n. The number of linearly uncorrelated rates is determined experimentally. The matrix equation S · v = 0 has a unique solution as long as it is greater than or equal to the degree of freedom F. Then, the undetermined metabolic reaction rates can be calculated, and thus, the flux distribution of the metabolic network can be determined.

Based on the intracellular composition and biosynthetic analysis of the various metabolites in the metabolic network, the following proposed steady-state equilibrium equation of the branch point intermediates was obtained:

G6P: r1-r2-r10-r11-0.0241r21 = 0.

F6P: r2-r3 + 0.6667r13 + r14 = 0.

Ru5P: r11-0.5288r12-r13-2r14-0.063r15 = 0.

GA3P: 2r3-r4 + 0.3333r13 + 0.063r15-10r16-0.313r21 = 0.

PG: r4-r5-0.1984r12-0.8987r20-0.1155r21 = 0.

PEP: r5-r6-2r14-0.063r15 = 0.

Pyr: r6-r7 + 0.063r15-9r16-2.0895r17-r18-0.4575r19 = 0.

AcCoA: r7-r8-0.5511r17-8.0862r21 = 0.

OAA: -r8 + r9-0.1941r12 + r18-1.4061r19 = 0.

Chr: r14-0.3878r15-r16 = 0.

AKG: r8-r9-1.3165r22 = 0.

NADPH: r8 + 2r11-0.911r12-r14-0.3878r15-10r16-2.0895r17-4.8055r19-1.1843r20-7.3028r21-3.2579r22 = 0.

When these equilibrium equations were constructed into a matrix equation, the rank of the matrix n = 12 was calculated. If the unknown rate variable J = 22, then the matrix has degrees of freedom F = J - n = 10. If the matrix equation is to have a unique solution, at least 10 rate variables must be determined. Of these 22 rate variables, r10, r12, r15, r17, r19, r20, r21, and r22 were calculated by measuring the synthesis rate of cell dry weight. Through the sugar consumption rate, r1 can be obtained. Through the synthesis rate of CoQ10, r16 can be obtained. From the data of these 10 rate variables, the remaining 12 unknown rate variables can be calculated using MATLAB software.

Analysis of the effects of two pH regulators on the synthetic metabolic flux of CoQ10

The synthesis rate of CoQ10 was significantly higher in the late logarithmic growth phase than in the mid and early growth periods when the pH was adjusted by CaCO3 or NH3·H2O. Then, metabolic flux analysis was performed to quantify the flux distribution of the CoQ10 synthesis pathway in the late logarithmic growth phase of these two fermentations. Flux analysis was performed using the 42nd hour of fermentation as the time point. The variable rates of metabolites (r1 to r22) were divided by the rate of glucose consumption (r1), and the result was the respective values of r1 to r22 (Table 2).

Table 2 The metabolic flux distribution of CoQ10 under different pH adjustment modes
$$C = \frac{{Flux\left( {CaC{O_3}adjustment} \right) - Flux(N{H_3} \cdot {H_2}Oadjustment)}}{{Flux(N{H_3} \cdot {H_2}Oadjustment)}}$$

It can be seen from Fig. 3 and Table 2 that the two pH regulators caused significant changes in the flux distribution of G6P, Ru5P, and PEP nodes. The values of r2 to r9 on the EMP pathway and TCA pathway were smaller, while the values of r11 and r13 on the HMP pathway were larger when CaCO3 was used to adjust pH compared to NH3·H2O. This indicates that CaCO3 adjustment strengthened the HMP pathway and weakened the EMP pathway and TCA pathway.

Fig. 3
figure 3

Metabolic network of CoQ10 synthesis in R. sphaeroides. AcCoA (Acetyl coenzyme A); AKG (a-Ketoglutarate); BIOM (biomass); Chr (Chorismate); F6P (Fructose-6-phosphate); G6P (Glucose-6-phosphate); GA3P (Glyceraldehyde-3-phosphate); GLC (Glucose); OAA (Oxaloacetate); PEP (Phosphoenolpyruvate); PG (3-Phosphoglycerate); PHB (poly-β-hydroxybutyrate); PSACH (Polysaccharides); PYR (Pyruvate); Ru5P (Ribulose-5-phosphate); NC (Nucleic Acids)

The change in r16 showed that the metabolic flux from Chr to CoQ10 under CaCO3 adjustment was 0.13, while NH3·H2O adjustment resulted in 0.03. The enhanced metabolic flux from Chr to CoQ10 increased the demand for NADPH, prompting more G6P to flow into the HMP pathway to form NADPH. However, CaCO3 adjustment has little effect on biosynthesis. Thus, the increased portion of G6P entering the HMP pathway was not used for the synthesis of bacterial cells, but formed G6P and GA3P and returned to the EMP pathway for oxidative catabolism in the TCA pathway.

The organic acids produced during fermentation react with CaCO3 to form Ca2+, which activates enzymes related to the CoQ10 synthesis pathway and stimulates the production of reactive oxygen species (ROS) that can damage cell membranes [27]. Since CoQ10 has antioxidant properties, it can alleviate the damage caused by ROS on cell membranes. With the increase in Ca2+ concentration, the formation of intracellular ROS was stimulated, which promoted the synthesis of CoQ10 and increased the yield of CoQ10.

A certain amount of Ca2+ has been proven to increase the yield of CoQ10. Whether different concentrations of Ca2+ can all achieve this effect is an aspect worthy of future research, and perhaps the best Ca2+ concentration can be found. Ca2+ is an activator of many enzymes, and the hypothesis that increased CoQ10 yield may also be due to its activation of the activity of related enzymes in the synthesis pathway can also be explored.

Changes in metabolic flux in G6P and Ru5P nodes

G6P is the common starting metabolite of both the EMP pathway and the HMP pathway. When CaCO3 was used to adjust the pH, the metabolic flux from G6P to Ru5P was 70.5, and the metabolic flux from G6P to F6P was 27.9 (Fig. 4). When NH3·H2O was used to adjust the pH, the metabolic flux from G6P to Ru5P was 68, and the metabolic flux from G6P to F6P was 30.4. The difference in metabolic flux caused by these two pH adjustment methods was equal between G6P to Ru5P and G6P to F6P. However, the increased metabolic flux of Ru5P had almost no effect on nucleic acid and histidine synthesis in the bacterium, which was mostly converted to F6P and GA3P (approximately 2.01). The metabolic flux from G6P to F6P migrates due to the addition of CaCO3 during fermentation or the effect of the reduction of NH4+ on fermentation. Some of the metabolic streams passed through the HMP pathway before entering the EMP pathway. This process has a completely different meaning than that of F6P, which is directly converted to F6P by the EMP pathway. Although the loss of carbon source material flux is small, the former provides more reducing power NADPH. Since the synthesis of isoprene pyrophosphate, the precursor of the CoQ10 side chain, requires NADPH to provide reducing power, 10 mol NADPH is directly involved for every 1 mol CoQ10 synthesis. Then, G6P enters the EMP pathway through the HMP pathway to provide NADPH for the synthesis of CoQ10.

Fig. 4
figure 4

The metabolic flux distribution of G6P and Ru5P nodes under different pH adjustment modes

Analysis of the effects of two pH regulators on the synthesis of CoQ10 by R. sphaeroides

The two pH regulators, CaCO3 and NH3·H2O, had almost no effect on the biosynthesis of the bacterium during the synthesis of CoQ10 by R. sphaeroides. The fermentation times used for conditioning with CaCO3 and NH3·H2O were 50 h and 45 h, respectively (Table 3). The difference in CoQ10 yield caused by the two pH regulators was significant (Fig. 5). When the fermentation endpoint was reached, the CoQ10 yield with CaCO3 conditioning was 274.43 mg·L-1 (8.71 mg·g-1 DCW), which was higher than that of NH3·H2O conditioning of 226.53 mg·L-1 (7.55 mg·g-1 DCW). The CoQ10 synthesis rates under the two pH adjustments were basically the same within 24 h from the beginning of fermentation. However, between 24 and 42 h, the rate of CoQ10 synthesis under NH3·H2O regulation was faster than that under CaCO3 regulation. The accelerated synthesis of CoQ10 in this process produces more organic acids. As an alkaline substance, NH3·H2O with high solubility can quickly resist the pH drop caused by organic acids, while CaCO3 with poor solubility reacts slowly. Therefore, NH3·H2O can better regulate the pH in this stage, which leads to a faster rate of CoQ10 synthesis. After 42 h of fermentation, the rate of NH3·H2O-regulated CoQ10 synthesis slowed down, whereas the CaCO3-regulated one was still able to maintain a high rate of growth. The reason was mainly that NH3·H2O was consumed in large quantities in the early stage, the concentration decreased dramatically, and the pH-regulating ability was weakened.

Table 3 Comparison of various parameters of the fermentation process under two pH adjustment methods
Fig. 5
figure 5

Effects of two pH regulators on the yield of CoQ10

Conclusions and future perspectives

In summary, pH adjustment with CaCO3 was more favorable for the synthesis of CoQ10 by R. sphaeroides compared with NH3·H2O adjustment. This study provides a reference for researchers in the selection of pH regulators. In addition, the main nodes of the metabolic network were identified as G6P and chorismic acid. This provides a theoretical basis for the modification of genes related to the CoQ10 synthesis pathway.

Due to the high value of CoQ10, an increasing number of researchers are investigating it. Its potential functions are constantly being developed and explored. Microbial fermentation as a production method will be the choice of more people in the future, and one of its most important goals is to maximize the yield of the desired fermentation product. The method is not only the well-known improvement and selection of strains, but also the analysis of the culture environment favorable to the synthesis of products according to the metabolic network of the microorganism, as used in this article. The application of metabolic engineering principles to improve the yield of microbial fermentation products has become a hot research topic.


  1. Liu J, Wu Q, He D, Ma T, Du L, Dui W, Guo X, Jiao R. Drosophila sbo regulates lifespan through its function in the synthesis of coenzyme Q in vivo. J Genet Genomics. 2011;38(6):225–34.

    Article  CAS  PubMed  Google Scholar 

  2. Raizner AE, Quiñones MA. Coenzyme Q10 for Patients with Cardiovascular Disease: JACC Focus Seminar. J Am Coll Cardiol. 2021;77(5):609–19.

    Article  CAS  PubMed  Google Scholar 

  3. Bentinger M, Tekle M, Dallner G. Coenzyme Q–biosynthesis and functions. Biochem Biophys Res Commun. 2010;396(1):74–9.

    Article  CAS  PubMed  Google Scholar 

  4. Kuriyama N, Nakamura T, Nakazawa H, Wen T, Berra L, Bittner EA, Goverman J, Kaneki M. Bioavailability of reduced Coenzyme Q10 (Ubiquinol-10) in burn patients. Metabolites. 2022;12(7):613.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Ruiz-Jiménez J, Priego-Capote F, Mata-Granados JM, Quesada JM, Luque de Castro MD. Determination of the ubiquinol-10 and ubiquinone-10 (coenzyme Q10) in human serum by liquid chromatography tandem mass spectrometry to evaluate the oxidative stress. J Chromatogr. 2007;1175(2):242–8.

    Article  Google Scholar 

  6. Vollmer DL, West VA, Lephart ED. Enhancing skin health: by oral administration of natural Compounds and Minerals with Implications to the dermal Microbiome. Int J Mol Sci. 2018;19(10):3059.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Žmitek K, Pogačnik T, Mervic L, Žmitek J, Pravst I. The effect of dietary intake of coenzyme Q10 on skin parameters and condition: results of a randomised, placebo-controlled, double-blind study. Biofactors. 2017;43(1):132–40.

    Article  PubMed  Google Scholar 

  8. Tian Y, Yue T, Yuan Y, Soma PK, Lo YM. Improvement of cultivation medium for enhanced production of coenzyme Q10 by photosynthetic Rhodospirillum rubrum. Biochem Eng J. 2010;51(3):160–6.

    Article  CAS  Google Scholar 

  9. Tian Y, Yue T, Yuan Y, Soma PK, Williams PD, Machado PA, Fu H, Kratochvil RJ, Wei CI, Lo YM. Tobacco biomass hydrolysate enhances coenzyme Q10 production using photosynthetic Rhodospirillum rubrum. Bioresour Technol. 2010;101(20):7877–81.

    Article  CAS  PubMed  Google Scholar 

  10. Wang Y, Chen S, Liu J, Lv P, Cai D, Zhao G. Efficient production of coenzyme Q10 from acid hydrolysate of sweet sorghum juice by Rhodobacter sphaeroides. RSC Adv. 2019;9(39):22336–42.

  11. Fan J, Xu W, Xu X, Wang Y. Production of Coenzyme Q10 by microbes: an update. World J Microbiol Biotechnol. 2022;38(11):194.

    Article  CAS  PubMed  Google Scholar 

  12. Lee SQ, Tan TS, Kawamukai M, Chen ES. Cellular factories for coenzyme Q10 production. Microb Cell Fact. 2017;16(1):39.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Yang Y, Li L, Sun H, Li Z, Qi Z, Liu X. Improving CoQ10 productivity by strengthening glucose transmembrane of Rhodobacter sphaeroides. Microb Cell Fact. 2021;20(1):207.

  14. Zhang J, Gao D, Jiajia C, Liu H, Qi Z. Improving coenzyme Q10 yield of Rhodobacter sphaeroides via modifying redox respiration chain. Biochem Eng J. 2018;135:98–104.

  15. Lee S, Yu J, Kim YH, Min J. Optimized Rhodobacter sphaeroides for the production of Antioxidants and the pigments with antioxidant activity. Mol Biotechnol. 2022;65(1):131–5.

    Article  PubMed  Google Scholar 

  16. Peng L, Lou W, Xu Y, Yu S, Liang C, Alloul A, Song K, Vlaeminck SE. Regulating light, oxygen and volatile fatty acids to boost the productivity of purple bacteria biomass, protein and co-enzyme Q10. Sci Total Environ. 2022;822:153489.

    Article  CAS  PubMed  Google Scholar 

  17. Yu J, Moon SK, Kim YH, Min J. Isoprene production by Rhodobacter sphaeroides and its antimicrobial activity. Res Microbiol. 2022;173:103938.

    Article  CAS  PubMed  Google Scholar 

  18. Cadirci BH. An electricity production study by Rhodobacter sphaeroides. Int J Hydrog Energy. 2018;43(38):18001–6.

    Article  CAS  Google Scholar 

  19. Giotta L, Agostiano A, Italiano F, Milano F, Trotta M. Heavy metal ion influence on the photosynthetic growth of Rhodobacter sphaeroides. Chemosphere. 2006;62(9):1490–9.

    Article  CAS  PubMed  Google Scholar 

  20. Imam S, Noguera DR, Donohue TJ. Global insights into energetic and metabolic networks in Rhodobacter sphaeroides. BMC Syst Biol. 2013;7(1):89.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Liu S, Zhang G, Li X, Wu P, Zhang J. Enhancement of Rhodobacter sphaeroides growth and carotenoid production through biostimulation. J Environ Sci. 2015;33:21–8.

    Article  CAS  Google Scholar 

  22. He S, Lu H, Zhang G, Ren Z. Production of coenzyme Q10 by purple non-sulfur bacteria: current development and future prospect. J Clean Prod. 2021;307:127326.

    Article  CAS  Google Scholar 

  23. Wang Y, Fan L, Huang J, Liang J, Wang X, Ren Y, Li H, Yue T, Gao Z. Evaluation of chemical composition, antioxidant activity, and gut microbiota associated with pumpkin juice fermented by Rhodobacter sphaeroides. Food Chem. 2022;401(4):134122.

    PubMed  Google Scholar 

  24. Zhang L, Liu L, Wang KF, Xu L, Zhou L, Wang W, Li C, Xu Z, Shi T, Chen H, Li Y, Xu H, Yang X, Zhu Z, Chen B, Li D, Zhan G, Zhang SL, Zhang LX, Tan GY. Phosphate limitation increases coenzyme Q10 production in industrial Rhodobacter sphaeroides HY01. Synth Syst Biotechnol. 2019;4(4):212–9.

  25. de Sena Aquino AC, Azevedo MS, Ribeiro DH, Costa AC, Amante ER. Validation of HPLC and CE methods for determination of organic acids in sour cassava starch wastewater. Food Chem. 2015;172:725–30.

    Article  PubMed  Google Scholar 

  26. Imam S, Yilmaz S, Sohmen U, Gorzalski AS, Reed JL, Noguera DR, Donohue TJ. iRsp1095: a genome-scale reconstruction of the Rhodobacter sphaeroides metabolic network. BMC Syst Biol. 2011;5(1):116.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Görlach A, Bertram K, Hudecova S, Krizanova O. Calcium and ROS: a mutual interplay. Redox Biol. 2015;6:260–71.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


This work was supported by the Natural Science Foundation of Fujian Province, China (Grant Number 2015J01127).

Author information

Authors and Affiliations



Yujun Xiao and Yi Zheng wrote the manuscript and analyzed the results. Yong Zhou performed the experiments and analyzed the results. Chaofan Yu and Ting-E Ye supervised the research. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Yi Zheng.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xiao, Y., Zheng, Y., Zhou, Y. et al. Metabolic flux analysis of coenzyme Q10 synthesized by Rhodobacter sphaeroides under the influence of different pH regulators. Microb Cell Fact 22, 206 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: