陳茹箐 王祖琪 馮繼濤 施琪浩 張英杰 余彬彬 金燕仙 秦利明 葉余原 蔣華江 章少陽 鐘愛國

摘? ? ? 要： 用密度泛函理論（DFT）和信息理論（IT）的方法，優(yōu)化了11種甘氨酸結(jié)構(gòu)并預(yù)測甘氨酸衍生物的 pKa 值。結(jié)果發(fā)現(xiàn)，甘氨酸及其10種衍生物的 H 原子的二級費歇爾熵與其酸度系數(shù)pKa值線性方程相關(guān)系數(shù)r值（Y =-34.960 + 9.751X ，r =0.868）最大。這一發(fā)現(xiàn)對研究甘氨酸衍生物結(jié)構(gòu)和相關(guān)性質(zhì)有著極其重要的意義。此外，還采用密度泛函理論模擬預(yù)測了甘氨酸在不同位置的一元取代紅外吸收、拉曼光譜、核磁共振和紫外吸收光譜，所得結(jié)果對制藥工程專業(yè)本科生開展拓展訓(xùn)練提供了好的材料。

關(guān)? 鍵? 詞：甘氨酸;酸度電離常數(shù);信息理論;分子光譜

中圖分類號：O621.25? ? ? ?文獻(xiàn)標(biāo)識碼： A? ? ? ?文章編號： 1671-0460（2020）09-1900-05

Abstract： Eleven glycine structures were optimized and pKa values of glycine derivatives were predicted by density functional theory and information theory. The results showed that the second-order Fisher entropy of H atom of glycine and its 10 derivatives had the maximum linear correlation coefficient（y=-34.960+9.751x， r=0.868） with the pKa. This discovery is very important for the study of the structure and related properties of glycine derivatives. In addition， the density functional theory was used to simulate and predict the infrared， Raman， nuclear magnetic and UV-Vis spectra of glycine at different positions. The results can provide good reference for the development training of pharmaceutical engineering undergraduates.

Key words： Glycine; Acidity coefficient; Information theory; Molecular spectrum

甘氨酸（Glycine，C2H5NO2），又名氨基乙酸，是結(jié)構(gòu)最簡單的α-氨基酸[1]，固態(tài)甘氨酸是白色結(jié)晶粉末或白色單斜晶系。在水中易溶，有特殊甜味，無毒無臭[2]，幾乎不溶于乙醇或乙醚。甘氨酸也是有機(jī)合成中一種重要的原料，廣泛應(yīng)用于農(nóng)藥、醫(yī)藥、化工、食品等領(lǐng)域。它不僅參與蛋白質(zhì)、卟啉、膽鹽等的代謝，具有保護(hù)細(xì)胞、調(diào)節(jié)免疫系統(tǒng)、光譜抗炎等作用，還在中樞神經(jīng)系統(tǒng)中作為神經(jīng)遞? 質(zhì)[3]。經(jīng)研究表明，甘氨酸具有復(fù)雜的生理功能并在許多疾病中發(fā)揮有效的治療作用。而甘氨酸可在人和哺乳動物的體內(nèi)降解，其主要是通過甘氨酸裂解系統(tǒng)的脫羧和脫氫[4]。甘氨酸的主要化學(xué)合成工藝有氯乙酸氨解法、海因法、施特雷克法三種方? ?法[5-6]。目前，國外已采用改進(jìn)的海因法和施特雷克法，但國內(nèi)仍舊采用已被國外淘汰的氯乙酸氨解法來生產(chǎn)甘氨酸[7]。甘氨酸在醫(yī)藥行業(yè)中作為抗胃酸沖劑、氨基酸制劑，也可用作治療重癥肌無力和進(jìn)行性肌肉萎縮，可治療兒童高脯氨酸血癥，可治療胃酯過多癥和慢性腸炎等。農(nóng)業(yè)上氨基酸主要作為家禽家畜等食用的飼料增加氨基酸的引誘劑與添加劑，可用作為水解蛋白添加劑和增效劑，合成除草劑固體草甘膦和殺菌劑異菌脲。本文采用信息理論方法，優(yōu)化甘氨酸及其衍生物的結(jié)構(gòu)，判斷化學(xué)結(jié)構(gòu)鍵對其酸度常數(shù)pKa是否影響[8-9]，通過計算得到相關(guān)性最高的關(guān)于預(yù)測酸度pKa的線性方程，從而起到預(yù)測甘氨酸其他衍生物的酸度pKa的作用。此外還用DFT理論模擬預(yù)測了甘氨酸一元取代紅外光譜、拉曼光譜、核磁和紫外吸收分子光譜[10-11]。

1? 實驗部分

1.1? 實驗與模型

裝有Chem3D、ChemDraw8.0、 GaussView5.0.9、 Gaussian09W、ACD-Labs 6.0、Origin 8.6、 Multiwfn 3.5 等軟件的 PC 電腦一臺。甘氨酸的結(jié)構(gòu)如圖1所示。

通過改變C和N上的取代基獲得其衍生物，來研究其4種熵（申農(nóng)熵、費歇爾熵、二級費歇爾熵和Parr熵）的變化與采用ACDLab6.0軟件測得的pKa的線性關(guān)系。本文用-Cl、-F、-Br、-CH3、-NO2作為取代基在不同位置上進(jìn)行一元取代，總共得到11種甘氨酸及其衍生物。使用Chem3D軟件，構(gòu)建11個甘氨酸及其衍生物的模型，并使用ChemDraw 8.0進(jìn)行其能量最低化;在GaussView5.0.9軟件中，對甘氨酸進(jìn)行構(gòu)型優(yōu)化，并保存mol格式和gif格式，在菜單欄中選Calculate-Gaussian Calculate Setup，分別計算11個甘氨酸及其衍生物的量化參數(shù)。記錄甘氨酸的紅外光譜、核磁光譜、紫外光譜和拉曼光譜的頻率值并用Origin 8.6 作出光譜圖;運用 ACD Labs 6.0 軟件，得到11個甘氨酸及其衍生物相對應(yīng)的酸度常數(shù)pKa預(yù)測值;再用Multiwfn 3.5分別計算11個甘氨酸及其衍生物的信息分子熵，記錄其所需O、H原子的熵值。最后分別將11個甘氨酸及其衍生物O、H原子的4種的熵值與得到的pKa用 Origin 8.6作圖，得到線性回歸方程，比較相關(guān)系數(shù)r值大小，確定r值最大的為預(yù)測方程。

1.2? 結(jié)果與討論

由軟件GaussView 5.0.9計算所得的甘氨酸及其10種衍生物O、H原子的申農(nóng)熵、費歇爾熵、二級費歇爾熵和Parr熵以及由軟件ACD-Labs6.0預(yù)測得到的相應(yīng)衍生物酸度pKa如表 1 和表2所示。

為了更明顯地表達(dá)出熵值變化與酸度pKa的線性關(guān)系，分別用甘氨酸及其10種衍生物的申農(nóng)熵、費歇爾熵、二級費歇爾熵和Parr熵的熵值與pKa值進(jìn)行Origin作圖分析，得出以下結(jié)果。甘氨酸及其10種衍生物O原子的申農(nóng)熵與pKa值關(guān)系相關(guān)系數(shù)r=0.565，線性方程為Y=-31.038+11.431X。甘氨酸及其10種衍生物H原子的申農(nóng)熵與pKa值關(guān)系，相關(guān)系數(shù)r為0.752，r2為0.518，線性方程為Y=-62.632+37.017X 。甘氨酸及其10種衍生物O原子的費歇爾熵與pKa 值關(guān)系相關(guān)系數(shù)r為0.614，r2為0.3079，線性方程為Y=-3550.172+7.964X。甘氨酸及其10種衍生物 H 原子的費歇爾熵與pKa值關(guān)系相關(guān)系數(shù)r為0.785，r2=0.574， 線性方程為Y=-62.959+18.648X。甘氨酸及其10種衍生物O原子的二級費歇爾熵與pKa值關(guān)系相關(guān)系數(shù)r為0.553， r2=0.229， 線性方程為Y=-997.994+2.250X。甘氨酸及其10種衍生物H原子的二級費歇爾熵與pKa值關(guān)系相關(guān)系數(shù)r為0.868，r2為0.727，線性方程為Y=-34.960+ 9.751X。甘氨酸及其10種衍生物H原子的Parr熵與pKa 值關(guān)系相關(guān)系數(shù)r為0.558，r2為 0.235，線性方程為Y=-288.967+5.109 28X。甘氨酸及其10種衍生物H原子的Parr熵與pKa值關(guān)系相關(guān)系數(shù)r為0.717， r2為0.460 6，線性方程為Y=-77.610+ 17.707X 。由以上Origin圖可知甘氨酸及其10種衍生物的H原子的二級費歇爾熵與酸度系數(shù)pKa值線性方程的r 值最大（r=0.868），即相關(guān)性最好。由此可知，可以運用H原子的二級費歇爾熵對酸度pKa值此條線性方程來預(yù)測其他甘氨酸衍生物的酸度系數(shù)，這一方法對研究甘氨酸衍生物結(jié)構(gòu)和相關(guān)性質(zhì)以及開發(fā)甘氨酸衍生物的化學(xué)性質(zhì)有著重要意義。

2? 分子光譜

2.1? 甘氨酸及其衍生物紅外光譜

紅外光譜法上是一種根據(jù)分子內(nèi)部原子間的分子轉(zhuǎn)動和相對振動等信息來確定分子結(jié)構(gòu)和鑒別物質(zhì)的分析方法。通過紅外光譜圖的分析可以得到分子結(jié)構(gòu)和化學(xué)鍵的特點，可以用通過與標(biāo)準(zhǔn)化合物的紅外光譜圖對比來鑒別未知物的結(jié)構(gòu)組成或確定其化學(xué)基團(tuán);而吸收譜帶的吸收強(qiáng)度與化學(xué)基團(tuán)的含量有關(guān)，可用于進(jìn)行定量分析和純度鑒定。本文使用 GaussView 5.0.9 軟件進(jìn)行建模、計算，采用 Frequency方法分別計算得到11種甘氨酸及其衍生物的譜圖頻率，使用 Origin 8.6 進(jìn)行譜圖處理。使用軟件GaussView 5.0.9 水平優(yōu)化結(jié)構(gòu)，采用 Frequency方法計算得到甘氨酸紅外光譜圖并使用Origin作圖分析如圖2（a）所示。甘氨酸在紅外光譜圖的吸收峰如下：在600 cm-1處的吸收峰為—COOH 的伸縮振動區(qū);在928 cm-1為O—H 的面外搖擺振動;在1 152 cm-1為C—C的伸縮振動;在1 568 cm-1為C=O 的伸縮振動。一元取代 R1=—F的甘氨酸衍生物在紅外光譜圖的吸收峰如下：在912 cm-1處的吸收峰為C—H的面外搖擺振動;在2 040 cm-1為—F的不對稱伸縮振動;在1 360 cm-1為 C—C 的伸縮振動。一元取代 R2=-F的甘氨酸衍生物在紅外光譜圖的吸收峰如下：在720 cm-1處的吸收峰為C—H 的面內(nèi)搖擺振動;在950 cm-1為C—H 的面外搖擺振動;在1 260 cm-1為C—O的面內(nèi)伸縮振動;在? ? ? 1 360 cm-1為C—C的伸縮振動;在2 040 cm-1為—F的不對稱伸縮振動;在3 740 cm-1為C—H的伸縮振動區(qū)。一元取代 R1=-Cl的甘氨酸衍生物在紅外光譜圖的吸收峰如下：570 cm-1處的吸收峰為O—H 的面內(nèi)搖擺振動;在950 cm-1為C—H 的面外搖擺振動;在1 230 cm-1為C—H 的面內(nèi)伸縮振動;在? ? ?1 360 cm-1為C—C的伸縮振動;在2 040 cm-1為-Cl的不對稱伸縮振動;在3 740 cm-1為C—H 的伸縮振動區(qū)。R1=-Cl的甘氨酸衍生物紅外光譜圖并使用 Origin作圖分析。一元取代 R1=-Cl的甘氨酸衍生物在紅外光譜圖的吸收峰如下：570 cm-1處的吸收峰為O—H的面內(nèi)搖擺振動;在 950 cm-1 C—H 的面外搖擺振動;在 1 230 cm-1 C—H的面內(nèi)伸縮振動;在 1 60 cm-1為C—C的伸縮振動;在2 040 cm-1為-Cl的不對稱伸縮振動;在3 740 cm-1峰為C—H 的伸縮振動區(qū)（未顯示）。

2.2? 甘氨酸及其衍生物拉曼光譜

拉曼光譜是印度物理學(xué)家Raman于1928年發(fā)現(xiàn)的，其檢測的原理是利用各種分子的振動能級不同，分子在接收到一束頻率為v的光后，分子與光相互作用，散射出頻率為v±v0 的光，而v0 則就稱為拉曼位移;對于同一種分子來說，入射光v 一定時，拉曼位移v0就一定，因此，可使用拉曼光譜來探究分子的結(jié)構(gòu)特征和主要官能團(tuán)。本文使用 GaussView5.0.9 軟件進(jìn)行建模、計算，采用 Frequency 方法分別計算得到11種甘氨酸及其衍生物的頻率，使用Origin 8.6 進(jìn)行譜圖處理。一元取代R2=-F的甘氨酸衍生物拉曼光譜圖并使用Origin作圖分析如圖2（b）所示。

一元取代 R2=-F的甘氨酸衍生物在拉曼光譜圖的吸收峰如下：在4 470 m-1處的吸收峰為OH 的對稱伸縮振動;在3 740 m-1為C—H對稱伸縮振動;在 2 700 m-1為-CN 伸縮振動;在1 800 cm-1為C—H的對稱伸縮振動;在1 740 cm-1為C-H的反對稱伸縮振動;在1 370 cm-1為C—C的伸縮振動;在10 300 cm-1為—OH 面外搖擺振動。拉曼光譜中所顯示的波峰與實驗數(shù)值有所浮動，這是由于 GaussView5.0.9 在建立模型時存在相應(yīng)的誤差，并不影響基本官能團(tuán)的判斷。一元取代 R1=-Cl的甘氨酸衍生物在拉曼光譜圖的吸收峰如下：在4 500 cm-1處的吸收峰為—OH的對稱伸縮振動;在3 740 cm-1為 C—H對稱伸縮振動;在1 740 cm-1為C—H 的反對稱伸縮振動;在1 360 cm-1處的吸收峰為C—C的伸縮振動;在1 230 cm-1處的吸收峰C—Cl的伸縮振動;在810 cm-1為碳骨架面內(nèi)搖擺振動。一元取代 R2=-Cl 的甘氨酸衍生物拉曼光譜圖并使用Origin 作圖分析如圖（未顯示）。一元取代 R2=-Cl的甘氨酸衍生物在拉曼光譜圖的吸收峰如下：4 500 cm-1處的吸收峰為—OH的對稱伸縮振動;在3 740 m-1 為C—H 對稱伸縮振動;在1 740 cm-1為C—H的反對稱伸縮振動;在1 360 cm-1為C—C的伸縮振動;在1 230 cm-1為C—Cl的伸縮振動;在810 cm-1為碳骨架面內(nèi)搖擺振動。拉曼光譜中所顯示的波峰與實驗數(shù)值有所浮動，這是由于GaussView 5.0.9在建立模型時存在相應(yīng)的誤差，并不影響基本官能團(tuán)的判斷。

2.3? 甘氨酸及其衍生物核磁共振吸收光譜

核磁共振是磁矩不為零的原子核，在外磁場作用下自旋能級發(fā)生塞曼分裂，共振吸收某一定頻率的射頻輻射的物理過程。它是對各種有機(jī)和無機(jī)物的成分、結(jié)構(gòu)進(jìn)行定性分析的最強(qiáng)有力的工具之一，有時亦可進(jìn)行定量分析。本文使用 GaussView 5.0.9軟件進(jìn)行建模、計算，采用 NMR 方法分別計算得到11種甘氨酸及其衍生物的頻率，使用 Origin 8.6進(jìn)行譜圖處理。使用軟件ChemDraw 8.0打開物質(zhì)模型，采用Predict1H-NMR shifts和Predict 13H-NMR shifts方法計算得到甘氨酸核磁光譜。表3顯示了甘氨酸C、H位移量。使用軟件 ChemDraw 8.0打開物質(zhì)模型，采用Predict 1 H-NMR shifts 和Predict 13H-NMR shifts 方法計算得到甘氨酸衍生物核磁光譜。表 3顯示了一元取代 R1=-F 的甘氨酸衍生物核磁 C、H 位移量。

2.4? 甘氨酸及其衍生物紫外吸收光譜

由于物質(zhì)分子內(nèi)部在進(jìn)行轉(zhuǎn)動、振動和電子運動，相應(yīng)狀態(tài)的能量即狀態(tài)的本征值是量子化的，因此分子具有轉(zhuǎn)動能級、振動能級和電子能。利用物質(zhì)分子對紫外光的吸收所產(chǎn)生的紫外可見光譜及吸收程度就可以對物質(zhì)的結(jié)構(gòu)、含量和組成進(jìn)行分，這就是紫外光譜法。本文使用GaussView 5.0.9 軟件進(jìn)行建模、計算，采用Energy方法分別計算得到11種甘氨酸及其衍生物的頻率，使用Origin 8.6進(jìn)行譜圖處理。使用軟件GaussView 5.0.9 水平優(yōu)化結(jié)構(gòu)，采用 Energy方法計算得到甘氨酸紫外光譜圖并使用 Origin8.6軟件作圖分析如圖 3（a）所示。

圖中顯示了甘氨酸在171.8 nm處有最大吸收波長，位于近紫外區(qū)。這是因為 n→σ* 的躍遷，吸收波長較長。一元取代 R1=-F 的甘氨酸衍生物紫外光譜圖并使用 Origin8.6 軟件作圖分析。顯示了一元取代 R1=-F 的甘氨酸衍生物在174.8 nm 處有最大吸收波長，位于近紫外區(qū)。這是因為π→π* 的躍遷，吸收波長較長。R2=-F的甘氨酸衍生物紫外光譜圖并使用Origin8.6軟件作圖顯示了一元取代R2=-F的甘氨酸衍生物在 299.2 nm處有最大吸收波長，位于近紫外區(qū)。這是因為n→σ*躍遷，吸收波長較長。使用軟件 GaussView 5.0.9 水平優(yōu)化結(jié)構(gòu)，一元取代 R1=-Cl的甘氨酸衍生物在192 nm 處有最大吸收波長，位于近紫外區(qū)。這是因為π→π* 的躍遷，吸收波長較短。R2=-Cl的甘氨酸衍生物紫外光譜圖并使用Origin8.6軟件作圖分析如圖3（b）所示。圖中顯示了一元取代R2=-Cl的甘氨酸衍生物在? 260.3 nm 處有最大吸收波長，位于近紫外區(qū)。這是因為n→σ*躍遷，吸收波長較長。使用軟件GaussView 5.0.9水平優(yōu)化結(jié)構(gòu)，采用Energy 方法計算得到一元取代 R1=-Br的甘氨酸衍生物紫外光譜圖并使用Origin8.6軟件作圖分析如圖 3（a）所示。圖中顯示了一元取代 R1=-Br的甘氨酸衍生物在 212.8 nm處有最大吸收波長，位于近紫外區(qū)。這是因為π→π*的躍遷，吸收波長較長。一元取代 R2=-Br甘氨酸衍生物紫外光譜圖并使用Origin8.6軟件作圖分析如圖3（b）所示。圖中顯示了一元取代R2=-Br 的甘氨酸衍生物在 271.5 nm 處有最大吸收波長，位于近紫外區(qū)。這是因為 n→σ* 的躍遷，吸收波長較長。

3? 結(jié) 論

本文將采用信息理論的方法，構(gòu)建并優(yōu)化甘氨酸及其10種衍生物模型，得出以下結(jié)論：

1）甘氨酸及其10種衍生物的H、O原子的費歇爾熵 、二級費歇爾熵存在較小的變化，而H、O原子的申農(nóng)熵和 Parr熵變化較為明顯，其中H原子的二級費歇爾熵與pKa 線性關(guān)系最好。

2）本文將采用信息理論的方法，構(gòu)建并優(yōu)化甘氨酸及其10種衍生物模型，通過計算得到相關(guān)性最高的關(guān)于預(yù)測酸度常數(shù)pKa的線性方程，其中H原子的二級費歇爾熵與pKa的線性方程r值最大（r=0.869，r2=0.727），即相關(guān)性最好。因此我們可以選用此條線性方程來預(yù)測其他甘氨酸衍生物的酸度常數(shù)，這對我院制藥工程專業(yè)學(xué)生今后學(xué)習(xí)氨基酸分子衍生物檢測提供參考數(shù)據(jù)。

3）本文對甘氨酸及部分衍生物進(jìn)行紅外吸收、拉曼光譜、核磁共振和紫外吸收光譜的進(jìn)行了模擬預(yù)測，并對譜圖中的特征峰進(jìn)行指認(rèn)和歸屬，為我院工科專業(yè)實施拓展訓(xùn)練材料及藥物光譜預(yù)測提供材料。

參考文獻(xiàn)：

[1]楊華. 甘氨酸生產(chǎn)應(yīng)用概況[J]. 四川化工， 2004（3）：31-34.

[2]黃平，李云政，張青山. 甘氨酸的應(yīng)用與合成工藝進(jìn)展[J]. 化學(xué)研究， 2004（2）：77-80.

[3]宋彥梅， 尹秋響， 王靜康. 甘氨酸的應(yīng)用及生產(chǎn)技術(shù)[J]. 氨基酸和生物資源， 2003（2）：55-60.

[4]徐澤輝，刁春霞， 黃亞茹， 等. 甘氨酸的生產(chǎn)現(xiàn)狀及發(fā)展趨勢[J]. 石油化工技術(shù)經(jīng)濟(jì)， 2004，20（5）：41-45.

[5]盧金龍. 5 000 t/a甘氨酸生產(chǎn)技術(shù)總結(jié)[J]. 氯堿工業(yè)， 2001（4）： 24-25.

[6]毛建衛(wèi)，崔艷麗. 甘氨酸合成與分離新工藝[J]. 北京化工大學(xué)學(xué)報（自然科學(xué)版），1997（4）：66-68.

[7]杜瑞平，張興夫，高民，等. 甘氨酸的免疫調(diào)節(jié)作用及其分子機(jī)制[J]. 動物營養(yǎng)學(xué)報，2015，27（3）：663-670.

[8]徐伯華，李來才，唐作華. 甘氨酸與水分子間相互作用的理論研究[J]. 原子與分子物理學(xué)報， 2003（2）：275-279.

[9]劉建慧，張愛萍. 國內(nèi)甘氨酸的生產(chǎn)與應(yīng)用[J]. 化學(xué)工業(yè)與工程技術(shù)， 1999（1）：28-30.

[10]徐洋，羅海霞，范強(qiáng)，等.新藥萘夫西林藥理性質(zhì)的密度泛函研究[J].當(dāng)代化工，2015，44（4）：699-701.

[11]陳凱浩，鐘愛國.甲基苯丙胺光譜性質(zhì)的密度泛函分析與指認(rèn)[J].當(dāng)代化工，2014，43（1）：29-31.