磺胺甲惡唑的振動光譜和密度泛函理論研究

呂少巖, 張 靜, 于 倩, 劉 春

(河北科技大學(xué) 環(huán)境科學(xué)與工程學(xué)院 河北省污染防治生物技術(shù)重點(diǎn)實(shí)驗(yàn)室, 石家莊 050018)

1 引 言

磺胺甲惡唑(Sulfamethoxazole，SMX)，又稱為新諾明，是一種常見的磺胺類抗生素. 作為一種廣譜類抗生素，SMX已被主要用于敏感細(xì)菌引起的尿路感染、呼吸系統(tǒng)感染、腸道感染和膽道感染等[1]. 過度使用SMX將會引起其在環(huán)境中的廣泛分布. 由于它的長期暴露會引起耐藥基因的出現(xiàn)，對人的身體健康造成威脅，從而可能會引起嚴(yán)重的公共衛(wèi)生問題[2-5]. 此外，SMX已被世界衛(wèi)生組織國際癌癥研究機(jī)構(gòu)列為3類致癌物，因此對SMX的檢測尤為重要.

目前抗生素的主要檢測方法是紫外分光光度法和高效液相色譜-質(zhì)譜法，但樣品需要預(yù)處理、檢測步驟繁瑣、成本較高. 紅外光譜和拉曼光譜法，具有樣品需求量少、無需特殊處理和檢測快速等優(yōu)點(diǎn)，被廣泛用于環(huán)境監(jiān)測和污染物結(jié)構(gòu)研究中[6，7]. 其中，紅外光譜屬于吸收光譜，是利用物質(zhì)對不同波長的紅外輻射的吸收特性來獲取分子結(jié)構(gòu)和化學(xué)組成；拉曼光譜屬于散射光譜，是利用測定散射光相對入射光頻率的變化來獲取分子內(nèi)部結(jié)構(gòu)信息，通過振動引起分子極化率變化而產(chǎn)生. 盡管它們的產(chǎn)生原理不同，但都是研究分子振動的重要手段，且在一定程度上可互相補(bǔ)充，可實(shí)現(xiàn)更好、更快捷的抗生素檢測. 密度泛函理論(Density functional theory，DFT)是一種量子化學(xué)計(jì)算方法，它的計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果有很好的吻合性，近年來在計(jì)算分子結(jié)構(gòu)性質(zhì)和振動光譜方面得到普遍應(yīng)用[8-10]. 當(dāng)前對SMX的研究主要集中在SMX殘留物的降解，而對分子結(jié)構(gòu)和振動光譜方面研究較少.

本文采用DFT計(jì)算方法，對SMX分子進(jìn)行幾何優(yōu)化及振動光譜的計(jì)算，將理論計(jì)算得到的光譜和實(shí)驗(yàn)光譜進(jìn)行對比，并采用GaussView 6.0對振動模式進(jìn)行歸屬分析，其結(jié)果可為SMX分子的檢測及其振動光譜研究提供參考. 同時對SMX分子結(jié)構(gòu)性質(zhì)進(jìn)行了分析，通過計(jì)算SMX分子表面靜電勢和前線分子軌道對反應(yīng)活性位點(diǎn)進(jìn)行預(yù)測，為以后SMX分子的降解過程研究提供理論參考.

2 實(shí)驗(yàn)和理論計(jì)算

磺胺甲惡唑購自上海易恩化學(xué)技術(shù)有限公司，粉末狀. 分子式為C10H11N3O3S，分子量為253.3，純度 ≥ 98 %. 紅外光譜實(shí)驗(yàn)選用BRUKER公司Tensor 27型傅里葉紅外光譜儀，KBr壓片處理法. 拉曼光譜實(shí)驗(yàn)選用HORIBA公司LabRAM HR Evolution型拉曼光譜儀，激光器波長選擇785 nm.

理論計(jì)算采用Gaussian 16軟件包[11]，分子初始構(gòu)型通過Chemoffice 18和GaussView 6.0軟件構(gòu)建. 在氣相環(huán)境下，采用DFT中雜化密度泛函B3LYP方法，選擇B3LYP/6-311 + + G (d，p)基組進(jìn)行結(jié)構(gòu)優(yōu)化和頻率計(jì)算. 頻率計(jì)算中添加了關(guān)鍵詞“freq = intmode”，計(jì)算每個振動模式的勢能分布[12]. 理論光譜的頻率修正選擇校正因子0.9630[13]，使用Multiwfn程序[14]將計(jì)算得到的拉曼活性轉(zhuǎn)換為拉曼光譜強(qiáng)度. 將計(jì)算得到的紅外和拉曼光譜與實(shí)驗(yàn)光譜進(jìn)行對比，通過GaussView 6.0對振動頻率歸屬進(jìn)行全面指認(rèn). 使用Multiwfn結(jié)合VMD程序[15]分析并繪制SMX分子表面靜電勢分布圖和高質(zhì)量分子軌道等值面圖.

3 結(jié)果與討論

3.1 SMX分子的幾何結(jié)構(gòu)

采用B3LYP/6-311 + + G (d，p)基組對初始結(jié)構(gòu)進(jìn)行幾何優(yōu)化，優(yōu)化后的SMX分子的幾何結(jié)構(gòu)如圖1所示，其中C、O、S、N、H分別分別采用黑、紅、黃、藍(lán)、灰進(jìn)行表示. 優(yōu)化結(jié)果收斂，無虛頻，表明優(yōu)化后的分子結(jié)構(gòu)為穩(wěn)定構(gòu)型，最終優(yōu)化后的分子能量為-1175.94 a.u..

圖1 優(yōu)化后的SMX分子結(jié)構(gòu)Fig. 1 Optimized structure of SMX molecule

由圖1可以看出，SMX分子結(jié)構(gòu)主要由磺酰胺基團(tuán)分別連接異惡唑環(huán)和苯環(huán)構(gòu)成. 1C連接一個甲基，14C連接一個氨基. SMX最終優(yōu)化后的部分幾何參數(shù)值如表1所示，包括鍵長、鍵角和二面角. 其中8S-7N和9C-8S鍵長分別為1.710 ?和1.783 ?. 其余主要鍵長基本約在1.400 ?左右，因此可以判斷，當(dāng)被其他試劑攻擊時發(fā)生化學(xué)反應(yīng)時，S-N和S-C鍵最可能發(fā)生斷裂. 7N-8S-9C鍵角為107.471°，8S-7N-3C-2C和8S-7N-3C-4N二面角分別為-50.517°和133.669°，12C-9C-8S-7N和16C-9C-8S-7N二面角分別為92.252°和-87.500°，9C-8S-7N-3C二面角為-64.193°，其他主要的二面角絕對值均在0°和180° 之間. 因此可以判斷與7N相連的異惡唑環(huán)甲基結(jié)構(gòu)以及8S相連的苯環(huán)氨基結(jié)構(gòu)均在一個平面上，但SMX分子整體為非平面結(jié)構(gòu).

表1 SMX分子優(yōu)化后的幾何參數(shù)

3.2 SMX的振動光譜分析

紅外光譜和拉曼光譜可以反應(yīng)分子官能團(tuán)和分子結(jié)構(gòu)等信息，對有機(jī)物的鑒定具有很好的參考價值[16-17]. 本文通過實(shí)驗(yàn)和理論相結(jié)合的方法，對SMX分子的紅外和拉曼光譜進(jìn)行了測定，得到相應(yīng)的振動頻率和強(qiáng)度，并將實(shí)驗(yàn)光譜和理論計(jì)算光譜進(jìn)行對比分析，繪制了4000-300 cm-1的光譜圖，SMX的實(shí)驗(yàn)和理論紅外光譜如圖2所示，SMX的實(shí)驗(yàn)和理論拉曼光譜如圖3所示. 通過對比發(fā)現(xiàn)實(shí)驗(yàn)光譜和理論計(jì)算光譜的峰型及振動頻率基本吻合，說明此計(jì)算方法的選擇是可取和可信的. 同時發(fā)現(xiàn)，存在部分峰的位置其實(shí)驗(yàn)和理論存在差異. 該原因可歸結(jié)為兩個方面，一是實(shí)驗(yàn)中SMX物質(zhì)以固體粉末形式存在，具有分子間作用力，而理論計(jì)算只是針對獨(dú)立單個分子進(jìn)行分析；二是實(shí)驗(yàn)設(shè)備本身會存在一定誤差[18].

圖2 SMX分子的理論和實(shí)驗(yàn)紅外光譜Fig. 2 Theoretical and experimental infrared spectra of SMX molecule

圖3 SMX分子理論和實(shí)驗(yàn)拉曼光譜Fig. 3 Theoretical and experimental raman spectra of SMX molecule

SMX分子在4000-300 cm-1范圍內(nèi)的振動頻率和歸屬如表2所示，包括理論計(jì)算未校正頻率、校正頻率、實(shí)驗(yàn)紅外頻率、實(shí)驗(yàn)拉曼頻率和振動模式的歸屬. 其中計(jì)算所得每個振動模式均有強(qiáng)度，但是一些振動并不是單一振動模式，而是幾個振動共同作用、相互影響的結(jié)果，彼此之間存在關(guān)聯(lián)，因此需要計(jì)算每個振動模式下的勢能分布(Potential energy distribution，PED). 通過Gauss View 6.0結(jié)合PED值(≥10%)對每個振動模式進(jìn)行了指認(rèn)歸屬，由表2可以看出在此頻率范圍內(nèi)共有66個簡正振動模式，根據(jù)振動集中程度，可將振動頻率分為3600-2800 cm-1和1700-300 cm-1兩個頻率區(qū)域.

表2 SMX分子的理論、實(shí)驗(yàn)振動頻率和歸屬

3.2.13600-2800cm-1區(qū)域

在3600-2800 cm-1振動區(qū)域內(nèi)存在11個振動模式，歸屬于N-H和C-H的伸縮振動. 在N-H伸縮振動中，3557 cm-1歸屬于氨基上17N-H的不對稱伸縮振動，實(shí)驗(yàn)紅外光譜中位于3467 cm-1，實(shí)驗(yàn)拉曼光譜中位于3465 cm-1；3456 cm-1歸屬于氨基上17N-H的對稱伸縮振動，實(shí)驗(yàn)紅外光譜中位于3378 cm-1，實(shí)驗(yàn)拉曼光譜中位于3379 cm-1；3420 cm-1歸屬于磺酰胺基團(tuán)上7N-H的伸縮振動，實(shí)驗(yàn)光譜中位于3298 cm-1，實(shí)驗(yàn)拉曼光譜中位于3300 cm-1. 在C-H伸縮振動中，3163 cm-1歸屬于異惡唑環(huán)上2C-H的伸縮振動，實(shí)驗(yàn)紅外光譜中位于3143 cm-1，實(shí)驗(yàn)拉曼光譜中位于3146 cm-1；該區(qū)域范圍內(nèi)，其他的振動模式均歸屬于C-H的伸縮振動，但在實(shí)驗(yàn)光譜中中均未檢測到，這可能是因?yàn)閷?yīng)的紅外強(qiáng)度或拉曼強(qiáng)度較低，加之固體粉末存在分子間作用力的影響，導(dǎo)致振動變?nèi)? 如，3080 cm-1歸屬于苯環(huán)上16C-H的伸縮振動，實(shí)驗(yàn)拉曼光譜中位于3072 cm-1，而實(shí)驗(yàn)紅外光譜中未檢測到；2976 cm-1歸屬于異惡唑環(huán)上6C-H的不對稱伸縮振動，但在實(shí)驗(yàn)紅外和拉曼光譜中均未檢測到. 研究可以發(fā)現(xiàn)，C-H伸縮振動的理論計(jì)算和實(shí)驗(yàn)光譜的頻率基本吻合，而N-H伸縮振動的理論計(jì)算和實(shí)驗(yàn)光譜的頻率存在一定差異，這主要是因?yàn)镹-H鍵極性較大，容易與其他分子產(chǎn)生牽制，導(dǎo)致頻率的偏移.

3.2.21700-300cm-1區(qū)域

在1700-300 cm-1振動頻率區(qū)域內(nèi)存在55個振動模式，并且理論計(jì)算光譜和實(shí)驗(yàn)光譜的振動頻率吻合度較高. 存在14個振動強(qiáng)度較大，實(shí)驗(yàn)紅外光譜和實(shí)驗(yàn)拉曼光譜中均具有較好的吻合. 如，1603 cm-1歸屬于氨基上17N-H的剪式彎曲振動和17N-14C的伸縮振動，實(shí)驗(yàn)紅外光譜中位于1611 cm-1，實(shí)驗(yàn)拉曼光譜中位于1609 cm-1；1587 cm-1歸屬于異惡唑環(huán)上1C = 2C的伸縮振動和2C-H的面內(nèi)搖擺振動，實(shí)驗(yàn)紅外光譜中位于1596 cm-1，實(shí)驗(yàn)拉曼光譜中位于1591 cm-1；1363 cm-1歸屬于甲基上6C-H的對稱彎曲振動，實(shí)驗(yàn)紅外光譜中位于1366 cm-1，實(shí)驗(yàn)拉曼光譜中位于1369 cm-1；1029 cm-1歸屬于氨基上17N-H的扭曲振動，實(shí)驗(yàn)紅外光譜中位于1027 cm-1，實(shí)驗(yàn)拉曼光譜中位于1026 cm-1；928 cm-1歸屬于苯環(huán)上C-H的扭曲振動，實(shí)驗(yàn)紅外光譜中位于927 cm-1，實(shí)驗(yàn)拉曼光譜中位于922 cm-1；684 cm-1歸屬于異惡唑環(huán)上2C-H的面外搖擺振動和C-C的面外搖擺振動，實(shí)驗(yàn)紅外光譜中位于685 cm-1，實(shí)驗(yàn)拉曼光譜中位于691 cm-1；其他振動模式由于對應(yīng)的振動強(qiáng)度較小和分子間作用力的影響或只具有單一光譜對應(yīng)或著均無對應(yīng). 如，1472 cm-1歸屬于苯環(huán)上C-C的伸縮振動和C-H的面內(nèi)搖擺振動，實(shí)驗(yàn)紅外光譜中位于1471 cm-1，實(shí)驗(yàn)拉曼光譜無對應(yīng)峰；624 cm-1歸屬于苯環(huán)上C-C的彎曲振動，拉曼實(shí)驗(yàn)光譜中位于629 cm-1，實(shí)驗(yàn)紅外光譜無對應(yīng)峰；1420 cm-1歸屬于甲基上6C-H的不對稱彎曲振動，實(shí)驗(yàn)光譜中均無對應(yīng)峰.

3.3 SMX分子表面靜電勢分析

空間中某點(diǎn)的靜電勢(Electrostatic potential，ESP)指的就是將單位正電荷從無窮遠(yuǎn)移動到該點(diǎn)所做的功，描述的是電子之間以及電子與核電荷間的相互作用. 分子在化學(xué)反應(yīng)初始階段大都通過靜電作用縮短距離，所以靜電勢在研究分子之間靜電作用、預(yù)測反應(yīng)活性位點(diǎn)等方面具有重要意義[19]. 通常，分子范德華表面上的靜電勢用來預(yù)測親電和親核反應(yīng)活性位點(diǎn). 可以通過定量分子表面分析確定靜電勢在分子范德華表面上的極小點(diǎn)和極大點(diǎn)，距離分子表面上靜電勢最小點(diǎn)(最大點(diǎn))位置最近的原子即為最有可能發(fā)生親電(親核)反應(yīng)的位點(diǎn). 采用Multiwfn結(jié)合VMD分析和繪制了SMX分子的表面靜電勢分布如圖4所示，從藍(lán)色到紅色代表著靜電勢由負(fù)值到正值的變化，藍(lán)色區(qū)域容易被親電試劑攻擊，紅色區(qū)域容易被親核試劑攻擊. 通過定量分析分子表面電勢得出，SMX分子在當(dāng)前計(jì)算模型下存在11個極小值點(diǎn)和11個極大值點(diǎn). 其中，ESP最小值為-38.32 kal/mol，位于磺酰胺基團(tuán)區(qū)域；ESP最大值為 + 47.07 kal/mol，位于氨基區(qū)域. 由此可以得出，在化學(xué)反應(yīng)初始階段，親電試劑和親核試劑通過靜電作用將攻擊磺酰胺基團(tuán)區(qū)域和氨基區(qū)域.

圖4 SMX分子的靜電勢分布Fig. 4 The electrostatic potential distribution of SMX molecule

3.4 SMX的前線分子軌道分析

最高占據(jù)分子軌道(Highest occupied molecular orbital，HOMO)和最低未占據(jù)分子軌道(Lowest unoccupied molecular orbital，LUMO)稱為前線分子軌道[20]. 前線分子軌道在反應(yīng)電子的能級躍遷和預(yù)測化學(xué)反應(yīng)活性位點(diǎn)等方面有重要的意義. 分子軌道成分的計(jì)算實(shí)質(zhì)上是考察某軌道上電子密度在各個原子的分布情況. 采用Multiwfn結(jié)合VMD分析和繪制了SMX的高質(zhì)量分子軌道等值面圖，如圖5所示. SMX分子的HOMO為軌道66，能量為-0.24 a.u.；LUMO為軌道67，能量為-0.05 a.u.. 因此HOMO和LUMO之間的能隙差(Egap)為0.19 a.u.，表明電子有從HOMO向LUMO躍遷的能力. 親電反應(yīng)和親核反應(yīng)也可以通過HOMO和LUMO的軌道成分進(jìn)行預(yù)測，而HOMO和LUMO的軌道成分主要分布在苯環(huán)位置，所以苯環(huán)相較于異惡唑環(huán)更容易發(fā)生親電反應(yīng)和親核反應(yīng).

圖5 SMX分子的前線分子軌道Fig. 5 Frontier molecular orbitals of SMX molecule

4 結(jié) 論

采用DFT中雜化密度泛函B3LYP方法，在B3LYP/6-311 + + G (d，p)基組水平上對SMX分子進(jìn)行了結(jié)構(gòu)優(yōu)化，對鍵長、鍵角和二面角進(jìn)行了計(jì)算. 在相同計(jì)算水平下對SMX分子的紅外光譜和拉曼光譜進(jìn)行了分析，并通過GaussView 6.0對不同頻率峰的振動模式進(jìn)行了指認(rèn)歸屬. 將理論計(jì)算得到的紅外和拉曼光譜與實(shí)驗(yàn)獲得的光譜進(jìn)行了對比分析，發(fā)現(xiàn)兩者較為吻合. 在結(jié)構(gòu)優(yōu)化基礎(chǔ)上，繪制了SMX分子的表面靜電勢分布圖，得到SMX分子靜電勢最小值位于磺酰胺基團(tuán)區(qū)域，其值為-38.32 kal/mol；ESP最大值位于氨基區(qū)域，其值為 + 47.07 kal/mol. 通過前線分子軌道計(jì)算，得出HOMO和LUMO之間能隙差為0.19 a.u.，軌道成分的分布結(jié)果分析得出苯環(huán)部分的反應(yīng)活性要高于異惡唑環(huán)部分. 該研究為SMX的振動光譜檢測儲備了基礎(chǔ)數(shù)據(jù)，為其在抗生素藥品殘留檢測領(lǐng)域的應(yīng)用及其結(jié)構(gòu)鑒定奠定了理論基礎(chǔ).