Engineering a Bifunctional Phr60-Rap60-Spo0A Quorum-Sensing

Jun 19, 2019 - Our dynamic pathway regulation led to a 40-fold improvement of MK-7 ..... pentose phosphate, shikimic acid, MEP and MK-7 pathways...
0 downloads 0 Views 4MB Size
Download PDF
Research Article Cite This: ACS Synth. Biol. XXXX, XXX, XXX−XXX

pubs.acs.org/synthbio

Engineering a Bifunctional Phr60-Rap60-Spo0A Quorum-Sensing Molecular Switch for Dynamic Fine-Tuning of Menaquinone‑7 Synthesis in Bacillus subtilis Shixiu Cui,†,‡ Xueqin Lv,†,‡ Yaokang Wu,†,‡ Jianghua Li,†,‡ Guocheng Du,†,‡ Rodrigo Ledesma-Amaro,§ and Long Liu†,‡,* †

Key Laboratory of Carbohydrate Chemistry and Biotechnology, Ministry of Education, Jiangnan University, Wuxi 214122, China Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi 214122, China § Department of Bioengineering, Imperial College London, London SW7 2AZ, U.K.

Downloaded via BUFFALO STATE on July 18, 2019 at 08:34:25 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

‡

S Supporting Information *

ABSTRACT: Quorum sensing (QS)-based dynamic regulation has been widely used as basic tool for fine-tuning gene expression in response to cell density changes without adding expensive inducers. However, most reported QS systems primarily relied on down-regulation rather than up-regulation of gene expression, significantly limiting its potential as a molecular switch to control metabolic flux. To solve this challenge, we developed a bifunctional and modular Phr60Rap60-Spo0A QS system, based on two native promoters, PabrB (down-regulation by Spo0A-P) and PspoiiA (up-regulation by Spo0A-P). We constructed a library of promoters with different capacities to implement down-regulation and upregulation by changing the location, number, and sequences of the binding sites for Spo0A-P. The QS system can dynamically balance the relationship between efficient synthesis of the target product and cell growth. Finally, we validated the usefulness of this strategy by dynamic control of menaquinone-7 (MK-7) synthesis in Bacillus subtilis 168, a model Gram-positive bacterium, with the bifunctional Phr60-Rap60-Spo0A quorum sensing system. Our dynamic pathway regulation led to a 40-fold improvement of MK-7 production from 9 to 360 mg/L in shake flasks and 200 mg/L in 15-L bioreactor. Taken together, our bilayer QS system has been successfully integrated with biocatalytic functions to achieve dynamic pathway regulation in B. subtilis 168, which may be extended for use in other microbes to fine-tune gene expression and improve metabolites production. KEYWORDS: menaquinone-7, Spo0A, Rap60-Phr60, quorum sensing, Bacillus subtilis

M

sensors, and quorum sensing (QS) systems.2,5 Metabolite responsive promoters are typically identified by screening for native promoters that are responsive to a specific metabolite. For example, using the farnesyl pyrophosphate-responsive promoter, the amorphadiene was increased by more than 2fold compared to that obtained from inducible or constitutive promoters.6 Protein-based biosensors, which are typically derived from transcription factors, are capable of binding to DNA and up-regulate/down-regulate gene expression upon interacting with metabolites.7−9 For instance, FadR, an acylCoA sensor, was used to repress the expression of adhb and pdc and prevent the accumulation of toxic intermediates, ultimately improving fatty acid ethyl ester production by 3fold.10 By using metabolite-responsive promoters or bio-

etabolic engineering aims to redistribute the metabolic flux of desired pathways by rewiring metabolic networks to achieve maximum production of valuable pharmaceuticals, functional nutraceuticals, and fine chemicals.1 Because of the complexity of microbial metabolic networks, commonly used static metabolic engineering strategies (i.e., gene overexpression and knockout) are prone to overload or disrupt the native metabolic network and cause metabolic imbalances, resulting in accumulation of toxic intermediates, reduced substrate conversion, and compromised product yield. Compared to static regulation, dynamic pathway regulation has been successfully applied to balance the trade-off between cell growth and synthesis of desired products. This is particularly useful when the engineered system needs to adapt to the changing intracellular or extracellular environments,2 which has been proven effective for fine-tuning of metabolic fluxes.3,4 Currently, dynamic regulation systems could be categorized into three groups: metabolite responsive promoters, bio© XXXX American Chemical Society

Received: March 31, 2019 Published: June 19, 2019 A

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

Figure 1. Schematic diagram of quorum sensing regulation system based on Spo0A. The activity of Spo0A is governed by KinA and KinB, two intermediate protein Spo0F and Spo0B. KinA and KinB phosphorylate the intermediate regulator Spo0F, the phosphoryl group is transferred from Spo0F−P to Spo0B, then to the response regulator Spo0A. Moreover, Rap60 influence Spo0A activity by modulating the phosphoryl system. As the cell grows, the Phr60 gradually increases and enters into the cell by the Opp transporter. Phr60 inhibits the activity of Rap60, and so Spo0A is activated to regulate gene expression.

Figure 2. Analysis of spore formation rate. (a) Domain architectures of KinA. PAS domains of the amino terminal signal input domain. H405 is catalytic histidine and ATP-binding site of the kinase domain are indicated. (b) Comparison of spore germination rates of different strains. CFU was measured by the production of heat-resistant (80 °C for 20 min) in DSM and LB medium All experiments were performed in triplicate and error bars show standard deviation (SD).

sensors, the titer of target products can be greatly improved;8 however, these promoters belong to pathway-specific regulation systems,11 severely limiting their application in other metabolic pathways. The QS system overcomes the limitations of metabolite-responsive promoters and transcriptional factorbased biosensors.12 A QS system could be rewired to control cell density-dependent processes such as biofilm formation, cellular motility, and sporulation,13,14 which has the benefit of not depending on a metabolic pathway or requiring inducers or other interventions.15 For instance, the Esa QS from Pantoea stewartia has been engineered for transcriptional repression of phosphofructokinase in Escherichia coli, resulting in a 5.5-fold increase in the myo-inositol titer.15 However, most reported QS systems tend to inhibit gene expression, while there are few reports of the simultaneous up-regulation and down-regulation of gene expression by using a bifunctional QS system.15,16 Such a system would be helpful for fine-tuning complex metabolic networks. Bacillus subtilis 168 is a Gram-positive model microorganism that is widely used to produce nutraceuticals. The transcription factor Spo0A is responsible for the initiation of spore formation in B. subtilis and can regulate the expression of 121 genes.16 Spo0A could be activated by histidine autokinases (KinA−E) and two phosphorelay proteins Spo0F and Spo0B.17,18 In addition, Rap60 can inhibit Spo0A-P activity and KinA phosphorylation.14,19 The activity of Rap60 is inhibited by the signal molecule Phr60, which depends on the cell density. Therefore, we hypothesize that a QS system can be developed by controlling the phosphorylation of Spo0A

with the regulation of protein Rap60 and signal molecule Phr60 (Figure 1). When Phr60 accumulates above a certain density during cell growth, Rap60 activity is greatly inhibited, after which phosphorylated Spo0A (Spo0A-P) binds the target promoter and performs its function. Notably, Spo0A can only be phosphorylated under carbon starvation conditions, while most bioprocesses are conducted in nutrient-rich medium; therefore, the Spo0A-based QS system should be preferentially engineered to function in nutrient-rich medium. In this study, we developed a bifunctional Phr60-Rap60Spo0A-based QS system for dynamic fine-tuning of complex pathways in B. subtilis 168. First, we reconstituted the Spo0A regulatory system and Rap60-Phr60 QS system to effectively down-regulate and up-regulate gene expression in nutrient-rich medium. Next, we systematically optimized the biosynthesis of the nutraceutical menaquinone-7 (MK-7) using the Rap60Phr60-Spo0A QS system. MK-7 has been used in prevention of osteoporosis, arterial calcification, cardiovascular disease and Parkinson’s disease.20,21 Applying this dynamic control led to a 40-fold improvement of MK-7 production in shaker flask. Our Phr60-Rap60-Spo0A QS system highlights the importance of using quorum-sensing based genetic circuits for fine-tuning of complex metabolic pathways in industrially relevant species. This strategy could also be generalized to improve metabolite production and construct efficient microbial cell factories for various bioproduction applications. B

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

Figure 3. Spo0A regulates gene expression. (a) The different concentrations of IPTG effect on PspoiiA-GFP expression. Pgrac100-Spo0A was constructed to obtain different concentrations of Spo0A-P. The GFP fluorescence signals increased significantly with increasing IPTG concentration. (b) The schematic of PspoiiA, delete the “0A box” in different locations of PspoiiA. Replace the sequence of “0A box” in the promoter with a conserved sequence. (c) Different promoters are regulated by Spo0A to varying degrees. The indicated strains SP1 (spo0A::Pgrac100-spo0A) were grown in LB medium. T0 represents the time at which IPTG (20 μM) was added to induce expression of Spo0A. (d) The different concentrations of IPTG effect on PabrB-GFP expression. The GFP fluorescence signals decreased significantly with increasing IPTG concentration. (e) The schematic of PabrB delete one “0A box” of PabrB. Replace the sequence of “0A box” in the promoter with a conserved sequence. (f) Different promoters are regulated by Spo0A to varying degrees. Induction conditions are the same as (b). The blue rectangle represents “0A box” in promoter. The “+1” represent the transcription initiation site. The purple rectangle represents the conserved sequence of “0A box”. All experiments were performed in triplicate and error bars show standard deviation (SD).

■

RESULTS AND DISCUSSION Reconstitution of Spo0A-Based Gene Expression Regulation System in B. subtilis. Spo0A, a member of the response regulator family of DNA-binding protein, is activated at the start of sporulation and can directly regulate the expression of approximately 121 genes.17 KinA and KinB are mainly responsible for phosphorylating Spo0A to activate sporulation pathways under carbon starvation conditions.22 The KinA amino terminus contains sensor domains with three PAS motifs (PAS-A, -B, and -C), and PAS-A senses carbon starvation and changes the conformation of KinA to cause selfphosphorylation of KinA23 (Figure 2a). To control Spo0A phosphorylation mainly by KinA, we first deleted the KinB gene and then deleted the PAS-A domain from KinA (marked as KinAΔPAS-A) to enable phosphorylation of Spo0A in nutrient-rich medium. We then replaced the native KinA promoter with the constitutive promoter Pveg to control the expression of KinAΔPAS-A and yield the strain BS1. When the mutant strains BS1 was cultured in DSM medium (favorable for sporulation) and LB medium (nutrient-rich medium), which results in the sporulation process regardless of whether the media was nutrient-rich (Figure 2b), indicating that Spo0A was phosphorylated, although the number of spores was significantly decreased in LB medium. The spore early genes spoIIA and spo0IIE in the BS1 strain were knocked out to prevent or decrease spore formation in nutrient-rich medium and the strain BS2 was obtained (Figure 2b). The results indicated that regardless of the nutrient conditions, signaling in

strain BS2 resulted in Spo0A phosphorylation, but spore synthesis was significantly inhibited. To verify whether different concentrations of Spo0A-P would regulate gene expression to different extents, spo0A was expressed under control of the IPTG inducible promoter Pgrac100, yielding the strain SP1. A transcriptional unit with spoIIA promoter droves the expression of GFP (Spo0A-P can bind to the spoIIA promoter and activate the transcription of spoIIA) was inserted into the pHT01 vector (Figure S1) and transformed into SP1 for transcriptional activity detection. As shown in Figure 3a, as the IPTG concentration was increased, GFP fluorescence intensity forms a quantitative dose−response curve with IPTG. Additionally, after 6−7 h of IPTG induction, strong GFP fluorescence signals were detected, indicating that Spo0A-P can effectively regulate the transcription of PspoIIA. Considering that the native PspoIIA promoter contains four spo0A binding sites named as the “0A box” sequence, we mutated or deleted the “0A box” to obtain a series of promoters with different affinities for Spo0A-P. Figure 3b shows that the activity of Pspoiia(cs‑3) was 3.5-fold higher than that of native PspoIIA when the “0A box” sequence was changed to conserved binding sites, surprisingly, an 8.0-fold transcriptional activity was observed when the “0A box” sequence at 1 and 3 positions was changed to conserved binding sites. When the binding site at 3 location was knocked out to create new promoter Pspoiia(OA‑3), the activity of the Pspoiia(OA‑3) promoter decreased to 64% of the native PspoIIA (Figure 3c). Additionally, the binding site sequence of Spo0A-P has a significant effect on C

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

Figure 4. Phr60-Rap60 affects the activity of protein Spo0A. (a) Effects of different concentrations of Rap60 on the activity of Spo0A. Expression of Rap60 with Pgrac‑100, and spo0A was expressed in genome with native promoter. The IPTG added in LB medium at first. (b) The relationship between Phr60 and Rap60. Expression of Rap60 and Phr60 with Phag promoter in genome, two Phr60 tandem expressions with two operons. (c) The schematic of Spo0A-P regulated PspoiiA and PabrB promoter. Phosphorylation level of Spo0A is regulated by Rap-Phr quorum sensing gene pair. Promoter spoiiA and abrB fused GFP to monitor Spo0A-P activity, Spo0A-P activate the transcription of PspoiiA-GFP, Spo0A-P inhibit transcription of PabrB-GFP. (d) Relationship between Phr60-Rap60-Spo0A QS and cell density, promoter spoiiA was fused with GFP to monitor Spo0A-P activity. (e) Relationship between Phr60-Rap60-Spo0A QS and cell density, promoter abrB was fused with GFP to monitor Spo0A-P activity. (f) The relationship between the transcription level of promoter Pspoiia (cs‑1,3) and the dry weight of cells. (g) The relationship between the transcription level of promoter PabrB(cs‑1) and the dry weight of cells.

and obtained the promoter PabrB(cs‑2), whose activity was only 46% of the native PabrB (Figure 3f). Next, a binding site was knocked out to produce PabrB(OA‑2), which significantly reduced the inhibition intensity to a greater extent than native PabrB. These results demonstrate that increasing the affinity of Spo0A-P to the promoter PabrB leads to stronger inhibition. Design and Construction of Phr60-Rap60-Spo0A QS System. After constructing a series of promoters with different transcriptional activity under the control of Spo0A-P, we further controlled the expression of Spo0A-P using the regulatory protein Rap60 and cell density-responsive molecule Phr60. Rap60 is known to inhibit Spo0A-P activity19 and KinA phosphorylation.14 To verify the inhibitory effect of Rap60 on Spo0A-P, we first knocked out the other two Spo0A-P activity inhibitors, RapA and RapB, and then fused the PspoIIA promoter with GFP to monitor Spo0A-P activity. Next, by constructing a

the affinity of Spo0A-P to promoter. By altering the Spo0A-P binding site sequence in Pspo0iia promoter, we obtained a series of promoters with different intensity of regulation by Spo0A-P. In addition to activating gene expression, Spo0A-P can repress gene expression.24 For quantitative detection, we created a transcriptional unit where the abrB promoter drove the expression of GFP (Spo0A can bind to the abrB promoter and inhibit the transcription of abrB). As expected, the expression level of GFP driven by the PabrB promoter decreased in the presence of IPTG (Figure 3d). There are two binding sites in the native PabrB promoter, and we obtained a series of promoters with different repression intensities by changing the sequence and number of the “0A box” in PabrB (Figure 3e). First, we changed the “0A box” sequence of site 1 and obtained the promoter PabrB(cs‑1), whose activity was only 15% of the native PabrB, and then changed the “0A box” sequence of site 2 D

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

Figure 5. Pathway of MK-7 in B. subtilis 168. MK-7 synthesis is a complex process involving the synergy of multiple metabolic pathways, including glycolysis, pentose phosphate, shikimic acid, MEP and MK-7 pathways. Enzymes: MenF, isochorismate synthase; MenD, 2-succinyl-5-enolpyruvyl6-hydroxy-3-cyclohexene-1-carboxylate synthase; YtxM, 2-succinyl-6-hydroxy-2,4-cyclohexadiene-1-carboxylate synthase; MenC, o-succinylbenzoate synthase; MenE, o-succinylbenzoic acid-CoA ligase; MenB, 1,4-dihydroxy-2-naphthoyl-CoA synthase; YuxO, 1,4-dihydroxy-2-naphthoyl-CoA hydrolase; MenA, 1,4-dihydroxy-2-naphthoate heptaprenyltransferase; MenH, demethylmenaquinone methyltransferase; EntC, isochorismate synthase from E. coli K12; PpsA, phosphoenolpyruvate synthase; AroGfbr, feedback-inhibition-resistant aroG; IspA, farnesyl diphosphate synthase; HepS/HepT, heptaprenyl diphosphate synthase component I/II; Dxs, 1-deoxyxylulose-5-phosphate synthase; Dxr, 1-deoxyxylulose-5-phosphate reductoisomerase; IspD, 2-C-methylerythritol 4-phosphate cytidylyltransferase; IspE, 4-diphosphocytidyl-2-C-methylerythritol kinase; IspF, 2-Cmethylerythritol 2,4-cyclodiphosphate synthase; IspG, 4-hydroxy-3-methylbut-2-enyl diphosphate synthase; IspH, 4-hydroxy-3-methylbut-2-enyl diphosphate reductase; DhbB, isochorismatase; Abbreviation of metabolites: Glu, glucose; G-6P, glucose-6- phosphate; F6P, fructose-6phosphate; F-1,2P, fructose-1,6-bisphosphate; Gly, glycerol; G3P, glyceraldehyde-3-phosphate; PEP, phosphoenolpyruvate; Pyr, pyruvate; E4P, erythrose 4-phosphate; DAHP, 3-deoxy-arabino-heptulonate 7-phosphate; DHQ, 3-dehydroquinate; DHS, 3-dehydroshikimate; SA, shikimate; S3P, shikimate 3-phosphate; EPSP, 5-O-(1-carboxyvinyl)-3-phosphoshikimate; CHA, Chorismite; DXP, 1-deoxyxylulose-5-phosphate; MEP, methyl-erythritol-4-diphosphate; HMBPP, 1-hydroxy-2-methyl-2-butenyl 4-diphosphate; DMAPP, dimethylallyl diphosphate; IPP, isopentenyl diphosphate; GPP, Geranyl diphosphate; FPP, farnesyl diphosphate; UPP, undecaprenyl phosphate; HDP, heptaprenyl diphosphate; ICHA, isochorismate; SEPHCHC, 2-succinyl-5-enolpyruvyl-6-hydroxy-3-cyclohexene-1-carboxylate; SHCHC, 2-succinyl-6-hydroxy-2,4-cyclohexadiene-1carboxylate; OSB, 2-succinylbenzoate; OSB-CoA, 2-succinyl benzoyl-CoA; DHNA-CoA, 1,4-dihydroxy-2-naphthoyl-CoA; DHNA, 1,4-dihydroxy2-naphthoate; DMK, 2-demethylmenaquinone; MK-7, menaquinone-7.

reached an OD600 of 1.289, the GFP fluorescence signal of PabrB-GFP decreased as cell density increased (Figure 4e). A QS system controls cell density-dependent processes and does not depend on a metabolic pathway or require inducers or other interventions.15 In this study, we have constructed promoter Pspoiia(cs‑1,3) and PabrB(cs‑1) whose activities were dependent on the cell density. When the cell dry weight reached 0.36 g/L, the transcription of Pspoiia (cs‑1,3) was activated. When the cell dry weight reached 0.25 g/L, transcription of PabrB(cs‑1) was inhibited (Figure 4f,g). These results indicate that the Phr60-Rap60-Spo0A QS system can be used as a bifunctional switch to increase or inhibit the transcription of target genes in response to cell density. Boosting the Synthesis of MK-7 Precursors in B. subtilis. MK-7 is a highly valuable vitamin K2 with wide applications in the prevention of osteoporosis, arterial calcification, cardiovascular disease, and Parkinson’s disease.20,21 The biosynthetic pathway of MK-7 in B. subtilis 168 is summarized in Figure 5. The quinone skeleton of MK-7 is synthesized from chorismate via the men gene cluster (including menF, menD, ytxM, menC, menE, and menB),25 and the variable side chain (C35) is derived from isoprenoid pathway with ispA and HepS/T.26 To improve the biosynthesis of MK-7, menF, and menE were overexpressed by the strong constitutive promoter P43, yielding strain BS3; menB was

Pgrac100 promoter-rap60 DNA fragment and integrating it into the genome, we obtained strain SP2; when the concentration of IPTG reached 40 μM, the fluorescence intensity of GFP decreased significantly to only 17% of the initial intensity, indicating that Rap60 severely repressed Spo0A-P activity (Figure 4a). To investigate the relationship between Phr60 and Rap60, rap60 and phr60 were respectively expressed by the constitutive promoter Phag, yielding strain SP3, and the Phr60 peptide bound to and inhibited Rap60, thereby restoring the ability of Spo0A-P to regulate gene expression. We also observed that the expression of Phr60-Rap60 caused the expression level of PspoiiA-GFP to be slightly higher than that of Rap60 alone, which may be because Phr60 did not completely restore the inhibition of Rap60 on Spo0A-P activity. To maximize the regulatory capacity of Spo0A-P, we evaluated the ratio of the copy number of Phr60 and Rap60 at 2:1 and found that Spo0A-P activity was restored to 85% of the native level when the cells entered the stationary phase, which supported our hypothesis (Figure 4b). To examine the effect of cell density on Spo0A-P activity, Phr60 and spo0A were expressed by their native promoters. Promoters spoiiA and abrB were fused with GFP to monitor Spo0A-P activity (Figure 4c). When the cell density reached an OD600 of 1.433, the GFP fluorescence signal of PspoiiA-GFP increased with increasing cell density in the shaking flask (Figure 4d). When cell density E

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

Figure 6. Strengthening precursor supply to push carbon flux to the MK-7 synthesis pathway. Panel (a) indicates the schematic of overexpression gene with promoter P43 or Phbs and gene knockout. Panels (b) and (c) indicate the production of MK-7 and cell growth of MK-7 after modular optimization men cluster by overexpression gene with P43-menF P43-menB Phbs-menE P43-entC and ΔdhbB. Panels (d) and (e) indicate the production of MK-7 and cell growth of MK-7 after gene overexpression with P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs- aroGf br:: lox72 P43-aroK. Panels (f) and (g) indicate overexpression of MEP pathway (P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs- aroGfbr::lox72 P43-aroK Phbs-ispA P43- hepS/T Phbs-kdpG P43-dxr P43-dxs P43-fni) impact on MK-7 production and cell growth. Panel (h) indicates the schematic of increasing the copy numbers of the menA gene at lacA, thrC and dacA locus. Panels (i) and (j) indicate the different copies of menA in strain BS12 effect on MK-7 production and cell growth. Integration sites are lacA, thrC and dacA. Slashes indicate cell precipitate, rectangular indicate of fermentation supernatant All experiments were performed in triplicate and error bars show standard deviation (SD).

overexpressed by the constitutive promoter Phbs in the genome, yielding strain BS4 (Figure 6a). Additionally, we used the P43 promoter to overexpress the heterogeneous gene entC (Km = 11.93 μM, kcat = 2.12s−1),27 which encodes isochorismate synthase in E. coli K12 to replace the isochorismatase gene dhbB in BS4, strengthening the synthesis of isochorismate and yielding the strain BS5. Strain BS5 yields MK-7 titer at 9.3 mg/ L in shaking flask after 6 days of cultivation, which was nearly the same as that of wild-type B. subtilis 168 (9 mg/L) (Figure 6b), and the growth of BS5 cells was similar to that of BS3 and BS4 (Figure 6c), indicating that the men gene cluster is not a rate-limiting step in MK-7 synthesis. Additionally, the glucose and glycerin consumption rate of BS5 was significantly lower than that of B. subtilis 168 (Figure S2a,b).

On the basis of the above results, we attempted to increase the supply of the precursors chorismate acid and heptaprenyl diphosphate (HDP) to improve MK-7 synthesis. To increase the supply of chorismate, we first enhanced the synthesis of erythrose-4-phosphate and phosphoenolpyruvate (PEP) by overexpressing transketolase encoded by tkt through in situ replacement of the native promoter with Phbs, and then overexpressed phosphoenolpyruvate synthase encoded by ppsA with P43 and deleted ptsG to obtain strains BS6 and BS7, respectively. Next, we overexpressed aroA and the feedbackinhibition-resistant aroG encoded by aroGf br,28 yielding strain BS8. Moreover, shikimate kinase aroK was expressed under the P43 promoter, yielding strain BS9 to eliminate the feedbackinhibition of phenylalanine, tyrosine, and tryptophan. The titer of MK-7 was increased to 19.1 ± 0.7 mg/L in strain BS9, a 2F

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

Figure 7. Spo0A dynamically regulates key genes in the MK-7 synthesis pathway. (a) The titers of MK-7 synthesized by strains BS19 and BS20. Spo0A inhibits the expression of pyk and uppS at appropriate time points to yield strain BS19. Spo0A promotes the expression of ispH and hepS/T basis on strain BS20. Slashes indicate cell precipitate, rectangular indicate of fermentation supernatant. (b) the cell growth of BS19 and BS20. (c) The relative transcriptional levels of genes which were overexpressed. (d) The relative transcriptional levels of gene pyk (in strain BS19), uppS (in strain BS19), hepS/T (in strain 20), ispH (in strain 20), compared with that of the corresponding genes in the control strain BS18. The transcriptional levels of pyk, uppS, hepS/T, ispH were regarded as 1 in strain BS18, respectively.

fold increase compared to in strain B. subtilis 168 (Figure 6d). The growth of BS9 was delayed at an early stage of culture, but the cell concentration was consistent with that of B. subtilis 168 after 30 h of cultivation (Figure 6e). We then integrated ispA and hepS/T in the genome to promote HDP synthesis, yielding strains BS10 and BS11, respectively. To generate additional metabolites such as PYR and G3P for the MEP pathway, we also integrated the heterologous gene kdpG (encode 2dehydro-3-deoxy-phosphogluconate aldolase in Zymomonas mobiliz) in the genome, yielding the recombinant strain BS12. dxr and dxs were overexpressed under the P43 promoter in situ to eliminate the bottlenecks of the MEP pathway, yielding BS13, and we further overexpressed f ni using the P43 promoter in situ in BS13, yielding BS14. The MK-7 titer of BS14 reached 33.1 ± 0.6 mg/L (Figure 6f), which was approximately 4-fold than that in wild-type B. subtilis 168. The cell growth of BS11, BS13, and BS14 was slow before 30 h but is comparable to wild-type B. subtilis 168 after 30 h (Figure 6g). Additionally, by increasing the copy number of menA at the lacA, thrC, and dacA locus in BS14, we obtained the recombinant strains BS15, BS16, and BS17, respectively (Figure 6h). The titer of MK-7 was increased to 75 mg/L in strain BS17 (Figure 6i) without cell growth inhibition (Figure 6j) presumably due to the consumption of toxic metabolite 1hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate (HMBPP). In addition, the glucose and glycerol consumption rates of strains BS15 and BS17 were significantly higher than those of BS14, but lower than those of the wild-type (Figure S2). These results demonstrate that menA is a critical player determining the efficiency of MK-7 synthesis pathway. Dynamic Fine-Tuning of MK-7 Synthesis in B. subtilis by Bifunctional Phr60-Rap60-Spo0A QS System.

Although there are lots of progress in recent decades to improve MK-7 production,29−33 the titer of MK-7 remained low because of the inherent constraints of the metabolic pathway. For example, PEP is an important precursor for the synthesis of chorismate acid,34 but most PEP forms pyruvate in the TCA cycle,35 and thus balancing the flux distribution at the PEP branch will be a critical intervention. In contrast, UppS is a prenyltransferase that catalyzes the condensation of farnesyl diphosphate with eight molecules of isopentenyl diphosphate to generate undecaprenyl pyrophosphate (UPP, C55), which plays an essential role in cell wall formation and serves as a carrier in the synthesis of peptidoglycan and another cell wall component; if this step is blocked, cell growth will be significantly inhibited.36 Additionally, the MK-7 precursors HMBPP and DMAPP in the MEP pathway are toxic to cells,37 and thus controlling the concentrations of HMBPP and DMAPP are also critical for efficient synthesis of MK-7 without significantly compromising growth fitness. Therefore, it posts a significant challenge to achieve the maximal synthesis of MK-7 using traditional metabolic engineering strategies such as direct gene overexpression and knockout. To maximize carbon flux toward MK-7 synthesis pathway, we use the bifunctional Phr60-Rap60-Spo0A QS system to fine-tune key genes in the MK-7 biosynthesis pathway. We first expressed Rap60 with Phag and expressed Phr60 under its native promoter in BS17, yielding strain BS18. By in situ replacement of pyk and uppS promoters with the PabrB(cs‑1) promoter in strain BS18, we obtained strain BS19, where we expect that the metabolic flux from PEP to the TCA pathway would be limited and consumption of IPP would be reduced. Compared to BS17, the synthesis of MK-7 was increased by approximately 110% to 169 ± 3.8 mg/L (Figure 7a). In the G

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

Figure 8. Production of MK-7 by BS20 in a 15-L bioreactor. (a) The production of MK-7 in fed-batch fermentation, including cell precipitate and fermentation supernatant. (b) Trends of cell growth and glucose concentration during fed-culture in a 15 L bioreactor.

Table 1. Strain and Plasmids Used in This Study names Strains BS168 BS1 BS2 BS3 BS4 BS5 BS6 BS7 BS8 BS9 BS10 BS11 BS12 BS13 BS14 BS15 BS16 BS17 BS18 BS19 BS20 SP1 SP2 SP3 Plasmids P7Z6 PDG148 pHT01 pHT01-Pspoiia-gfp

characteristics

reference

Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 Phbs-aroA Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG P43-dxr P43-dxs Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG P43-dxr P43-dxs P43-fni Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG P43-dxr P43-dxs P43-fni Pmena-menA:: lacA Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG P43-dxr P43-dxs P43-fni Pmena-menA:: lacA Pmena-menA:: thrC Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG P43-dxr P43-dxs P43-fni Pmena-menA:: lacA Pmena-menA:: thrC Pmena-menA:: dacA Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG P43-dxr P43-dxs P43-fni Pmena-menA:: lacA Pmena-menA:: thrC Pmena-menA:: dacA Phag-Rap60 Pnative-Phr60:: hag Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG P43-dxr P43-dxs P43-fni Pmena-menA:: lacA Pmena-menA:: thrC Pmena-menA:: dacA Phag-Rap60 Pnative-Phr60:: hag Pabrb-pyk::pyk Pabrb-uppS::uppS Bacillus subtilis 168 Pveg-kinA-ΔPAS-A ΔkinB ΔspoIIA Δspo0IIE, P43-menF P43-menB Phbs-menE P43-entC ΔdhbB Phbs-tkt P43-ppsA ΔptsG Phbs-aroGf br:: lox72 P43-aroK Phbs-ispA P43-hepS/T Phbs-kdpG P43-dxr P43-dxs P43-fni Pmena-menA:: lacA Pmena-menA:: thrC Pmena-menA:: dacA Phag-Rap60 Pnative-Phr60:: hag Pabrb-pyk::pyk Pabrb-uppS::uppS PspoiiA-ispH::spoIIE PspoiiA-HepS/T::spoIIA BS168 derivate, expressing of spo0A with Pgrac100Bacillus subtilis 168 Pgrac100-spo0A:: spo0A BS168 derivate, expressing of Rap60 with Pgrac100Bacillus subtilis 168 Pgrac100-Rap60:: amye BS168 derivate, expressing of Rap60 and Phr60 with PhagBacillus subtilis 168 Phag-rap60, Phag-phr60:: hag Pmd18-T containing lox71-zeo-lox66 cassette Amp, Km, E. coli−B. subtilis shuttle vector, containing cre under the control of Pspac Amp, Cm, E. coli−B. subtilis shuttle vector pHT01 ferivate with GFP gene under the control of Pspoiia

H

Lab stock This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study This study

This study This study This study

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

dynamically fine-tune metabolic networks in B. subtilis to produce other chemicals.

late stage, cell growth decreased rapidly because of the limited pyruvate synthesis (Figure 7b). Considering that the MEP pathway is the only route to synthesize IPP in B. subtilis and that the intermediate of HMBPP and dimethylallyl pyrophosphate (DMAPP) are toxic to cells, we expressed ispH and hepS/T driven by the Pspoiia(cs‑1,3) promoter, yielding strain BS20, in which the yield of MK-7 reached 360 mg/L (Figure 7a). The cell growth of BS20 decreased and the maximum OD600 was 60% of B. subtilis 168 (Figure 7b), the specific MK7 production in BS20 is 85 mg/L/gDCW. MK-7 was predicted to play a crucial role in shuttling electrons between membranes bound to protein complexes in the electron transport chain. Whether a high concentration of MK-7 can switch the physiology of cells and lead to cell lysis remains unclear. Because the titer of MK-7 changed significantly in 3−4 days, we detected the relative expression level of genes on the third day. Figure 7c shows that the transcription level of genes was ranging from 2 to 24-fold; particularly, the changes in aroA and menE were more than 20-fold. The transcriptional level of hepS/T and ispH gradually increased over the culture time and the transcription level of hepS/T increased by 5.0-fold compared to in BS18. As expected, the transcription of pyk and uppS was inhibited by Spo0A-P (Figure 7d). These results indicated that the Phr60-Rap60-Spo0A QS system can effectively fine-tune gene expression and improve MK-7 production in B. subtilis 168. On the basis of the flask culture results, the recombinant strain BS20 was used to produce MK-7 in a 15-L bioreactor via fed-batch culture. Figure 8a shows the production of MK-7 in fed-batch culture with control of glucose concentration, and the maximum production of MK-7 reached 200 mg/L on the third day. On the first day of fermentation, no MK-7 was found in the fermentation supernatant, while it appeared in the fermentation supernatant on the second day. The MK-7 titer in fermentation supernatant is almost as much as that in the cell precipitate. In fed-batch cultivation with control of glucose concentration between 3 and 6 g/L (Figure 8b), a total of 1500 mL of feeding solution was added to the bioreactor during the culture period. Cell density reached OD600 = 30 at 17 h, then decreased at 50 h and stabilized at OD600 = 20, the specific MK7 production in the engineered cell is 30 mg/L/ gDCW. These results demonstrated that fed-batch cultivation with control of glucose concentration was favorable for MK-7 production by strain BS20 and provided a basis for large-scale production of MK-7. Other reported QS systems including LuxR system from Vibrio f ischeri and Esa system from Pantoea stewartia have been well-documented. Using the LuxR system, Liu et al. achieved delayed induction of key enzymes in the production of 1,4butanediol from xylose,38 and Soma et al. has constructed a cell concentration-dependent toggle switch that could divert carbon flux from TCA cycle to maximize isopropanol production.39 Gupta et al. used the Esa system to inhibit the transcription of phosphofructokinase in E. coli, which improved glucaric acid production by dynamically shifting carbon flux away from glycolysis.13 Compared to other reported QS regulation systems, our Phr60-Rap60-Spo0A system can be used as a bifunctional switch for simultaneous upregulation and down-regulation of gene expression. We used the Phr60Rap60-Spo0A system to inhibit the expression of pyk and uppS, which prevented the adverse effect on cell growth. This bifunctional Phr60-Rap60-Spo0A system can be used to

■

CONCLUSIONS In this study, we developed a bifunctional Phr60-Rap60-Spo0A QS switch to simultaneously increase and reduce the expression of multiple genes in B. subtilis. With this QSbased genetic switch, we also systematically optimized MK-7 synthesis pathways in B. subtilis and dynamically tuned the expression of key enzymes with the Phr60-Rap60-Spo0A system. The MK-7 titer in shaking flask was improved to 360 mg/L, which is the highest level ever reported, demonstrating the effectiveness of the Phr60-Rap60-Spo0A system and bringing the biotechnological production of MK-7 one step closer to industrially relevant titers.

■

METHODS

Microorganism and Medium. All strains with the genotypes of the constructed and engineered B. subtilis and plasmids used in this study are listed in Table 1. All strains were cultivated in Luria−Bertani (LB) liquid culture or on LB agar plates at 37 °C for genetic experiments. The fermentation medium consisted of 5% (w/v) glucose, 5% (w/v) glycerol, 5% (w/v) soy peptone, and 0.06% (w/v) KH2PO4. Media were sterilized for 30 min at 121 °C. Appropriate antibiotics were added to the medium: kanamycin (50 μg/mL), zeocin (20 μg/ mL), ampicillin (100 μg/mL), spectinomycin (50 μg/mL), and chloramphenicol (5 μg/mL). All chemicals were purchased from Sangon Biotech Co., Ltd. (Shanghai, China). Pure MK-7 was purchased from ChromaDex (Irvine, CA, USA). Soy peptone, glycerol, and glucose were purchased from Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). Methanol, dichloromethane, 2-propanol, and n-hexane were obtained from Sigma-Aldrich (St. Louis, MO, USA). Shake Flask Culture of Engineered B. subtilis. The engineered B. subtilis strains were grown in 3 mL LB at 37 °C with shaking at 220 rpm for 10 h in a 15 mL tube. A seed culture (1 mL) was inoculated into a 250 mL baffled flask containing 20 mL of fermentation medium and grown at 40 °C with shaking at 220 rpm for 7 days. A volume of 0.5 mL of fermentation liquid was sampled every 24 h for MK-7 concentration analysis. B. subtilis Transformation and DNA Assembly. DNA constructs were generated using previously described20,40 assembly methods and transformed into B. subtilis by electroporation. The modified DNA fragments and plasmids were sequenced by GENEWIZ (South Plainfield, NJ, USA). Primers used for genetic modification are listed in Supplementary Table S1. In Situ Replacement of Promoter, Genomic Integration of Heterogeneous Genes and Gene Knockout. A marker-free genetic modification strategy was used as previously reported.20,40 The expression of homologous genes was strengthened by replacing the native promoter with the strong constitutive promoter P43 or Phbs. The upstream, target sequences, and P43 or Phbs promoter sequences were amplified from B. subtilis. These three fragments and the lox71−zeo−lox66 cassettes were joined by triple-fusion PCR. The purified PCR products were used to transform competent B. subtilis cells. The resistance marker from the host strain was removed by using the Cre/lox system. The PDG148 plasmid was introduced into zeor clones to promote the recombination I

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

Biolabs, Ipswich, MA, USA) was used to synthesize cDNA from total RNA according to the manufacturer’s instructions. qRT-PCR amplifications were carried out using the Bio-Rad IQ5 Real-Time PCR system (Hercules, CA, USA). The PCR samples were prepared in a total volume of 20 μL using 2 × Power SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA) with 1.5 μL of cDNA and a final primer concentration of 200 nM. The 16 S rRNA gene was used as an internal standard. PCR conditions were 95 °C for 3 min, followed by 50 three-step cycles at 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s. Fed-Batch Culture in 15-L Bioreactor. The production of MK-7 by fed-batch culture of BS20 was performed at an initial glucose concentration of 10 g/L. The fermentation medium used for fed-batch culture consisted of 50 g/L soy peptone, 5 g/L glycerol and 0.6 g/L KH2PO4. The feeding solution contained 500 g/L glucose and 500 g/L glycerol. Seed culture was carried out in 2 L shake flasks containing 500 mL of seed medium with shaking at 220 rpm and 37 °C for 12 h on rotary shakers. The seed culture was inoculated into 15-L fermenter containing 10 L fermentation medium. The pH was natural, and the temperature was maintained at 40 °C. The aeration rate and agitation speed were 10 vvm and 500 rpm, respectively. In fed-batch culture, whenever the residual glucose concentration fell to blew 5 g/L, the feeding solution was pumped into the fermenter to restore the glucose concentration to 3−6 g/L. The feeding rates were adjusted every 2 h based on the concentration of residual glucose in the fermentation medium. Statistical Analysis. All experiments were independently carried out at least three times and the results were expressed as mean ± standard deviation (SD).

of lox71 and lox66, which removed the resistance marker cassette. Next, the intracellular plasmid PDG148 was removed by incubation at 50 °C for 12 h. Heterogeneous genes were expressed under control of the strong constitutive promoter P43 or Phbs. Five fragments, including the upstream and downstream sequences (1000 base pairs, bp) flanking the integration locus, lox71−zeo−lox66 cassettes, P43 or Phbs strong constitutive promoter, and target gene sequences, were joined by triple-fusion PCR, and the purified PCR products were used to transform competent B. subtilis cells. Heterologous genes were synthesized with codon optimization by GenScript Biotech (Piscataway, NJ, USA). The method used for gene knockout was similar to that used for gene overexpression. Briefly, the upstream and downstream sequences (both approximately 1000 bp) flanking the deletion targets were amplified from B. subtilis. These two fragments and the lox71−zeo−lox66 cassettes, which were amplified from the plasmid p7Z6, were joined by triple-fusion PCR. Purified PCR products were used to transform competent B. subtilis cells. The promoter sequence used in this study was the same as that used by Yang et al.21 MK-7 Extraction.41 Centrifugal fermentation media collects supernatant as fermentation supernatant, then resuspend the cell pellet as cell precipitate using an equal volume of sterile water. MK-7 was extracted from the fermentation supernatant and cell precipitate using a mixture of 2-propanol and n-hexane (1:2, v/v) in a 4:1 ratio (organic:liquid, v/v). The mixture was vigorously shaken with a vortex mixer for 10 min and then centrifuged at 5000g for 5 min to separate the two phases. The organic phase was then separated to recover the extracted MK-7. Sporulation Assays. Overnight cultures were inoculated into Difco nutrient broth sporulation medium (DSM; Detroit, MI, USA) and fresh LB cultures were grown at 37 °C with vigorous aeration for approximately 20 h after entry into the stationary phase and heated to 80 °C for 20 min. Serial dilutions were prepared with each culture using sterile water, and 100 μL of the culture was spread onto LB plates. Viable colony-forming units were evaluated after overnight incubation at 37 °C. Analytical Methods. Cell density was determined from the optical density at 600 nm with a spectrophotometer after suitable dilution with deionized water. The concentration of glycerol was assayed with free glycerol determination Kit (Jiancheng, Nanjing, China). The glucose concentration in the medium was measured with a glucose−glutamate analyzer (SBA-40C; Biology Institute of Shandong Academy of Sciences, Shandong, China). High-performance liquid chromatography (Agilent 1260, Santa Clara, CA, USA) equipped with a photon diode array UV detector and C18 ODS column (5 μm, 250 × 4.6 mm, Thermo Fisher Scientific, Waltham, MA, USA) was used at 40 °C to determine the MK-7 concentration. Methanol: dichloromethane (9:1, v/v) was used as the mobile phase at a flow rate of 1 mL/min. A wavelength of 254 nm was used for calibration and analysis. The MK-7 calibration curve was linear between 1 mg/L and 100 mg/L (R2 = 0.998). Quantitative Reverse-Transcriptase PCR (qRT-PCR) Analysis. Cells were harvested during different growth phases from the shaking flask cultures. Total RNA was purified with the RNeasy Mini Kit (Qiagen, Hilden, Germany). A ProtoScript II First Strand cDNA Synthesis Kit (New England

■

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acssynbio.9b00140.

■

Table S1, Figures S1−S2, Supplementary Sequence (PDF)

AUTHOR INFORMATION

Corresponding Author

*Tel: +86-510-85918312. Fax: +86-510 85918309. E-mail: [email protected]. ORCID

Long Liu: 0000-0002-9679-9130 Author Contributions

L.L. designed the research. S.X.C., X.Q.L., and Y.K.W. performed the experiments. L.L., S.X.C., J.H.L., G.C.D., and L.A.R analyzed the data and wrote the manuscript. Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS This work was financially supported by the National Natural Science Foundation of China (31871784, 31870069, 21676119, and 31671845). Postgraduate Research & Practice Innovation Program of Jiangsu Provence J

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology

(20) El Asmar, M. S., Naoum, J. J., and Arbid, E. J. (2014) Vitamin k dependent proteins and the role of vitamin k2 in the modulation of vascular calcification: a review. Oman Med. J. 29, 172−177. (21) Scheiber, D., Veulemans, V., Horn, P., Chatrou, M. L., Potthoff, S. A., Kelm, M., Schurgers, L. J., and Westenfeld, R. (2015) High-dose menaquinone-7 supplementation reduces cardiovascular calcification in a murine model of extraosseous calcification. Nutrients 7, 6991− 7011. (22) LeDeaux, J. R., Yu, N., and Grossman, A. D. (1995) Different roles for KinA, KinB, and KinC in the initiation of sporulation in Bacillus subtilis. J. Bacteriol. 177, 861−863. (23) Stephenson, K., and Hoch, J. A. (2001) PAS-A domain of phosphorelay sensor kinase A: A catalytic ATP-binding domain involved in the initiation of development in Bacillus subtilis. Proc. Natl. Acad. Sci. U. S. A. 98, 15251−15256. (24) Hamon, M. A., Stanley, N. R., Britton, R. A., Grossman, A. D., and Lazazzera, B. A. (2004) Identification of AbrB-regulated genes involved in biofilm formation by Bacillus subtilis. Mol. Microbiol. 52, 847−860. (25) Dairi, T. (2012) Menaquinone biosyntheses in microorganisms. Methods Enzymol. 515, 107−122. (26) Guan, Z., Xue, D., Abdallah, I., Dijkshoorn, L., Setroikromo, R., Lv, G., and Quax, W. J. (2015) Metabolic engineering of Bacillus subtilis for terpenoid production. Appl. Microbiol. Biotechnol. 99, 9395−9406. (27) Lin, Y., Shen, X., Yuan, Q., and Yan, Y. (2013) Microbial biosynthesis of the anticoagulant precursor 4-hydroxycoumarin. Nat. Commun. 4, 2603. (28) Kikuchi, Y., Tsujimoto, K., and Kurahashi, O. (1997) Mutational analysis of the feedback sites of phenylalanine-sensitive 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase of Escherichia coli. Appl. Environ. Microbiol. 63, 761−762. (29) Yang, S., Cao, Y., Sun, L., Li, C., Lin, X., Cai, Z., Zhang, G., and Song, H. (2019) Modular pathway engineering of Bacillus subtilis to promote de novo biosynthesis of menaquinone-7. ACS Synth. Biol. 8, 70−81. (30) Ebrahiminezhad, A., Varma, V., Yang, S., and Berenjian, A. (2016) Magnetic immobilization of Bacillus subtilis natto cells for menaquinone-7 fermentation. Appl. Microbiol. Biotechnol. 100, 173− 180. (31) Mahanama, R., Berenjian, A., Valtchev, P., Talbot, A., Biffin, R., Regtop, H., Dehghani, F., and Kavanagh, J. M. (2011) Enhanced production of menaquinone 7 via solid substrate fermentation from Bacillus subtilis. Int. J. Food Eng., DOI: 10.2202/1556-3758.2314. (32) Mahdinia, E., Demirci, A., and Berenjian, A. (2017) Strain and plastic composite support (PCS) selection for vitamin K (Menaquinone-7) production in biofilm reactors. Bioprocess Biosyst. Eng. 40, 1507−1517. (33) Singh, R., Puri, A., and Panda, B. P. (2015) Development of menaquinone-7 enriched nutraceutical: inside into medium engineering and process modeling. J. Food Sci. Technol. 52, 5212−5219. (34) Licona-Cassani, C., Lara, A. R., Cabrera-Valladares, N., Escalante, A., Hernandez-Chavez, G., Martinez, A., Bolivar, F., and Gosset, G. (2014) Inactivation of pyruvate kinase or the phosphoenolpyruvate: sugar phosphotransferase system increases shikimic and dehydroshikimic acid yields from glucose in Bacillus subtilis. J. Mol. Microb Biote. 24, 37−45. (35) Liu, D. F., Ai, G. M., Zheng, Q. X., Liu, C., Jiang, C. Y., Liu, L. X., Zhang, B., Liu, Y. M., Yang, C., and Liu, S. J. (2014) Metabolic flux responses to genetic modification for shikimic acid production by Bacillus subtilis strains. Microb. Cell Fact. 13, 40. (36) Wang, Y., Desai, J., Zhang, Y., Malwal, S. R., Shin, C. J., Feng, X., Sun, H., Liu, G., Guo, R. T., and Oldfield, E. (2016) Bacterial cell growth Inhibitors targeting undecaprenyl diphosphate synthase and undecaprenyl diphosphate Phosphatase. ChemMedChem 11, 2311− 2319. (37) Li, Q., Fan, F., Gao, X., Yang, C., Bi, C., Tang, J., Liu, T., and Zhang, X. (2017) Balanced activation of IspG and IspH to eliminate

(KYCX18_1786). The Fundamental Research Funds for the Central Universities (JUSRP51713B).

■

REFERENCES

(1) Furusawa, C., Horinouchi, T., Hirasawa, T., and Shimizu, H. (2012) Systems metabolic engineering: the creation of microbial cell factories by rational metabolic design and evolution. Adv. Biochem. Eng./Biotechnol. 131, 1−23. (2) Xu, P. (2018) Production of chemicals using dynamic control of metabolic fluxes. Curr. Opin. Biotechnol. 53, 12−19. (3) Anesiadis, N., Cluett, W. R., and Mahadevan, R. (2008) Dynamic metabolic engineering for increasing bioprocess productivity. Metab. Eng. 10, 255−266. (4) Gadkar, K. G., Doyle, F. J., Iii, Edwards, J. S., and Mahadevan, R. (2005) Estimating optimal profiles of genetic alterations using constraint-based models. Biotechnol. Bioeng. 89, 243−251. (5) Tan, S. Z., and Prather, K. L. (2017) Dynamic pathway regulation: recent advances and methods of construction. Curr. Opin. Chem. Biol. 41, 28−35. (6) Dahl, R. H., Zhang, F., Alonso-Gutierrez, J., Baidoo, E., Batth, T. S., Redding-Johanson, A. M., Petzold, C. J., Mukhopadhyay, A., Lee, T. S., Adams, P. D., and Keasling, J. D. (2013) Engineering dynamic pathway regulation using stress-response promoters. Nat. Biotechnol. 31, 1039−1046. (7) Zhang, F., and Keasling, J. (2011) Biosensors and their applications in microbial metabolic engineering. Trends Microbiol. 19, 323−329. (8) Xu, P., Li, L., Zhang, F., Stephanopoulos, G., and Koffas, M. (2014) Improving fatty acids production by engineering dynamic pathway regulation and metabolic control. Proc. Natl. Acad. Sci. U. S. A. 111, 11299−11304. (9) Xu, P., Wang, W., Li, L., Bhan, N., Zhang, F., and Koffas, M. (2014) Design and kinetic analysis of a hybrid promoter-regulator system for malonyl-CoA sensing in Escherichia coli. ACS Chem. Biol. 9, 451−458. (10) Zhang, F., Carothers, J. M., and Keasling, J. D. (2012) Design of a dynamic sensor-regulator system for production of chemicals and fuels derived from fatty acids. Nat. Biotechnol. 30, 354−359. (11) Wan, X., Marsafari, M., and Xu, P. (2019) Engineering metabolite-responsive transcriptional factors to sense small molecules in eukaryotes: current state and perspectives. Microb. Cell Fact. 18, 61. (12) Lv, Y., Qian, S., Du, G., Chen, J., Zhou, J., and Xu, P. (2019) Coupling feedback genetic circuits with growth phenotype for dynamic population control and intelligent bioproduction. Metab. Eng. 54, 109−116. (13) Lyon, G. J., and Novick, R. P. (2004) Peptide signaling in Staphylococcus aureus and other Gram-positive bacteria. Peptides 25, 1389−1403. (14) Boguslawski, K. M., Hill, P. A., and Griffith, K. L. (2015) Novel mechanisms of controlling the activities of the transcription factors Spo0A and ComA by the plasmid-encoded quorum sensing regulators Rap60-Phr60 in Bacillus subtilis. Mol. Microbiol. 96, 325−348. (15) Gupta, A., Reizman, I. M., Reisch, C. R., and Prather, K. L. (2017) Dynamic regulation of metabolic flux in engineered bacteria using a pathway-independent quorum-sensing circuit. Nat. Biotechnol. 35, 273−279. (16) Anesiadis, N., Kobayashi, H., Cluett, W. R., and Mahadevan, R. (2013) Analysis and design of a genetic circuit for dynamic metabolic engineering. ACS Synth. Biol. 2, 442−452. (17) Piggot, P. J., and Hilbert, D. W. (2004) Sporulation of Bacillus subtilis. Curr. Opin. Microbiol. 7, 579−586. (18) Molle, V., Fujita, M., Jensen, S. T., Eichenberger, P., GonzálezPastor, J. E., Liu, J. S., and Losick, R. (2003) The Spo0A regulon of Bacillus subtilis. Mol. Microbiol. 50, 1683−1701. (19) Koetje, E. J., Hajdo-Milasinovic, A., Kiewiet, R., Bron, S., and Tjalsma, H. (2003) A plasmid-borne Rap-Phr system of Bacillus subtilis can mediate cell-density controlled production of extracellular proteases. Microbiology 149, 19−28. K

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX

Research Article

ACS Synthetic Biology MEP intermediate accumulation and improve isoprenoids production in Escherichia coli. Metab. Eng. 44, 13−21. (38) Liu, H., and Lu, T. (2015) Autonomous production of 1, 4butanediol via a de novo biosynthesis pathway in engineered Escherichia coli. Metab. Eng. 29, 135−141. (39) Soma, Y., and Hanai, T. (2015) Self-induced metabolic state switching by a tunable cell density sensor for microbial isopropanol production. Metab. Eng. 30, 7−15. (40) Yan, X., Yu, H.-J., Hong, Q., and Li, S.-P. (2008) Cre/lox system and PCR-based genome engineering in Bacillus subtilis. Appl. Environ. Microb. 74, 5556−5562. (41) Yang, S., Du, G., Chen, J., and Kang, Z. (2017) Characterization and application of endogenous phase-dependent promoters in Bacillus subtilis. Appl. Microbiol. Biotechnol. 101, 4151−4161.

L

DOI: 10.1021/acssynbio.9b00140 ACS Synth. Biol. XXXX, XXX, XXX−XXX