趙小水，李成偉

（1.河南省項城市第一高級中學(xué)，河南 項城 466200；2.廈門大學(xué) 信息科學(xué)與技術(shù)學(xué)院，福建 廈門 361005）

0 引言

葡萄糖氧化酶（簡稱GOx）屬于脫氫酶[1]類，具有專一性，高效性等特點(diǎn)，在有氧環(huán)境中，能夠?qū)Ｒ恍源呋?D-葡萄糖，使其生成葡萄糖酸以及過氧化氫。GOx主要分布在在微生物體內(nèi)[2]，作用繁多，在食品、醫(yī)療、生物燃料等方面均有著良好的發(fā)展前景。

GOx主要是由青霉[3]和黑曲霉產(chǎn)生，許多學(xué)者已具體研究了其結(jié)構(gòu)。GOx由兩個亞基構(gòu)成，兩個亞基的結(jié)構(gòu)性質(zhì)完全相同，每個亞基的序列長度約為580，每個亞基中都有一種輔因子，名為黃素腺嘌呤二核苷酸，簡稱為FAD。在青霉中組成葡萄糖氧化酶的每個亞基中都有587個氨基酸殘基，1個FAD,3個甘露糖以及6個N-乙酰葡萄糖糖胺。而在黑曲霉中提取的葡萄糖氧化酶除了個別殘基[4]之外，肽鏈結(jié)構(gòu)與青霉中的GOx大致相同。

GOx是一種氧化還原酶，它能夠與過氧化氫酶形成一個氧化還原系統(tǒng)，在有氧環(huán)境中，消耗氧氣得到葡萄糖酸。酶的催化效率和反應(yīng)速率一般而言受到酸堿性、溫度、氧氣、動力學(xué)參數(shù)等各種因素影響。對于GOx催化過程來說，不同的底物分子與GOx反應(yīng)酶活力是不同的，其中，當(dāng)β-D-葡萄糖作為底物分子時，GOx對其具有高度專一性，相對活力強(qiáng)，而底物分子6-脫氧-D-葡萄糖（簡稱 6-deoxy-D-glucose）的相對活力就弱的多，催化效率較低。

目前，盡管國內(nèi)外學(xué)者針對GOx的研究已取得許多進(jìn)展[5-7]。但尚無法通過實(shí)驗(yàn)的方法得到GOx與底物分子反應(yīng)的中間產(chǎn)物，只能通過分子模擬的方法進(jìn)行分析。本文基于生物分子動力學(xué)模擬軟件，對GOx與底物分子6-脫氧-D-葡萄糖的復(fù)合體系進(jìn)行了模擬分析，建立了完整的復(fù)合物體系模型，計算了體系能量指標(biāo)等相關(guān)數(shù)據(jù)，探討了GOx與底物分子6-脫氧-D-葡萄糖的作用機(jī)制，可為后續(xù)有關(guān)GOx的實(shí)驗(yàn)提供一定的理論依據(jù)。

2 模擬過程

基于Amber軟件進(jìn)行分子動力學(xué)模擬，具體過程包括：①晶體結(jié)構(gòu)的獲??；②晶體結(jié)構(gòu)文件的預(yù)處理；③小分子模板和缺失參數(shù)的生成；④tleap計算坐標(biāo)文件以及拓?fù)湮募纳?；⑤體系能量最小化處理；⑥體系升溫過程；⑦體系平衡過程；⑧模擬體系平衡的判斷；⑨體系分子動力學(xué)模擬：分析能量、氫鍵、殘基質(zhì)心距離分析、原子距離分析等數(shù)據(jù)指標(biāo)。

3 結(jié)果與討論

3.1 相互作用能分析

圖1為GOx與其底物6-deoxy-D-glucose之間的靜電作用勢能和范德華作用勢能隨時間的變化。從圖中可看出靜電作用勢能與范德華作用勢能存在差距。在0-50 ns的時間范圍內(nèi)，靜電作用勢能波動明顯，在50-80ns的時間范圍內(nèi)，靜電作用勢能發(fā)生大幅度減小，而在80ns-200ns的時間范圍內(nèi)波動較小。范德華作用勢能在0-200ns的時間范圍內(nèi)基本穩(wěn)定。該結(jié)果表明靜電作用和范德華作用共同促進(jìn)GOx與底物6-deoxy-D-glucose的結(jié)合與穩(wěn)定。

圖1 GOx與6-deoxy-D-glucose間相互作用的能量隨時間的變化

3.2 RMSD分析

圖2 GOx的 RMSD值隨時間的變化圖

圖3 GOx的回轉(zhuǎn)半徑隨時間的變化

圖4 GOx活性口袋的回轉(zhuǎn)半徑隨時間的變化

圖2為GOx及含有底物分子6-deoxy-D-glucose的GOx的RMSD值?？煽闯鯣Ox的RMSD值在0-80ns的時間范圍內(nèi)波動明顯，在80-130 ns的時間范圍內(nèi)，GOx的RMSD值較為平穩(wěn)，130-200ns之間，GOx的RMSD值先增大至3 ?，隨后保持穩(wěn)定狀態(tài)。而添加底物后GOx的RMSD值在0-200 ns的時間范圍內(nèi)小幅度波動，逐步達(dá)到平衡。

3.3 回轉(zhuǎn)半徑分析

圖3為添加底物6-deoxy-D-glucose前，GOx回轉(zhuǎn)半徑隨時間的變化有較明顯的波動。GOx在0-65 ns時間范圍內(nèi)，回轉(zhuǎn)半徑波動幅度較小，在23.9 ?上下。70-200 ns的時間范圍內(nèi)波動較大，回轉(zhuǎn)半徑先增大后減小。該結(jié)果說明GOx前期體積先發(fā)生膨脹，最終基本維持穩(wěn)定。圖4表示為GOx添加底物前酶活性口袋的回轉(zhuǎn)半徑隨時間的變化情況?？煽闯鯣Ox活性口袋的回轉(zhuǎn)半徑在前期有明顯波動，40 ns之后時間范圍內(nèi)保持穩(wěn)定，維持在11.3 ?左右。這表明添加底物前酶活性口袋部位體積先發(fā)生膨脹隨后緊縮并維持穩(wěn)定。

4 結(jié)論

基于分子動力學(xué)軟件Amber對含有底物6-deoxy-D-glucose的GOx進(jìn)行了模擬分析。結(jié)果顯示，靜電作用和范德華作用共同促進(jìn)了GOx與6-deoxy-D-glucose的結(jié)合與穩(wěn)定；通過分析RMSD和回轉(zhuǎn)半徑，發(fā)現(xiàn)含有底物分子的GOx活性口袋的RMSD值較低，6-deoxy-D-glucose底物在GOx活性口袋中柔性大，該底物的添加造成GOx活性口袋不夠穩(wěn)定。

[1] 康輝,李紅飛,郝曉霞.葡萄糖氧化酶簡介及其應(yīng)用[J].農(nóng)業(yè)工程技術(shù)·農(nóng)產(chǎn)品加工業(yè),2008,(1):16-19.

[2] 李艷,李靜.葡萄糖氧化酶及其應(yīng)用[J].食品工程,2006,(3):9-11.

[3] 劉超,袁建國,王元秀,等.葡萄糖氧化酶的研究進(jìn)展[J].食品與藥品,2010,(7):285-289.

[4] H.J.Bright,M.Appleby.The pH dependence of the individual steps in the glucose oxidase reaction[J]. The Journal of biological chemistry,1969,244(13):3625-34.

[5] 李艷,李靜.葡萄糖氧化酶及其應(yīng)用[J].食品工程,2006,(3):9-11.

[6] 熊云霞,馬現(xiàn)永,胡友軍,等.葡萄糖氧化酶在畜牧業(yè)中的應(yīng)用及檢測方法研究進(jìn)展[J].飼料工業(yè),2016,(4):15-20.

[7] 安文亭,劉樹棟,孫展英,等.葡萄糖氧化酶對仔豬生長性能、健康、飼料養(yǎng)分利用率及相關(guān)理化指標(biāo)的影響[J].黑龍江畜牧獸醫(yī),2014,(19):77-79.