Skip to main content

Systems metabolic engineering of Corynebacterium glutamicum for the bioproduction of biliverdin via protoporphyrin independent pathway



Biliverdin, a prospective recyclable antioxidant and one of the most important precursors for optogenetics, has received growing attention. Biliverdin is currently produced by oxidation of bilirubin from mammalian bile using chemicals. However, unsustainable procedures of extraction, chemical oxidation, and isomer separation have prompted bio-based production using a microbial cell factory.


In vitro thermodynamic analysis was performed to show potential candidates of bottleneck enzymes in the pathway to produce biliverdin. Among the candidates, hemA and hemL were overexpressed in Corynebacterium glutamicum to produce heme, precursor of biliverdin. To increase precursor supply, we suggested a novel hemQ-mediated coproporphyrin dependent pathway rather than noted hemN-mediated protoporphyrin dependent pathway in C. glutamicum. After securing precursors, hmuO was overexpressed to pull the carbon flow to produce biliverdin. Through modular optimization using gene rearrangements of hemA, hemL, hemQ, and hmuO, engineered C. glutamicum BV004 produced 11.38 ± 0.47 mg/L of biliverdin at flask scale. Fed-batch fermentations performed in 5 L bioreactor with minimal medium using glucose as a sole carbon source resulted in the accumulation of 68.74 ± 4.97 mg/L of biliverdin, the highest titer to date to the best of our knowledge.


We developed an eco-friendly microbial cell factory to produce biliverdin using C. glutamicum as a biosystem. Moreover, we suggested that C. glutamicum has the thermodynamically favorable coproporphyrin dependent pathway. This study indicated that C. glutamicum can work as a powerful platform to produce biliverdin as well as heme-related products based on the rational design with in vitro thermodynamic analysis.


Biliverdin, a tetrapyrrolic pigment, is a product of heme catabolism that consists of several isomers such as biliverdin IX alpha (α), beta (β), gamma (γ) and delta (δ) [1]. Especially, biliverdin IXα, which is typically known as biliverdin, has been considered to be a novel recyclable antioxidant for humans [2]. In human bodies, glutathione is recycled in the glutathione cycle to protect water-soluble proteins from reactive oxygen species (ROS), and biliverdin protects lipids from ROS as a component of the biliverdin-bilirubin cycle [3]. Biliverdin is gaining interest in many medical fields including research on antiviral agents [4], inflammation regulators [5], and medicine for lung graft injury [6] or liver ischemia reperfusion injury [7]. Moreover, it has been recognized as an important precursor of various chromophores for phytochromes used in material sciences [8], optogenetics [9] and synthetic biology [10].

Current methods to produce biliverdin are related to the chemical oxidation of bilirubin that was extracted from mammalian bile [11, 12]. These methods raised some problems such as environmental disruption, impurity, and isomer production. In contrast to chemical processing, bio-based production using microbial cell factories has been considered as a rational source for the production of medicines, materials, fuels, and a broad range of chemicals [13]. Biliverdin is biologically produced via the heme biosynthesis pathway. Heme, a direct precursor of biliverdin, is an important material for living organisms; therefore, the heme biosynthesis pathway is active in most organisms. The heme biosynthesis pathway is composed of almost ten enzymes in two forks [14]. The first fork aims to produce 5-aminolevulinic acid from amino acids such as glycine or glutamate. One route of the first fork is the C4 pathway that directly condensates glycine and succinyl-CoA in mammals, birds, and fungi [15]. Another route is the C5 pathway, which starts from glutamate using glutamyl-tRNA synthetase, glutamyl-tRNA reductase, and glutamate-1-semialdehyde aminotransferase in most bacteria, archaea and, plants [16]. The second fork aims to produce heme from coproporphyrinogen III via two different pathways. The protoporphyrin dependent pathway uses protoporphyrin IX as a frame for heme, but the recently discovered coproporphyrin dependent pathway utilizes a coproporphyrin III chassis for heme [17].

A few biological production methods of biliverdin using E. coli have been suggested, but the yield was low despite the constant addition of nitrogen sources such as yeast extract or peptone [18]. Furthermore, there have been some efforts to produce phycobilins using mammalian cell or E. coli but titers were too low due to the low accumulation of biliverdin, the precursor of phycobilins [9, 19, 20]. Thus, it is necessary to develop another strategy for efficient bio-based production of biliverdin. One of the attractive biosystems in synthetic biology is Corynebacterium glutamicum which has been used for the industrial production of amino acids such as glutamate or lysine [21]. It has already been recognized as an efficient biosystem to produce various compounds such as chemicals, fuels, and materials for healthcare [22]. Even though C. glutamicum had the limelight for production of 5-aminolevulinic acid which is intermediate in biliverdin biosynthesis pathways, studies on the production of other tetrapyrroles were limited [23]. There are some advantages to use C. glutamicum which is known to have the C5 pathway and protoporphyrin-dependent pathway for biliverdin production. It naturally produces large amounts of glutamate, the first intermediate of the C5 pathway. Moreover, C. glutamicum contains native heme oxygenase which can produce biliverdin from heme [24]. However, there are no enzymes that can produce isomers of biliverdin such that additional separation methods are unnecessary, unlike current approaches [11, 12]. Additionally, it is generally regarded as a safe (GRAS) gram-positive bacteria which is appropriate for producing pharmaceutical products such as biliverdin.

Recently, the rational design of microbial cell factories has become increasingly important due to the abstruseness of the control of carbon flow in microorganisms [25]. Push and pull engineering is one of the methods for the reasonable design of microbial cell factories [26]. Push engineering overexpresses enzymes that act as a bottleneck for pushing carbon flow to produce plenty of precursors, and pull engineering aims to pull carbon flow to targeted final products by strengthening terminal enzymes. This method is commonly based on previous studies about each enzyme related to the pathway, but in vitro thermodynamic analysis has considerable potential as an alternative way to find bottlenecks in the pathway [27, 28]. In vitro thermodynamic analysis calculates the Gibbs free energy of each enzyme in the path and represents hypothetical rate-limiting steps [29].

In this study, we developed C. glutamicum platform strains for the efficient production of biliverdin. A push and pull strategy with in vitro thermodynamic analysis was used because of a convoluted and unrevealed heme biosynthesis pathway (Fig. 1a). Push and pull modules were constructed after the investigation of candidate genes. Moreover, a combination of modules through modular optimization and the addition of biotin enabled the production of biliverdin in C. glutamicum. Finally, fed-batch fermentations were performed using 5 L bioreactor for the scalable production of biliverdin.

Fig. 1

Schematic representation and thermodynamic data of biosynthetic pathway for production of biliverdin. a Metabolic pathway and main push and pull strategies for biliverdin biosynthesis in C. glutamicum. The red arrows represent the push steps and the blue arrows represent the pull steps. b Comparing the changes in Gibbs free energy of reaction as a function of pathway progress from glutamate to biliverdin. Closed circle, protoporphyrin dependent pathway; Open circle, coproporphyrin dependent pathway. Abbreviations: UPG, uroporphyrinogen III; UP, uroporphyrin III; CPG, coproporphyrinogen III; CP, coproporphyrin III, Fe-CP, Coproheme III; PPG, protoporphyrinogen IX; PP, protoporphyrin IX; GltX, glutamyl-tRNA synthetase; HemA, glutamyl-tRNA reductase; HemL, glutamate-1-semialdehyde 2,1-aminomutase; HemB, porphobilinogen synthase; HemC, hydroxymethylbilane synthase; HemD, uroporphyrinogen III synthase; HemE, uroporphyrinogen decarboxylase; HemY, prtoporphyrinogen/coproporphyrinogen III oxidase; HemH, protoporphyrin/coproporphyrin ferrochelatase; HemQ, Fe-coproporphyrin III decarboxylase; HemN, coproporphyrinogen III oxidase; HmuO, heme oxygenase


Identification of the genes for a push module using in vitro thermodynamic analysis

The heme biosynthesis pathway is known as a tightly regulated pathway and its regulation mechanism is still ambiguous [30]. Moreover, it is hard to accurately quantify the intermediates of this pathway because porphyrinogens easily change to porphyrins due to their ability to autoxidize [31]. To efficiently identify the rate-limiting steps, levels of change in Gibbs free energies in the heme biosynthesis pathway were calculated (Fig. 1b and Table 1). Thermodynamically favorable reactions are exergonic and have negative ΔG0′. On the contrary, reactions that have a positive ΔG0′ are endergonic and thermodynamically unfavorable. The calculations showed that only hemA (glutamyl-tRNA reductase) and hemH (protoporphyrin ferrochelatase) had positive ΔG0′ in the pathway, and the values were 372.69 kJ/mol and 244.49 kJ/mol, respectively. hemL (glutamate-1-semialdehyde aminotransferase) had nearly zero ΔG0′ even though other enzymes had highly negative values. The overall pathway of C. glutamicum that is known to have a protoporphyrin-dependent pathway using hemN (coproporphyrinogen III oxidase) can be written as the reaction below [24]:

Table 1 Calculated ΔG0′ values for the reaction to produce biliverdin starting from l-glutamate

8 l-Glutamate + 2 AdoMet + 8 ATP + 8 NADPH + 6 Fdred + 12 H+ + 6 O2 → Biliverdin + 2 l-Methionine + 2 5'-dAdo + 8 AMP + 8 PPi + 8 NADP+ + 6 Fdox + 3 H2O2 + 11 H2O + 4 NH3 + CO2 + CO.

The calculated ΔG0′ of the overall reaction was − 2661.39 kJ/mol and it was expected to be a highly thermodynamically favorable pathway when enzymes related to bottlenecks were overexpressed (Fig. 1b). Among candidates HemA, HemL, and HemH analyzed through in vitro thermodynamic analysis, HemA and HemL are used eight times to produce biliverdin but HemH used once (Table 1). Based on these data, hemA from Salmonella typhimurium and hemL from Escherichia coli were selected as rational components for push modules. Especially, hemAM (mutated hemA) that contained additional lysine near the N-terminus was used for the prevention of the degradation of hemA by heme [32]. To identify the effect of overexpression of hemAM and hemL, HM002 strain that contains pEKEx2-hemAL and pMTZ vectors in C. glutamicum ATCC 13826, known as a glutamate overproducing strain, was constructed. The HM002 strain showed similar optical density at 60 h, but the specific growth rate was decreased to 0.14 ± 0.01 h− 1 compared with 0.24 ± 0.03 h− 1 of control strain, HM001 (Fig. 2a-b). Even though HM001 produced little or no porphyrins, HM002 produced 6.29 ± 0.78 mg/L and 2.91 ± 0.11 mg/L of uroporphyrin III (UP) and coproporphyrin III (CP), respectively (Fig. 2d-e). Moreover, the HM002 strain produced 15.59 ± 0.79 mg/L of heme (HM), a 5.48-fold increase compared with control strain HM001. According to the thermodynamic data, we anticipated accumulation of protoporphyrin IX because of the highly positive ΔG0′ of hemH, which converts protoporphyrin IX to heme (Table 1). However, protoporphyrinogen IX and protoporphyrin IX were not detected.

Fig. 2

Effect of hemA, hemL, and hemN overexpression as push parts for heme production in C. glutamicum. a, b, and c Bacterial cell growth and glucose consumption rate of HM001, HM002, and HM003, respectively. d, e, and f Porphyrin intermediate analysis at 60 h cultivation of HM001, HM002, and HM003, respectively. Closed circles, bacterial cell growth; open circle, residual glucose concentration

Therefore, we speculated that overexpression of hemN was required to strengthen carbon metabolic flow to protoporphyrin IX for heme production rather than coproporphyrin III, which might not act as a precursor of heme in C. glutamicum. Since the amino acid sequences of the hemN from C. glutamicum ATCC 13826 are the same as that from C. glutamicum ATCC 13032, we used the hemN from C. glutamicum ATCC 13032. As another part for the push module, hemN was cloned into a pMTZ vector and expressed in C. glutamicum ATCC 13826 with hemAM and hemL. This HM003 strain showed a slightly increased specific growth rate (0.16 ± 0.01 h− 1) and optical density at 60 h (54.03 ± 2.73) compared with HM002 (Fig. 2c). However, it could still not produce any protoporphyrinogen IX or protoporphyrin IX. In addition, even heme, uroporphyrin III, and coproporphyrin III were decreased to 11.81 ± 0.83 mg/L, 2.08 ± 0.49 mg/L and 1.30 ± 0.02 mg/L, respectively (Fig. 2f).

Discovery of hemQ as a part for the unrevealed coproporphyrin dependent pathway in C. glutamicum

It is generally known that C. glutamicum has a protoporphyrin dependent pathway using hemN, hemY (protoporphyrinogen III oxidase), and hemH (Fig. 1a). However, a recent hypothesis that some Actinobacteria and Firmicutes were unable to produce protoporphyrin and used coproporphyrin III as a precursor for heme production has been raised [17]. The coproporphyrin dependent pathway is composed of hemY, hemH, and hemQ (Fe-coproporphyrin III decarboxylase). Differently from their roles in the protoporphyrin-dependent pathway, hemY and hemH act as coproporphyrinogen III oxidase and coproporphyrin ferrochelatase, respectively, in this pathway. The overall heme biosynthesis pathway that uses coproporphyrin can be written as the reaction below:

8 l-Glutamate + 8 ATP + 8 NADPH + 6 Fdred + 12 H+ + 6. O2 Biliverdin + 8 AMP + 8 PPi + 8 NADP+ + 6 Fdox + H2O2 + 15 H2O + 4 NH3 + 6CO2 + CO.

The ΔG0′ of the whole reaction was − 3108.60 kJ/mol and this value was 447.22 kJ/mol less than the value from the protoporphyrin-dependent heme biosynthesis pathway (Fig. 1b and Table 1). Therefore, it was a more thermodynamically favorable reaction and used fewer products such as S-adenosyl-l-methionine. Based on this information, Cgl1899 in C. glutamicum was selected as a putative hemQ gene as part of the pull module using BLAST. Crude enzyme extract assays of HemH and HemQ were performed to identify the functional activity of HemQ using HM004, HM005, and HM006 that harboring pMTZ, pMTZ-hemH, and pMTZ-hemQNAT, respectively (Fig. 3). When HemH and HemQ of C. glutamicum were assayed with coproporphyrin III, 8.55 ± 0.45 μM of heme, 3.08 times higher than no enzyme expression, was accumulated, confirming that HemQ (Cgl1899) has functional activity. Furthermore, crude extract assays with protoporphyrin IX were performed to compare the protoporphyrin dependent pathway and coproporphyrin dependent pathway. When HemH alone or HemH and HemQ were assayed with protoporphyrin IX, the amount of heme was not significantly increased compared with the results using coproporphyrin III as a substrate. Based on the results, we confirmed the possibility of hemQ-mediated coproporphyrin dependent pathway in C. glutamicum and hemQ was selected as a pull part to produce biliverdin.

Fig. 3

In vitro heme synthesis from protoporphyrin IX or coproporphyrin III using HemH and HemQ from C. glutamicum. Plus symbols represent the addition of enzymes (HemH or HemQ) or substrates (protoporphyrin IX or coproporphyrin III)

To identify the effect of combination of push and pull parts, pMTZ-hemQNAT carrying the native hemQ (Cgl1899) gene that uses GTG as a start codon was constructed and cotransformed with pEKEx2-hemAL (named HM007). The HM007 showed a similar growth pattern but the specific growth rate was increased to 0.18 ± 0.01 h− 1 compared with HM002 (Fig. 4a). As the hemQ gene was added, heme production was increased to 17.66 ± 0.42 mg/L and coproporphyrin III was decreased to 1.14 ± 0.01 mg/L in the HM007 strain (Fig. 4d). To increase HemQ activity, pMTZ-hemQ carrying Cgl1899 with a change of the GTG start codon to ATG was constructed (named HM008). The HM008 represented a similar growth pattern but produced 27.22 ± 0.65 mg/L of heme, 9.56-fold and 1.74-fold increase compared with HM001 and HM002, respectively (Fig. 4b and e). It produced 1.42 ± 0.07 mg/L of coproporphyrin III and 6.78 ± 0.54 mg/L of uroporphyrin III. However, biliverdin was not detected even though precursors were supplied.

Fig. 4

Effect of hemQNAT, hemQ, and hmuO overexpression as pull parts for biliverdin production in C. glutamicum. a, b, and c Bacterial cell growth and glucose consumption rate of HM007, HM008, and BV001, respectively. d, e, and f Porphyrin intermediate analysis at 60 h cultivation of HM007, HM008, and BV001, respectively. Closed circles, bacterial cell growth; open circle, residual glucose concentration

Enhanced metabolic flux for biliverdin biosynthesis using hmuO as a part of the pull module

Because all recombinant strains only produced negligible amounts of biliverdin, the hmuO gene from C. glutamicum was selected as another part of the pull module and cloned into pMTZ, generating pMTZ-hmuO. The BV001 strain was developed by coexpressing pEKEx2-hemAL and pMTZ-hmuO. During cultivation, the cell growth pattern of the BV001 strain was similar to that of other strains (Fig. 4c). However, it consumed all of the glucose at 28 h while other strains did this at 24 h. During 60 h batch cultivation in a 500 mL baffled flask, 7.37 ± 0.40 mg/L of biliverdin was produced (Fig. 4f). Compared to HM002 that had only a push module, heme which is a direct precursor of biliverdin was decreased to 7.93 ± 0.29 mg/L. Decreased amounts of uroporphyrin III (4.25 ± 0.62 mg/L) and coproporphyrin III (3.27 ± 0.37 mg/L) were also detected. These results showed that the BV001 strain has the potential for production of biliverdin, but more genetic engineering was needed to supply precursors.

On the basis of the results, we speculated that the utilization of both hemQ and hmuO as part of the pull module could produce biliverdin more efficiently (Fig. 5a). BV002 harboring pEKEx2-hemALQ and pMTZ-hmuO was constructed and produced 9.39 ± 0.22 mg/L of biliverdin, which was 1.27-fold higher than the BV001 strain (Fig. 5b). Furthermore, another three strains that had both push parts (hemAM and hemL) and pull parts (hemQ and hmuO) with different genetic arrangements were developed. The BV004 strain harboring pEKEx2-hemAL and pMTZ-hemQO gave better production of biliverdin than the other strains. At 60 h, the BV004 strain produced 11.38 ± 0.47 mg/L of biliverdin, while 9.14 ± 1.79 mg/L of heme remained (Fig. 5b). Meanwhile, 10.93 ± 0.23 mg/L and 6.00 ± 0.27 mg/L of biliverdin were produced in BV003 and BV005, respectively, while 11.04 ± 0.16 mg/L and 15.09 ± 1.32 mg/L of heme had accumulated, respectively. As shown in Additional file 1: Figure S1, four different modules produced different concentrations of uroporphyrin III and coproporphyrin III, but protoporphyrin IX was still not detected.

Fig. 5

Modular optimization of biliverdin production by combining push module and pull module. a Schematic representation of different combinations. The genes hemAM, hemL, hemQ, hmuO are assembled and rearranged with the plasmids that have different origins. b Biliverdin and heme production of recombinant strains with different gene combinations. White box, concentration of heme; black box, concentration of biliverdin

In addition, the biotin concentration was optimized in BV004 culture for enhanced biliverdin production. Biotin limitation is one of the storage tactics in C. glutamicum for the overproduction of glutamate which is an important precursor for biliverdin production. Nevertheless, the concentration of biotin should be optimized because a low level of biotin could inhibit the growth of C. glutamicum which is a biotin auxotroph [33]. For the optimization of biotin concentrations, C. glutamicum BV004 was cultivated in modified CGX medium with eight different biotin concentrations: 1 μg/L, 5 μg/L, 20 μg/L, 50 μg/L, 100 μg/L, 200 μg/L, 300 μg/L and 500 μg/L (Table 2). BV004 had similar growth patterns when cultivated with biotin concentrations from 20 μg/L to 200 μg/L, and the final OD600 reached almost 50. A total of 1 μg/L and 5 μg/L of biotin had a lower final OD600 (40.12 ± 0.22 and 49.80 ± 0.63), but 300 μg/L and 500 μg/L of biotin had a final OD600 up to almost 59. Among the different concentrations of biotin, BV004 cultivated with 200 μg/L produced the highest biliverdin titer, which was 11.38 ± 0.47 mg/L.

Table 2 Biotin optimization for biliverdin production in C. glutamicum BV004

Biliverdin production by C. glutamicum BV004 in fed-batch fermentation

To evaluate the ability of the BV004 strain that produced the highest concentration of biliverdin among all strains, scaled-up fed-batch cultures were carried out in 5 L bioreactor to produce biliverdin (Fig. 6 and Additional file 1: Figure S2). Based on the optimization of the biotin concentration, 200 μg/L of biotin was used in fed-batch fermentation (Table 2). Starting with 80 g/L glucose as the sole carbon source, total 183.02 ± 7.53 g of glucose was consumed. During fed-batch fermentations, 68.74 ± 5.00 mg/L of biliverdin was produced, which is the highest titer using microbial cell factories to the best of our knowledge. The yield and productivity were 0.84 ± 0.07 mg/g glucose and 0.95 ± 0.07 mg/L/h.

Fig. 6

Averaged fed-batch fermentation profile of C. glutamicum BV004. The error bars represent the means and standard deviations of triplicate fermentation profiles. Closed circles, bacterial cell growth; open circle, residual glucose concentration; green square, biliverdin concentration


Biliverdin has attracted increasing interest due to its wide application in various fields. In particular, it has become apparent that biliverdin is one of the recyclable antioxidants such as glutathione and it is expected that a lower amount of biliverdin was needed for reducing ROS compared to other antioxidants because of its recyclability [2, 3]. In this study, we developed recombinant C. glutamicum strains for the efficient production of biliverdin via novel coproporphyrin dependent pathway. Even though C. glutamicum is known to have heme oxygenase, biliverdin was not detected during the culture of the wild type strain. Thus, we introduced a rational push and pull strategy with in vitro thermodynamic analysis to produce biliverdin. The strategies used in this study were as follows: (i) overexpression of hemAM and hemL as push parts based on thermodynamic data; (ii) discovery of unnoted hemQ as a pull part after comparison of known hemN and hemQ; (iii) overexpression of hmuO as a pull part; (iv) modular optimization of push and pull modules; and (v) fed-batch fermentation to a scalable production of biliverdin using a proper biotin concentration.

The push and pull strategy is one of the logical tactics to produce targeted materials and it has been used in many biosystems such as Yarrowia lipolytica [26], Saccharomyces cerevisiae [34] and E. coli [35]. Briefly, ‘push’ means the overexpression of enzymes that are recognized as bottlenecks in pathways and ‘pull’ means the overexpression of terminal enzymes. This strategy minimizes target genes overexpression for the efficient production of materials. Overexpression of essential enzymes, rather than all of the related enzymes, avoids metabolic burdens and excessive use of antibiotics [36]. Moreover, it is more important in the heme biosynthesis pathway because the enzymes related to this intricately regulated pathway influence the mRNA expression levels of one another [30]. In the case of E. coli, the overexpression of hemD upregulated almost all of the enzymes in the heme biosynthesis pathway. Conversely, the overexpression of enzymes such as hemB or hemG reduced the relative mRNA expression levels of other enzymes. Candidates for the components of the push and pull strategy are commonly selected based on many studies on the related enzymes but there were a few articles on the biliverdin biosynthesis pathway in C. glutamicum. Therefore, we primally introduced in vitro thermodynamic analysis of all enzymes that were related to the pathway [37, 38]. This analysis provided information on the thermodynamic stability of enzymes using the calculation of Gibbs free energy. Enzymes with higher Gibbs free energies are thermodynamically unfavorable and have greater possibilities to act as rate-limiting steps in pathways [39, 40]. The success of the overexpression of genes based on the calculations in this study corresponded with previous studies that overexpressed hemA and hemL or hemH for heme production [41,42,43]. Therefore, hemH can be a tempting candidate for further studies to enhance heme production in C. glutamicum.

The result that the overexpression of hemN negatively effects heme production is noteworthy (Fig. 2f). We speculated on the existence of a hemQ-mediated pathway because protoporphyrin IX is not detected in C. glutamicum despite hemN overexpression. The result that C. glutamicum does not use protoporphyrin IX but instead uses coproporphyrin III for heme makes C. glutamicum a reasonable microbial cell factory to produce materials related to heme. There are some different enzymes related to the protoporphyrin-dependent pathway that use different cofactors [14]. For example, in the case of E. coli, hemF and hemG can replace hemN and hemY, respectively (Additional file 1: Table S1). There are five methods for producing heme from coproporphyrinogen III. Based on the calculation of the Gibbs free energies, the reaction using hemY, hemH, and hemQ, as with C. glutamicum, has the greatest potential among the tested reactions. The ΔG0′ of this reaction was − 940.33 kJ/mol, and this value was at least 190.18 kJ/mol less than the other reactions. Moreover, it utilizes the least number of substrates for heme production. It only uses 1 mol coproporphyrinogen III, 3 mol oxygen, and 1 mol ferrous ion to synthesize 1 mol of heme, even though other methods need additional cofactors such as 5′-deoxyadenosines or menaquinones.

Balancing the expression of gene rearrangements may be an important factor in producing materials in synthetic biology [44]. The utilization of substrates or the production of targeted chemicals can be diverse even though the same genes were overexpressed [45, 46]. A pEKEx2 vector that contains a pBL1 origin and a pMTZ vector that carries a pCG1 origin produce moderate copy numbers of plasmids, but a pBL1 origin produces a slightly higher copy number per cell than a pCG1 origin [47]. Four different strains were constructed with different gene arrangements using the pEKEx2 and pMTZ vectors (Fig. 5a). All strains, except for BV005, produced more biliverdin than the BV001 strain. Interestingly, BV005 showed similar production of biliverdin compared to the BV001 strain. It seems that the order of genes is one of the important factors because of various translation levels in a synthetic operon [48]. Modular optimization may be important in the heme biosynthesis pathway because it is tightly regulated and related genes affect each other. Despite using the same genes, such as hemAM, hemL, hemD, and hemF, 5-aminolevulinic acid production was 180 times different depending on modular optimization in E. coli [30]. Apart from the biliverdin concentration, uroporphyrin III also seemed to be influenced by gene arrangements. In the heme biosynthesis pathway, uroporphyrinogen III is converted to coproporphyrinogen III by hemE or is spontaneously changed to uroporphyrin III without an enzymatic reaction. The amount of accumulated uroporphyrin III is highly affected by the overexpression of genes related to the heme biosynthesis pathway in E. coli [49]. When hemB, hemC, and hemD genes were overexpressed in E. coli, uroporphyrin III was increased to 4.67-, 2.17-, and 17.5-fold compared to the control strains, respectively. Likewise, the overexpression of hemQ with hmuO using a moderate copy number vector seems to highly affect the concentration of uroporphyrin III. BV004 and BV005 in which hemQ and hmuO were overexpressed with the pMTZ vector accumulated almost 20 mg/L of uroporphyrin III, while other strains produced nearly 5 mg/L of uroporphyrin III (Additional file 1: Figure S1). Therefore, hemE that strengthens the carbon flow to coproporphyrinogen III and reduces accumulation of uroporphyrin III can be a reasonable candidate for further studies.


In summary, we suggested the possibility of C. glutamicum as a reasonable biosystem to produce biliverdin. We demonstrated the existence of the coproporphyrin dependent pathway in C. glutamicum and produced the highest titer of biliverdin using this pathway. Based on this study, further improvements in upstream bioprocessing, such as rerouting the carbon flux of the TCA cycle and introducing a novel pathway to produce more precursors, and downstream bioprocessing, such as developing an efficient purification method, will be necessary for the industrial utilization of biliverdin.


Strains, plasmids, and primers

All bacterial strains and plasmids used in this study are listed in Table 3. All primers used in this study are presented in Table 4. E. coli DH5α was used as the host for gene cloning and plasmid construction, and C. glutamicum ATCC 13826 was used as the host strain to produce biliverdin. For the two-vector system of C. glutamicum, pMTZ vector was constructed by replacing the kanamycin resistance gene of pMT-tac with a zeocin resistance gene from pPsADHα by digestion with SalI and SnaBI. The hemN, hemH, hemQNAT, hemQ, and hmuO genes were amplified from the genomic DNA of C. glutamicum. To construct plasmids pMTZ-hemN, pMTZ-hemH, pMTZ-hemQNAT, pMTZ-hemQ, and pMTZ-hmuO, the hemN, hemH, hemQNAT, hemQ, and hmuO genes were amplified using primers N-F/N-R, H-F/H-R, QNAT-F/QNAT-R, Q-F/Q-R, and O-F/O-R, respectively. The amplified genes were digested with ClaI/NotI (hemN, hemH, hemQNAT, and hemQ) and ClaI/BamHI (hmuO) and cloned into plasmid pMTZ. To construct plasmids pEKEx2-hemALQ and pEKEx2-hemALO, hemQ and hmuO genes were amplified using primers ALQ-F/ALQ-R and ALO-F/ALO-R, respectively. The amplified genes were digested with BamHI/KpnI and cloned into plasmid pEKEx2-hemAL. To construct plasmids pMTZ-hemQO and pMTZ-hemOQ, hmuO and hemQ genes were amplified using primers QO-F/QO-R and OQ-F/OQ-R, respectively. The amplified genes were digested with NotI/NotI (hmuO) and BamHI/NotI (hemQ) and cloned into plasmid pMTZ-hemQ and pMTZ-hmuO, respectively.

Table 3 Strains and plasmids in this study
Table 4 Primers used in this study

Thermodynamic analysis

The Gibbs free energies for the compounds involved in the pathway to produce biliverdin were obtained from MetaCyC [53]. Net changes in the Gibbs free energy for each reaction were calculated using the equation ΔG0′ = Σ(G0products) − Σ(G0substrates). The values given in kcal/mol were converted to kJ/mol (1 kcal/mol = 4.184 kJ/mol).

Media composition

LB medium (10 g/L NaCl, 10 g/L tryptone, and 5 g/L yeast extract) was used for molecular genetic procedures. BHIS medium (37 g/L brain heart infusion and 91 g/L sorbitol) was used for the preculture of C. glutamicum. Modified CGXII medium consisting of 42 g 3-morpholinopropane-1-sulfonic acid (MOPS), 20 g ammonium sulfate, 5 g urea, 1 g potassium dihydrogen phosphate, 1 g potassium phosphate dibasic anhydrous, 1 g calcium chloride anhydrous, 0.25 g MgSO4·7H2O, 10 mg CaCl2, 10 mg FeSO4·7H2O, 1 mg ZnSO4·7H2O, 0.2 mg biotin, 0.31 mg CuSO4·5H2O, 0.1 mg MnSO4·H2O and 0.02 mg NiCl2·6H2O per liter was used for the culture of C. glutamicum with glucose as a sole carbon source. When appropriate, kanamycin (50 μg/mL for E. coli or 25 μg/mL for C. glutamicum), ampicillin (100 μg/mL) or zeocin (25 μg/mL) was added to the medium.

Flask cultivation

Cells were inoculated into 100 mL baffled flasks containing 20 mL BHIS medium and were cultivated at 30 °C for 12 h with 150 rpm. The precultures were transferred to 500 mL baffled flasks containing 100 mL modified CGXII medium with 40 g/L glucose and 1 mM Isopropyl β-d-1-thiogalactopyranoside (IPTG) at an initial OD600 of 1. Then, cells were cultivated at 30 °C for 60 h with 150 rpm.

Fed-batch cultivation

Cells were inoculated into 100 mL baffled flasks containing 20 mL BHIS medium and were cultivated at 30 °C for 12 h with 150 rpm. The first precultures were transferred to 1 L baffled flasks containing 200 mL BHIS medium and were cultivated at 30 °C for 12 h with 150 rpm. The second precultures were washed once with modified CGXII medium and were transferred to 5 L bioreactor containing 1.8 L modified CGXII medium with 80 g/L glucose and 1 mM IPTG. Temperature, agitation, and aeration rate were maintained at 30 °C, 600 rpm, and 2 vvm, respectively. The pH was maintained at 7.0 by automatically adding 10% H3PO4 and 4 N KOH. When the glucose concentration decreased to approximately 20 g/L, 50% glucose was added to maintain the level of glucose at approximately 20 g/L. Foam was removed using 10% antifoam 204 (Sigma-Aldrich, USA).

Analytical procedures

Bacterial cell growth was followed by measuring the optical density at 600 nm using a UV/Vis spectrophotometer (Mecasys Co., Ltd., Korea). The glucose concentration in the culture medium was measured with glucose assay kit (Sigma-Aldrich, USA). Porphyrin intermediates and heme concentrations were measured with a high-performance liquid chromatography (HPLC) system (Waters Corporation, USA) with a UV detector (Waters 2487) at 400 nm [41]. The initial mobile phase composition was 20% solvent A (1:9 methanol:acetonitrile) and 80% solvent B (0.5% trifluoroacetic acid in water, pH 2.60). The solvent gradient consisted in 40 min linear change at 95% solvent A and 5% solvent B with 1.0 ml/min flow velocity at 40 °C. Biliverdin concentration was analyzed by HPLC system with a UV detector at 376 nm [54]. The initial mobile phase composition was 25% solvent C (methanol) and 75% (1 M ammonium acetate, pH 5.16). The solvent gradient consisted of 8 min linear change at 95% solvent C and 5% solvent D, and then 2 min at these conditions and 8 min at 25% solvent C and 75% solvent D with 1.5 ml/min flow velocity at 70 °C.

Enzymatic assays

Cells were cultivated at 30 °C for 16 h with 150 rpm, harvested by centrifugation (4000 rpm at 4 °C for 20 min), washed twice in 66 mM Tris-HCl (pH 8.0). The cells were resuspended in same buffer and sonicated for 20 min. Cell debris was removed by centrifugation (4000 rpm at 4 °C for 20 min) and the supernatant was used to crude extract enzyme assays. The protein concentration was measured by Bradford method [55]. The coupled assays of HemH and HemQ were performed as previously reported [56]. The reactions were performed at 30 °C for 15 min in 66 mM Tris-HCl (pH 8.0) containing 3.3% Tween 20, 5 mM glutathione, 100 μM 2-mercaptoethanol, 100 μM ferrous ammonium sulfate, 25 μM substrate (protoporphyrin IX or coproporphyrin III), and 1 mg/mL crude enzyme extracts.









Coproporphyrin III


Coproporphyrinogen III

Fdox :

Oxidized ferredoxin

Fdred :

Reduced ferredoxin


Coproheme III




Protoporphyrin IX


Protoporphyrinogen IX


Uroporphyrin III


Uroporphyrinogen III


  1. 1.

    O’Carra P, Colleran E. Separation and identification of biliverdin isomers and isomer analysis of phycobilins and bilirubin. J Chromatogr A. 1970;50:458–68.

    Article  Google Scholar 

  2. 2.

    Ziberna L, Martelanc M, Franko M, Passamonti S. Bilirubin is an endogenous antioxidant in human vascular endothelial cells. Sci Rep. 2016;6:1–6.

    Article  Google Scholar 

  3. 3.

    Sedlak TW, Saleh M, Higginson DS, Paul BD, Juluri KR, Snyder SH. Bilirubin and glutathione have complementary antioxidant and cytoprotective roles. Proc Natl Acad Sci. 2009;106:5171–6.

    Article  Google Scholar 

  4. 4.

    Zhang A, Duan H, Li N, Zhao L, Pu F, Huang B, et al. Heme oxygenase-1 metabolite biliverdin, not iron, inhibits porcine reproductive and respiratory syndrome virus replication. Free Radic Biol Med. 2017;102:149–61.

    Article  Google Scholar 

  5. 5.

    Wegiel B, Gallo D, Csizmadia E, Roger T, Kaczmarek E, Harris C, et al. Biliverdin inhibits toll-like receptor-4 (TLR4) expression through nitric oxide-dependent nuclear translocation of biliverdin reductase. Proc Natl Acad Sci. 2011;108:18849–54.

    Article  Google Scholar 

  6. 6.

    Zhou H, Qian H, Liu J, Zhu D, Ding W, Pan P, et al. Protection against lung graft injury from brain-dead donors with carbon monoxide, biliverdin, or both. J Hear Lung Transplant. 2011;30:460–6.

    Article  Google Scholar 

  7. 7.

    Andria B, Bracco A, Attanasio C, Castaldo S, Cerrito MG, Cozzolino S, et al. Biliverdin protects against liver ischemia reperfusion injury in swine. PLoS One. 2013;8:1–8.

    Article  Google Scholar 

  8. 8.

    Watermann T, Elgabarty H, Sebastiani D. Phycocyanobilin in solution-a solvent triggered molecular switch. Phys Chem Chem Phys. 2014;16:6146–52.

    Article  Google Scholar 

  9. 9.

    Uda Y, Goto Y, Oda S, Kohchi T, Matsuda M, Aoki K. Efficient synthesis of phycocyanobilin in mammalian cells for optogenetic control of cell signaling. Proc Natl Acad Sci. 2017;114:11962–7.

    Article  Google Scholar 

  10. 10.

    Fernandez-Rodriguez J, Moser F, Song M, Voigt CA. Engineering RGB color vision into Escherichia coli. Nat Chem Biol. 2017;13:706–8.

    Article  Google Scholar 

  11. 11.

    Lin GL, Himes JA, Cornelius CE. Bilirubin and biliverdin excretion by the chicken. Am J Phys. 1974;226:881–5.

    Article  Google Scholar 

  12. 12.

    Ding ZK, Xu YQ. Purification and characterization of biliverdin IXalpha from Atlantic salmon (Salmo salar) bile. Biokhimiya/Biochemistry. 2002;67:927–32.

    Google Scholar 

  13. 13.

    Jullesson D, David F, Pfleger B, Nielsen J. Impact of synthetic biology and metabolic engineering on industrial production of fine chemicals. Biotechnol Adv. 2015;33:1395–402.

    Article  Google Scholar 

  14. 14.

    Layer G, Reichelt J, Jahn D, Heinz DW. Structure and function of enzymes in heme biosynthesis. Protein Sci. 2010;19:1137–61.

    Article  Google Scholar 

  15. 15.

    Burnham BF. δ-Aminolevulinic acid synthase (Rhodopseudomonas sphaeroides). Methods Enzymol. 1970;17A:195–204.

    Article  Google Scholar 

  16. 16.

    Woodard SI, Dailey HA. Regulation of Heme biosynthesis in Escherichia coli. Arch Biochem Biophys. 1995;316:110–5.

  17. 17.

    Dailey HA, Gerdes S, Dailey TA, Burch JS, Phillips JD. Noncanonical coproporphyrin-dependent bacterial heme biosynthesis pathway that does not use protoporphyrin. Proc Natl Acad Sci U S A. 2015;112:2210–5.

    Article  Google Scholar 

  18. 18.

    Chen D, Brown JD, Kawasaki Y, Bommer J, Takemoto JY. Scalable production of biliverdin IX alpha by Escherichia coli. BMC Biotechnol. 2012;12:89.

    Article  Google Scholar 

  19. 19.

    Stiefelmaier J, Ledermann B, Sorg M, Banek A, Geib D, Ulber R, et al. Pink bacteria—production of the pink chromophore phycoerythrobilin with Escherichia coli. J Biotechnol. 2018;274:47–53.

    Article  Google Scholar 

  20. 20.

    Ge B, Chen Y, Yu Q, Lin X, Li J, Qin S. Regulation of the heme biosynthetic pathway for combinational biosynthesis of phycocyanobilin in Escherichia coli. Process Biochem. 2018;71:23–30.

    Article  Google Scholar 

  21. 21.

    Lee JH, Wendisch VF. Production of amino acids – genetic and metabolic engineering approaches. Bioresour Technol. 2017;245:1575–87.

    Article  Google Scholar 

  22. 22.

    Becker J, Rohles CM, Wittmann C. Metabolically engineered Corynebacterium glutamicum for bio-based production of chemicals, fuels, materials, and healthcare products. Metab Eng. 2018;50:122–41.

  23. 23.

    Ramzi AB, Hyeon JE, Kim SW, Park C, Han SO. 5-Aminolevulinic acid production in engineered Corynebacterium glutamicum via C5 biosynthesis pathway. Enzym Microb Technol. 2015;81:1–7.

    Article  Google Scholar 

  24. 24.

    Yu X, Jin H, Cheng X, Wang Q, Qi Q. Transcriptomic analysis for elucidating the physiological effects of 5-aminolevulinic acid accumulation on Corynebacterium glutamicum. Microbiol Res. 2016;192:292–9.

    Article  Google Scholar 

  25. 25.

    Fisher AK, Freedman BG, Bevan DR, Senger RS. A review of metabolic and enzymatic engineering strategies for designing and optimizing performance of microbial cell factories. Comput Struct Biotechnol J. 2014;11:91–9.

    Article  Google Scholar 

  26. 26.

    Tai M, Stephanopoulos G. Engineering the push and pull of lipid biosynthesis in oleaginous yeast Yarrowia lipolytica for biofuel production. Metab Eng. 2013;15:1–9.

    Article  Google Scholar 

  27. 27.

    Ataman M, Hatzimanikatis V. Heading in the right direction: thermodynamics-based network analysis and pathway engineering. Curr Opin Biotechnol. 2015;36:176–82.

    Article  Google Scholar 

  28. 28.

    Shi T, Han P, You C, Zhang YHPJ. An in vitro synthetic biology platform for emerging industrial biomanufacturing: bottom-up pathway design. Synth Syst Biotechnol. 2018;3:186–95.

    Article  Google Scholar 

  29. 29.

    Kallscheuer N, Vogt M, Marienhagen J. A novel synthetic pathway enables microbial production of polyphenols independent from the endogenous aromatic amino acid metabolism. ACS Synth Biol. 2017;6:410–5.

    Article  Google Scholar 

  30. 30.

    Zhang J, Kang Z, Chen J, Du G. Optimization of the heme biosynthesis pathway for the production of 5-aminolevulinic acid in Escherichia coli. Sci Rep. 2015;5:8584.

    Article  Google Scholar 

  31. 31.

    Mauzerall D. The porphyrinogens. In: Dolphin D, editors.The porphyrins V2: Structure and Synthesis, Part B. New York: Academic Press; 1978. p. 96–8.

  32. 32.

    Wang L, Wilson S, Elliott T. A mutant HemA protein with positive charge close to the N terminus is stabilized against heme-regulated proteolysis in Salmonella typhimurium. J Bacteriol. 1999;181:6033–41.

    Google Scholar 

  33. 33.

    Sato H, Orishimo K, Shirai T, Hirasawa T, Nagahisa K, Shimizu H, et al. Distinct roles of two anaplerotic pathways in glutamate production induced by biotin limitation in Corynebacterium glutamicum. J Biosci Bioeng. 2008;106:51–8.

    Article  Google Scholar 

  34. 34.

    Lv X, Xie W, Lu W, Guo F, Gu J, Yu H, et al. Enhanced isoprene biosynthesis in Saccharomyces cerevisiae by engineering of the native acetyl-CoA and mevalonic acid pathways with a push-pull-restrain strategy. J Biotechnol. 2014;186:128–36.

    Article  Google Scholar 

  35. 35.

    He L, Xiao Y, Gebreselassie N, Zhang F, Antoniewicz MR, Tang YJ, et al. Central metabolic responses to the overproduction of fatty acids in Escherichia coli based on 13C-metabolic flux analysis. Biotechnol Bioeng. 2014;111:575–85.

    Article  Google Scholar 

  36. 36.

    Bentley WE, Mirjalili N, Andersen DC, Davis RH, Kompala DS. Plasmid-encoded protein : the principal factor in the “metabolic burden” associated with recombinant bacteria. Biotechnol Bioeng. 1990;35:668–81.

    Article  Google Scholar 

  37. 37.

    Nagai H, Masuda A, Toya Y, Matsuda F, Shimizu H. Metabolic engineering of mevalonate-producing Escherichia coli strains based on thermodynamic analysis. Metab Eng. 2018;47:1–9.

    Article  Google Scholar 

  38. 38.

    You C, Shi T, Li Y, Han P, Zhou X, Zhang YHP. An in vitro synthetic biology platform for the industrial biomanufacturing of myo-inositol from starch. Biotechnol Bioeng. 2017;114:1855–64.

    Article  Google Scholar 

  39. 39.

    Wang W, Liu M, You C, Li Z, Zhang YHP. ATP-free biosynthesis of a high-energy phosphate metabolite fructose 1,6-diphosphate by in vitro metabolic engineering. Metab Eng. 2017;42:168–74.

    Article  Google Scholar 

  40. 40.

    Thompson B, Pugh S, Machas M, Nielsen DR. Muconic acid production via alternative pathways and an synthetic “metabolic funnel”. ACS Synth Biol. 2017;7:565–75.

    Article  Google Scholar 

  41. 41.

    Ko YJ, Joo Y-C, Hyeon JE, Lee E, Lee M-E, Seok J, et al. Biosynthesis of organic photosensitizer Zn-porphyrin by diphtheria toxin repressor (DtxR)-mediated global upregulation of engineered heme biosynthesis pathway in Corynebacterium glutamicum. Sci Rep. 2018;8:14460.

    Article  Google Scholar 

  42. 42.

    Ramzi AB, Hyeon JE, Han SO. Improved catalytic activities of a dye-decolorizing peroxidase (DyP) by overexpression of ALA and heme biosynthesis genes in Escherichia coli. Process Biochem. 2015;50:1272–6.

    Article  Google Scholar 

  43. 43.

    Franken ACW, Werner ER, Haas H, Lokman BC, Van Den Hondel CAMJJ, Ram AFJ, et al. The role of coproporphyrinogen III oxidase and ferrochelatase genes in heme biosynthesis and regulation in Aspergillus niger. Appl Microbiol Biotechnol. 2013;97:9773–85.

    Article  Google Scholar 

  44. 44.

    Ajikumar PK, Xiao WH, Tyo KEJ, Wang Y, Simeon F, Leonard E, et al. Isoprenoid pathway optimization for Taxol precursor overproduction in Escherichia coli. Science (80-). 2010;330:70–4.

    Article  Google Scholar 

  45. 45.

    Kang MK, Lee J, Um Y, Lee TS, Bott M, Park SJ, et al. Synthetic biology platform of CoryneBrick vectors for gene expression in Corynebacterium glutamicum and its application to xylose utilization. Appl Microbiol Biotechnol. 2014;98:5991–6002.

    Article  Google Scholar 

  46. 46.

    Huang Q, Lin Y, Yan Y. Caffeic acid production enhancement by engineering a phenylalanine over-producing Escherichia coli strain. Biotechnol Bioeng. 2013;110:3188–96.

    Article  Google Scholar 

  47. 47.

    Baritugo K-A, Kim HT, David Y, Choi J, Hong SH, Jeong KJ, et al. Metabolic engineering of Corynebacterium glutamicum for fermentative production of chemicals in biorefinery. Appl Microbiol Biotechnol. 2018;102:3915–37.

    Article  Google Scholar 

  48. 48.

    Hiroe A, Tsuge K, Nomura CT, Itaya M, Tsuge T. Rearrangement of gene order in the phaCAB operon leads to effective production of ultrahigh-molecular-weight poly[(R)-3-hydroxybutyrate] in genetically engineered Escherichia coli. Appl Environ Microbiol. 2012;78:3177–84.

    Article  Google Scholar 

  49. 49.

    Ju LM, Kim HJ, Lee JY, Kwon AS, Jun SY, Kang SH, et al. Effect of gene amplifications in porphyrin pathway on heme biosynthesis in a recombinant Escherichia coli. J Microbiol Biotechnol. 2013;23:668–73.

    Article  Google Scholar 

  50. 50.

    Shin SK, Hyeon JE, Kim YI, Kang DH, Kim SW, Park C, et al. Enhanced hydrolysis of lignocellulosic biomass: Bi-functional enzyme complexes expressed in Pichia pastoris improve bioethanol production from Miscanthus sinensis. Biotechnol J. 2015;10:1912–9.

  51. 51.

    Eikmanns BJ, Kleinertz E, Liehl W, Sahm H. A family of Corynebacterium glutamicum/Escherichia coli shuttle vectors for cloning, controlled gene expression, and promoter probing. Plasmid. 1991;102:93–8.

  52. 52.

    Joo YC, Hyeon JE, Han SO. Metabolic Design of Corynebacterium glutamicum for Production of l -Cysteine with Consideration of Sulfur-Supplemented Animal Feed. J Agric Food Chem. 2017;65:4698–707.

  53. 53.

    Caspi R, Billington R, Fulcher CA, Keseler IM, Kothari A, Krummenacker M, et al. The MetaCyc database of metabolic pathways and enzymes. Nucleic Acids Res. 2018;46:D633–9.

  54. 54.

    Butler MW, Ligon RA. Interactions between biliverdin, oxidative damage, and spleen morphology after simulated aggressive encounters in veiled chameleons. PLoS One. 2015;10:1–16.

  55. 55.

    Bradford MM. A Rapid and Sensitive Method for the Quantification of Microgram Quantities of Protein Utilizing the Principle of Protein-Dye Binding. Anal Biochem. 1976;72:248–54.

  56. 56.

    Dailey TA, Boynton TO, Albetel AN, Gerdes S, Johnson MK, Dailey HA. Discovery and characterization of HemQ: An essential heme biosynthetic pathway component. J Biol Chem. 2010;285:25978–86.

Download references


We would like to express gratitude to School of Life Sciences and Biotechnology for BK21 PLUS, Korea University for processing charges in open access journals.


This work was supported by the National Research Foundation of Korea (N.R.F.) with a grant funded by the Korean government (MSIP) (No. 2018R1A2B2003704).

Availability of data and materials

Gene sequences used in this project are from GenBank ( and the material and data supporting their findings can be found in the main paper and Additional file.

Author information




SOH supervised and coordinated this study. JS and SOH designed, performed most of the experiments and mainly wrote the manuscript. YJK and MEL participated in the research. YJK, MEL, JEH and SOH revised and edited the manuscript. All authors reviewed and approved the final manuscript.

Corresponding author

Correspondence to Sung Ok Han.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional file

Additional file 1:

Figure S1. Production of porphyrin intermediates of recombinant strains with different gene combinations. Black, white and gray bar represent concentration of uroporphyrin III, protoporphyrin IX and coproporphyrin III, respectively. Protoporphyrin IX was not detected in any recombinant strains. Figure S2. Non-averaged fed-batch fermentation profiles of C. glutamicum BV004 related to the Fig. 6. a, b, and c The first, second, and third rounds of fermentation, respectively. Closed circle, bacterial cell growth; open circle, residual glucose concentration; green square, biliverdin concentration. Table S1. Calculated ΔG0′ values for the reaction related to heme biosynthesis pathway. (DOCX 419 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Seok, J., Ko, Y.J., Lee, M. et al. Systems metabolic engineering of Corynebacterium glutamicum for the bioproduction of biliverdin via protoporphyrin independent pathway. J Biol Eng 13, 28 (2019).

Download citation


  • Biliverdin
  • Corynebacterium glutamicum
  • In vitro thermodynamic analysis
  • Coproporphyrin dependent pathway
  • Synthetic biology
  • Metabolic engineering