樊占青，劉曉婷，魏 貞，王哲人，王亞南，高 欣，閔偉紅

（吉林農(nóng)業(yè)大學(xué)食品科學(xué)與工程學(xué)院，小麥和玉米深加工國家工程實(shí)驗(yàn)室，吉林 長春 130118）

甲硫氨酸（Met）是一種不能在體內(nèi)合成而必須由外界獲取的人體必需氨基酸[1]，參與生物體內(nèi)的新陳代謝和蛋白質(zhì)合成，是多種代謝物質(zhì)的合成前體[2]。目前微生物發(fā)酵法是最有前景生產(chǎn)Met的方法[3]，但是由于生物體內(nèi)代謝路徑多樣復(fù)雜，各種代謝物質(zhì)彼此調(diào)控，所以該方法尚未完全實(shí)現(xiàn)工業(yè)化生產(chǎn)[4]。在生物體內(nèi)，天冬氨酸經(jīng)過一系列代謝合成Met，而天冬氨酸激酶（aspartate kinase，AK）是該代謝途徑中的第一關(guān)鍵限速酶[5]，其活性同時(shí)受賴氨酸（Lys）和蘇氨酸（Thr）的協(xié)同反饋抑制[6]，導(dǎo)致碳流競爭激烈，Met合成供應(yīng)減少，使得最終代謝產(chǎn)物Met不能大量積累[7-8]。目前研究發(fā)現(xiàn)，不同生物體內(nèi)AK具有不同的存在形式，擬南芥（Arabidopsis）AK具有5 種形式，包括3種單功能AKs（AKI、AKII和AKIII）和2 種雙功能AK-高絲氨酸脫氫酶（homoserine dehydrogenase，HSDH）[9-10]；大腸桿菌中含有1 種單功能AKIII和2 種雙功能AK-HSDHI、AK-HSDHII[11-12]，此外，在一些生物體中也存在著更簡單的抑制機(jī)制，詹氏酵母和嗜熱菌中只含有一種形式的AK[13]。本課題組首次發(fā)現(xiàn)北京棒桿菌天冬氨酸激酶（Corynebacterium pekinenseaspartate kinase，CpAK）是一種新型單體別構(gòu)酶，且仍受下游產(chǎn)物L(fēng)ys和Thr的協(xié)同反饋抑制[5]。

由于CpAK在代謝途徑中被下游產(chǎn)物協(xié)同反饋抑制，導(dǎo)致下游產(chǎn)物天冬氨酸族氨基酸產(chǎn)量不能大量積累，因此通過削弱或解除反饋抑制提高最終代謝產(chǎn)物產(chǎn)量成為研究熱點(diǎn)。目前，基因敲除[14]、基因過表達(dá)[15-17]、定點(diǎn)突變[18-23]等基因工程技術(shù)是改造AK的主要技術(shù)手段，其中，定點(diǎn)突變技術(shù)的應(yīng)用比較成熟廣泛。Han Caijing等[5]利用定點(diǎn)突變技術(shù)削弱了下游產(chǎn)物L(fēng)ys對(duì)CpAK的反饋抑制，獲得酶活力提高11.32 倍的突變體A380I，為Met產(chǎn)量積累提供了理論依據(jù)。也有學(xué)者對(duì)抑制劑Lys和Thr周圍結(jié)合位點(diǎn)和底物結(jié)合位點(diǎn)的周圍關(guān)鍵殘基進(jìn)行突變研究，得到了AK活力顯著提高且反饋抑制作用明顯減弱的突變株R169P、R169Y、T379L和M372I/T379W，為強(qiáng)化AK代謝途徑中下游產(chǎn)物碳流量的積累、構(gòu)建高產(chǎn)Met工程菌株提供了理論依據(jù)和技術(shù)支持[16-17,24]。

目前文獻(xiàn)報(bào)道對(duì)AK的突變改造主要圍繞Lys和Thr抑制劑周圍關(guān)鍵位點(diǎn)[20-26]，而對(duì)催化活性中心ATP以及底物結(jié)合位點(diǎn)的研究相對(duì)較少。AK的一般抑制機(jī)制已從先前研究中得出結(jié)論：效應(yīng)子結(jié)合誘導(dǎo)從底物-ATP結(jié)合的松弛狀態(tài)（R狀態(tài)，緊密構(gòu)象）過渡到無效的緊張狀態(tài)（T狀態(tài)，開放構(gòu)象），因此底物和ATP結(jié)合位點(diǎn)的構(gòu)象變化均會(huì)阻礙催化作用[27]。Li Changcheng等[28]對(duì)銅綠假單胞菌AK催化域中ATP結(jié)合位點(diǎn)D182和R184進(jìn)行雙突變改造，獲得酶活力降低60%的雙突變株D182/R184A，證明ATP結(jié)合口袋的改變對(duì)AK活力的重要性。本研究選擇AK催化中心ATP周圍5 ?范圍內(nèi)關(guān)鍵殘基絲氨酸（Ser227）位點(diǎn)，進(jìn)行四突變株的構(gòu)建，篩選獲得高酶活力、性狀良好的突變株，并結(jié)合分子動(dòng)力學(xué)模擬，對(duì)其機(jī)制進(jìn)行分析，旨在為優(yōu)化AK代謝途徑提供參考。

1 材料與方法

1.1 材料與試劑

1.1.1 菌株、質(zhì)粒與培養(yǎng)基

北京棒桿菌E31（吉林農(nóng)業(yè)大學(xué)發(fā)酵工程實(shí)驗(yàn)室保藏）提供本實(shí)驗(yàn)AK基因；重組大腸桿菌突變株（pET-28a-T379N/A380C/G171I-AK）由吉林農(nóng)業(yè)大學(xué)發(fā)酵工程實(shí)驗(yàn)室提供；表達(dá)宿主大腸桿菌BL21（DE3）購自北京鼎國昌盛生物技術(shù)有限責(zé)任公司。

本實(shí)驗(yàn)選用的培養(yǎng)基均為LB培養(yǎng)基。液體LB培養(yǎng)基：0.5%酵母浸粉、1% NaCl、1%蛋白胨。固體培養(yǎng)基則在液體培養(yǎng)基中加入2%的瓊脂粉，溶于蒸餾水中，適宜濃度的NaOH溶液將其pH值調(diào)至7.0，121 ℃高壓蒸汽滅菌20 min。

1.1.2 試劑

十二烷基硫酸鈉-聚丙烯酰氨凝膠電泳（sodium dodecyl sulfate-polyacrylamide gelelectrophoresis，SDSPAGE）凝膠試劑盒、卡那霉素（kana）、丙基硫代-β-D-半乳糖苷（isopropyl-β-D-thiogalactopyranoside，IPTG）北京鼎國昌盛生物技術(shù)有限公司；DpnI消化酶、LaTaq酶TaKaRa（大連）公司；質(zhì)粒DNA抽提試劑盒 生工生物工程（上海）股份有限公司；非變性鎳柱及柱料美國GE公司；本實(shí)驗(yàn)所用引物由吉林省庫美生物科技有限公司合成，基因測序由生工生物工程（上海）股份有限公司完成。

1.2 儀器與設(shè)備

梯度聚合酶鏈?zhǔn)椒磻?yīng)（polymerase chain reaction，PCR）儀 德國Eppendorf AG公司；DYY-6C核酸電泳儀北京六一儀器廠；瓊脂糖凝膠成像儀 德國耶拿分析儀器股份公司；Z36HK臺(tái)式低溫高速冷凍離心機(jī) 德國Hermle公司；SE260蛋白電泳儀 美國GE公司；SPECTRA Max酶標(biāo)儀 美國Molecular Devices公司；WFH-201BJ蛋白印跡半干轉(zhuǎn)膜儀 美國Bio-Rad公司；Scientz超聲細(xì)胞粉碎儀 寧波新芝生物公司。

1.3 方法

1.3.1 定點(diǎn)飽和突變

首先通過Primer 5.0軟件進(jìn)行突變位點(diǎn)的引物設(shè)計(jì)，上游引物：5’-CTGGAACTTGCTGCTGTTGGCNNNAAG ATTTTGG-3’；下游引物：5’-GCGCAGCACCAAAATCT TNNNGCCAACAGC-3’。由吉林省庫美生物科技有限公司合成。

根據(jù)質(zhì)粒DNA抽提試劑盒說明書提取T379N/A380C/G171I-AK的質(zhì)粒，以其為模板進(jìn)行梯度PCR擴(kuò)增全質(zhì)?；颍磻?yīng)條件為：94 ℃預(yù)變性1 min；94 ℃變性1 min，56 ℃退火1 min，72 ℃延伸10 min（18個(gè)循環(huán)）；72 ℃延伸10 min；4 ℃保存。

1.3.2 轉(zhuǎn)宿主大腸桿菌BL21

1.3.1節(jié)完成PCR擴(kuò)增，將產(chǎn)物用1%瓊脂糖核酸電泳驗(yàn)證擴(kuò)增成功。PCR突變產(chǎn)物用DpnI酶消化，反應(yīng)體系為20 μL，依次加入10×T Buffer 2 μL、PCR產(chǎn)物2 μL、DpnI酶0.3 μL，最后用ddH2O補(bǔ)足，金屬浴37 ℃恒溫2 h。

消化結(jié)束，將2～4 μL產(chǎn)物轉(zhuǎn)入120 μL大腸桿菌BL21感受態(tài)中，冰浴30 min，42 ℃金屬浴熱激90 s，繼續(xù)冰浴5 min后加入900 μL不含kana的LB液體培養(yǎng)基，37 ℃、170 r/min搖床擴(kuò)大培養(yǎng)2 h，5 000 r/min離心5 min后棄掉800 μL上清液，用剩余上清重懸菌體，將其均勻涂布于LB固體培養(yǎng)基（含kana），37 ℃倒置培養(yǎng)12～16 h。

1.3.3 高通量篩選及測序驗(yàn)證

挑取成功轉(zhuǎn)入的單菌落于含kana的LB液體培養(yǎng)基中，180 r/min、37 ℃擴(kuò)大培養(yǎng)12 h，加入1 mmol/L IPTG，適宜溫度、130 r/min誘導(dǎo)表達(dá)8 h后進(jìn)行高通量篩選，方法見參考文獻(xiàn)[5]，挑取酶活力顯著提高的突變株加入10 mL試管培養(yǎng)基（含kana）于搖床180 r/min、37 ℃過夜擴(kuò)大培養(yǎng)。將擴(kuò)大培養(yǎng)后的突變株以菌液為模板，用克隆引物（AK引物，上游引物5’-GAATTCCATATGGC CCTGGTCGTACAGAA-3’，下游引物5’-GGAATTCTTA GCGTCCGGTGCCTGCAT-3’）進(jìn)行PCR擴(kuò)增驗(yàn)證，反應(yīng)條件為：94 ℃預(yù)變性1 min，94 ℃變性30 s，58 ℃退火1 min，72 ℃延伸90 s，循環(huán)30 次，4 ℃保存[20]。將PCR產(chǎn)物進(jìn)行1%瓊脂糖核酸電泳驗(yàn)證，挑取目的條帶明亮的菌株送到生工生物工程（上海）股份有限公司，進(jìn)行全質(zhì)粒基因測序驗(yàn)證。

1.3.4 AK分離純化及蛋白電泳驗(yàn)證

測序驗(yàn)證成功的突變株以2%接種量于試管培養(yǎng)基活化12 h，吸取2 mL活化的菌液加入到100 mL LB液體培養(yǎng)基培養(yǎng)1.5～2 h后，加入終濃度為1 mmol/L IPTG誘導(dǎo)表達(dá)8 h，于低溫高速冷凍離心機(jī)8 000 r/min離心，將上清液舍棄后收集菌體。參照朱云明等[29]實(shí)驗(yàn)方法制備粗酶液及純化酶樣。用SDS-PAGE和Western Blot對(duì)粗酶樣及純化樣進(jìn)行驗(yàn)證，方法參考文獻(xiàn)[30]。

1.3.5 AK活力測定及動(dòng)力學(xué)分析

AK活力以單位時(shí)間內(nèi)單位質(zhì)量的AK催化反應(yīng)速率表示，單位為U/（mg·min），酶活力表示為1 000×A540nm。以L-天冬氨酸（L-Asp）為底物，在1 mL反應(yīng)體系中加入酶樣純化液，1 mL反應(yīng)體系為800 mmol/L KCl、10 mmol/LL-Asp、100 mmol/L Tris-HCl、800 mmol/L NH4OH、10 mmol/Lβ-巰基乙醇、10.4 mmol/L ATP、1.6 mmol/L MgSO4[25]，保持反應(yīng)體系以及其他測定條件不變，改變底物濃度（0.5、1、3、5、7、9、10、12、14、16 mmol/L）反應(yīng)30 min后，于540 nm波長處測其吸光度A540nm，通過Hill方程進(jìn)行非線性擬合計(jì)算AK活力。

1.3.6 AK酶學(xué)性質(zhì)分析

1.3.6.1 最適溫度

反應(yīng)條件及反應(yīng)體系不變，改變反應(yīng)溫度（15、20、25、26、28、30、35、40、45、50 ℃），反應(yīng)30 min后終止反應(yīng)，測定酶活力，每個(gè)溫度檢測3 組平行，最高酶活力定義為100%。

1.3.6.2 最適pH值

反應(yīng)條件及反應(yīng)體系不變，改變反應(yīng)體系中Tris-HCl緩沖液pH值（6.0、6.5、7.0、7.5、8.0、8.5、9.0、9.5、10.0），反應(yīng)30 min后終止反應(yīng)，測定酶活力，每個(gè)pH值檢測3 組平行，最高酶活力定義為100%。

1.3.6.3 穩(wěn)定性

反應(yīng)體系不變，選擇最適溫度和最適pH值，將酶樣純化液加入體系中，130 r/min搖床反應(yīng)，間隔1 h取樣一次測酶活力，每組樣品檢測3 組平行，將0 h酶活力定義為100%。

1.3.6.4 底物抑制劑對(duì)AK活力的影響

反應(yīng)體系不變，分別加入不同濃度（終濃度分別為0.2、1.0、5.0、10.0 mmol/L）的底物抑制劑，分別為Lys、Thr、Met、Lys+Thr、Lys+Met、Thr+Met、Lys+Thr+Met，測定抑制劑對(duì)AK活力的影響，對(duì)照組加入等量蒸餾水每組樣品檢測3 組平行，將對(duì)照組相對(duì)酶活力定義為100%。

1.3.7 分子動(dòng)力學(xué)模擬分析

以3aaw和2hmf為結(jié)構(gòu)模板，同源模建獲得AK三維結(jié)構(gòu)，建模后丟失的原子通過VMD方法補(bǔ)充，同時(shí)通過局部的結(jié)構(gòu)檢查確定His的質(zhì)子態(tài)；借助VMD的solvate plugin程序添加水溶液環(huán)境，并添加Na+或Cl-平衡體系，力場參數(shù)使用charmm36。將篩選出的高酶活力突變株T379N/A380C/G171I/S227D進(jìn)行不少于110 ns的NPT自由模擬（分子動(dòng)力學(xué)），分析其酶活力提高機(jī)制。

在前期實(shí)驗(yàn)室構(gòu)建的PDB結(jié)構(gòu)文件基礎(chǔ)上，通過Discovery Studio軟件分析Ser227位點(diǎn)突變前后與ATP及周圍殘基的作用力及殘基間距離的變化。

1.4 數(shù)據(jù)分析

2 結(jié)果與分析

2.1 ATP結(jié)合位點(diǎn)的篩選

通過Discovery Studio軟件對(duì)AK模型分析，Ser227位點(diǎn)是ATP附近5 ?范圍內(nèi)的關(guān)鍵殘基；并且利用Clustal X軟件進(jìn)行同源序列比對(duì)，Ser227位點(diǎn)是保守位點(diǎn)，如圖1所示。Han Caijing等[5]對(duì)CpAK突變體A380I進(jìn)行100 ns分子動(dòng)力學(xué)模擬，表明突變使A380I提高了ATP結(jié)合口袋處殘基Ieu220～Glu270的均方根偏差（root mean square deviation，RMSD）值，與ATP交互作用變強(qiáng)，同時(shí)突變使ATP（N5）- Ser227（OG）中氫鍵占用率從31.74%增加到39.14%，ATP和CpAK之間氫鍵的占用比增多，增加了ATP和CpAK之間的結(jié)合，從而增強(qiáng)了CpAK的催化活性。因此本研究選擇Ser227位點(diǎn)，設(shè)想通過飽和突變使其與ATP形成氫鍵，增強(qiáng)ATP與CpAK的結(jié)合，達(dá)到提高AK活力的預(yù)期效果。

圖1 ATP結(jié)合位點(diǎn)的篩選Fig. 1 Screening of ATP binding sites

2.2 突變株構(gòu)建

以T379N/A380C/G171I（含有質(zhì)粒pET-28a-AK）為模板，在Ser227引物作用下PCR擴(kuò)增進(jìn)行定點(diǎn)飽和突變。PCR產(chǎn)物經(jīng)1%核酸電泳驗(yàn)證結(jié)果如圖2A所示，在7 000 bp處有明亮條帶，初步證明成功突變。取定量PCR產(chǎn)物經(jīng)DpnI酶消化后轉(zhuǎn)入大腸桿菌感受態(tài)細(xì)胞（BL21）中，經(jīng)過高通量篩選后獲得高酶活力突變株。將突變株活化后進(jìn)行菌液PCR驗(yàn)證，如圖2B所示，在1 000～2 000 bp處有明亮單一條帶，與目的基因CpAK長度1 266 bp相符，進(jìn)一步證明目的基因成功大量復(fù)制。將篩選的高酶活力突變株送去測序，結(jié)果如圖2C所示，Ser227位點(diǎn)由殘基Ser突變?yōu)闅埢鵄sp，最終證明突變株T379N/A380C/G171I/S227D構(gòu)建成功。

2.3 SDS-PAGE和Western Blot驗(yàn)證

圖3 SDS-PAGE（A）及Western Blot（B）結(jié)果Fig. 3 Results of SDS-PAGE (A) and Western Blot (B)

WT菌株和突變株經(jīng)IPTG誘導(dǎo)表達(dá)目的蛋白，分別進(jìn)行非變性鎳柱純化，將得到的粗酶樣和純化樣進(jìn)行SDS-PAGE和Western Blot驗(yàn)證，結(jié)果如圖3所示。在48 kDa處有明顯條帶，證明AK蛋白表達(dá)成功。Western Blot驗(yàn)證結(jié)果表明，AK蛋白與抗體（His-tag27E8 Mouse mAb）特異性結(jié)合，純化樣是單一AK蛋白。

2.4 酶動(dòng)力學(xué)分析

以不同濃度L-Asp為底物，對(duì)WT和突變株進(jìn)行酶動(dòng)力學(xué)分析，測定其酶活力。測定結(jié)果用Origin 9.0軟件進(jìn)行非線性擬合，根據(jù)Hill方程V=VmaxSn/（Kn+Sn）得到WT和突變株的最大催化速率Vmax。動(dòng)力學(xué)結(jié)果見圖4、表1，該突變株T379N/A380C/G171I/S227D的Km值減小為1.35 mmol/L，底物親和力增大，利于AK與底物的結(jié)合，從而提高其催化效率；Vmax值為242.05 U/（mg·min），較T379N/A380C/G171I（187.88 U/（mg·min））提高到1.28 倍，較WT AK（ 3.01 U/（mg·min））提高到80.41 倍，表明突變后AK的催化活力顯著提高。

圖4 WT和突變株的動(dòng)力學(xué)分析Fig. 4 Kinetic analysis of WT and mutant strains

表1 WT AK和突變株的動(dòng)力學(xué)參數(shù)Table 1 inetic parameters of WT AK and mutant strains

2.5 酶學(xué)性質(zhì)表征

2.5.1 最適溫度

圖5 WT和突變株的酶學(xué)性質(zhì)表征Fig. 5 Characterization of enzymatic properties of WT and mutant strains

由圖5A可知，WT AK最適催化溫度25 ℃，T379N/A380C/G171I AK最適催化溫度28 ℃，該突變株T379N/A380C/G171I/S227D最適催化溫度30 ℃，較WT和T379N/A380C/G171I分別提高5 ℃和2 ℃，表明突變?cè)鰪?qiáng)了AK的耐熱性，有效提高了AK的最適溫度。當(dāng)溫度高于30 ℃后，該突變株的酶活力仍然較高于WT和T379N/A380C/G171I，說明多次突變有利于AK的耐高溫性能的優(yōu)化。

2.5.2 最適pH值

由圖5B可知，WT AK最適pH值為8.0，T379N/A380C/G171I AK最適pH值為9.0，該突變株T379N/A380C/G171I/S227D最適pH值為8.5。與T379N/A380C/G171I相比，該突變株最適pH值降低0.5，耐酸性提高，同時(shí)，pH值耐受范圍變大，這可能是Ser227位點(diǎn)由中性氨基酸Ser突變?yōu)樗嵝园被酇sp的緣故。在pH值為10時(shí)，該突變株的酶活力仍然保持在32%，較高于WT和T379N/A380C/G171I。

2.5.3 AK穩(wěn)定性

相對(duì)酶活力下降到最初活力50%所用時(shí)間為酶的半衰期。選擇最適溫度和最適pH值，分別研究WT和突變株的穩(wěn)定性。如圖5C所示，WT AK的半衰期為4.9 h，而T379N/A380C/G171I AK的半衰期為3.1 h，該突變株T379N/A380C/G171I/S227D的半衰期為3.9 h。與WT AK相比，該突變體AK在反應(yīng)初（0～3 h）和反應(yīng)末（7～10 h）相對(duì)酶活力顯著高于WT AK。與T379N/A380C/G171I AK相比，該突變株AK半衰期延長0.8 h，且相對(duì)酶活力幾乎始終高于T379N/A380C/G171I AK。

2.5.4 底物抑制劑對(duì)AK活力的影響

由表2可知，WT AK受Thr和Lys兩種抑制劑的協(xié)同反饋抑制，而Met對(duì)其抑制作用相對(duì)較弱。隨著Thr和Lys濃度的增加，抑制作用明顯增強(qiáng)，當(dāng)Lys濃度增加到10 mmol/L時(shí)，對(duì)WT AK的抑制率大約達(dá)到39%。當(dāng)兩種抑制劑組合濃度為10 mmol/L時(shí)，WT AK抑制率最高為45.04%。當(dāng)3種抑制劑組合濃度為10 mmol/L時(shí)，WT AK 抑制率僅比Lys+Thr組合時(shí)高4.15%，再次說明AK受Met抑制作用相對(duì)較弱，而主要受Thr和Lys的抑制，這與先前Yoshida等[31]研究一致。突變后，AK的被抑制作顯著減弱，甚至被明顯激活。在10 mmol/L Thr+Met、Lys+Met和Thr+Lys+Met不同組合的情況下，AK的催化活性被顯著激活，T379N/A380C/G171I AK的相對(duì)酶活力分別為98.81%、100.14%和108.46%，而該突變體相對(duì)活力提高更為顯著，分別為109.37%、121.25%和143.35%，且呈現(xiàn)明顯的濃度計(jì)量依賴，這與2.4節(jié)動(dòng)力學(xué)分析結(jié)果（Km減小，Vmax增大）相對(duì)應(yīng)，進(jìn)一步證明本研究位點(diǎn)Ser227是影響AK活力的關(guān)鍵殘基位點(diǎn)，為后期工程菌的構(gòu)建提供理論依據(jù)。

表2 不同底物抑制劑條件下WT和突變株的相對(duì)酶活力Table 2 Relative enzyme activity of WT and mutant strains under different substrate inhibitor conditions%

2.6 突變體的空間結(jié)構(gòu)分析

以構(gòu)建的AK模型為基礎(chǔ)，通過Discovery Studio軟件對(duì)Ser227位點(diǎn)突變前后與ATP和其周圍殘基的作用力以及距離變化進(jìn)行分析。如圖6A、B所示，突變前Ser227殘基與ATP及周圍氨基沒有直接作用，突變后，S227D殘基側(cè)鏈向ATP延伸，與ATP和I229殘基產(chǎn)生氫鍵，氫鍵占有率增加，連接更加緊密，更容易促進(jìn)催化反應(yīng)發(fā)生，這與Han Caijing等[5]研究報(bào)道一致；如圖6C、D所示，突變前Ser227殘基與ATP距離為2.953 ?，突變后距離縮短為2.562 ?。突變后，距離縮短，作用力增強(qiáng)，從而提高酶活力。

圖6 Ser227位點(diǎn)突變前后成鍵變化Fig. 6 Bonding changes before and after mutation at Ser227

2.7 分子動(dòng)力學(xué)模擬

對(duì)該突變體T379N/A380C/G171I/S227D AK+Asp+ATP+Lys+Thr復(fù)合物進(jìn)行分子動(dòng)力學(xué)模擬研究，揭示其活性提高機(jī)制。如圖7A所示，復(fù)合物主鏈C骨架RMSD平均值為0.90 nm，體系在19 ns處達(dá)到平衡界點(diǎn)，之后一直平衡穩(wěn)定于0.9 nm左右，說明體系在整個(gè)模擬過程中始終保持穩(wěn)定狀態(tài)[32]。Rg是在模擬過程中根據(jù)整個(gè)蛋白體積和性狀計(jì)算出的蛋白回旋半徑，反映了體系在模擬過程中的柔韌性，如圖7B所示，復(fù)合物在整個(gè)模擬過程中基本平穩(wěn)，平均Rg收斂于26.66 ?，整個(gè)體系柔韌性良好，這可能是由于突變使得蛋白結(jié)構(gòu)發(fā)生重排趨于穩(wěn)定狀態(tài)。蛋白表面溶劑可及表面積（solvent accessible surface area，SASA）是描述蛋白質(zhì)親疏水性的重要參數(shù)[33]，對(duì)蛋白的空間構(gòu)象及其穩(wěn)定性具有重要意義，由圖7C可知，復(fù)合物的模擬軌跡沒有較大波動(dòng)，一直保持平穩(wěn)，整個(gè)模擬過程平均SASA值為222.56 nm2，表明突變使得T379N/A380C/G171I/S227D 的親水性增強(qiáng)，蛋白表面溶劑的親和能力增強(qiáng)，從而影響其催化活性，這可能是由于突變后Asp殘基側(cè)鏈像ATP延伸，同時(shí)形成氫鍵。以上3個(gè)參數(shù)都是在蛋白質(zhì)整體上進(jìn)行分析研究，而均方根波動(dòng)（root mean square fluction，RMSF）是從整個(gè)模擬催化過程單個(gè)氨基酸殘基的相對(duì)波動(dòng)情況進(jìn)行分析[34]。由圖7D可知，整個(gè)模擬過程中，復(fù)合物體系都表現(xiàn)出很好的靈活性，利于催化反應(yīng)。

圖7 T379N/A380C/G171I/S227D AK復(fù)合物分子動(dòng)力學(xué)模擬分析Fig. 7 MDS analysis of T379N/A380C/G171I/S227D AK complex

3 討 論

大量研究報(bào)道對(duì)AK的抑制劑結(jié)合位點(diǎn)進(jìn)行突變，可以有效減弱或解除抑制作用，顯著提高酶活力。Gao Yunna等[24]對(duì)抑制劑Lys周圍結(jié)合位點(diǎn)M372和T379關(guān)鍵殘基進(jìn)行定點(diǎn)突變，獲得酶活力提高16.51 倍的雙突變體M372I/T379W AK，并通過分子動(dòng)力學(xué)模擬對(duì)其機(jī)制進(jìn)行分析。Min Weihong等[19]等通過虛擬篩選獲得一個(gè)底物結(jié)合和催化反應(yīng)的關(guān)鍵殘基位點(diǎn)R169，通過定點(diǎn)突變技術(shù)，改造其空間結(jié)構(gòu)，得到了酶活力提高4.71 倍的突變株R169Y。表明底物結(jié)合和催化反應(yīng)位點(diǎn)的改造對(duì)AK活力提高有重要影響。魏貞等[35]研究通過對(duì)ATP周圍結(jié)合位點(diǎn)Y198和D201進(jìn)行定點(diǎn)飽和突變，獲得性能優(yōu)良、酶活力提高18.26 倍的雙突變株Y198N/D201M AK。因此本研究在前期研究的基礎(chǔ)上，繼續(xù)對(duì)ATP結(jié)合位點(diǎn)進(jìn)行突變改造，獲得酶活力顯著提高的突變體，進(jìn)一步表明ATP結(jié)合位點(diǎn)的重要性。

本研究推測T379N/A380C/G171I/S227D酶活力提高的原因是T379N/A380C/G171I突變株基礎(chǔ)上繼續(xù)突變，引起ATP結(jié)合口袋重排，同時(shí)Ser227位點(diǎn)與ATP由非鍵作用變?yōu)橹苯託滏I作用力，致使AK與ATP和底物結(jié)合更緊密，反應(yīng)體系更加穩(wěn)定，導(dǎo)致在催化過程中競爭性抑制劑不足以與底物ASP以及ATP競爭，從而減弱或解除抑制作用，顯著提高AK催化活力，使得碳流競爭減弱，更多的流向下游代謝途徑，進(jìn)而達(dá)到高產(chǎn)下游代謝產(chǎn)物的目的，為后期構(gòu)建高產(chǎn)蛋氨酸工程菌提供理論支持。

4 結(jié) 論

本研究對(duì)CpAK中催化域ATP周圍關(guān)鍵殘基Ser227位點(diǎn)進(jìn)行定點(diǎn)飽和突變，以達(dá)到減弱或解除抑制，酶活力顯著提高并獲得優(yōu)良多突變株的目的。通過軟件Primer Premier 5 進(jìn)行突變引物設(shè)計(jì)，以T379N/A380C/G171I AK為模板，經(jīng)1%核酸電泳驗(yàn)證，成功構(gòu)建突變株T379N/A380C/G171I/S227D。并將其轉(zhuǎn)入宿主大腸桿菌BL21表達(dá)，經(jīng)SDS-PAGE和Western Blot驗(yàn)證AK成功表達(dá)。動(dòng)力學(xué)分析結(jié)果表明，該突變株的Vmax較T379N/A380C/G171I提高到1.28 倍，較WT AK提高到80.41 倍，酶活力顯著提高，并且與底物的親和力也明顯增大，有利于與底物的結(jié)合。酶學(xué)性質(zhì)研究表明：該突變株T379N/A380C/G171I/S227D的最適反應(yīng)溫度增至30 ℃，較WT和T379N/A380C/G171I分別提高5 ℃和2 ℃；最適pH值為8.5，較T379N/A380C/G171I（pH 9.0）減小0.5；半衰期為3.9 h，較T379N/A380C/G171I延長0.8 h；且在不同濃度、不同組合抑制劑存在時(shí)，該突變株表現(xiàn)出明顯的激活作用，抑制效果減弱顯著，相對(duì)酶活力最高達(dá)143.35%，這與動(dòng)力學(xué)結(jié)果相符。同時(shí)分子動(dòng)力學(xué)模擬從個(gè)體和單個(gè)氨基酸殘基的變化，其中包括氫鍵的形成以及復(fù)合物系統(tǒng)中的分子間作用力，對(duì)酶活力提高進(jìn)一步說明，為構(gòu)建高產(chǎn)工程菌提供參考。