γ-谷氨酰转肽酶(γ-glutamyltranspeptidase,GGT)在临床诊断[1-2]、催化合成氨基酸衍生物[3]、药物或药物前体[4]以及改善苦味氨基酸的苦味等方面具有广泛应用[5]。研究表明氨基酸L-苯丙氨酸(L-phenylalanine,L-Phe)在食物中可以呈现出苦味[6],又是人体必需氨基酸[7]。SUZUKI等为改善其苦味,利用Escherichia coli K-12来源的GGT 0.5 U/mL,以200 mmol/L L-谷氨酰胺(L-glutamine,L-Gln)为供体底物,200 mmol/L L-Phe为受体底物,最终γ-谷氨酰苯丙氨酸(γ-glutamylphenylalanine,γ-Glu-Phe)产率达到70%[6]。
Bacillus subtilis 168(B.subtilis 168)作为GGT表达的宿主细胞,存在分泌效率低、表达量不高的问题[8-9]。研究团队前期引入分子伴侣脂蛋白PrsA以及优化Poly(A/T)尾巴的长度,酶活力达到30.4 U/mL[10]。转录水平的高低对于基因表达量具有重要影响[11],启动子对于调节基因的转录水平发挥着重要的作用[12]。通过对-10区上游的保守区域进行非理性改造,可以提高目的蛋白表达量[13-14]。然而,目的基因在相同启动子的调控下,表达量也会受到起始翻译速率的影响[15]。核糖体结合位点(ribosome binding site,RBS)是影响起始翻译速率的重要生物学元件。通过优化编码目的基因的RBS序列,可使目标酶表达量提高10倍[16]。根据文献报道,RBS Calculator(https://www.denovodna.com/software/)可以计算得到具有不同起始翻译速率的RBS序列[17-18]。
传统方法获得高表达量的GGT菌株耗时且困难。本研究建立了一种基于显色反应的高通量筛选方法,能快速且准确筛选GGT高表达量菌株。通过对于启动子HpaII的-10区上游关键区域进行两轮迭代突变建立了启动子突变文库,经高通量筛选得到最优启动子。在此基础上,构建人工RBS序列文库,经筛选得到最优人工RBS序列,最大程度提高GGT表达量。另外,优化了GGT催化L-Phe生成γ-Glu-Phe的转化条件,提高了γ-Glu-Phe产率,从而为GGT应用于食品行业提供实验基础。
1.1.1 菌株与质粒
枯草芽孢杆菌B.subtilis 168、pMA5-ggt-prsA/B.subtilis 168、大肠杆菌Escherichia coli JM109、质粒pMA5由本实验室保藏。
1.1.2 酶和试剂
BamH I、Mlu I、Dpn I限制性核酸内切酶、DNA聚合酶,TaKaRa公司;高保真酶、同源重组酶克隆试剂盒,南京诺维赞生物科技有限公司;琼脂糖凝胶DNA回收试剂盒、小量质粒提取试剂盒、细菌DNA基因组提取试剂盒,上海捷瑞生物工程有限公司;卡那霉素、氨苄霉素、γ-谷氨酰对硝基苯胺一水合物、双甘二肽、L-Phe、L-Gln,生工生物工程(上海)股份有限公司;γ-Glu-Phe标样、其他分析纯试剂,国药集团化学试剂有限公司。
1.1.3 培养基
LB液体培养基(g/L):酵母粉5.0,胰蛋白10.0,NaCl 10.0[19]。
LB固体培养基(g/L):酵母粉5.0,胰蛋白栋10.0,NaCl 10.0,琼脂粉20。
SpI以及SpII培养基:
SpA(g/L):K2HPO4·3H2O 28.0,KH2PO4 12.0,(NH4)2SO4 4.0,柠檬酸三钠2.0。
SpB(g/L):MgSO4·7H2O 0.4。
100×CAYE(g/L):酵母粉100.0,水解酪蛋白20.0, MgCl2·6H2O 5.08,CaCl2 5.5。
发酵培养基(g/L):葡萄糖10.0,酵母粉20.0,胰蛋白胨15.0,K2HPO4·3H2O 2.0,MgSO4·7H2O 1.0,pH 7.2。
1.2.1 质粒提取和引物的设计
pMA5-ggt-prsA质粒的提取按照质粒提取试剂盒的方法进行。在本文中涉及到的引物如表1所示。
表1 本研究中使用的引物
Table 1 Primers used in this study
引物名称引物序列(5′-3′)长度/bpF1ATATAAAGTATAGTGNNNTATACTTTACTTGGA33R1TCCAAGTAAAGTATANNNCACTATACTTTATAT33F2ATATAAAGTATNNNNTGATATACTTTACTTGGA33R2TCCAAGTAAAGTATATCANNNNATACTTTATAT33F3CGAAAATGCCTCACATTTGTACTGTATTCACTAAAAAAGGTAGGTCGTATTTCCATATGAGTTAT65F4CGAAAATGCCTCACATTTGTCCCGTCCGCGGCCAGGAGCCCTTACCATATGAGTTA56R3GACCTCTAGAACGCGTTTATTTAGAATTGCTT32F5ACGCGTTCTAGAGGTCGAAATTCAC25R4AAATGCCTCACATTTGTGCCACCTAA25
注;表中下划线代表突变序列:斜体代表筛选得到的RBS361和RBS362序列
1.2.2 构建启动子突变文库
以提取的质粒pMA5-ggt-prsA作为模板,以引物F1和R1进行PCR扩增,PCR产物在Dpn I的作用下37 ℃消化1 h,将消化产物转化入B.subtilis 168感受态细胞中,将重组细胞涂布在LB固体培养基(含50 μg/mL卡那霉素)上培养,菌落PCR筛选鉴定后,建立第一批启动子突变文库。在第一轮得到的正向启动子菌株中提取质粒,送至苏州金唯智生物科技有限公司测序,再以此质粒作为模板,以引物F2以及R2进行扩增,PCR产物在Dpn I的作用下,37 ℃消化1 h,将消化产物转化入B.subtilis 168感受态细胞中,将重组细胞涂布在LB固体培养基(含50 μg/mL卡那霉素)上培养,菌落PCR筛选鉴定后,建立起第二批启动子突变文库。
1.2.3 设计RBS序列文库
在得到最优启动子的基础上,以相应的质粒作为模板,以引物F3、F4和R3进行PCR,将人工RBS序列替换原始RBS,另外,以质粒pMA5为模板,利用引物F5和R4进行PCR,PCR产物在Dpn I 的作用下,37 ℃消化1 h。将消化产物、改造后的ggt基因经核酸胶回收后,利用同源重组试剂盒进行连接,将连接产物转化入B.subtilis 168感受态细胞中,涂布在LB固体培养基(含50 μg/mL卡那霉素)上培养,菌落PCR筛选鉴定后,建立起RBS文库。
1.2.4 高通量筛选操作流程
首先挑取1.2.2、1.2.3中的突变体文库中的阳性克隆子,在含有3 mL发酵培养基的96孔板中,37 ℃,200 r/min条件下培养48 h。离心获得含有GGT的发酵上清液,在含有50 mmol/L L-Phe、50 mmol/L L-Gln 的3 mL的发酵培养基中,在反应温度为37 ℃,用0.1 mol/L Tri-HCl控制反应pH为10.0,转速为200 r/min的条件下,反应10 min。取反应液10 μL与100 μL 20 mmol/L邻苯二甲醛(O-phthalaldehyde,OPA)工作液充分混合,振荡混匀20 min,然后先加入100 μL二甲基亚砜(dimethyl sulfoxide,DMSO),再加入50 μL 10%三氯乙酸(trichloroacetic acid,TCA)充分混匀10 min,利用酶标仪在600 nm处记录吸光值。
1.2.5 实时荧光定量PCR进行基因转录丰度分析
取培养好的细胞发酵液2 mL置于离心管中,12 000 r/min的条件下离心5 min,弃去上清液,菌体可利用诺唯赞试剂盒进行RNA提取。利用实时荧光定量PCR分析不同菌株GGT转录水平的具体步骤参考文献所述[10]。
1.2.6 酶的表达与纯化
取-40 ℃冰箱保存的菌种于固体LB平板划线活化,挑取单菌落接入含有卡那霉素(50 μg/mL)抗性的10 mL LB培养液中。37 ℃过夜培养后按照6%的接种量转接入含相同浓度抗生素50 mL发酵培养基中继续培养48 h,离心后获得含有目的蛋白的发酵液上清液。
首先配制pH 10.0的饱和硫酸铵溶液,以0%~20%,20%~40%,40%~60%,60%~80%,80%~100%饱和度硫酸铵梯度分级沉淀。取体积分数为40%的饱和硫酸铵沉淀的活性蛋白上清液,再以0%~20%,20%~40%,40%~60%,60%~80%,80%~100%饱和硫酸铵梯度沉淀,从而得到含有活性蛋白的沉淀,以Tri-HCl缓冲液(pH 8.0,0.1 mol/L)复溶,经透析后以Hitrap DEAE FF柱层析分离,上样缓冲液为Tri-HCl(pH 8.0,0.05 mol/L),洗脱液为Tri-HCl(pH 8.0,0.05 mol/L)+0.5 mol/L NaCl,以0%~100%B,30 min梯度洗脱,收集洗脱峰[20]。将得到的纯酶通过SDS-PAGE方法进行分析,以及采用Bradford法测定蛋白浓度。
1.2.7 酶活力的测定
取培养48 h的菌液50 mL,于4 ℃条件下10 000 r/min离心30 min收集含有活性目的蛋白的上清液。酶活力测定方法如下:1 mL的标准反应体系中包含100 mmol/L Tris-HCl(pH 10.0),2.5 mmol/L γ-L-谷氨酰-对硝基苯胺一水合物,60 mmol/L双甘二肽和5 μL适当稀释的酶液。37 ℃反应15 min,加入200 μL 3.5 mol/L的醋酸终止反应,以不添加双甘二肽,添加等体积的水的反应体系作为对照,利用酶标仪在410 nm处测定吸光度值。通过实验组和对照组的吸光度差值计算转肽酶活力。酶活力定义为:1 min 生成1 μmol的对硝基苯胺所需的酶量为1个酶活力单位。
1.2.8 γ-Glu-Phe的检测
以GGT作为催化剂,以L-Gln、L-Phe为底物,在适宜催化反应条件下进行转化,定时取样,每次取样过程中加入等体积的10%三氟乙酸终止转肽反应。得到的反应产物,通过0.22 μm过滤膜过滤。利用高效液相色谱仪1200 Infinity系列进行检测,检测条件为:紫外检测器,X BridgetM C18柱(Agilent 5 μm,4.6 mm×250 mm),流动相A为50 mmol/L乙酸钠(含1%的N, N-二甲基甲酰胺),流动相B为50%乙腈,流速1 mL/min,柱温40 ℃。
和OPA反应生成无色化合物,与TCA和DMSO混合时,可通过紫外光谱法定量测定[21]。GGT催化L-Gln和L-Phe生成γ-Glu-Phe的过程中,会释放在底物L-Gln、L-Phe和产物γ-Glu-Phe不会影响高通量筛选过程的情况下(表2),释放的浓度与600 nm的吸光值呈正相关(图1)。基于此,建立了一种高通量筛选方法,用于筛选获得GGT高表达量的菌株。具体操作流程见方法1.2.4。
表2 OPA与反应体系中化合物之间的显色反应
Table 2 Chromogenic reactions between OPA and
compounds present in the reaction system
反应物(10 mmol/L)OPAL-GlnL-Pheγ-Glu-PheNH+4H2OA600nm0.060.060.070.060.510.06
注:H2O为空白对照,为阳性对照,OPA为阴性对照
图1 释放的浓度和A600nm的关系
Fig.1 The relationship between the concentration of released
HpaII是强组成型启动子,其-10区上游7个碱基是调控基因转录水平的关键区域。通过对关键区域进行非理性改造,提高基因的转录水平以及表达量。在得到的目的蛋白表达量提高的正向启动子突变中,通过实时荧光定量PCR、SDS-PAGE与酶活力检测等方式分析了基因的转录水平、蛋白表达以及酶活力情况。结果如图2、图3所示,较原始启动子HpaII相比,突变体EA13HpaII、EB6HpaII、ED11HpaII、EF3HpaII、EG2HpaII(基因序列见表3)启动子调控的基因转录水平明显提高,其中EF3HpaII的转录丰度较原始启动子HpaII提高了5倍。但是,检测pMA5-EF3HpaII-ggt-prsA/B.subtilis 168酶活力较pMA5-HpaII-ggt-prsA/B.subtilis 168(WT)仅提高了63%。这也就表明GGT的表达量不仅仅和基因的转录水平相关。
图2 不同启动子菌株中GGT的表达量以及转录丰度
Fig.2 Transcription level and GGT production of recombinant
strains with different promoters
注:GGT 以HpaII(WT)、EA13HpaII、 EB6HpaII、
ED11HpaII、EF3HpaII及EG2HpaII为启动子进行表达
M-Marker;1-阴性对照,B.subtilis 168/pMA5;2~7-GGT在
B.subtilis 168中以 HpaII、ED11HpaII、EB6HpaII、EG2HpaII、
EA13HpaII及EF3HpaII为启动子进行表达
图3 SDS-PAGE分析不同启动子重组菌株中GGT的表达量
Fig.3 SDS-PAGE analysis of GGT production in recombinant
strains with different promoters.
表3 不同启动子-10区上游的7个碱基序列
Table 3 The 7 bases upstream of the -10 region in
different promoters
启动子-10区上游的7个碱基序列WTAGTGTGTEA13AGCGTGAEB6AATCTGAED11TGGGTGAEF3TGTGTGAEG2ATTATGA
为了探究2.2中ggt基因转录水平与蛋白表达提高程度存在重大差异的原因,首先,我们用DNAMAN软件分析了基因起始密码子前40 bp左右的EF3HpaIIRBS序列,在EF3HpaIIRBS的SD序列中,存在着多个明显的颈环二级结构,这严重影响基因的起始翻译速率,更有甚者导致基因无法表达[21-22]。
另外,通过RBS Calculator分析EA13HpaII、EB6HpaII、ED11HpaII、EF3HpaII、EG2HpaII以及HpaII起始翻译速率(表4),HpaII的起始翻译速率明显低于EF3HpaII等其他正向启动子。通过以上结果,表明基因的起始翻译速率可能就是造成基因转录水平与蛋白表达提高程度不一致的原因。因此,找到与正向启动子更加适配的RBS序列是进一步提高目的基因表达量的关键。
表4 不同RBS序列的起始翻译速率
Table 4 The initial translation speeds of different
RBS sequences
RBS起始翻译速率/auHpaIIRBS353 243.79EA13HpaIIRBS1 054 346.62EB6HpaIIRBS754 881.51ED11HpaIIRBS1 553 243.46EF3HpaIIRBS1 734 332.37EG2HpaIIRBS98 525.41
为了进一步提高GGT的表达量,我们利用RBS Calculator设计与启动子正向突变菌株pMA5-EF3HpaII-ggt-prsA/B.subtilis 168更加适配的RBS序列,从而得到具有梯度起始翻译速率(1 000~4 986 270.7 au)的500个RBS序列文库。按照其起始翻译速率不同,将其共分为10组,1×103~2×103、2×103~5×103、5×103~1×104、1×104~5×104、5×104~1×105、1×105~2×105、2×105~5×105、5×105~1×106、1×106~2×106、2×106~4 986 270.7 au,每组50个RBS。
最终,我们得到362个较EF3HpaII原始RBS酶活力提高的人工RBS序列。根据不同RBS序列调控表达的GGT的转肽酶活力差异,从低到高排列为RBS1~RBS362,其中RBS361、RBS362起始翻译速率为2 986 410.3和3 974 351.5 au。结果如图4所示,重组菌株RBS361和RBS362经培养表达后,转肽酶活力为56.95和60.25 U/mL,是菌株pMA5-HpaII-ggt-prsA/B.subtilis 168(WT)的1.91和2.02倍。基于此,我们又利用生物分析软件DNAMAN分析了RBS361和RBS362的SD序列,结果如图5所示,经优化后的人工RBS361和RBS362的颈环二级结构得到明显改善。以上研究结果表明,在菌株pMA5-EF3HpaII-ggt-prsA/B.subtilis 168中,人工RBS序列比原始RBS序列更适配,从而能更大程度上提高GGT的表达量。
图4 不同RBS序列重组菌株GGT的表达量
Fig.4 GGT production of recombinant strains with different
RBS sequences
GGT在HpaII(WT)、EF3HpaII、EF3HpaIIRBS361及
EF3HpaIIRBS362调控下进行表达
a-EF3HpaII RBS;b-EF3HpaII RBS361;c-EF3HpaII RBS362
图5 不同RBS中SD序列的二级结构
Fig.5 The secondary structure of SD sequences in
different RBS
在前期的研究中,我们已经对于GGT的酶学性质进行了探索[10]。本研究在GGT催化合成γ-Glu-Phe的催化反应过程中,优化了反应体系中底物浓度配比、反应温度、初始反应pH以及酶量等催化条件,提高了γ-Glu-Phe产率。
2.5.1 底物浓度对GGT催化合成γ-Glu-Phe的影响
由于提高受体浓度有利于降低GGT水解副反应,优化底物浓度配比可以提高γ-Glu-Phe产率。控制反应温度为40 ℃,反应体系pH 10.0,其他催化条件不变,控制供体底物50 mmol/L L-Gln与受体底物L-Phe浓度配比为1∶1、1∶2、1∶3、1∶4、1∶5、1∶6,进行体外转化。其结果如图6-a所示,当供体底物50 mmol/L L-Gln与受体底物L-Phe浓度配比为1∶4时,产物产量达到最高,产量为10.62 g/L,产率达到70%。
2.5.2 温度对GGT催化合成γ-Glu-Phe的影响
温度可以影响酶促反应速率,从而影响酶催化合成γ-Glu-Phe。控制供体底物50 mmol/L L-Gln与受体底物L-Phe浓度配比为1∶4,反应体系pH 10.0,其他催化条件一致的情况下,分别在25、30、35、40、45 ℃ 下进行酶催化反应合成γ-Glu-Phe产物。结果如图6-b所示,当反应温度为35 ℃时,产物产量最高,达到11.98 g/L,产率达到79%。
2.5.3 初始pH对GGT催化合成γ-Glu-Phe的影响
控制供体底物50 mmol/L L-Gln与受体底物L-Phe 浓度配比为1∶4,反应温度为35 ℃,其他催化条件一致的情况下,设置转化体系初始pH 6.0、7.0、8.0、9.0、10.0、11.0、12.0条件下进行酶催化反应合成γ-Glu-Phe产物。结果如图6-c所示,在初始反应pH<10.0时,产物产率随pH的升高而升高。当pH为10.0左右时,产物产量达到最高,产量为12.14 g/L,产率为80%,在初始反应pH>10.0时,产物产量大幅度下降。这是由于GGT最适转肽酶活力pH=10.0,当pH>10.0时,GGT转肽活性降低[10],从而其催化合成产物γ-glu-Phe的能力下降,导致产物产率下降。
2.5.4 酶浓度对于GGT催化合成γ-Glu-Phe的影响
实际工业化应用中,为了最大程度提高生产效率和降低生产成本,将催化反应体系中酶浓度控制在适宜的水平也是至关重要的。控制其他催化条件相同,采用供体底物50 mmol/L L-Gln与受体底物200 mmol/L L-Phe,将反应温度设为35 ℃,反应初始pH 10.0,酶浓度分别为0.2、0.6、1、1.4、1.8 U/mL 条件下进行酶催化反应合成γ-Glu-Phe产物。结果如图6-d所示,当反应体系中酶浓度<1 U/mL时,底物转化率以及产物产率随酶浓度增加而逐渐升高,当酶量在1 U/mL左右时,产物产量达到最高,产量为14.35 g/L,产率达到94%。当GGT浓度过高时,会发生自转肽反应,从而生成谷氨酸等副产物[3]。所以,当酶量>1 U/mL时,γ-Glu-Phe产物产率大幅度下降。
a-底物浓度配比对GGT催化合成γ-Glu-Phe的影响;b-温度对GGT催化合成γ-Glu-Phe的影响;c-pH对GGT催化合成γ-Glu-Phe的影响;
d-酶浓度对GGT催化合成γ-Glu-Phe的影响
图6 催化条件对GGT催化合成γ-Glu-Phe的影响
Fig.6 Influence of catalytic conditions on the synthesis of γ-Glu-Phe by GGT
利用传统的HPLC方法表征单位体积酶的产量费时费力,基于高通量筛选方法能够快速准确初步筛选酶的高产量菌株。本研究基于GGT催化合成γ-Glu-Phe会释放的特性,建立了基于定量显色的高通量筛选方法,该方法能够准确、快速以及高效定量体系中浓度,从而表征GGT高产菌株。启动子作为一种生物学元件对于调节基因的转录水平发挥着重要的作用。基于高通量筛选技术,对启动子HpaII的-10区上游7个碱基进行非理性改造,获得了突变菌株,较原始菌株的GGT表达量提高了63%。在对启动子进行突变后,找到与该启动子更加适配的RBS序列可以进一步提高目的基因表达量。我们利用RBS Calculator建立了由500个RBS序列组成的文库,通过高通量筛选获得了与启动子EF3HpaII更加适配的RBS362,获得的突变菌株GGT的表达量是原始菌株的2.02倍。因此,本文建立的高通量筛选方法可以准确高效获得GGT高表达量菌株。最后,我们对GGT体外催化合成γ-Glu-Phe的反应条件进行了探索优化。最终,γ-Glu-Phe产率达到94%。以上研究成果为GGT在食品行业进行工业化生产应用奠定了理论和实验基础。但在研究中发现,GGT转肽酶最适pH为10.0,同时GGT比酶活力仍需要提高,为了进一步提高其工业化应用价值,我们后续拟基于该高通量筛选方法,尝试改造GGT最适pH和比酶活力等特性,进一步改善GGT工业属性。
[1] LI Z L, ZHU R T, LIU Y Q, et al.γ-Glutamyltranspeptidase from Bacillus amyloliquefaciens:Transpeptidation activity enhancement and L-theanine production[J].Enzyme and Microbial Techology, 2020, 140:109644.
[2] YAO T T, QIAN J D, WANG G Q.Efficacy of ursodeoxycholic acid combined with prednisolone and immunosuppressant triple therapy in the treatment of refractory primary biliary cholangitis[J].Medicina Clínica, 2020, 155(4):165-170.
[3] 帅玉英, 张涛, 江波, 等.γ-谷氨酰转肽酶高产菌株的筛选、鉴定及其产酶条件优化[J].食品与发酵工业, 2010, 36(3):20-24.
SHUAI Y Y, ZHANG T, JIANG B, et al.Isolating of high GGT-producing strain and optimizing of enzyme formation conditions[J].Food and Fermentation Industries, 2010, 36(3):20-24.
[4] FUJII H, DOI H, KO T, et al.Frequently abnormal serum gamma-glutamyl transferase activity is associated with future development of fatty liver:A retrospective cohort study[J].BMC Gastroenterology, 2020, 20(1):213-217.
[5] SUZUKI H, KATO K, KUMAGAI H.Development of an efficient enzymatic production of gamma-D-glutamyl-L-tryptophan (SCV-07), a prospective medicine for tuberculosis, with bacterial gamma-glutamyltranspeptidase[J].Journal of Biotechnology, 2004, 111(3):291-295.
[6] SUZUKI H, KAJIMOTO Y, KUMAGAI H.Improvement of the bitter taste of amino acids through the transpeptidation reaction of bacterial gamma-glutamyltranspeptidase[J].Journal of Agricultural and Food Chemistry, 2002, 50(2):313-318.
[7] ZHU X P, SUN-WATERHOUSE D X, CHEN J H, et al.Comparative study on the novel umami-active peptides of the whole soybeans and the defatted soybeans fermented soy sauce[J].Journal of the Science of Food and Agriculture, 2021, 101(1):158-166.
[8] PLOSS T N, REILMAN E, MONTEFERRANTE C G, et al.Homogeneity and heterogeneity in amylase production by Bacillus subtilis under different growth conditions[J].Microbial Cell Factories, 2016, 15:57.
[9] MU W M, ZHANG T, JIANG B.An overview of biological production of L-theanine[J].Biotechnology Advances, 2015, 33(3-4):335-342.
[10] YANG T W, IRENE K, LIU H L, et al.Enhanced extracellular gamma glutamyl transpeptidase production by overexpressing of PrsA lipoproteins and improving its mRNA stability in Bacillus subtilis and application in biosynthesis of L-theanine[J].Journal of Biotechnology, 2019, 302:85-91.
[11] LI X, XU S Q, ZHANG X, et al.Design of a high-efficiency synthetic system for L-asparaginase production in Bacillus subtilis[J].Engineering in Life Sciences, 2019, 19(3):229-239.
[12] LIEBETON K, LENGEFELD J, ECK J.The nucleotide composition of the spacer sequence influences the expression yield of heterologously expressed genes in Bacillus subtilis[J].Journal of Biotechnology, 2014, 191(2):214-220.
[13] LIU M F, TOLSTORUKOV M, ZHURKIN V, et al.A mutant spacer sequence between-35 and-10 elements makes the P-lac promoter hyperactive and cAMP receptor protein-independent[J].Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(18):6 911-6 916.
[14] HAN L C, CUI W J, SUO F Y, et al.Development of a novel strategy for robust synthetic bacterial promoters based on a stepwise evolution targeting the spacer region of the core promoter in Bacillus subtilis[J].Microbial Cell Factories, 2019, 18(1):96-101.
[15] BENTELE K, SAFFERT P, RAUSCHER R, et al.Efficient translation initiation dictates codon usage at gene start[J].Molecular Systems Biology, 2013, 9(2):675-682.
[16] ZHANG B, ZHOU N, LIU Y M, et al.Ribosome binding site libraries and pathway modules for shikimic acid synthesis with Corynebacterium glutamicum[J].Microbial Cell Factories, 2015, 14(1):1-14.
[17] CHAPPELL J, WESTBROOK A, VEROSLOFF M, et al.Computational design of small transcription activating RNAs for versatile and dynamic gene regulation[J].Nature Communications, 2017, 8:1051.
[18] NA D, LEE S, LEE D.Mathematical modeling of translation initiation for the estimation of its efficiency to computationally design mRNA sequences with desired expression levels in prokaryotes[J].BMC Systems Biology, 2010, 4(1):71.
[19] 吴傲, 张显, 徐美娟, 等.天蓝色链霉菌海藻糖合酶的异源表达、活性分析及重组菌全细胞转化合成海藻糖的条件优化[J].生物工程学报, 2019, 35(7):1 348-1 358.
WU A, ZHANG X, XU M J, et al. Heterologous expression and activity analysis of trehalose synthase from streptomyces coelicolor and production of trehalose by a recombinant E.coli BL21 whole-cell biocatalyst[J].Chinese Journal of Biotechnology, 2019, 35(7):1 348-1 358.
[20] 汪前, 姚忠, 荀志金, 等.γ-谷氨酰转肽酶的酶学性质及其转肽反应机制[J].高校化学工程学报, 2008,22(2):288-293.
WANG Q, YAO Z, XUN Z J, et al.Properties and catalytic mechanism of γ-Glutamyltranspeptidase from B.subtilis NX-2[J].Journal of Chemical Engineering of Chinese Universities, 2008,22(2):288-293.
[21] VELLANOWETH R L, RABINOWITZ J C.The influence of ribosome-binding-site elements on translational efficiency in Bacillus subtilis and Escherichia coli in vivo[J].Molecular Microbiology, 2010, 6(9):1 105-1 114.
[22] ISAACS F J, DWYER D J, DING C, et al.Engineered riboregulators enable post-transcriptional control of gene expression[J].Nature Biotechnology, 2004, 22(7):841-847.