郭麗，徐剛，吳堅平，楊立榮

（浙江大學化學工程與生物工程學系生物工程研究所，浙江 杭州 310027）

脂肪酶對底物具有良好的立體選擇性，并且可以適應包括疏水性有機溶劑在內的非天然反應環(huán) 境[1-4]，因此成為應用最為廣泛的生物催化劑之一。長期以來，提高脂肪酶在有機溶劑中的催化性能是酶學研究的熱點。Kazlauskas 等[5]提出了脂肪酶對手性底物立體識別的分子機理，認為脂肪酶的催化中心附近區(qū)域是一個容納底物的空腔結構，這一空腔可以看作是由大小不同的兩個子空腔組成，即大、小口袋。當手性底物靠近酶的催化中心時，一種構型的手性中心上較大的基團位于“大口袋”中、較小基團位于“小口袋”中，則該構型的底物受到的空間位阻較小，因而容易與催化中心結合而被催化；相反，如果較大基團位于小口袋中，受到的空間位阻較大，因此與催化中心結合較為困難而不易被催化。脂肪酶正是依靠這種催化中心空腔大小口袋與底物手性基團的性狀匹配來識別底物的手性構型的[6-7]。

根據上述機理，很多研究者成功地通過基因修飾[8-12]或者化學修飾[13-16]等方法縮小脂肪酶的“小口袋”，使其不能容納原本勉強能夠容納的手性基團，從而進一步提高脂肪酶的立體選擇性。如 Chen 等[17]通過基因突變在產堿假單胞菌脂肪酶（Pseudomonas alcaligenes lipase，PaL）活性中心附近引入一個 Cys272，并使用 5,5-二硫代雙(2-硝基苯甲酸) 對半胱氨酸進行特異性修飾，將其活性中心的小口袋減小，從而將該酶催化水解 L-薄荷醇丙酸酯的選擇性提高了3 倍。Tuomi 等[18]使用四甲基甲烷硝基化 PcL 的酪氨酸，將其活性中心的小口袋減小，從而將該酶催化水解 2-苯氧基-1-丙醇乙酸酯的選擇性由 E = 17 提高到了 36。

然而，上述修飾改造在提高選擇性的同時，卻嚴重降低了脂肪酶的催化活力（如 Chen 等的修飾脂肪酶殘余酶活僅有原始酶的 1/35）。此外，盡管脂肪酶能夠很好的適應疏水有機溶劑，現(xiàn)有關于脂肪酶的選擇性強化，卻主要是針對水溶液中的反應體系。對有機溶劑中脂肪酶立體選擇性的計算模擬，以及在此指導下的理性改造鮮有報道。本文基于分子動力學（molecular dynamics，MD）技術，開發(fā)了模擬有機溶劑中脂肪酶的立體識別的新方法。并運用該方法，成功將 Kazlauskas 規(guī)則推廣到疏水有機溶劑體系中。以PcL 在正己烷中催化多種手性仲醇的選擇性轉酯化為模型，在模擬計算指導下對PcL 的催化空腔進行特異性化學修飾，在不損失酶活的條件下，從而成功提高了其對多種手性仲醇的選擇性，最好E 值提高到了 48.1。

1 實驗部分

1.1 實驗材料

洋蔥假單胞菌脂肪酶（Pseudomonas cepacia lipase，縮寫為PcL，商品名Lipase PS）購自Amano Enzyme Inc.。胰蛋白酶（trypsin）購自北京 Promega 公司。手性仲醇，如仲丁醇、仲己醇、仲辛醇及其光學純對映體標準品購自 Tokyo Chemical Industry Co. Ltd.。?；w如乙酸乙烯酯和化學修飾劑 N-乙酰咪唑購于 Sigma-Aldrich Inc.。其余試劑均為分析純。

1.2 分子模擬

首先在 Discovery Studio 3.0 軟件中處理PcL蛋白結構（PDB：2NW6）[19]。使用Calculate Protein Ionization and Residue pK 工具設定蛋白表面可電離氨基酸的質子化狀態(tài)，并補全晶體結構中缺失的氫原子，刪除晶體結構中原有的配體小分子（POT612），保留全部結晶水（HOH615～893）和離子（CA613，NA614），再將催化中心三聯(lián)體 Ser87-His286-Asp264中組氨酸的質子化狀態(tài)設為 HIP，從而得到受體蛋白結構。使用 Discovery Studio 的 Sketching 工具生成底物（手性仲醇和?；w）同時與催化中心結合時的四面體過渡態(tài)結構，并使用 LigandFit 工具[20]將其對接到PcL 活性中心，再以共價單鍵連接底物羰基碳原子與PcL的催化中心 Ser87-Oγ。

將得到的酶-底物過渡態(tài)復合物導入到 AMBER11 軟件[21]中，給復合物周圍加上六面體形的TIP3P 水溶劑盒子[22]，溶劑邊界到蛋白表面最遠距離 10?（1?=0.1nm），并以Na+和Cl-平衡蛋白所帶電荷，進行2.5ns 恒溫恒壓MD 以獲得接近真實情況的酶-底物復合物構象。接下來進行5ns 正己烷中的MD 模擬。從水溶液MD 結果中剔除后來加入的TIP3P 溶劑水，將提取出的帶有結晶水、離子的復合物導入Packmol 軟件，在其周圍加上六面體形的正己烷溶劑盒子，溶劑邊界到蛋白表面最遠距離10?，再導入到AMBER11 軟件中進行5ns 有機相MD 模擬。模擬結束后使用AMBER PTRAJ 工具分析結果。

1.3 手性仲醇的動力學拆分

手性仲醇（100μmol）、?；w（100μmol）、正己烷（1mL）與PcL 酶粉（10mg，約含 0.2mg 蛋白）充分混合，200 r/min 攪拌下 25℃溫浴一段時間后離心取上清液。反應式如圖式1 所示。反應液用氣相色譜分析[Hewlett-Packard 1890 型氣相色譜儀，Varian CP-Cyclodextrin-β-2，3，6-M-19 毛細管手性柱（50m × 0.25mm × 0.25μm）]。參照文獻[23]，計算反應的對映體選擇性：E = ln[1 - c(1 + eep)]/ln[1 - c(1 - eep)]，其中eep為產物對映體過剩值，c 為底物轉化率，eep= {[(R)-酯] - [(S)-酯] } / { [(R)-酯] + [(S)-酯] }，方括號表示濃度。

圖1 PcL 催化拆分仲醇

1.4 PcL 的化學修飾

N-乙酰咪唑(N-acetylimidazole，NAI) 修飾酪氨酸殘基(Tyr)的方法參照文獻[24]，10mg 脂肪酶PcL溶于1mL 磷酸緩沖液中（50mmol/L，pH = 7.5），于2mL 平底離心管中配成酶液，加入5mg NAI，并加入200μL 界面激活劑正己烷，形成兩相界面，激活脂肪酶打開“蓋子”，從而利于口袋內部氨基酸修飾[25]。冰浴下200r/min 攪拌反應3h，使用脫鹽柱（HiTrap Desalting，GE Healthcare Corp.）將修飾酶液中的目標蛋白與各種鹽離子和反應后的修飾劑 NAI 進行分離，最終溶解在去離子水中，立即置于 -80℃冰箱預凍 6～8h，然后真空冷凍干燥。冷凍干燥方法參照文獻[26]，凍干完成后將樣品取出，儲存于 -20℃保存待用。PcL 的剩余水解活力對修飾時間作圖得到修飾失活動力學曲線，水解酶活測定采用 p-NPP 法[27]。

1.5 修飾PcL 蛋白質譜分析

化學修飾的PcL 蛋白通過十二烷基磺酸鈉-聚丙烯酰胺凝膠電泳（SDS-PAGE）膠內酶切方法獲取肽段[28]，采用基質輔助激光解析電離飛行時間質譜（MALDI-TOF-MS）檢測方法，在 UltraFleXtreme 質譜儀（Bruker Corp.）上進行肽段質譜分析，從而獲得修飾位點的相關信息。本實驗采取基質 α-氰-4-羥基肉桂酸（HCCA）對肽段進行質譜分析。PcL蛋白經胰蛋白酶（trypsin）酶解的肽譜由 ExPASy PeptideMass（http：//web.expasy.org/peptide_mass/）在線工具檢索產生。搜索時漏切數選擇1 并顯示質荷比 Z/A 大于 500 Da 的肽段序列。

2 結果與討論

2.1 PcL 在正己烷中對手性仲醇的立體識別的模擬計算

分別對PcL 在正己烷中催化(R)-及(S)-2-丁醇、2-己醇、2-辛醇與乙酸乙烯酯的轉酯化反應進行了MD 模擬。以(R)-2-丁醇與乙酸乙烯酯的轉酯化為例，圖2顯示了2-丁醇與乙酸乙烯酯的乙?；赑cL催化中心結合形成的過渡態(tài)構象。圖中青色的部分為底物分子，其醇部較大的手性基團（乙基）朝向催化空腔開口一側，其附近空間開闊，即為PcL 的“大口袋”，較小的手性基團（甲基）朝向催化空腔內側，附近被構成空腔壁的氨基酸殘基所包圍，較為擁擠，即為“小口袋”。Tyr29殘基位于底物的大小手性基團之間，是大小口袋的分界點。Tyr29的酚羥基朝向催化空腔內側（即“小口袋”），在這一位置通過化學修飾連接合適的基團，會使“小口袋”更為擁擠，根據Kaslauskaz 規(guī)則，這將提高PcL 對底物的選擇性。表1 列出了對Tyr29酚羥基進行乙酰化修飾的PcL 以及原始的PcL 催化一系列手性仲醇的酶-底物過渡態(tài)結合能的計算數據。對于一對手性底物，酶與兩者的結合能差的絕對值越大，選擇性越好。表1 數據顯示，Tyr29-酚羥基乙?；腜cL結合自由能差的絕對值明顯大于原始PcL，表明這一修飾能夠顯著提高選擇性。

2.2 對PcL 的 Tyr29 的化學修飾

NAI 是一種修飾酪氨酸酚羥基的特異性化學試劑。在水-疏水界面上用NAI 可以修飾PcL 靠近活性中心的Tyr 殘基。如圖3 所示，NAI 修飾PcL，修飾反應持續(xù)6hPcL 的本征活力(p-NPP 水解活力) 已降至最低，表明NAI 修飾影響到了PcL 的活性中心附近區(qū)域。將原始PcL 酶液和脫鹽后的NAI 修飾PcL 酶液，進行SDS-PAGE 分析，如圖4 所示。2 號泳道（原始PcL）和3 號泳道（修飾PcL）都僅在33 kDa 左右出現(xiàn)的單一條帶，與文獻報道的PcL 分子量一致[19]。這一電泳結果結合圖3 的修飾-失活過程，表明PcL 分子上僅有少數Tyr 被修飾。

用特異性切開 N-端為 Arg/Lys 殘基的胰蛋白酶水解PcL，可以得到特異性的定酶切肽指紋圖譜，譜中的肽段分子量可以通過軟件計算，通過MALDI-TOF-MS 蛋白質譜表征，也可以精確測定這些肽段的分子量。通過對比計算預測的肽譜和測得的肽譜中特定肽段的質量遷移，即可定位發(fā)生修飾反應的氨基酸殘基。表2 給出了預測和測定的肽指紋圖譜數據，可以看到10 號肽段的分子量相比理論值（9 號肽段質量）高出 42 Da，即乙?；姆肿恿浚?Tyr29剛好位于這一肽段中，從而證明修飾特異性發(fā)生在 Tyr29位點。圖4 是原始PcL 的 MALDI-TOF-MS 原始圖譜，可以看到其中9 號肽段的分子量與表2 中的預測值一致，而圖5 中9 號肽段的響應值幾乎消失，10 號肽段則與9 號肽段相差乙?；姆肿恿?。對比圖4 和圖5，也證明了修飾發(fā)生在位于10 號肽段上的 Tyr29殘基。

圖2 2-丁醇與乙酸乙烯酯的乙?；赑cL 催化中心結合形成的過渡態(tài)構象

表1 分子動力學模擬計算的PcL 與不同手性醇過渡態(tài)復合物的結合能

圖3 NAI 修飾PcL 失活動力學

圖4 脂肪酶PcL 蛋白的 SDS-PAGE 結果

表2 NAI-修飾的PcL 的MALDI-TOF 質譜表征

圖5 原始脂肪酶PcL 的胰蛋白酶酶切蛋白樣品指紋圖譜

2.3 修飾 Tyr29 對選擇性的影響

圖6 NAI 修飾脂肪酶的胰蛋白酶酶切蛋白樣品指紋圖譜

前文表1 已經對PcL 在正己烷中識別手性底物進行了模擬計算，考察了PcL 催化3 種不同結構仲 醇的不對稱轉酯化，Tyr29-酚羥基乙?；腜cL 結合自由能差的絕對值明顯大于原始PcL，表明這一修飾能夠顯著提高拆分仲醇的選擇性，見圖6。進而對PcL 進行了 NAI 特異性修飾實驗研究，對表3 中考察的反應而言，Tyr29-酚羥基乙?；疨cL 的選擇性較原始PcL 均有明顯提高，其中修飾PcL 催化 2-己醇與乙酸乙烯酯的轉酯化反應的選擇性提高最多，由原來的 E = 12.6 提高到 48.1，這與表1 的模擬計算是一致的。此外，在選擇性提高的同時，修飾PcL 催化這些仲醇底物的活力還有小幅提高。

表3 對PcL 氨基酸殘基側鏈的特異性化學修飾

3 結 論

根據 Kazlauskas 規(guī)則，通過有機相中MD 模擬的新方法，提出了對PcL 催化空腔中的關鍵殘基進行特異性修飾來提高酶的立體選擇性的研究策略。并通過MALDI-TOF-MS 表征和PcL 催化3 種不對稱轉酯化反應的實驗數據證明這一修飾的有效性。這一方法成功地在不損失酶的催化活力的情況下提高了其立體選擇性，并且其研究方法借助獨立于分子生物學之外的化學修飾方法，結合蛋白質譜解析和分子動力學模擬等現(xiàn)代技術手段，開創(chuàng)了一條嶄新的靈活、高效、可靠的酶構效關系研究。

[1] 崔麗娟，徐剛，孟梟，等. 多元電解質對脂肪酶有機相拆分炔戊醇的激活[J]. 化工進展，2014，33（8）：2150-2154.

[2] 戴曉庭，孟梟，徐剛，等. 酰基供體對動態(tài)動力學拆分1-四氫萘胺的影響[J]. 化工進展，2014，33（9）：2421-2424.

[3] Bajaj A，Lohan P，Jha P N，et al. Biodiesel production through lipase catalyzed transesterification：An overview[J]. J. Mol. Catal. B：Enzym.，2010，62（1）：9-14.

[4] Fan Y X，Qian J Q. Lipase catalysis in ionic liquids/supercritical carbon dioxide and its applications[J]. J. Mol. Catal. B：Enzym.，2010，66（1）：1-7.

[5] Laane C，Boeren S，Vos K，et al. Rules for optimization of biocatalysis in organic solvents[J]. Biotechnol. Bioeng.，1987，30（1）：81-87.

[6] Yadav G D，Lathi P S. Kinetics and mechanism of synthesis of butyl isobutyrate over immobilized lipases[J]. Biochem. Eng. J.，2003，16（3）：245-252.

[7] Lemke K，Lemke M，Theil F. A three-dimension predictive active site model for lipase from Pseudomonas cepacia[J]. J. Org. Chem.，1997，62（18）：6268-6273.

[8] Fernandez L R. Lipase from Thermomyces lanuginosus：Uses and prospects as an industrial biocatalyst[J]. J. Mol. Catal. B：Enzym.，2010，62（3）：197-212.

[9] Contesini F J，Lopes D B，Macedo G A，et al. Aspergillus sp. lipase：Potential biocatalyst for industrial use[J]. J. Mol. Catal. B：Enzym.，2010，67（3）：163-171.

[10] Li N，Zong M H. Lipases from the genus Penicillium：Production，purification，characterization and applications[J]. J. Mol. Catal. B：Enzym.，2010，66（1）：43-54.

[11] Rodrigues R C，F(xiàn)ernandez L R. Lipase from Rhizomucor miehei as a biocatalyst in fats and oils modification[J]. J. Mol. Catal. B：Enzym.，2010，66（1）：15-32.

[12] Rodrigues R C，F(xiàn)ernandez L R. Lipase from Rhizomucor miehei as an industrial biocatalyst in chemical process[J]. J. Mol. Catal. B：Enzym.，2010，64（1）：1-22.

[13] Cabrera Z，F(xiàn)ernandez L G，F(xiàn)ernandez L R，et al. Enhancement of Novozym-435 catalytic properties by physical or chemical modification[J]. Proc. Biochem.，2009，44：226-231.

[14] Palomo J M，F(xiàn)ernandez L G，Guisan J M，et al. Modulation of immobilized lipase enantioselectivity via chemical amination[J]. Adv. Synth. Catal.，2007，349（7）：1119-1127.

[15] Rodrigues R C，Berenguer M A，F(xiàn)ernandez L R. Coupling chemical modification and immobilization to improve the catalytic performance of enzymes[J]. Adv. Synth. Catal.，2011，353（13）：2216-2238.

[16] Mateo C，Palomo J M，F(xiàn)ernandez L G，et al. Improvement of enzyme activity，stability and selectivity via immobilization techniques[J]. Enzym Microb. Technol.，2007，40（6）：1451-1463.

[17] Chen H，Wu J P，Yang L R，et al. A combination of site-directed mutagenesis and chemical modification to improve diastereopreference of Pseudomonas alcaligenes lipase[J]. Biochim. Biophys Acta，2013，1834（12）：2494-2501.

[18] Tuomi W V，Kazlauskas R J. Molecular basis for enantioselectivity of lipase from Pseudomonas cepacia toward primary alcohols. Modeling，kinetics，and chemical modification of Tyr29 to increase or decrease enantioselectivity[J]. J. Org. Chem.，1999，64（8）：2638-2647.

[19] Kim K K，Song H K，Shin D H，et al. The crystal structure of a triacylglycerol lipase from Pseudomonas cepacia reveals a highly open conformation in the absence of a bound inhibitor[J]. Sturcture，1997，5（2）：173-185.

[20] Peretz A，Pell L，Gofman Y，et al. Targeting the voltage sensor of Kv7.2 voltage-gated K+channels with a new gating-modifier[J]. PANS，2010，107（35）： 15637-15642.

[21] Case D A，Darden T A，Cheatham T E ，et al. AMBER 11[C]. San Francisco：University of California，2011，

[22] Jorgensen W L，Chandrasekhar J，Madura J D，et al. Comparison of simple potential functions for simulating liquid water[J]. J. Chem. Phys.，1983，79（2）：926-935.

[23] Chen C S，F(xiàn)ujimoto Y，Girdaukas G，et al. Quantitative analyses of biochemical kinetic resolutions of enantiomers[J]. J. Am. Chem. Soc.，1982，104（25）：7294-7299.

[24] Furth A J，Hope D B. Studies on the chemical modification of the tyrosine residue in bovine neurophysin-II[J]. Biochem. J.，1970，116：545-553.

[25] Norin M，Olsen O，Svendsen A，et al. Theoretical studies of Rhizomucor miehei lipase activation[J]. Protein Eng.，1993，6（8）：855-863.

[26] Ru M T，Dordick J S，Reimer J A，et al. Optimizing the salt-induced activation of enzymes in organic solvents：Effects of lyophilization time and water content[J]. Biotechnol. Bioeng.，1999，63（2）：233-241.

[27] Gupta N，Rathi P，Gupta R. Simplified para-nitrophenyl palmitate assay for lipases and esterases[J]. J. Med. Chem.，2003，46（4）：499-511.

[28] Zehl M，Lescic L，Abramic M，et al. Characterization of covalently inhibited extracellular lipase from Streptomyces rimosus by matrix-assisted laser desorption/ionization time-of-flight and matrix-assisted laser desorption/ionization quadrupole ion trap reflectron time-of-flight mass spectrometry：Localization of the active site serine[J]. J. Mass Spectrom.，2004，39（12）：1474-1483.