钩吻镇痛活性成分与吗啡的计算研究

Journal of Advances in Physical Chemistry 物理化学进展, 2016, 5(1), 9-17 Published Online February 2016 in Hans.

Computational Study on the Active Analgesic Components of Gelsemium and Morphine

Xinxin Lv, Jiazhen Xu, Zhengxue Li, Binjun Chen, Guang’an Lin, Ruixia Wang, Chaojie Wang*

School of Pharmaceutical Science, Wenzhou Medical University, Wenzhou Zhejiang

th th th

Received: Jan. 28, 2016; accepted: Feb. 14, 2016; published: Feb. 18, 2016

Copyright 2016 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Abstract

The geometrical and electronic structures, spectral properties of 16 active analgesia components of Gelsemium and Morphine were calculated at the X3LYP/6-31++G(d,p) level and Concept-DFT was also employed. The computed results indicate that: Structure of the active analgesia compo-nents of Gelsemium is similar to the partial structure of morphine. The infrared absorption peak of active analgesia components is mainly distributed in the three regions in wave number of 1700 - 1832 cm−1, 3500 - 3600 cm−1 and 3730 - 3860 cm−1. To characterize the solvent effect, the polarity of all compounds is larger than that of gaseous states. The values of electronegativeity, hardness, softness and electrophilicity indexes of gw7, gw8, gw12 and gw15 are in good agreement with Morphine. Among all of the compounds, the pharmacokinetic evaluation of gw7 is the optimal.

Keywords

Gelsemium, Morphine, Density Functional Theory, Pharmacokinetic

钩吻镇痛活性成分与吗啡的计算研究

吕鑫鑫，徐佳禛，李郑雪，陈炳君，林广安，王瑞霞，王朝杰*

温州医科大学药学院，浙江 温州

收稿日期：2016年1月28日；录用日期：2016年2月14日；发布日期：2016年2月18日

*

通讯作者。

文章引用: 吕鑫鑫, 徐佳禛, 李郑雪, 陈炳君, 林广安, 王瑞霞, 王朝杰. 钩吻镇痛活性成分与吗啡的计算研究[J]. 物

吕鑫鑫 等

摘 要

在X3LYP/6-31++G(d,p)水平，对钩吻16种镇痛活性成分和吗啡的几何结构、谱学性质和电子结构进行了详细的计算研究，并根据概念密度泛函理论对所有化合物的化学性质进行考察。研究结果表明：钩吻镇痛活性成分的几何结构与吗啡的部分结构类似；钩吻镇痛活性成分的红外光谱吸收峰主要分布于波数1700~1832 cm−1，3500~3600 cm−1和3730~3860 cm−1的三个区域；乙醇溶剂中钩吻和吗啡的极性增加，电负性、硬度、软度和亲电指数显示gw7、gw12和gw15与吗啡的性质最为相似。药代动力学数据表明gw7的药代动力学评价最佳。

关键词

钩吻，吗啡，密度泛函理论，药代动力学

1. 引言

如今，全球每年新增加癌症患者约1000万人，在我国，癌症患者人数也逐年高速增长。因此，有效规范的镇痛治疗是提高患者生存质量的重要措施。癌症疼痛治疗主要使用效果显著的强阿片类药物吗啡等，但该类药易引起一系列的副作用，比如：便秘、恶心、尿潴留等[1]。更多的是长期使用吗啡可产生耐受性而使得用药剂量增加，可产生痛觉过敏以及成瘾性，从而限制了其在临床上的应用。目前临床上陷入了对癌症患者“用药荒”的困境，因此急需开发出新型的不易产生耐受性的镇痛药物代替吗啡。

钩吻又名断肠草，味辛、苦、性温、有剧毒，故多以外用为主[2]。根据目前对钩吻的药用研究，发现其药理作用主要有抗肿瘤、镇痛、恢复造血功能、增强免疫、降压等；临床调查及基础研究表明，钩吻生物总碱具有显著的镇痛作用，调节疼痛阈值增大的能力仅次于吗啡，明显高于其他镇痛药物，且无成瘾性、耐受性和累积性毒性等副作用，但其毒性较大，治疗剂量与中毒剂量相近[3]-[6]。黄志毅首次报道钩吻素子具有显著的治疗外周神经损伤所致神经病理性疼痛的作用，其镇痛作用的有效剂量远离LD 50，安全性强，作用持续时间长[7]。相关研究人员采用大鼠骨癌疼痛模型，深入研究了钩吻素甲对癌症疼痛的镇痛作用及作用机制，证实了“钩吻可以有效减轻癌症等病患所带来的疼痛”，这为癌症患者带来了新的希望。

本研究将钩吻的镇痛活性成分的几何结构与吗啡的几何结构进行量子化学分析，比较两者之间的异同性[8]-[12]，具体化合物分子式见表1。利用钩吻(断肠草) 的有效镇痛的成分，研究其分子结构与药代动力学参数，为研制出不易产生耐药性的镇痛药提供理论基础。

2. 计算方法

本文采用X3LYP 密度泛函理论(DFT)在6-31++G(d,p)基组水平上，对已报道的16种钩吻镇痛活性成分和吗啡的几何结构，电子结构性质和振动光谱进行了计算研究，并在相同计算水平上采用极化连续模型(PCM)计算了16种钩吻镇痛活性成分和吗啡在乙醇环境中的性质变化，所有计算工作用G09程序包完成[13]。运用ACD/Percepta软件平台统一模拟计算其脂水分布系数(LogP)，药物与血浆蛋白结合能力(PPB)和血脑屏障通过性(CNS)等基础数据，由此预测药代动力学性质。

3. 结果与讨论

3.1. 几何结构

如图1所示吗啡为五环结构，其中含有部分氢化的菲环。吗啡分子中存在C5、C6、C9、C13和C14

吕鑫鑫 等

Table 1. Active analgesia components of Gelsemium and Morphine 表1. 钩吻镇痛活性成分和吗啡的分子式

编号 gw1 gw2 gw3 gw4 gw5 gw6 gw7 gw8 gw9 gw10 gw11 gw12 gw13 gw14 gw15 gw16 morphine

化学名

Gelsemine (钩吻素甲) Gelsemine N-oxide

Gelsevirine

19-(R)-Hydroxydihydro gelsemine

Koumine (钩吻素子) (19R)-Kouminol Koumine N-oxide Humantenine Gelsenicine (钩吻素己) Gelselegine (钩吻精碱) Humantenmine (葫蔓藤碱甲)

Gelsemoxonine 16-epi-voacarpine Gelseamydine 19-(Z)-Taberpsychine Gelsemamide (钩吻内酰胺)

Morphine (吗啡)

分子式 C 20H 22N 2O 2 C 20H 22N 2O 3 C 21H 24N 2O 3 C 20H 24N 2O 3 C 20H 22N 2O C 20H 24N 2O 2 C 20H 22N 2O 2 C 21H 26N 2O 3 C 19H 22N 2O 3 C 20H 26N 2O 4 C 19H 22N 2O 3 C 19H 22N 2O 5 C 21H 24N 2O 4 C 29H 36N 2O 6 C 20H 24N 2O C 20H 24N 2O 4 C 17H 19NO 3

分子量 322 338 352 340 306 324 322 344 326 358 326 358 368 508 308 356 285

可旋转键

1 1 2 1 1 1 1 3 2 3 2 3 3 5 0 1 0

共5个手性碳，左旋体吗啡作用于µ受体产生药理作用，而右旋体吗啡则无生理活性。C5、C6和C14上的氢与乙胺链成顺式，A 环与B 、C 、E 环均成顺式，D 环(乙胺链) 突出于其他平面之上。吗啡环A 和环C 上的酚羟基的酰化和醚化，6位羟基的氧化和氮上取代基的改变以及引入新的基团都会影响吗啡的活性。与吗啡的几何结构相比，钩吻的镇痛活性成分结构中也具有一平坦的芳环结构，苯环可以通过范德华力与作用位点结合产生药理效应；除gw3、gw5、gw6、gw7、gw14和gw15这6种活性成分，其余钩吻镇痛活性成分均有一碱性中心，能在生理pH 条件下大部分电离为阳离子，与受体的阴离子部位以静电结合[14]。吗啡中N 原子被不同取代基取代，可从激动剂转变为拮抗剂，如吗啡17位的氮原子上的甲基换成烯丙基得到纳洛酮，是吗啡受体的纯拮抗剂，但在钩吻镇痛活性成分结构中没有类似构型。

3.2. 谱学性质

优化得到的稳定钩吻镇痛活性成分和吗啡的IR 光谱数据如图2和表2所示。主要存在3个较强的红外活性特征振动：C=O、N-H 和O-H 键的伸缩振动。除gw5、gw6、gw7、gw15和吗啡外，其他钩吻镇痛活性成分的红外光谱数据显示在1700~1832cm-1范围内均存在的一个尖峰，对应C=O伸缩振动，与文献[15]中的数据吻合。在3500~3600 cm−1范围出现较弱吸收带，对应N-H 伸缩振动。gw4、gw6、gw10、gw12、gw13和gw14都在3730~3860 cm−1范围内出现醇羟基的伸缩振动，ν(N-H)较ν(O-H)谱带弱、尖。吗啡分子中存在一个酚羟基和一个六元环上的醇羟基，吗啡六元环C=C伸缩振动为1664 cm−1，环中-CH=CH-的反式CH 面外弯曲振动为937 cm−1，与文献[16]中的C=C和-CH=CH-振动伸缩频率红移了约20 cm−1。

吕鑫鑫 等

Figure 1. Structure of active analgesia components of Gelsemium and morphine 图1. 钩吻镇痛活性成分和吗啡的结构

(a) (b)

吕鑫鑫 等

(c) (d)

(e) (f)

(a), (b), (c)

Infared spectrum of gw1-gw16 and morphine; (d), (e), (f)UV-Vis spectra of gw1-gw16 and morphine

Figure 2. Infared spectrum and UV-Vis spectrum of active analgesia components of Gelsemium and Morphine 图2. 钩吻镇痛活性成分与吗啡的红外谱图和紫外–可见光谱

采用密度泛函理论中的X3LYP/6-31++(d,p)方法计算了钩吻16种镇痛活性成分和吗啡的UV-Vis 的光谱性质(如图2所示) 。钩吻镇痛活性成分中gw1、gw3、gw4、gw5和gw9在220~240 nm和260~300 nm范围内存在两个吸收峰，构型中存在和n →σ*和n →π*跃迁。图2(f)显示gw16的C=O吸收峰与吗啡相比红移了15 nm。gw2和gw13在300 nm以上存在吸收峰，说明分子中存在较大的共轭体系[17]。考虑溶剂化效应后，吗啡在233 nm处的吸收峰发生蓝移，在224 nm处出现吸收峰；对于钩吻镇痛活性成分中C=O的吸收峰在乙醇环境中也出现蓝移。

3.3. 电子结构分析

钩吻镇痛活性成分和吗啡的偶极矩如图3所示，所有化合物的偶极矩均大于1，其中结构gw6的偶极矩最大，约为16.07 × 10−30 C ∙m 。gw16的偶极矩仅次与gw6，gw2的偶极矩最小，极性最小。从图3可以看出gw3、gw6、gw11和gw16的偶极矩大于吗啡的偶极矩，表示这四种镇痛活性成分的极性较大。采用极化连续模型(PCM)考虑溶剂效应后，钩吻镇痛活性成分和吗啡的偶极矩均大于气相中的偶极矩，

因此在乙醇环境中，钩吻镇痛活性成分和吗啡的极性均增大。

吕鑫鑫 等

Table 2. The stretching vibrational frequencies (ν(C=O), ν(N-H) and ν(O-H)) of active analgesia components of Gelsemium and Morphine

表2. 钩吻镇痛活性成分与吗啡C=O、N-H 和O-H 键的伸缩振动频率

编号 gw1 gw2 gw3 gw4 gw5 gw6 gw7 gw8 gw9 gw10 gw11 gw12 gw13 gw14 gw15 gw16 morphine

ν(C=O) 1710.8 1712.8 1729.1 1700.6 - - - 1780.0 1789.1 1795.6 1795.5 1793.1 1799.0 1793.2 1785.0 1831.1 - 1747.0 -

ν(N-H) 3560.6 3558.8 - 3561.4 - - - - - 3500.6 - 3561.6 3683.4 - 3677.4 3476.7 -

ν(O-H) - - - 3841.8 - 3735.9 - - - 3852.7 - 3757.8 3777.5 3843.8 3843.2 - - 3806.5 3836.4

Figure 3. Dipole moment (μ) of active analgesia components of Gelsemium and Morphine

图3. 钩吻镇痛活性成分和吗啡结构的偶极矩(μ)

3.4. 概念密度泛函分析

采用概念密度泛函方法按照下式计算了钩吻16种镇痛活性成分和吗啡的电负性(χ) 、硬度(η) 、软度(S ) 和亲电指数(ω) [18] [19]：

吕鑫鑫 等

+εLUMO εI +A

=−HOMO 22η=I −A =εLUMO −εHOMO

χ=

S =

1

η

χ2

ω=

2η

具体数值见表3。应用概念密度泛函理论定义的电负性和硬度，从而在理论上探讨化合物的化学性质。gw7、gw12和吗啡的电负性较大，电负性越大，得电子能力越强，与受体结合作用越强。由η值可知，吗啡和gw15的硬度大于钩吻其余镇痛活性成分，其分子中可旋转键最少，而gw14的硬度最小。gw15的亲电指数最小，亲电性最弱，则其对应的酸性则最强。因此，从电负性、硬度、软度和亲电指数来看gw7、gw12和gw15与吗啡的性质最为相似。

3.5. 药代动力学分析

钩吻的镇痛活性成分和吗啡的脂水分布系数评价都较好。吗啡与不同脑区的阿片受体结合，模拟内源性阿片肽的作用，抑制P-糖蛋白的释放，从而干扰痛觉冲动传入中枢发挥镇痛作用，从表4中可以看出吗啡进入中枢神经系统的比例是35%，而gw7进入中枢神经系统的比例高达97%；gw7更容易作用于

Table 3. The values of electronegativity (χ), hardness (η), softness (S) and electrophilicity (ω) of active analgesia compo-nents of Gelsemium and Morphine calculated at the X3LYP/6-31++G(d,p) level

表3. 在X3LYP/6-31++G(d,p)水平计算得到钩吻镇痛活性成分和吗啡的电负性(χ) 、硬度(η) 、软度(S)和亲电指数(ω)

钩吻 gw1 gw2 gw3 gw4 gw5 gw6 gw7 gw8 gw9 gw10 gw11 gw12 gw13 gw14 gw15 gw16 morphine

Gas 0.12 0.13 0.13 0.12 0.13 0.11 0.14 0.12 0.11 0.13 0.13 0.14 0.11 0.13 0.11 0.11 0.14

χ/eV

Ethanol 0.12 0.14 0.13 0.12 0.13 0.13 0.15 0.12 0.13 0.13 0.13 0.14 0.11 0.13 0.11 0.13 0.11

Gas 0.19 0.19 0.17 0.18 0.17 0.16 0.19 0.19 0.17 0.19 0.19 0.20 0.15 0.20 0.18 0.18 0.10

η/eV

Ethanol 0.17 0.19 0.16 0.17 0.16 0.16 0.19 0.18 0.20 0.20 0.20 0.20 0.17 0.20 0.18 0.22 0.20

Gas 5.31 5.14 5.93 5.56 5.94 6.12 5.25 5.35 5.98 5.29 5.16 5.10 6.66 5.05 5.65 5.62 10.45

S/eV−1

Ethanol 5.72 5.19 6.14 5.75 6.25 6.13 5.22 5.44 5.11 5.03 5.04 4.98 5.78 5.04 5.43 4.60 4.94

Gas 0.04 0.04 0.05 0.04 0.05 0.04 0.05 0.04 0.04 0.04 0.04 0.05 0.04 0.04 0.03 0.04 0.10

ω/eV

Ethanol 0.04 0.05 0.05 0.04 0.05 0.05 0.06 0.04 0.04 0.04 0.04 0.05 0.03 0.04 0.04 0.04 0.03

吕鑫鑫 等

Table 4. The pharmacokinetic parameters of active analgesia components of Gelsemium and Morphine 表4. 钩吻镇痛活性成分和吗啡的药代动力学参数

gw1 gw2 gw3 gw4 gw5 gw6 gw7 gw8 gw9 gw10 gw11 gw12 gw13 gw14 gw15 gw16 morphine

CYP1A2 抑制剂 0.31 0.35 0.30 0.39 0.39 0.29 0.48 0.27 0.24 0.31 0.26 0.32 0.34 0.28 0.49 0.38 0.15

CYP2C9抑制剂 0.33 0.31 0.31 0.37 0.37 0.34 0.39 0.35 0.30 0.37 0.31 0.36 0.34 0.33 0.39 0.37 0.14

CYP2C19抑制剂 0.42 0.41 0.43 0.44 0.44 0.41 0.51 0.37 0.39 0.40 0.39 0.45 0.42 0.38 0.44 0.44 0.22

CYP2D6抑制剂 0.40 0.28 0.37 0.51 0.51 0.50 0.47 0.37 0.40 0.38 0.39 0.37 0.48 0.37 0.51 0.38 0.50

CYP3A4抑制剂 0.47 0.48 0.45 0.50 0.50 0.46 0.51 0.48 0.54 0.49 0.46 0.41 0.51 0.50 0.49 0.51 0.33

LogP 2.23 1.85 2.19 3.23 2.23 2.11 2.72 3.39 2.19 1.17 2.00 1.67 1.97 1.59 4.09 1.48 0.69

PPB −2.37 −2.65 −2.78 −2.31 −2.31 −3.03 −2.80 −2.56 −2.02 −3.49 −2.18 −2.28 −2.38 −3.61 −1.72 −2.72 −2.63

CNS 80% 90% 80% 92% 92% 91% 97% 87% 27% 51% 39% 18% 43% 58% 86% 82% 35%

中枢神经系统。药物和血浆蛋白的结合率是影响药物的分布和代谢的重要参数，是药物在体内动力学的一个重要过程，游离的血药浓度能够影响药物的排泄和代谢速率，更为重要的是它密切关系到药物的作用强度，计算结果显示，gw13和gw16有较强的血清蛋白结合率。

影响药物的代谢主要是指药物对药物代谢酶的影响[20]，根据报道90%以上的药物是通过肝微粒体的细胞色素P450酶完成[21]。人类参与药物代谢的主要细胞色素氧化酶CYP1A2、CYP2C9、CYP2C19、CYP2D6和CYP3A4在此激活过程中具有决定性的催化作用。计算结果显示，吗啡对CYP1A2、CYP2C9、CYP2C19酶没有其抑制作用，相比之下钩吻的镇痛活性成分中gw11的抑制作用较小，CYP2D6酶的抑制作用较强，在钩吻镇痛活性成分中gw4、gw5、gw6、gw7和gw15对CYP2D6酶抑制作用与吗啡相似。综上所述，gw7的药代动力学评价与其他镇痛活性成分相比最佳。

4. 结论

采用X3LYP 密度泛函理论(DFT)在6-31++G(d,p)基组水平上计算钩吻16种镇痛活性成分和吗啡的几何结构，同时采用极化连续模型(PCM)模拟乙醇溶液环境表明溶剂化效应对几何结构的影响。计算结果显示：1) 钩吻的镇痛活性成分分子中也具有一平坦的芳环结构，苯环可以通过范德华力与作用位点结合；钩吻镇痛活性成分中的N 原子被不同取代基取代，可从激动剂转变为拮抗剂。2) 红外谱图中钩吻的镇痛活性成分存在三个较强的吸收峰即：C=O、N-H 和O-H 键；gw16的C=O吸收峰与吗啡相比红移了15 nm。考虑溶剂化效应后，紫外-可见光谱中钩吻镇痛活性成分和吗啡的羰基吸收峰发生蓝移。3) 从电负性、硬度、软度和亲电指数来看gw7、gw12和gw15与吗啡的性质最为相似。4) 通过计算机模拟平台，预测钩吻镇痛活性成分的药代动力学性质发现gw7的药代动力学的评价与其他钩吻镇痛活性成分相比最佳。

吕鑫鑫 等

致 谢

感谢浙江省大学生科技创新活动计划(新苗人才计划) 项目(2014R413020)的资助。

参考文献 (References)

[1] 李桂英, 蔡卫梅. 饮食干预对阿片类药物治疗癌痛致便秘的影响[J]. 中国实用医药, 2013, 8(27): 1-3. [2] 全国中草药汇编编写组. 全国中草药汇编: 下册[M]. 北京: 人民卫生出版社, 1983: 118-121. [3] 王坤, 杨柯. 钩吻的药理研究及临床应用现状[J]. 广西中医学院学报, 2001, 4(4): 124-126. [4] 周名璐, 黄冲, 杨小平. 钩吻总碱的镇痛、镇静及安全性研究[J]. 中成药, 1998, 20(1): 35-36. [5] 谭建权, 邱成之, 郑林忠. 钩吻碱的镇痛作用和无依赖性[J]. 中药药理与临床, 1988, 4(1): 24-28. [6] 黄茜, 王春. 上海交大揭示古老中药“钩吻”镇痛奥秘[J]. 科技日报, 2013: 7.

[7] 黄志毅. 钩吻素子治疗神经病理疼痛的作用[D]: [硕士学位论文]. 福建: 福建医科大学, 2010: 6. [8] 吴达荣, 秦瑞, 蔡晶. 钩吻素子抗肿瘤作用研究[J]. 中药药理与临床, 2006, 22(5).

[9] Hua, W., Guo, T., Zhang, L., Wu, L.J. and Zhao, Q.C. (2007) Chemical Constituents of Gelsemium elegans Benth.

Chinese Journal of Medicinal Chemistry, 76, 108-110. [10] Hua, W., Zhao, Q.C., Yang, J., Shi, G.B., Wu, L.J. and Guo, T. (2008) Two New Benzofuran Lignan Glycosides from

Gelsemium elegans. Chinese Chemical Letters, 19, 1327-1329. [11] Kogure, N., Ishii, N., Kitajima, M., Wongseripipatana, S. and Takayama, H. (2006) Four Novel Gelsenicine-Related

Oxindole Alkaloids from the Leaves of Gelsemium elegans Benth. Organic Letters, 8, 3085-3088. [12] 刘霞, 胡应喜, 尹玉英. 吗啡的结构和旋光方向[J]. 北京石油化工学院学报, 1995, 3(2): 40-43.

[13] Frisch, M.J., Trucks, G.W., Schlegel, H.B., et al. (2010) Gaussian 09, Revision B.01. Gaussian Inc., Wallingford. [14] 郑虎, 翁玲玲. 药物化学: 第五版[M]. 北京: 人民卫生出版社, 2003.

[15] Rychly, J., Rychla, L., Stloukal, P. and Fiedlerova, A. (2013) UV Initiated Oxidation and Chemiluminescence from

Aromatic-Aliphatic Co-Polyesters and Polylactic Acid. Polymer Degradation and Stability, 98, 2556-2563. [16] 张敏敏, 金芩, 杭睿, 曾狄勤, 郑 武, 薛雅楠, 何碧如, 王朝杰. α-螺旋型寡聚L-乳酸的计算研究[J]. 物理化学

进展, 2015, 4(3): 103-110. [17] 李贺珍. 盐酸吗啡、咖啡因FTIR 检验与谱图解释分析[J]. 中国中医药咨询, 2012, 4(1): 400-402. [18] 孟令芝, 龚淑玲, 何永炳. 有机波谱分析: 第三版[M]. 武汉: 武汉大学出版社, 2010.

[19] Alparone, A. (2013) Response Electric Properties of α-Helix Polyglycines: A CAM-B3LYP DFT Investigation. Chem-ical Physics Letters, 563, 88-92. [20] 赵浩, 朱桂丽, 赵喆. 中西药物间相互作用的分析及研究[J]. 医药论坛杂志, 2006, 27(16): 116-117. [21] Park, B.K. (2000) Cytochrome P450 Enzymes in the Heart. Lancet , 355, 945-946.