拓展實驗阿莫西林的譜學(xué)分析與指認

洪文燕　施琪浩　張英杰　蔣鑫　陳詩雨　章少陽　陳茹箐　秦利明　鐘愛國

摘????? 要：用密度泛函理論（DFT）的方法，在B3YLYP/6-31+G水平下優(yōu)化阿莫西林的結(jié)構(gòu)，通過改變其椅式構(gòu)型上1-H、3-H、4-H的直立鍵和平伏鍵相互轉(zhuǎn)變，發(fā)現(xiàn)-OH化學(xué)位結(jié)構(gòu)鍵的轉(zhuǎn)變對其分子溶解度（logS）值沒有明顯影響;采用密度泛函理論的方法，在B3YLYP/6-31+G水平下優(yōu)化了7種阿莫西林衍生物的結(jié)構(gòu)，通過計算得到相關(guān)性最高ESP-S5電荷值與logS值的線性方程（y=13.51x-0.72， R²=0.912）。還采用密度泛函理論的方法預(yù)測了7種阿莫西林衍生物紅外吸收、拉曼吸收、紫外吸收、熒光發(fā)射、磷光發(fā)射和核磁共振譜吸收光譜圖，并對其特征峰進行指認和歸屬。為制藥專業(yè)學(xué)生拓展訓(xùn)練實驗識別阿莫西林及其衍生物以及構(gòu)效關(guān)系提供了基礎(chǔ)數(shù)據(jù)。

關(guān)? 鍵? 詞：制藥專業(yè)實驗;密度泛函理論; 阿莫西林;溶解度;譜學(xué)分析

中圖分類號：G822.8?????? 文獻標(biāo)識碼： A?????? 文章編號： 1671-0460（2020）07-1321-05

Spectroscopic Analysis and Identification of Amoxicillin

for Pharmaceutical Professional Expansion Experiments

HONG Wen-yan， SHI Qi-hao， ZHANG Ying-jie， JIANG Xin， CHEN Shi-yu，

ZHANG Shao-yang， CHEN Ru-qian， QIN Li-min， ZHONG Ai-guo*

（College of Pharmaceutical， Chemical and Materials Engineering， Taizhou University， Taizhou Zhejiang 318000， China）

Abstract： The structure of amoxicillion was optimized by the density functional theory（DFT） at B3LYP/6-31+G level. By changing the structure of 2-H， 3-H， 4-H， the results showed that the change of chemical bond had no significant effect on the solubility （logS） value. The method of density functional theory was used to optimize the structures of 7 amoxicillin derivatives at B3YLYP/6-31+G level. And the linear equation （y=13.51x-0.72） with the highest correlation between the ESP-S5 charge value and the logS value was obtained by calculating the structure of the amoxicillion derivative. At the same time， the method of density functional theory was used to predict the infrared， Raman， UV-Vis and 1H-NMR spectra of seven kinds of amoxicillion derivatives， and the characteristic peaks were identified and assigned. This paper can provide powerful reference for the study of the spectral prediction and structure activity relationship of other derivatives.

Key words： Pharmaceutical engineering experiments; DFT; Amoxicillion; Solubility; Spectral enalysis

阿莫西林（AMXL）為人工半合成的耐酸和堿的廣譜殺菌青霉素類抗生素類藥物。該分子藥抗菌抗病毒譜廣，殺菌力較強^[1]。阿莫西林也是醫(yī)院臨床用藥中一種常見的具有氨基側(cè)鏈的青霉素分子藥物。它在常溫常壓態(tài)下也是一種結(jié)晶狀的白色粉末藥品，能微溶于溫水。它在醫(yī)院臨床上應(yīng)用范圍非常廣泛和用料大^[2]。如在病人消化不良、病人慢性支氣管炎、病人膀胱炎、病人皮疹以及病人消化性潰瘍等疾病的治療等方面，皆需要應(yīng)用到該種粉末狀藥物。但是，一般的病人患者服用該種藥物后，極有可能會出現(xiàn)病人惡心、病人嘔吐、病人腹瀉、病人皮疹以及蕁麻疹等多種不良反應(yīng)發(fā)生^[3]。阿莫西林的臨床抗菌活性和氨芐西林非常相似，但阿莫西林的耐酸性更強，殺菌效果及殺菌速度也強于氨芐西林^[4]，但是阿莫西林用于治療腦膜炎患者會引起嚴重不良反應(yīng)，因而阿莫西林禁用于治療腦膜炎患者。阿莫西林在體內(nèi)的半衰期為62 min^[5]。阿莫西林的可吸收性較好，在酸性條件下可穩(wěn)定吸收，且腸道吸收率高達90%，吸收速度遠遠大于氨芐西林。阿莫西林能夠穿透細胞壁， 和大分子量青霉素結(jié)合蛋白結(jié)合， 與分子內(nèi)酚胺可以水解， 與菌體轉(zhuǎn)膚酶結(jié)合， 讓菌體失去活性， 切斷菌體合成糖膚，細菌細胞壁被破壞后， 由于菌體失去細胞壁而脹破滅亡^[6]。哺乳動物細胞無細胞壁， 所以阿莫西林對人體毒性小， 但仍會有惡心、嘔吐等不良反應(yīng)存在。王晨等在《中國新藥雜志》利用計算機模擬技術(shù)對阿莫西林膠囊溶出度方法的研究，利用計算機模擬軟件建立并驗證阿莫西林在人體內(nèi)的生理藥動學(xué)（physiologically based pharmacokinetic， PBPK）模型，用于阿莫西林膠囊溶出度方法的評價研究。方通過溶出特征曲線優(yōu)化，建立并優(yōu)化PBPK模型。通過模型計算阿莫西林膠囊藥時曲線對比人口服實測藥時曲線計算結(jié)果的準確性，基于該PBPK模型，以去卷積分算法模擬計算出阿莫西林膠囊在人體內(nèi)的釋放曲線，與現(xiàn)行版《中華人民共和國藥典》和《美國藥典》中阿莫西林膠囊溶出度限度進行比較與評價，證明上述藥典中各溶出度限度均高于制劑體內(nèi)累計吸收量要求?，F(xiàn)行方法是基于阿莫西林膠囊制劑的不同處方，用于進行工藝控制，保證制劑生產(chǎn)一致性的指標(biāo)，目的在于降低制劑工藝對產(chǎn)品質(zhì)量波動的影響風(fēng)險，而非緊密關(guān)聯(lián)體內(nèi)有效性水平的限度指標(biāo)。高光武等在《世界最新醫(yī)學(xué)信息文摘》發(fā)表論文，針對阿莫西林與莫西沙星治療幽門螺桿菌胃病的有效性觀察，評價阿莫西林與莫西沙星治療幽門螺桿菌胃病的有效性，分析其藥理作用。他們選取本院2017年4月至2019年3月收治的146例幽門螺桿菌胃病患者，其中73例使用莫西沙星進行治療的患者作為研究組，另外使用使用阿莫西林進行治療的73例患者作為對照組。將兩組患者的臨床治療效果、幽門螺桿菌（HP）根除率、藥物不良反應(yīng)等指標(biāo)，進行數(shù)據(jù)的統(tǒng)計記錄，并使用統(tǒng)計學(xué)軟件進行統(tǒng)計數(shù)據(jù)的處理和比較分析。結(jié)果使用莫西沙星進行治療的研究組患者，其臨床治療總有效率、幽門螺桿菌根除率指標(biāo)均顯著高于對照組患者;治療后胃痛癥狀評分、不良反應(yīng)發(fā)生率兩項指標(biāo)，研究組患者均顯著低于對照組。兩組患者組間數(shù)據(jù)比較差異有統(tǒng)計學(xué)意義（P<0.05）。結(jié)論本次研究數(shù)據(jù)的結(jié)果表明，對于幽門螺桿菌胃病患者，使用口服莫西沙星進行治療相比口服阿莫西林，臨床效果更加顯著。高曉華等在《中國現(xiàn)代藥物應(yīng)用》發(fā)表了阿莫西林與莫西沙星治療幽門螺桿菌胃病的臨床效果探討，他們探討阿莫西林與莫西沙星治療幽門螺桿菌胃病的臨床效果。用100例幽門螺桿菌胃病患者，隨機分為阿莫西林組與莫西沙星組， 每組50例。阿莫西林組采取阿莫西林治療， 莫西沙星組采取莫西沙星治療。比較兩組患者的治療效果; 根除幽門螺桿菌時間、幽門螺桿菌胃病愈合時間、住院時間; 治療前后生活質(zhì)量評分、胃腸道癥狀積分;不良反應(yīng)發(fā)生情況。結(jié)果莫西沙星組的治療總有效率100.00%高于阿莫西林組的76.00%，差異具有統(tǒng)計學(xué)意義（P<0.05）。治療后，莫西沙星組的生活質(zhì)量評分（98.81±1.55）分高于阿莫西林組的（85.68±2.99）分，胃腸道癥狀積分（0.45±0.01）分低于阿莫西林組的（1.42±0.21）分，差異均具有統(tǒng)計學(xué)意義（P<0.05）。莫西沙星組的根除幽門螺桿菌時間、幽門螺桿菌胃病愈合時間、住院時間分別為（6.50±1.02）、（6.50±1.09）、（8.50±2.01）d，均短于阿莫西林組的（7.35±2.28）、（8.50±2.99）、（10.54±2.08） d， 差異均具有統(tǒng)計學(xué)意義（P<0.05）。兩組不良反應(yīng)發(fā)生率比較，差異無統(tǒng)計學(xué)意義（P>0.05）。阿莫西沙星治療幽門螺桿菌胃病的臨床效果顯著，能夠縮短治療時間，提高患者的生活質(zhì)量，改善臨床癥狀，值得推廣應(yīng)用。吳春華等人在《抗感染藥學(xué)》雜志發(fā)表論文，探討了阿莫西林和替硝唑三聯(lián)療法對幽門螺桿菌陽性的慢性胃炎患者的療效及其對Hp清除的影響，評價了阿莫西林、替硝唑、奧美拉唑三聯(lián)療法對幽門螺桿菌（Hp）陽性的慢性胃炎患者的療效及其對炎癥因子和Hp清除的影響。他們選取2017年9月—2018年9月間收治的Hp陽性慢性胃炎患者100例資料，按治療方法的不同將其隨機分為對照組（50例）和觀察組（50例）;對照組患者給予奧美拉唑和阿莫西林治療，觀察組患者在對照組基礎(chǔ)上加用替硝唑治療;比較兩組患者治療前后血清白細胞介素-6（IL-6）、腫瘤壞死因子-α（TNF-α）水平測得值的變化情況，以及治療后的總有效率及Hp清除率的差異。結(jié)果表明，兩組患者治療前血清IL-6、TNF-α水平測得值經(jīng)組間比較其差異均無統(tǒng)計學(xué)意義（P>0.05），治療后優(yōu)于治療前（P<0.05）;觀察組患者治療后血清IL-6、TNF-α水平測得值均低于對照組（P<0.05）;觀察組患者治療后的總有效率及Hp清除率均高于對照組（P<0.05）。采用奧美拉唑、阿莫西林、替硝唑三聯(lián)療法治療Hp陽性慢性胃炎患者的療效較為確切，有利于促進炎癥的改善，能明顯提高了其對Hp的清除率。王德宜在《臨床醫(yī)藥文獻電子雜志》上，針對阿莫西林所致不良反應(yīng)及藥學(xué)臨床分析，他們探討了阿莫西林所致不良反應(yīng)及藥學(xué)臨床分析。他們選取2018年1月-2019年1月院收治的阿莫西林治療不良反應(yīng)患者50例作為研究對象，回顧性分析所有患者阿莫西林的使用量、用藥情況，以及患者不良反應(yīng)發(fā)生時間、表現(xiàn)等藥學(xué)臨床反應(yīng)。發(fā)現(xiàn)本次研究患者服藥的目的都是為了預(yù)防感染，而用藥途徑則是以口服為主。用藥1～7 d出現(xiàn)不良反應(yīng)現(xiàn)象的患者明顯要高于其他時間段，不良反應(yīng)表現(xiàn)則是過敏為主，相較于它癥狀的不良反應(yīng)情況發(fā)生數(shù)量而言，數(shù)據(jù)比較，差異有統(tǒng)計學(xué)意義（P<0.05）。結(jié)論患者在使用阿莫西林的時候，一定要確保用藥的規(guī)范與合理性，結(jié)合患者實際情況來合理用藥，同時在用藥過程中做好不良反應(yīng)測試與規(guī)避工作，這樣才能降低這一藥物對于患者所造成的不良影響。

本文采用密度泛函理論的方法，B3YLYP/ 6-31+g水平下優(yōu)化阿莫西林的結(jié)構(gòu)，通過改變其椅式構(gòu)型上1-H、3-H、4-H的直立鍵和平伏鍵相互轉(zhuǎn)變來判斷化學(xué)結(jié)構(gòu)鍵對溶解度（logS）值是否存在影響。鑒于先前已經(jīng)研究采用密度泛函理論來預(yù)報與測量分子藥物衍生物的pK_b值，并且已經(jīng)得到了良好的線性相關(guān)系數(shù)關(guān)于分子藥物pK_b預(yù)測的線性方程式。因此本論文將采用密度泛函理論和密度泛函活性理論（DFT and DFRT）的方法與方式，B3YLYP/6-31+g水平下優(yōu)化阿莫西林衍生物的結(jié)構(gòu)，通過模擬計算來預(yù)測阿莫西林其他衍生物的溶解度的作用。此外本文還采用密度泛函理論模擬預(yù)測了阿莫西林一元取代物的紅外、拉曼、紫外和核磁共振譜圖，對譜圖中的特征峰進行指認和歸屬，這對阿莫西林的化工合成、光譜與理論性質(zhì)預(yù)測與構(gòu)效關(guān)系研究提供理論支持。本文試圖從波譜學(xué)和分子溶解性兩個角度分析其特征，所得結(jié)果可為AMXL檢測提供基礎(chǔ)數(shù)據(jù)參考^[7-10]。

1? 計算方法

在P4-DELL電腦上安裝Chem Draw 8.0軟件、Chem 3D軟件、Gauss View 5.0.9軟件、Gaussian 09W軟件、Origin 8.5等插件和軟件。用 Chem3D 模塊和 ChemDraw 8.0 軟件模塊，構(gòu)建了7個AMXL極其衍生物的模型分子，首先使其分子優(yōu)化結(jié)構(gòu)，使其能量最低化;然后在GaussView 5.0.9 軟件模塊中，先使其能量構(gòu)型優(yōu)化到最低，然后使用密度泛函理論DFT/ B3LYP 與方法，借助6-31+G 水平下優(yōu)化分子模型結(jié)構(gòu)，其中計算關(guān)鍵詞分別為iop（5/33）=1， pop=（nbo， chelpg， hirshfeld） 的約束條件下，計算出7個AMXL及其衍生物分子的原子電荷。我們通過變換AMXL結(jié)構(gòu)模型位置上的取代基團以及其取代基位置來獲得相關(guān)的衍生物分子，進一步來研究其4種電荷：Mulligen、ESP、NBO和HIR電荷變化;采用ACD Lab6.0軟件測得的logS的線性關(guān)系。本文使用了—NH₂、—CH₃作為取代基團在R₁、R₂、R₃、R₄位置進行一元取代基、二元取代基，總共得到7種阿莫西林的衍生物。阿莫西林的結(jié)構(gòu)如圖1所示，使用分別將椅式結(jié)構(gòu)1-H、3-H、4-H的平伏鍵改成直立鍵，總共獲得7種阿莫西林的化學(xué)結(jié)構(gòu)鍵模型來研究在不同化學(xué)結(jié)構(gòu)鍵的條件下，5-S的4種電荷——Mulligen、ESP、NBO和HIR電荷值變化與采用ACD Lab6.0軟件測得logS值之間的關(guān)系。在此基礎(chǔ)上，借助量子化學(xué)計算可以提高譜峰歸屬準確性。

2? 結(jié)果與討論

2.1? AMXL及其衍生物的原子電荷與溶解度關(guān)系

本論文首先使用Chem Draw8.0模塊和Chem Draw 3D軟件模塊，構(gòu)建了8個阿莫西林及其衍生物分子的模型，運行軟件并使其結(jié)構(gòu)和能量最低化;在Gauss View 5.0.9軟件模塊中，使用B3LYP泛函方法在6-31+G（d，p）基組水平下優(yōu)化其分子結(jié)構(gòu)，我們的計算關(guān)鍵詞為pop=（nbo，chelpg，hirshfeld）的約束條件下，分別計算出8個阿莫西林衍生物的DFT和DFRT量化參數(shù)。記錄5-S的Mulligen、ESP、NBO和HIR的電荷值;運用ACD Labs6.0軟件，得到8種阿莫西林衍生物相應(yīng)水溶液中的溶解度（logS）預(yù)報值;然后分別將8種阿莫西林及其衍生物的Mulligen電荷、ESP電荷、NBO電荷和HIR的電荷值擬合，用Origin 8.5作圖，得到阿莫西林及其衍生物溶解度（logS）相應(yīng)的線性回歸方程，比較R值（相關(guān)系數(shù)）大小，來確定R值最大的就是預(yù)測方程。由表1可知，阿莫西林椅式結(jié)構(gòu)上H的平伏鍵變成直立鍵時，logS保持不變。這說明化學(xué)結(jié)構(gòu)鍵的變化對于阿莫西林的溶解度沒有影響;阿莫西林椅式結(jié)構(gòu)上H的平伏鍵變成直立鍵時，5-S的Mulligen、NBO和HIR的電荷值存在較小的變化，而5-S的ESP電荷值變化明顯（y=13.51x-0.72， R²=0.912）。

2.2? 紫外-可見吸收光譜

準確測定有機化合物的分子結(jié)構(gòu)，對從分子水平去認識物質(zhì)世界，推動近代有機化學(xué)的發(fā)展是十分重要的。采用現(xiàn)代儀器分析方法，可以快速、準確地測定有機化合物的分子結(jié)構(gòu)。在有機化學(xué)中應(yīng)用最廣泛的測定分子結(jié)構(gòu)的方法是四大光譜法：紫外光譜、紅外光譜、核磁共振和質(zhì)譜。紫外和可見光譜（ultraviolet and visible spectrum）簡寫為UV。在TD DFT B3LYP/6-31+G//OPT DFT/6-31+G水平下優(yōu)化AMXL分子結(jié)構(gòu)，計算得到阿莫西林紫外光譜。有機化合物分子中主要有三種電子：形成單鍵的σ電子、形成雙鍵的π電子、未成鍵的孤對電子，也稱n電子?；鶓B(tài)時σ電子和π電子分別處在σ成鍵軌道和π成鍵軌道上，n電子處于非鍵軌道上。僅從能量的角度看，處于低能態(tài)的電子吸收合適的能量后，都可以躍遷到任一個較高能級的反鍵軌道上。模擬結(jié)果顯示了阿莫西林在303.3 nm處有最大吸收波長，且位于近紫外區(qū)，這是由苯環(huán) π→π^* 躍遷電子所致（圖省去）。

2.3? 核磁共振吸收光譜

在密度泛函理論DFT B3LYP/6-31+G水平下優(yōu)化結(jié)構(gòu)，采用NMR-GIAO方法計算得到阿莫西林氫核共振譜圖（圖2，¹H-NMR）。苯環(huán)上4個氫（弱，（6.6～6.89）×10^-6）;主鏈上—CH（中，4.85×10^-6）;氨基上的H（中，2.0×10^-6）;甲基上的H（中，1.35×10^-6）;受旁邊的羰基、S原子等強吸引電子作用，其主鏈上的CH的化學(xué)位移發(fā)了變化（中，（4.72～4.85）×10^-6）。羧基中的H原子（弱，11.0×10^-6）。

2.4? 紅外和拉曼吸收光譜

拉曼光譜（Raman spectra），是一種散射光譜。拉曼光譜分析法是基于印度科學(xué)家C.V.拉曼（Raman）所發(fā)現(xiàn)的拉曼散射效應(yīng)，對與入射光頻率不同的散射光譜進行分析以得到分子振動、轉(zhuǎn)動方面信息，并應(yīng)用于分子結(jié)構(gòu)研究的一種分析方法。拉曼效應(yīng)起源于分子振動（和點陣振動）與轉(zhuǎn)動，因此從拉曼光譜中可以得到分子振動能級（點陣振動能級）與轉(zhuǎn)動能級結(jié)構(gòu)的知識。在密度泛函理論DFT B3LYP、6-31+G水平下優(yōu)化結(jié)構(gòu)，采用Frequency方法計算得到圖3結(jié)構(gòu)的阿莫西林紅外光譜圖3（a）。3 940 cm^-1為平面內(nèi)—OH伸縮振動;1 800 cm^-1為平面內(nèi)C=O伸縮振動;1 230 cm^-1平面內(nèi)—CH擺動，環(huán)上H對稱擺動;1 100 cm^-1平面外環(huán)H不對稱擺動。有些吸收峰在譜圖中均有出現(xiàn)，但其頻率較小，Origin 8.5軟件中不做波峰處理。

在密度泛函理論B3LYP水平下優(yōu)化結(jié)構(gòu)，采用Frequency方法計算得到圖3（b）結(jié)構(gòu)的阿莫西林拉曼光譜圖。3 770 cm^-1為平面內(nèi)-OH伸縮振動;1 800 cm^-1為平面內(nèi)C=O伸縮振動;1 540 cm^-1為平面內(nèi)環(huán)上H對稱擺動;1 060 cm^-1為平面外-CH₃對稱擺動;980 cm^-1為平面外環(huán)上H不對稱擺動。拉曼光譜中所有的波峰比實驗數(shù)值有所浮動，這是由于GaussView 5.0.9在模型建立中存在相對應(yīng)的誤差，但并不影響基本官能團的判斷。

2.5? 熒光和磷光發(fā)射光譜

通過 OPT/TD DFT方法，我們實現(xiàn)了 AMXL 建模和結(jié)構(gòu)優(yōu)化，熒光和磷光發(fā)射光譜圖頻率則通過Energy 方法計算得到，再借助 Origin 8.5 作圖分析（圖省去）顯示了AMXL 在 343.9 nm 處有最大熒光發(fā)射波長，且在近紫外區(qū)范圍內(nèi)。而其磷光圖發(fā)射峰在 323.8 nm 處，有較短發(fā)射波長原因是發(fā)生苯環(huán)分子內(nèi)發(fā)生電子 n→π* 躍遷。

3? 結(jié) 論

1）在DFT B3YLYP水平下優(yōu)化阿莫西林衍生物結(jié)構(gòu)，改變其1-H、3-H、4-H的化學(xué)結(jié)構(gòu)鍵位置，結(jié)果顯示直立/平伏鍵的轉(zhuǎn)變對阿莫西林的溶解度沒有影響;阿莫西林椅式結(jié)構(gòu)上H的平伏鍵變成直立鍵時，5-S的Mulligen、NBO和HIR的電荷值存在較小的變化，而5-S的ESP電荷值變化明顯（y=13.51x-0.72，R²=0.912）。

2）本文對7種不同的阿莫西林衍生物分別進行紅外、紫外、拉曼光譜、熒光和磷光和核磁共振進行了模擬預(yù)測，對以上光譜圖中的特征峰進行指認和歸屬，對AMXL及部分衍生物進行了圖譜模擬預(yù)測和描述，這對我院制藥專業(yè)學(xué)生學(xué)習(xí)AMXL衍生物檢測認知提供參考數(shù)據(jù)。值得在我院其他專業(yè)拓展訓(xùn)練中推廣與應(yīng)用。實踐證明該方法，對于培養(yǎng)優(yōu)秀的制藥工程專業(yè)的復(fù)合型人才，有十分積極的影響。

參考文獻：

[1]萬蕾. 阿莫西林的特性， 藥理和臨床應(yīng)用[J]. 黑龍江科技信息， 2016 （14）： 134-134.

[2]邸小雨. 關(guān)于阿莫西林的研究綜述[J]. 赤子（上中旬）， 2015 （5）： 299-301.

[3]尹彬. 康復(fù)新聯(lián)合奧美拉唑、阿莫西林膠囊三聯(lián)療法治療消化性潰瘍療效分析[J]. 中外醫(yī)學(xué)研究， 2010， 08 （26）： 17-18.

[4]陳格， 余瓊娥. 不同晶型阿莫西林原料對制劑溶出行為的影響[J]. 中外醫(yī)學(xué)究， 2013， 11 （13）： 136-137.

[5]LIU G Q. The pharmacological mechanism and clinical application of amoxicillin in Western Medicine[J]. World Latest Medicne Information （Electronic Version）， 2016， 65 （16）： 122-124.

[6]再曼， 王虎. 西藥阿莫西林的藥理作用機制以及臨床應(yīng)用[J]. 世界最新醫(yī)學(xué)信息文摘： 連續(xù)型電子期刊， 2015 （81）： 112-113.

[7]古麗努爾·蘇力坦汗. 西藥阿莫西林的藥理學(xué)機制及臨床應(yīng)用[J]. 廣東化工， 2017， 44 （9）： 178-178.

[8]李余， 裴詩恩， 喬駿馳， 等. 抗埃博拉病毒藥物法匹拉韋及其衍生物pKb值預(yù)測[J]. 當(dāng)代化工， 2016， 45 （7）： 1503-1504.

[9]陳虹， 丁俊杰， 丁曉琴， 等. 從頭計算研究水溶液中有機化合物pKa值[J]. 計算機與應(yīng)用化學(xué)， 2007， 24 （5）： 580-582.

[10]焉炳飛， 張明霞， 方圣濤， 等. 麥角甾醇及其衍生物1H-和13C-NMR光譜的密度泛函理論研究[J]. 煙臺大學(xué)學(xué)報， 2016， 29 （2）： 83-88.

基金項目： 2019 年度臺州市大學(xué)生科技創(chuàng)新項目資助（項目編號：2019121）。

收稿日期： 2019-09-28

作者簡介： 洪文燕（1997-），女，浙江省三門縣人，研究方向：藥物合成。E-mail：997202351@qq.com。

通訊作者： 鐘愛國（1964-），男，教授，碩士，研究方向：計算化學(xué)。E-mail：xg2268@163.com。