新藥萘夫西林藥理性質的密度泛函研究

徐 洋，羅海霞，范 強, 鐘愛國

新藥萘夫西林藥理性質的密度泛函研究

徐 洋，羅海霞，范 強, 鐘愛國

（臺州學院醫藥化工學院， 浙江 臺州 316000）

使用密度泛函理論（DFT），在B3LYP方法和6-31+G(d)基組上，對抗菌消炎新藥萘夫西林分子進行結構優化，對其分子構型、偶極距、紅外光譜、紫外光譜、NBO電荷分布進行了分析。結果表明，萘夫西林發揮其藥理和毒理作用的最大可能部位在酰胺基的N原子上。

萘夫西林；密度泛函理論；定量計算

萘夫西林英文名為Nafcilline、Unipen 、Naftopen等，為以 6-氨基青霉素烷酸（6－APA）化學半合成的一種既口服也可注射的新型青霉素[1-3]。為微黃白色或白色粉末，有引濕性，易溶于水、乙醇、氯仿，微臭。在酸和青霉素酶中穩定。主要試用于治療心內膜炎、急性化膿性的腦膜炎、肺炎、敗血癥（金黃色葡萄球菌引起）等。與其他青霉素類藥物一樣的是，它也存在過敏反應。抗耐藥金葡菌的活性與苯唑青霉素比較接近,比甲氧苯青霉素要強,對溶血性鏈球菌和肺炎鏈球菌的活性比苯唑西林強。特別的是，對溶血性鏈球菌和肺炎球菌的效果非常好。口服的吸收效果不是很好，尿液排泄量少，血藥濃度也較低。通過肌肉注射萘夫西林的血藥濃度高于口服，但是仍然比較低，膽汁中的藥的濃度高于血藥濃度，且90 % 都經膽汁排泄排出。血漿蛋白結合率高,約為88 %。阿拉賓度公司（Aurobindo）在2012年12月29日宣布了萘夫西林注射液（nafcillin injections）已獲得了FDA的批準， 該公司也已馬上啟動了該藥物的裝運工作。目前，用密度泛函研理論究萘夫西林方面的文獻還比較缺乏[4,5]。

本文是采用密度泛函理論方法, 對萘夫西林的分子光譜如紅外吸收，拉曼吸收，紫外-可見吸收，核磁共振吸收以及自然鍵軌道理論(NBO)等進行了理論模擬和指認，通過先采用GaussView 3.7構造出萘夫西林的分子結構，再提交給Gaussian 09W 用相應的方法和基組計算光譜數據，取得了與實驗相吻合的結果， 為萘夫西林的分析與檢測提供幫助。并結合所得數據，討論分析了萘夫西林的藥理及毒理性質。

1 實驗部分

采用GaussView看圖軟件，構建了萘夫西林分子結構模型，采用 DFT 理論的 B3LYP 方法，在6-31+G(d)基組水平上，進行了結構優化計算。在優化結構的基礎上，采用頻率分析方法，進行了頻率分析。表明所有振動頻率全部為正，表明其計算結果是可以相信的。全部計算工作在 Gaussian 09W程序及 PC 機上完成。圖 1顯示了萘夫西林（C21H22N2O5S）分子結構和原子編號。

2 結果與討論

2.1 分子的幾何結構分析

用Gaussian 09W程序在B3LYP/6-31+G(d)基組上進行優化，模擬得到了萘夫西林分子的鍵長, 鍵角等數據。正常的C-O單鍵的鍵長（（1.20～1.43）×10-10m），C-C單鍵（1.54×10-10m左右），C=O雙鍵（1.20×10-10m左右），C-N單鍵（1.48×10-10m），C-S單鍵（1.82×10-10m），O-H（0.98×10-10m），N-H（1.01×10-10m）表中列出了優化后的分子部分鍵長數據表明：除C4-C1的鍵長偏小一點外，其他鍵長均屬于正常范圍。表中列出的重要鍵角數據表明，稠雜環的大環(N1-C2-S6-C7-C3原子組成的五元雜環)和小環(N1-C2-C5-C4原子組成的四元雜環)。

圖1 萘夫西林分子結構及編號Fig.1 Structure and atomic of nafcilline

2.2 紅外吸收光譜

物質的IR光譜是其分子骨架結構的反映，圖2譜圖中的吸收峰與分子中各基團振動形式相對應。

圖2 萘夫西林分子紅外吸收光譜Fig.2 IR of nafcilline

使用 DFT/B3LYP/6-31+G(d) 軟件模擬顯示, 萘夫西林藥物分子在指紋區 616 cm-1(弱，2S, 5R, 6R結構)處存在紅外吸收峰，在特征區1 318 cm-1(中，甲基C-H的彎曲振動吸收峰)，1 516 cm-1(弱，萘環上的C-C骨架振動吸收峰), 1 777 cm-1（強，羰基C=O鍵伸縮振動吸收峰）, 3 064 cm-1（弱，萘環上的C-H伸縮振動吸收峰）, 3 586 cm-1（弱，萘酰胺基上的N-H伸縮振動吸收峰）, 3 676 cm-1（弱，羧基上的羥基O-H的伸縮振動吸收峰）處顯示較強的IR吸收峰，這些峰與萘夫西林小分子的實驗IR譜（分別在608，1 305，1 556，1 787，3 101，3 576，3 656 cm-1處有紅外吸收峰）基本是吻合的。

2.3 拉曼吸收光譜

分子的拉曼吸收效應起源于分子振動與轉動,從拉曼吸收光譜中可得到與分子振動互補的信息。采用DFT/B3LYP 方法，在 6-31+G(d)基組上進行分子優化和頻率分析，得到7個拉曼吸收特征峰(如圖3)。我們對分子的7個特征峰的分子振-轉吸收模式進行了歸屬指認。萘夫西林（（2S,5R,6R)-6-(2-乙氧基-1-萘酰氨基)-3,3-二甲基-7-羰基-4-硫雜-1-氮雜-二環[3.2.0]庚烷-2-羧酸）是一個稠雜環的化合物，也可以看作是萘環的β位上氫原子被一個乙氧基所取代，α位上的氫原子被一個由酰胺基和稠雜環構成的支鏈（青霉素母核－6－氨基青霉烷酸（6－APA））所取代而形成的雙取代萘類分子。圖 3是理論模擬萘夫西林分子的拉曼散射譜特征吸收峰，其主要位于3 676 cm-1（弱，羧基上的游離羥基O-H的伸縮振動），3 586 cm-1（弱，萘酰胺基上的N-H伸縮振動）,3 064 cm-1（強，二取代萘環上的C-H對稱呼吸振動）,1 777 cm-1（弱，羰基C=O鍵伸縮振動）, 1 516 cm-1(中，萘環上的C-C骨架振動)，1 318 cm-1(中，甲基C-H的彎曲振動)，616 cm-1(弱，環內振動) 。對比實驗和理論拉曼吸收光譜，它們的拉曼吸收峰基本是吻合的。

圖3 模擬萘夫西林分子拉曼吸收光譜Fig.3 Laman of nafcillin

2.4 紫外吸收光譜

采用 TD DFT/B3LYP/6-31+G(d)理論方法, 模擬顯示了萘夫西林分子在 247 nm (強吸收)處有強的紫外吸收峰，這是源于電子從最高占據軌道(HOMO)躍遷到最低空軌道(LUMO)所顯示的吸收峰，在367 nm（弱吸收）處有弱的紫外吸收峰，這是源于電子從次高占據軌道(HOMO-1)躍遷到最低空軌道(LUMO)所顯示的峰。這與其標準圖譜顯示萘夫西林分子在247.5 nm和357.5 nm 處有較強的紫外吸收峰基本是相吻合(見圖4)。

圖4 萘夫西林分子紫外-可見吸收光譜Fig.4 UV-Vis of nafcilin

2.5 核磁共振光譜

用TMS B3LYP/6-311+G(2d,p) GIAO的方法理論模擬了萘夫西林分子的氫核磁共振譜（見圖5所示）。

其萘環上六個氫(弱，7.218 9×10-6)，乙氧基上的CH2(弱，3.581 4, 3.917 9×10-6); CH3(中,1.328 5，1.356 7，0.786 7×10-6); 酰胺鍵上的氫(弱，5.085 5 ×10-6) ，稠雜環6號位上的H（弱，4.216 9×10-6），5號位上的H(5.267 5×10-6),與3號位相連的2個CH3（強，1.521 3，1.291 1 ，1.199 7，1.292 8，1.269 7，1.157 8×10-6)。這些吸收峰較好地與其核磁共振實驗吸收峰值吻合起來，驗證了模擬計算的準確性[6-7]。

圖5 萘夫西林的核磁共振氫譜Fig.5 1-HNMR of nafcilline

2.6 偶極矩和原子的自然電荷

分子的偶極矩是衡量其分子極性的重要物理量。采用 Gaussian 09W，我們計算得到標題分子的偶極矩為5.5418 Debye，這表明萘夫西林的極性比較大，血藥濃度也較低，通過肌肉注射萘夫西林的血藥濃度高于口服，但是仍然比較低。血漿蛋白結合率高，膽汁中的藥的濃度高于血藥濃度，所以尿液排泄量少，主要經膽汁排泄排出。從文獻上查到的疏水參數計算參考值（logP）為2.9，這可能是因為萘夫西林的羧基的存在，增大了萘夫西林的親水性和溶解度。

原子的自然電荷值是研究其化學活性點位及確定親核或親電特性強弱的一種有效方法。采用B3LYP/6-31+G(d)基組，我們分析了萘夫西林分子的自然鍵軌道值(NBO)。結果顯示電荷最正的是羧基上的碳原子（8C，0.831e）；其次為β-內酰胺環上的羰基上的碳原子（4C，0.711e）；第三為酰胺鍵上的碳原子（15C，0.703e）；第四為酰胺鍵上的氫原子（20H，0.433e），電荷最負的為羧基上的羥基氧（14O，-0.708e），其次為酰胺鍵上的氮原子（10N，-0.697e）；第三為酰胺鍵上的氧原子（16O，-0.687e）；第四為β-內酰胺環上的羰基上的氧原子（9O，-0.686e）。所以β-內酰胺環及與其相連的酰胺鍵（-CO-NH-）上的原子很可能是其發揮藥理和藥理活性的親電和親核反應中心。這也符合了青霉素類抗生素的藥理是由于β-內酰胺類作用于細菌的細胞壁而產生的。

3 結 論

對萘夫西林分子(（2S,5R,6R)-6-(2-乙氧基-1-萘酰氨基)-3,3-二甲基-7-羰基-4-硫雜-1-氮雜-二環[3.2.0]庚烷-2-羧酸），我們基于DFT理論和方法，在 6-31+G(d)基組水平上進行了分子結構優化。對其藥物分子光譜(UV-Vis，IR, 1HNMR )等理化性質進行了理論模擬和指認，取得了與實驗光譜基本吻合的結果，表明了本研究方法是可信的。模擬的幾何優化結構表明, 萘夫西林分子具有穩定的分子構型和良好的疏水參數，在極性溶液中溶解能力將會較強。自然電荷（NBO）值計算結果表明，β-內酰胺環及與其相連的酰胺鍵（-CO-NH-）上的N原子很可能是其發揮藥理和藥理活性的親電和親核反應中心，符合青霉素類抗生素的藥理是由于β-內酰胺類作用于細菌的細胞壁而產生。對其分子的模擬光譜及其指認，將有助于萘夫西林類藥物的分析與檢測。

［1］昝林森, 邱懷, 袁志發. 環核苷酸(CNT)在豚鼠體內代謝規律的研究[J]. 畜牧獸醫學報,2000，12(2):39-45.

［2］張坐省, 傅建熙. β—受體興奮劑在畜牧業生產上的意義[J]. 畜牧獸醫雜志, 1996,14(4):34-38.

［3］涂秀英, 曹智勇, 聶復禮. 絞股藍抗衰老藥理實驗研究[J]. 江西中醫學院學報, 1999, 25(2):12-16.

［4］曾淑君，沈寶蓮，文良珍. 云芝糖肽對裸鼠人肺腺癌抗癌作用的研究[J]. 中國藥理學通報, 1995, 3(12):12-15.

［5］李敬芝, 姚鋒. 老藥新用[J]. 中國醫刊, 1980, 12(1):22-25.

［6］彭聚寶.急性痛風的非秋水仙堿治療[J]. 中國醫院藥學雜志, 1992, 31(2):21-24.

［7］黃慶國, 韓朔睽, 王連生. 利用前線分子軌道能預測有機物生物毒性的方法[J].環境化學, 1994, 12(1):23-26.

Density Functional Theory Study on Pharmacological Properties of Nafcilline

XV Yang，LUO Hai-xia,FAN Qiang,ZHONG Ai-guo
（College of Chemical Engineering and Pharmacy, Taizhou University, Zhejiang Taizhou 316000，China）

Using the density functional theory B3LYP method with 6-31+G(d) basis set,nafcilline molecules were optimized, and their molecular configuration, molecular dipole moment, IR spectra, UV spectra, NBO charge distribution were analyzed. The results show that, the maximum possible parts of pharmacological and toxicological effects of nafcilline are N atoms on the amide.

Nafcilline; Density functional theory; Quantitative calculation

TQ 624.2

： A

： 1671-0460（2015）04-0699-03

浙江省大學生科技創新項目（新苗人才計劃），項目號：2014R428015。

2015-03-15

徐洋（1991-），男，浙江紹興人，臺州學院高分子專業，研究方向：從事高分子研究工作。E-mail：451872652@qq.com。