將 epTtTS 基因使用兩種限制酶進行切割(SalI 和 NdeI,Biolab) 以產生兩端不同的 sticky end,另一方面,為了減少背景值干擾,使 用載有不同基因之 pET-23b 質體進行相同限制酵素的切割,並操作 Gel elution kit 取得質體片段與 epTtTS 基因。將兩種 DNA 片段以莫 爾數 1 : 3 進行混合,並添加 T4 DNA Ligase (總反應體積的 1/20)與 10X Ligase buffer,放於 4℃ 隔夜。之後取 Ligation 樣本 5μL 進行
- 27 - 71.2%,麥芽糖生成仍只有 28.3%,海藻酮糖生成則提升至 87.8%,
但並無生成麥芽酮糖(圖二)。 催化效率(catalytic efficiency),即 Kcat 與 KM比值,顯示麥芽糖為 TtTS 最適作用受質,有最高的催化效率為 1.83 (圖三 H),額外添加 15mM
- 28 -
三、GI-coupled TtTS 反應分析
在GI-coupled TtTS 之反應,依方法七.加入 GI 的量為 TtTS量 (50μg)的2000倍(100mg),會使 TtTS 異構化作用明顯偏向水解反應 生成葡萄糖,比起沒有GI的組別,TtTS反應 30 分鐘產生的海藻糖
- 29 - 為已經解出結構的海藻糖合成酶,分別來自 Mycobacterium
smegmatis (PDB entry: 3ZO9,MsTS)、Mycobacterium tuberculosis (PDB entry: 4LXF,MtTS)、Deinococcus radiodurans(PDB entry:
4TVU,DrTS)其胺基酸序列與 TtTS 有 50%以上的一致性,另外依照 催化性質還有兩類酵素,一類是只有水解活性的水解酶(hydrolase),
如: Geobacillus HTA-462 α-glucosidase、Halothermothrix orenii α-amylase、Bacillus subtilis isomaltase;另一類是具有異構化 (isomerization)活性的蔗糖異構酶(sucrose isomerase,SI),如:
Pseudomonas mesoacidophila MX-45 trehalulose synthase(PDB entry :1ZJA)或 Protaminobacter rubrum isomaltulose synthase(PDB entry :3GBD)(表三)。其中與 TtTS 相似最高的蛋白 MsTS,因胺基酸 Leu344 佔據,阻擋了受質結合位,相似度次之的 MtTS 也因為其結 構缺少 S7 區域而使活性位呈現開放狀態,因此我們選擇相似度也極 高而活性位置呈現完整封閉的 DrTS 作為 TtTS 結構模擬模板,方便 進行酵素與受質作用關係的探討(附錄六)。利用 Swiss-Model 進行 TtTS 蛋白質三級結構模擬,以 DrTS 為模板,模擬片段為 3 到 534 個胺基酸,胺基酸序列之間的一致性為 57%。模擬的 TtTS 結構中央 由(β/α)8的桶狀結構所組成,胺基酸包含 3-102、175-303、351-372
- 30 -
與 440-451(紅色),末端 C 區域由反向平行的 β 摺疊所組成,胺基酸 包含 452-534(淡藍色),DrTS 參與受質進出開關的 S7 區域以及 B 區 域,對應 TtTS 為 304-350(S7 區域,綠色)、103-174(B 區域,深藍色),
還有富含環狀但無明顯二級結構的 S8 區域,胺基酸包含 373-439(黃 色)(圖七)。
五、模擬 TtTS 與麥芽糖的複合物立體結構
利用 AutoDock Vina 軟體將麥芽糖嵌合入 TtTS 蛋白質分子三級 結構,嵌合範圍選定催化性胺基酸 ASP199,ASP308,GLU242 位置,
進行能量最適化運算,挑選結合能最低的結構(圖八)。結果顯示,麥 芽糖所結合的位置與 DrTS 活性中心模擬的麥芽糖結合位置相似,其 中非還原性端葡萄糖在-1 subsite 結構相同,而+1subsite 有些微不同 (圖九),選定 TtTS 與麥芽糖距離 5Å 之內共有 18 個胺基酸(圖十),
依據空間位置分類,在-1subsite 胺基酸分別為 Asp60、Tyr63、His103、
Gln167、Arg197、Asp199、Glu242、His307、Ap308、Arg387,而 +1subsite 胺基酸為 Ile140、Phe141、Phe163、Ala200、Tyr203、Asn244、
Glu309、Glu313,利用 Pymol 軟體分析麥芽糖分子與胺基酸作用力 整理(表四),其中+1subsite 和麥芽糖無明顯作用力的胺基酸有 6 個分 別為 Ile140、Phe141、Phe163、Ala200、Tyr203,Asn244。比對與 TtTS 三級模擬結構相似的 GH13 family 蛋白(表五),分析位於+1 subsite 與麥芽糖無明顯作用力的 6 個胺基酸,發現 Ile140、Phe141、
Phe163、Asn244 於水解酵素中有相同的胺基酸,推測可能因此參與 水解作用,能透過突變進行改造;而 Ala200、Tyr203 屬於 TS unique 的胺基酸,對 TS 異構化作用可能扮演重要的腳色,因此不納入突變
- 31 - Serine,分別為 : Phe141Tyr(TTC→TAC)、Phr141Ser(TTC→TCC)、
Phe163Tyr(TTC→TAC)、Phe163Ser(TTC→TCC)。依據方法十.3.進行 為 9.6%,逆反應水解為 4.6%,蔗糖水解為 2.1% (圖十四),Phe141Ser 突變株正反應水解為 6.2%,逆反應水解為 2.5%,蔗糖水解為 1.4% (圖 十五);而 Phe163Tyr 突變株仍保有異構化能力,但正反應、逆反應 以及蔗糖異構化效率降為野生型 TtTS 的 1/3,分別為 21.3%、13.2%
- 32 -
變,分別為 : Phe141X(TTC→NNK)、Phe163X(TTC→NNK)、
Ile140(ATC→NNK)、Asn244(AAC→NNK),依據方法十.3.進行定位
- 33 - 此胺基酸仍維持為 Proline(CCC→CCT)而非 Leucine(CTT),突變率依 據判讀核苷酸數量中所含突變之比例,平均 5 端與 3 端定序突變率, 合為 I140、F140L、F163V、N244,其餘皆為 WT;失去活性之突變
- 34 -
株則突變成不同胺基酸(表八)。
6. 定序之突變株轉化蔗糖產率分析
確認序列之突變株依方法五.1.反應蔗糖之海藻酮糖最大轉化率 與反應麥芽糖所得產物組成結果合併(表九),可發現能異構化麥芽糖 之突變株,能力與異構化蔗糖相當,而完全無法異構化或水解麥芽 糖之突變株,對於蔗糖同樣沒有反應。如多點飽和突變株 I140、
F140L、F163V、N244,海藻糖最大轉化率為 50.4%,海藻酮糖最大 轉化率 58.2%,與 WT 相當,而突變株 I140A 對麥芽糖與蔗糖皆無 反應,這樣的結果顯示 TtTS 具有相同活性中心,是透過類似的催化 方式異構化麥芽糖與蔗糖;其中,N244V 突變株異構化蔗糖生成產 物偏向異麥芽酮糖,從 WT 海藻酮糖與異麥芽酮糖組成比為 30 : 1 變為 1: 1.5,顯示 N244 位點可決定產物專一性(product specificity)。
- 35 -
- 36 - 得最後製得麥芽糖晶體以 β-構型居多(Xiaoling, 1998)。雖然當麥芽糖
晶體溶於水中會發生自發性變旋作用產生 α-麥芽糖,但最終平衡時 結果來看(圖四、B),只有一小時內 TtTS 反應速率有因為 α-amylase 處理較久而提升,再長時間的海藻糖產量不論任何處理組皆相近,
- 37 -
三、GI-coupled TtTS 反應分析
加入 GI 至 TtTS 反應麥芽糖溶液中,可將水解副產物葡萄糖異構
- 38 - 實驗只使用單一引子進行 QuickChange 突變,即便後續有以 DpnI 處 理,但含有半甲基化(hemimethylated)之親代股 DNA 抵抗 DpnI 作用 能力較高,不易清除,且後續含有缺口(nick)之新合成 DNA 又較完 整親代股 DNA 難轉型至 E.coli,使得最後突變庫中包含高比例之親 代股基因(Vovis and Lacks, 1977)。這樣的現象可透過互補引子對的 設計,排除傳統 QuickChange 操作時無法以新合成股作為模板進行 放大的缺點,以提高基因合成效率、減少親代股基因相對數量(Liu and Naismith, 2008)。
而本實驗使用 NNK 進行多點飽和突變,所需篩選突變庫數量可 減少為 NNN 之一半(Patrick et al., 2003),但為了避免前述所產生之現 象,需要進行更多的篩選(over sampling),勢必會造成實驗負擔,且 NNK 仍有胺基酸冗贅現象(codon redundancy),如 Arginine、Leucine 包含了三種密碼子組成,會因此增加特定胺基酸出現頻率。若使用 NDT/VHG/TGG (D,G = no C、V = no T、H = no G) 搭配,只需要 22
- 39 - 株異構化能力不到 WT 的一半(圖二十五)。根據 Arnold and Georgiou 的統計分析結果,失活突變株約占突變庫之 40%~50%為最適合操作 定向演化的突變率,其基因(平均為含一千個鹼基)約含有 1 至 2 個胺 基酸突變(Arnold and Georgiou, 2003)。越高突變率則可能因為引入終 止密碼子或不利於結構摺疊之胺基酸而提高失活突變株的比例(Guo 證 96 孔突變株之反應結果,可發現利用 Maltase-coupled TS assay 作 為篩選改善異構化能力之突變株方式,分析所得突變株轉化率須高
- 40 -
20 倍體積培養後,分離 TtTS 蛋白再反應結果與高通量反應結果一 致,顯示此高通量系統大致上能反映酵素異構化程度,但顧及 Maltase-coupled TS assay 會有誤差存在,因此後續進行 I140X、N244X 以及多點飽和與隨機突變所得可能突變株,直接取反應結果(尚未進 定值,稱為溫度因素(temperature factor or B factor),若 B factor 越高 代表此胺基酸於整體酵素結構中受溫度影響的搖擺程度越大,不利 於酵素整體結構的穩定(Yuan et al., 2003),因此藉由突變高 B factor 胺基酸可能會使酵素耐熱性與穩定性上升,例如以 Bacillus subtilis 脂肪酶結構中高 B factor 胺基酸建立飽和突變庫,透過重複飽和突變
- 41 -
(iterative saturation mutagenesis),疊加具有改善耐熱性質之突變株,
最終篩選到耐熱性提升 490 倍之突變株(Reetz et al., 2006);或是以 Serratia plymuthica AS9 麥芽糖異構酶結構中高 B factor 胺基酸建立 定位突變,找到熱穩定性提升 7.6 倍且催化效率增加 38.2%的突變株 (Duan et al., 2016),雖然 B factor 存在於解析度較高之蛋白晶體時最 能反映胺基酸受溫度影響的狀態(Parthasarathy and Murthy, 2000),但 透過軟體計算,亦能直接分析預測一級序列胺基酸之 B factor
(Schlessinger et al., 2006)。
TS 結構中控制受質以及產物進出如同閥門的 subdomain 7 與 選數量(Bosley and Ostermeier, 2005),需要仰賴更高效率之篩選系 統。結合流式細胞儀與螢光激活訊號之超高通量篩選技術,可減少 所需樣本與反應試劑體積,且每小時可分析高達 107之突變株,非常 適合進行酵素定向演化之篩選工具(Agresti et al., 2010),參考 Colin
- 42 -
等人於 2015 發表的水解酶篩選方式,搭配本實驗可規畫初步流程為 附錄十三。將多點飽和、隨機或 B factor 突變策略所建立之突變基 因,轉型至 E. coli,蒐集轉型後細胞株搭配破菌液以及酵素反應受質 麥芽糖,通過油包水珠之生成器,因為油水不互融之特性可於體外 模擬細胞與外界環境的分隔狀態(in vitro compartmentalization),生成 的每一個油包水珠中包含單一個突變株細胞,經過一段反應時間,
使細胞破菌後,釋出的酵素與受質反應,並利用雙層包覆方式“水- 油 - 水”(water in oil in water),混合兩種油包水珠(紅色通道的酵素與綠 色通道之麥芽糖水解酶和 GOD/POD 檢測劑),通入另一個水包油珠 生成器,經過穩定油包水珠之脂溶性溶劑(黃色通道)與去離子水(藍 色通道)形成水包油珠,靜置一段時間使麥芽糖水解酶將未轉化為海 藻糖之殘餘麥芽糖分解成葡萄糖後,通入分析器。由於替換 POD 受 質,從本實驗使用的 o-Dianisidine dihydrochloride 變為能激發出螢光 訊號之 Amplex RedR,當突變株含有較弱螢光訊號則代表異構化能力 提升,蒐集此種突變株後,定序其 DNA,可獲得改善性質之突變酵 素。綜合本研究之結果,可以提供從事 TS 之蛋白質工程的參考。
- 43 -
陸、參考文獻
Agresti, J.J., Antipov, E., Abate, A.R., Ahn, K., Rowat, A.C., Baret, J.-C., Marquez, M., Klibanov, A.M., Griffiths, A.D., and Weitz, D.A.
(2010). Ultrahigh-throughput screening in drop-based microfluidics for directed evolution. Proceedings of the National Academy of Sciences 107, 4004-4009.
Arnold, F.H., and Georgiou, G. (2003). Directed enzyme evolution:
screening and selection methods (Springer Science & Business Media).
Bailey, J., Fishman, P.H., and Pentchev, P. (1967). Studies on
mutarotases I. Purification and properties of a mutarotase from higher plants. Journal of Biological Chemistry 242, 4263-4269.
Belton, P.S., and Gil, A.M. (1994). IR and Raman spectroscopic studies of the interaction of trehalose with hen egg white lysozyme.
Biopolymers 34, 957-961.
Biasini, M., Bienert, S., Waterhouse, A., Arnold, K., Studer, G., Schmidt, T., Kiefer, F., Cassarino, T.G., Bertoni, M., Bordoli, L., et al. (2014).
SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res 42, W252-258.
Bornscheuer, U.T., and Pohl, M. (2001). Improved biocatalysts by directed evolution and rational protein design. Current Opinion in Chemical Biology 5, 137-143.
Bosley, A.D., and Ostermeier, M. (2005). Mathematical expressions useful in the construction, description and evaluation of protein libraries. Biomolecular engineering 22, 57-61.
Caner, S., Nguyen, N., Aguda, A., Zhang, R., Pan, Y.T., Withers, S.G., and Brayer, G.D. (2013). The structure of the Mycobacterium smegmatis trehalose synthase reveals an unusual active site configuration and acarbose-binding mode. Glycobiology 23, 1075-1083.
Carpenter, J.F., and Crowe, J.H. (1989). An infrared spectroscopic study of the interactions of carbohydrates with dried proteins.
Biochemistry 28, 3916-3922.
- 44 -
Chen, F., Nakamura, T., and Wada, H. (2004). Development of New Organ Preservation Solutions in Kyoto University. Yonsei Med J 45, 1107-1114.
Chen, Y.-S., Lee, G.-C., and Shaw, J.-F. (2006). Gene Cloning,
Expression, and Biochemical Characterization of a Recombinant Trehalose Synthase from Picrophilus torridus in Escherichia coli.
Journal of Agricultural and Food Chemistry 54, 7098-7104.
Chica, R.A., Doucet, N., and Pelletier, J.N. (2005). Semi-rational
approaches to engineering enzyme activity: combining the benefits of directed evolution and rational design. Curr Opin Biotechnol 16, 378-384.
Chow, S.-Y., Wang, Y.-L., Ye, L.-C., and Liaw, S.-H. (2015). Crystal Structures of Trehalose Synthase from Deinococcus Radiodurans Reveal a Closed Conformation for Intramolecular Isomerization Catalysis and Mutant Induction of an Active-Site Aperture.
Biophysical Journal 108, 376a-377a.
Colin, P.-Y., Kintses, B., Gielen, F., Miton, C.M., Fischer, G., Mohamed, M.F., Hyvönen, M., Morgavi, D.P., Janssen, D.B., and Hollfelder, F.
(2015). Ultrahigh-throughput discovery of promiscuous enzymes by picodroplet functional metagenomics. Nature communications 6.
Crowe, J.H., Crowe, L.M., Oliver, A.E., Tsvetkova, N., Wolkers, W., and
Crowe, J.H., Crowe, L.M., Oliver, A.E., Tsvetkova, N., Wolkers, W., and