Enhanced biosynthesis of the pseudouridine precursor in E. coli
The starting substrates and the enzymes used for pseudouridine enzymatic synthesis were existed in E. coli MG1665 [22, 23]. Nevertheless, no pseudouridine was detected by HPLC in the shake flask fermentation broth of E. coli MG1655 (Fig. 2). When phoA and EcpsuG (from MG1665) were overexpressed in MG1665, the strain MG1655(phoAEcpsuG) can produce 29.5 mg/L pseudouridine after flask fermentation (Fig. 2, Additional file: Fig. S1/S2/S3), which suggested that phosphatase and ΨMP glycosidase played a crucial role in pseudouridine biosynthesis.
To increase pseudouridine production, the pseudouridine precursors biosynthesis pathway was engineered to increase metabolite flux to pseudouridine biosynthesis (Fig. 1). thrA and argF were successively knocked out to block the branching pathway [22, 24]. Subsequently the negative regulator pepA [25, 26]was deleted to enhance pyrimidine biosynthesis, which yielded strain ZM3. After that, plasmid pZH2 was transformed into ZM3, yielding strain ZM3(phoAEcpsuG). The pseudouridine yield of ZM3(phoAEcpsuG) increased to 102.7 mg/L (Fig. 2). It was more than threefold higher than that of MG1655(phoAEcpsuG) (Fig. 2), which illustrated that increasing precursors biosynthesis was effective in improving pseudouridine production.
-----------------------------Fig. 2----------------------------
Identification of phosphatase gene for efficient pseudouridine production
The ΨMP synthesis reaction involving ΨMP glycosidase was favoured for pseudouridine production and was promoted by ΨMP dephosphorylation. Therefore, ΨMP glycosidase and phosphatase genes were the focus of this study. First, phosphatase genes were screened to determine the most suitable genes for pseudouridine production. SDT1, YKL033W-A and PHM8 from Saccharomyces cerevisiae have been reported to show dephosphorylation activity of ΨMP [27]. pumD was involved in pseudouridine biosynthesis and was predicted to be a phosphatase gene [28], so pumD was selected to test whether it has ΨMP dephosphorylation activity. Then, the codons of the four abovementioned genes, namely, SDT1, YKL033W-A, PHM8 and pumD, were optimized and overexpressed with EcpsuG in ZM3, which yielded strains ZM3(SDT1EcpsuG), ZM3(YKL033W-AEcpsuG), ZM3(PHM8EcpsuG) and ZM3(pumDEcpsuG), respectively. All four genes led to improved pseudouridine yields after fermentation in shake flasks that were approximately threefold higher than that produced with phoA (Fig. 3A). Among these genes, pumD was the most efficient for pseudouridine production, and the yield of ZM3(pumDEcpsuG) was the highest, reaching 395.7 mg/L (Fig. 3A).
-----------------------------Fig. 3----------------------------
SDS-PAGE (Fig. 3B) analysis of the supernatants showed that all the phosphatases tested were correctly expressed, although there were some variations in expression levels. Enzymatic assays revealed that PhoA had the lowest activity among the five phosphatase (Table 1), which was consistent with fermentation results. PumD had the highest activity and substrate affinity (Table 1). Therefore, pumD was chosen for subsequent experiments.
Table 1
Kinetic parameters* of phosphatase
Enzyme
|
Substrate
|
Km(mM)
|
kcat (S − 1)
|
kcat/Km (mM − 1S − 1)
|
PhoA
|
ΨMP
|
86 ± 24.7
|
3 ± 0.4
|
0.035
|
SDT1
|
ΨMP
|
22.8 ± 3.9
|
7 ± 1.2
|
0.307
|
YKL033W-A
|
ΨMP
|
24.7 ± 4.6
|
6 ± 1.5
|
0.243
|
PHM8
|
ΨMP
|
29.5 ± 3.5
|
9 ± 1.7
|
0.305
|
PumD
|
ΨMP
|
9.8 ± 1.7
|
10 ± 1.4
|
1.02
|
*Data are presented as mean ± s.d. (n = 3).
It is worth noting that all the four aforementioned phosphatases belong to the haloacid dehalogenase protein family [27, 28], which suggests that haloacid dehalogenases are potential sources for screening more efficient phosphatases for use in pseudouridine production. Haloacid dehalogenase family phosphatases detoxify ΨTP and ΨMP in Saccharomyces cerevisiae and humans [27], and genome analysis revealed that Group I C-glycosynthases usually coexist with a haloacid dehalogenase family phosphatase gene located in the vicinity of the C-glycosynthase gene [29], which may support the speculation that haloacid dehalogenases are candidates for screening more efficient phosphatases against ΨMP.
Identification of ΨMP glycosidase gene for efficient pseudouridine production
psuG is widely presented in the genomes of prokaryotes [19]. In this study, EcpsuG and four other psuG genes from different sources were tested. KspsuG (from Klebsiella spallanzanii), RspsuG (from Rhizobium sp. CF142), SspsuG (from Saccharopolyspora spinosa) and SppsuG (from Streptomyces platensis) came from different bacteria belonged to different genera and had different genetic distance to EcpsuG (Fig. 4A). The plasmids carrying the four aforementioned codon-optimized psuG genes co-overexpressed with pumD were transformed into ZM3, yielding strains ZM3(pumDKspsuG), ZM3(pumDRspsuG), ZM3(pumDSspsuG), and ZM3(pumDSppsuG). Compared with that of ZM3(pumDEcpsuG), the titre of pseudouridine from ZM3(pumDRspsuG) was increased by 28% and reached 506.6 mg/L, which was the highest in all the aforementioned strains (Fig. 4B). Pseudouridine were detected in the fermentation broth of Streptomyces platensis [30], but SppsuG was not efficient for pseudouridine production, and the yield of ZM3(pumDSppsuG) was only 78.2 mg/L (Fig. 4B).
-----------------------------Fig. 4----------------------------
SDS-PAGE (Fig. 4C) analysis of the supernatants showed that all the PsuG tested were correctly expressed. Enzymatic assays revealed that RspsuG had the most activity and substrate affinity than the others, and it was 6-fold more efficient and showed 1.57-fold more substrate affinity than EcpsuG (Table 2), which was consistent with fermentation results.
Table 2
Kinetic parameters* of ΨMP glycosidase
Enzyme
|
Substrate
|
Km(µM)
|
kcat (S − 1)
|
kcat/Km (mM − 1S − 1)
|
EcpsuG
|
Uracil
|
223 ± 61
|
5.5 ± 0.52
|
24.7
|
KspsuG
|
Uracil
|
378 ± 72
|
2.5 ± 0.31
|
6.6
|
RspsuG
|
Uracil
|
127 ± 18
|
18.6 ± 2.1
|
146.4
|
SspsuG
|
Uracil
|
836 ± 185
|
1.2 ± 0.09
|
1.4
|
SppsuG
|
Uracil
|
766 ± 162
|
1.5 ± 0.17
|
1.9
|
*Data are presented as mean ± s.d. (n = 3).
Fine-tuning of gene expression
Various factors, such as plasmid copy number, promoter strength, gene copy number and order of genes, can affect gene expression level [31, 32]. To further improve the yield of pseudouridine, five recombinant E. coli strains were constructed according to different genetic strategies (Fig. 5). There were significant differences in pseudouridine yield under different genetic strategies. ZM3(pZH23) exhibited the highest yield, reached to 592 mg/L (Fig. 5). The only difference between ZM3(pZH23) and ZM3(pumDEcpsuG) was the cloning order of pumD and RspsuG, however, the yield of ZM3(pumDEcpsuG) was lower than that of ZM3(pZH23). When the T7 promoter of pZH23 was replaced by trc promoter, the yield of ZM3(pZH20) decreased. When pumD and RspsuG were cloned in pCDFDuet-1 under individual T7 promoter, although the copy number increased, the yield of ZM3(pZH23) was the lowest (Fig. 5).
-----------------------------Fig. 5----------------------------
Knockout of pseudouridine catabolism-related genes
The involvement of psuK has been identified in pseudouridine metabolism in Escherichia coli UTI89 and is specifically phosphorylated pseudouridine to be pseudouridine 5’-phosphate. psuT was predicted to function as pseudouridine uptake which coorperated with psuK and psuG to metabolite extracellular pseudouridine [19, 33]. Based on their functions, knockout of psuK and psuT was expected to contribute to pseudouridine accumulation. psuK and psuT were successively deleted in ZM3, yielding ZM4 and ZM5, respectively. Plasmid pZH23, which exhibited the most efficient pseudouridine production, was transformed into ZM4 and ZM5, yielding ZM4(pZH23) and ZM5(pZH23), respectively. The deletion of psuK exerted no positive effect on pseudouridine production (Fig. 6). After deleting psuT, the yield of ZM5(pZH23) increased by 8%, reaching 631 mg/L (Fig. 6). The result of fermentation of the psuT mutant suggested that blocking pseudouridine uptake was conducive to the accumulation of pseudouridine.
-----------------------------Fig. 6----------------------------
Fed-batch fermentation of pseudouridine in a 5-L bioreactor
The strain ZM5(pZH23), which exhibited the highest yeild of pseudouridine in shake flasks, was used to explore the production potential of a fed-batch fermentation strategy under the indicated cultivation conditions. As shown in Fig. 7, during the growth stage, the 10 g/L glucose concentration initially used was almost exhausted after 4 h, and glucose was then fed at the appropriate rate. In the initial 8 h, the cells grew exponentially but no product accumulated. IPTG was added to the medium after 8 h of cultivation, and the fermentation temperature was decreased to 30°C which was optimal for the reaction of ΨMP glycosidase synthesis [20]. Then, pseudouridine began to accumulate. The cell concentration reached maximum at 32 h, after that the OD600 value began to decrease steadily, but pseudouridine continued to accumulate, exhibiting a cell growth-independent production profile. At 52 h, the final pseudouridine titre reached 7.5 g/L, indicating a productivity rate of 0.14 g/L/h.
-----------------------------Fig. 7----------------------------