文章快速检索     高级检索
  山东大学学报(工学版)  2016, Vol. 46 Issue (3): 99-105  DOI: 10.6040/j.issn.1672-3961.0.2016.115
0

引用本文 

白艳艳, 曹月欣, 张宪玺, 张鲁格, 孙德志, 延辉, 张翀. SDS对香豆素在SB-16胶束增溶行为影响[J]. 山东大学学报(工学版), 2016, 46(3): 99-105. DOI: 10.6040/j.issn.1672-3961.0.2016.115.
BAI Yanyan, CAO Yuexin, ZHANG Xianxi, ZHANG Luge, SUN Dezhi, YAN Hui, ZHANG Chong. The influence of SDS on C-343 solubilizing behavior in SB-16 micelle[J]. Journal of Shandong University(Engineering Science), 2016, 46(3): 99-105. DOI: 10.6040/j.issn.1672-3961.0.2016.115.

作者简介

白艳艳(1990— ),女,山东莘县人,硕士研究生,主要研究方向为多孔碳基材料亲水性质提升与气体分离的理论研究.

通讯作者

张翀(1970— ),男,山东济宁人,副教授,博士,主要研究方向为多孔碳基材料亲水性质提升与气体分离的理论研究. E-mail:zhangchong@lcu.edu.cn

文章历史

收稿日期:2016-04-07
网络出版时间:2016-05-09
SDS对香豆素在SB-16胶束增溶行为影响
白艳艳1, 曹月欣2, 张宪玺1, 张鲁格1, 孙德志1, 延辉1, 张翀1     
1. 聊城大学化学与化工学院,山东 聊城 252059;
2. 山东省医学高等专科学校药学系,山东 济南 250002
摘要: 为了模拟药物香豆素在十二烷基硫酸钠(SDS)和十六烷基磺基甜菜碱(SB-16)复配体系中的分配规律,采用分子动力学的方法模拟了香豆素分子(C-343)在SB-16水溶液中的增溶行为,在上述溶液的基础上又继续添加了表面活性剂SDS,并研究了SDS对香豆素在SB-16胶束增溶行为影响。研究结果表明,在298K,101.3kPa下的平衡条件下,C-343可以自发地增溶到SB-16球形胶束中。SDS分子的烷基尾链部分插入到SB-16球形胶束内部,并在SB-16表面形成一层疏水基团在内,亲水基团在外的SDS分子薄层。C-343分子在SDS分子的疏水作用和静电作用下,从SB-16胶束内部逐步迁移至SDS胶束的外层表面上。从理论证实了添加表面活性剂可以对药物分子在胶束溶液中的增溶与释放过程进行有效控制,以期为该类体系的试验研究提供理论依据。
关键词: 分子动力学模拟    十六烷基磺基甜菜碱    十二烷基硫酸钠    香豆素分子    增溶    药物释放    
The influence of SDS on C-343 solubilizing behavior in SB-16 micelle
BAI Yanyan1, CAO Yuexin2, ZHANG Xianxi1, ZHANG Luge1, SUN Dezhi1, YAN Hui1, ZHANG Chong1     
1. School of Chemistry and Chemical Engineering, Liaocheng University, Liaocheng 252059, Shandong, China ;
2. Pharmacy Department, Shandong Medical College, Jinan 250002, Shandong, China
Abstract: In order to explore the dissolutive law of coumarin (C-343) in the SDS(sodium dodecyl sulfate)/SB-16(N-hexadecyl-N,N-dimethylammonio-1-propanesulfonate) synergistic system, the solvation of C-343 in pure SB-16 micelle solution and the additive of SDS surfactant influencing this solvation process were performed based on molecular dynamics method. The results showed that the C-343 could be spontaneously dissolved in the pure SB-16 micelle in equilibrium at 298 K and 101.3kPa. After addition of SDS, the SB-16 micelle appeared slightly swollen and its outer surface was inserted by the hydrophobic tails of SDS, forming a thinner SDS outshell covering the SB-16 micelle. Simultaneously, migration of C-343 from the inner of SB-16 micelle back to the outsurface of SDS shell took place due to the hydrohobic and electrostatic interactions. This work provided one convincing evidence that adding certain surfactants in some kinds of micelles could effectively control the solubilization or release behaviors of some drugs in these solutions, and thus give one helpful guidance for those corresponding experimental studies accordingly.
Key words: molecular dynamics simulation    N-hexadecyl-N,N-dimethylammonio-1-propanesulfonate (SB-16)    sodium dodecyl sulfate (SDS)    coumarin (C-343)    solubilization    drug release    
0 引言

香豆素主要来自于一些高等植物[1-4],具有多样化的结构[5]以及抗炎[6]、抗菌、抗癌、抗抑郁[7]、抗氧化[8]等丰富特性[5-10]。目前,香豆素药物的提取技术已较成熟,而其在人体内的运输机制尚不清楚。在香豆素的使用过程中,载体的选择至关重要,磺基甜菜碱型离子型表面活性剂就是其中候选载体之一[11-12]。试验证明,当阴离子表面活性剂、两性表面活性剂等混合形成复配体系时,往往呈现出单一体系所不具备的优点,如低刺激性、耐硬水性以及生物降解性等[13-17],因而受到人们越来越多地关注。

比如,Chiranjib Banerjee等人通过试验研究了阴离子表面活性剂十二烷基硫酸钠(sodium dodecyl sulfate,SDS)和两性表面活性剂N,N-二甲基-N-十六烷基丙磺酸(N-hexadecyl-N,N-dimethylammonio-1-propanesulfonate,SB-16)的复配体系与药物分子姜黄色素之间的作用机理和光学性质[18],当在溶有姜黄色素的SB-16球状胶束溶液中加入SDS之后,姜黄色素发生了激发态分子内质子转移。同时,他们还将阴离子表面活性剂SDS加入到嵌有香豆素分子(Coumarin-343,C-343)的SB-16两性表面活性剂胶束溶液中,发现SDS可使香豆素分子从SB-16胶束中释放出来[19]

然而,Chiranjib Banerjee等人所做的工作,仅仅描述了SDS(阴离子表面活性剂)和香豆素分子在SB-16(两性表面活性剂)胶束复配体系中所发生的宏观现象,而微观增溶机理则没有涉及。为此,本研究采用分子动力学方法,模拟分析了香豆素分子C-343在SB-16胶束溶液中的增溶行为,以及SDS的加入对这种增溶行为的影响,还解释了香豆素分子在SDS/SB-16复配体系和纯SB-16体系内发生不同增溶现象的微观原因,以期为某些药物分子在表面活性剂复配体系中运输机理的研究提供理论依据。

1 模拟体系及方法 1.1 模拟体系

图 1为模拟研究的两性表面活性剂(SB-16)、香豆素(C-343)以及阴离子表面活性剂(SDS) 3种单体分子的构型示意图。三种单体分子之间所构成的两种复配体系:即(Ⅰ) C-343/SB-16体系和(Ⅱ) C-343/SB-16/SDS体系,两体系的对应参数列于表 1

图 1 表面活性剂分子及香豆素C-343的分子构型示意图 Figure 1 The molecular strcutures of sufactant molecules SB-16,SDS and coumarin-343
表 1 模拟体系Ⅰ和体系Ⅱ的组分构成及其相关参数 Table 1 The components and corresponding parameters of system Ⅰ and system Ⅱ

体系Ⅰ是C-343/SB-16复配体系,用于研究一个香豆素分子C-343在SB-16胶束中自发增溶的行为。在构建过程中,首先,将一个由70个SB-16表面活性剂单体组成的球形胶束(其半径r~2.10 nm)置于一12×12×12 nm3的立方盒子中心,并在盒内其余部分充满水,此时SB-16浓度为6.25×10-2 mol/L,远大于文献所报道的临界胶束浓度(3.0×10-5mol/L)[20]。然后,再在胶束周围的水环境中随机添加一个香豆素分子,这就构成了体系Ⅰ,见图 2(a)。为了便于观察,图 2(a)中的周期性格子以及其内部所填充的水分子均未显示。经过30 ns的NPT动力学模拟后,体系Ⅰ便达到热力学平衡,见图 2(d)

图 2 体系Ⅰ(C-343/SB-16) 和Ⅱ(C-343/SB-16/SDS)初始构型和平衡构型示意图 Figure 2 The associated initial and equilibrium configurations of system Ⅰ (C-343/SB-16) and Ⅱ (C-343/SB-16/SDS)

体系Ⅱ是C-343/SB-16/SDS复配体系,用于模拟添加SDS后C-343分子在SB-16胶束溶液中增溶机理的变化过程。在构建体系Ⅱ的过程中,以体系Ⅰ所得到的C-343/SB-16胶束平衡构型(半径r~2.10 nm)为核心,利用Packmol程序包再将74个SDS分子以烷烃尾链为近心端[21],以硫酸根头基为远心端(SDS分子长度约为1.60 nm)的直链方式随机分布在SB-16胶束表面外2.20~3.80 nm的壳层里,形成一个类球形的复配C-343/SB-16/SDS胶束体系,然后将之置于一12×12×12 nm3的立方盒子中心并将其余部分充满水。此时,SDS添加剂的浓度为6.61×10-2 mol/L。由于体系Ⅱ(C-343/SB-16/SDS)有3种组分,若在同一幅图中同时显示,则各组分之间将难以分辨。所以,为了更加清晰,本研究拟以两幅示意图来说明,即图 2(b)图 2(c)(二者分别显示了C-343/SB-16组分和SB-16/SDS组分的初始分布)。为了便于观察,两图中周期性盒子以及盒内水分子都没显示。

1.2 模拟方法

所有模拟都采用GROMACS 4.0.5 程序包与GROMOS 53A6联合原子力场[22]。在模拟过程中,采用可很好描述水分子介电性质和热力学性质的SPC/E模型,使用最陡下降法进行能量最小化,NPT系综的条件设定为101.3 kPa和298 K,温度和压力的控制方法为Berendsen和Isotropic方法[23],采用LINCS[24]算法约束各分子的键长,采用PME方法描述各原子间的静电相互作用力[25],Lennard-Jones势的范德华力截断半径设定为1.20 nm,模拟步长设定为2 fs,每隔2 ps收集一次动力学轨迹信息,使用VMD 1.8.5[26]程序观察动力学轨迹。

2 结果与讨论 2.1 SB-16/C-343体系

为了研究香豆素分子(C-343)在SB-16胶束溶液中的增溶行为,本研究模拟了体系Ⅰ (C-343/SB-16),并把模拟时间设定为30 ns。通过VMD程序发现,当浓度为6.25×10-2 mol/L的SB-16水溶液内随机放入一个C-343分子后(图 2(a)),在疏水作用及静电作用下,体系会在约17 ns内(图 3)达到平衡,完成C-343分子从水溶液体系中逐渐向胶束接近,直到完全嵌入胶束的整个过程(平衡构型见图 2(d))。

图 3 模拟体系Ⅰ(C-343/SB-16)中香豆素分子质心与SB-16胶束质心之间的距离随模拟时间的变化示意图 Figure 3 The sketch of the distance between the mass centers of coumarin and SB-16 micelle along with simulation time of system Ⅰ (C-343/SB-16)

为了进一步探讨C-343分子在平衡胶束中的具体分布,绘制了C-343质心及SB-16分子的末端C16以及头基N、S原子相对胶束质心的径向分布概率图,详见图 4。分析图 4可知,SB-16分子的C16相对胶束质心的距离分布较广,几乎在0 ~ 2.10 nm的范围内都有分布,这表明,SB-16的疏水烷基长链在平衡时实际上已蜷曲成一团,形成了胶束的疏水核心。与此同时,SB-16亲水基团(-NO'+(CH3)2-(CH2)3-SO'-O3)的N和S原子的分布则相当集中,分别位于距胶束质心约2.06和2.08 nm处,组成了胶束的最外层。而香豆素的最大径向分布值为1.82 nm。这说明,在体系Ⅰ (C-343/SB-16)达到平衡时,C-343的质心已经包裹在胶束内部的疏水栅栏层内了。有趣的是,由图 2(f)可以看出,C-343分子真正完全嵌入胶束内的是其疏水部分,而其亲水端羧基部分则处在整个球形胶束的浅层表面上,与胶束外的水分子有直接接触的可能。

图 4 香豆素分子的质心、SB-16分子的末端碳(C16)、以及头基中N,S原子相对于SB-16胶束质心的径向分布概率示意图 Figure 4 The radial density distribution of the mass center of coumarin,C16,N and S atoms of SB-16 molecules relative to the mass center of the whole micelle
2.2 SB-16/C-343/SDS体系

研究体系Ⅰ可知,C-343分子在平衡时可自发地增溶到SB-16胶束内的疏水区域内。而实践上,在利用表面活性剂自组装体系进行药物输运时,不仅要考虑对药物摄取的控制,而且还要考虑对药物释放的控制。因此本研究配置了体系Ⅱ(C-343/SB-16/SDS),以观察SDS阴离子表面活性剂对C-343在SB-16胶束体系中分布位置的影响,具体参数见表 1图 2(b)(c)分别显示了体系Ⅱ中C-343/SB-16组分和SB-16/SDS组分的初始构型,图 2(e)(f)则分别为对应的平衡构型。

2.2.1 C-343/SB-16组分

图 5为体系Ⅱ中C-343质心相对于SB-16胶束质心的距离随模拟时间t的变化图。分析图 5可知,体系在3 ns后就达到了平衡。达到平衡后,C-343分子的质心相对于整个胶束质心的最大径向分布概率在2.50 nm附近,这个距离显然比体系Ⅰ(C-343/SB-16)的对应距离1.82 nm外移了差不多近0.68 nm,甚至还超出了SB-16胶束外径(~2.25 nm,此处用SB-16最外端的头基S原子相对于胶束质心的径向分布概率定义) 0.25 nm。这说明,SDS的添加不仅使SB-16的球形胶束发生了膨胀,而且还使原来增溶于SB-16胶束内部的C-343分子迁移出来(见图 2(e))。

图 5 香豆素分子(C-343)的质心相对SB-16胶束质心的距离随模拟时间的变化 Figure 5 The curve of the distance between the mass center of coumarin and that of SB-16 micelle depending on the simulation time

分析图 6可知,S原子的径向分布范围为1.00~3.00 nm,比体系Ⅰ中的1.50 ~ 2.70 nm(见图 4)也增大了很多。这说明,SDS的加入对SB-16的头基分布产生了较大影响,使原来较为集中的头基分布发生了很大的分散。从图 2(e)也可以直观地看出这些特点。图 7为体系Ⅱ中SB-16胶束的末端C原子以及头基中的N、S原子和SDS表面活性剂的C12及其头基中的S原子相对胶束质心的径向分布密度图。分析图 7还知,SB-16的末端C16的径向分布范围为0~3.00 nm,相比体系Ⅰ的0~2.40 nm范围也变广了,这证明了在SDS分子的作用下,胶束SB-16的疏水核心不仅发散,而且蜷曲程度也增加了。

图 6 香豆素分子质心和SB-16胶束表面的S原子相对胶束质心的径向分布概率 Figure 6 The radial density distribution of the mass center of coumarin and S of SB-16 molecules relative to the mass center of micelle
图 7 SB-16胶束的末端C原子以及头基中的N、S原子和SDS表面活性剂的C12及其头基中的S原子相对胶束质心的径向分布密度 Figure 7 The radial density distribution of C16,N,S atoms of SB-16 molecules and C12 and S atoms of SDS molecules relative to the mass center of micelle
2.2.2 SB-16/SDS组分

由于香豆素分子在SB-16胶束中分布的变化,以及SB-16胶束本身的形态变化都是由SDS的加入引起的,所以对体系Ⅱ中SB-16/SDS组分关系的研究对理解整个胶束体系的分子动力学行为非常重要。

分析图 7可知,在体系Ⅱ中,SDS的末端C12最大概率分布约在2.30 nm,与SB-16的头基S原子的最大概率分布 (~2.30 nm)十分相近。这说明,SDS表面活性剂分子和SB-16胶束由于都具有很长的疏水链而发生直接接触。然而,这二者的径向分布概率的范围延展却都很大,前者为0.90~3.40 nm,后者为1.00~3.20 nm,此时SDS不仅紧紧贴在SB-16表面,而且尾端烃基有很大一部分伸入到SB-16胶束内部,形成一种相互扦插的复配胶束体系。此外,分析图 7还可知,SDS的头基S原子的最大概率分布约为2.70 nm,仅比其尾基C12的最大概率(~ 2.30 nm)大0.40 nm,但却小于SDS的分子链长度1.60 nm,这说明有的SDS分子在平衡时扦插于SB-16胶束内,而有的则缠绕在SB-16胶束表面上(见图 2(f)),二者通过疏水作用和静电作用包合在一起,共同为香豆素分子C-343从SB-16胶束疏水栅栏层的迁出创造了空间条件。

3 结论

本研究采用分子动力学的方法详细探讨了C-343在SB-16胶束中的增溶行为,以及SDS对这种增溶行为的影响。通过体系平衡时分子构型直观图的观察以及C-343分子与SB-16胶束质心之间的径向密度分布和回转半径的分析,发现在纯SB-16表面活性剂体系中,C-343分子可以自发地增溶到胶束的疏水栅栏层。而当向SB-16胶束中添加SDS阴离子表面活性剂时,C-343分子又可从SB-16胶束的疏水栅栏层迁移至胶束表面之外。本研究可为后续研究者对其他种类的表面活性剂复配体系作为药物分子载体的应用机理探讨提供理论依据。

参考文献
[1] 路静静, 曹家庆, 李巍, 等. 甘草废渣中一个新的香豆素类化合物[J]. 中草药,2015, 46 (2) : 174-177.
LU Jingjing, CAO Jiaqing, LI Wei, et al. A new coumarins from residue of Glycyrrhiza uralensis[J]. Chinese Traditional and Herbal Drugs,2015, 46 (2) : 174-177. (0)
[2] 周堃, 王月德, 董伟, 等. 当归中一个新香豆素类化合物及其细胞毒活性[J]. 中草药,2015, 46 (18) : 2680-2682.
ZHOU Kun, WANG Yuede, DONG Wei, et al. A new coumarin from roots of Angelica sinensis and its cytotoxicity[J]. Chinese Traditional and Herbal Drugs,2015, 46 (18) : 2680-2682. (0)
[3] 张韶瑜, 孟林, 高文远, 等. 香豆素类化合物生物活性研究进展[J]. 中国中药杂志,2005, 30 (6) : 410-414.
ZHANG Shaoyu, MENG Lin, GAO Wenyuan, et al. Advances on biological activities of coumarins[J]. China Journal of Chinese Materia Medica,2005, 30 (6) : 410-414. (0)
[4] 齐晓丽, 田克立, 张瑞典, 等. 冬凌草甲素诱导胰腺癌PANC-1细胞凋亡及对G2/M细胞周期的阻滞作用[J]. 山东大学学报(医学版),2011, 49 (2) : 9-13.
QI Xiaoli, TIAN Keli, ZHANG Ruidian, et al. Induction of apoptosis and G2/M cell cycle arrest by oridonin on human pancreatic cancer PANC-1 cells[J]. Journal of Shandong University (Health Science),2011, 49 (2) : 9-13. (0)
[5] GARDELLY M, TRIMECH B, BELKACEM M A, et al. Synthesis of novel diazaphosphinanes coumarin derivatives with promoted cytotoxic and anti-tyrosinase activities[J]. Bioorganic & Medicinal Chemistry Letters,2016, 26 (10) : 2450-2454. (0)
[6] 赵春苗, 李亮亮. 白芷总香豆素对疮疡模型的影响[J]. 中药药理与临床,2014, 30 (1) : 61-64.
ZHAO Chunmiao, LI Liangliang. Effect of total coumarin of Angelicae dahuricae on sores model[J]. Pharmacology and Clinics of Chinese Materia Medica,2014, 30 (1) : 61-64. (0)
[7] ESCUDERO G E, LAINO C H, ECHEVERRIA G A, et al. Improving the antidepressant action and the bioavailability of sertraline by co-crystallization with coumarin 3-carboxylate. Structural determination[J]. Chemico-Biological Interactions,2016, 249 : 46-55. (0)
[8] 裴媛, 马贤德, 易杰, 等. 独活香豆素对帕金森病模型大鼠抗氧化功能及谷氨酸含量的影响[J]. 中国老年学杂志,2014, 34 (5) : 1272-1274.
PEI Yuan, MA Xiande, YI Jie, et al. Effect of total coumarins from Radix Angelicae on antioxidant function and the content of glutamate in rat model of Parkinson disease[J]. Chinese Journal of Gerontology,2014, 34 (5) : 1272-1274. (0)
[9] 梁新丽, 朱梦良, 赵国巍, 等. 白芷中呋喃香豆素类成分对葛根素、芍药苷、长春新碱肠道转运的影响[J]. 中草药,2015, 46 (7) : 1007-1011.
LIANG Xinli, ZHU Mengliang, ZHAO Guowei, et al. Effect of furan coumarins from Angelica Dahuricae Radix on intestinal transport absorption of puerarin, paeoniflorin, and vincristine[J]. Chinese Traditional and Herbal Drugs,2015, 46 (7) : 1007-1011. (0)
[10] 李娜, 孙捷, 孙敬勇, 等. 4-芳基香豆素类化合物的活性及合成方法研究进展[J]. 中国医药工业杂志,2015, 46 (12) : 1352-1365.
LI Na, SUN Jie, SUN Jingyong, et al. Progress in bioactivities and synthetic methods of 4-Arylcoumarins[J]. Chinese Journal of Pharmaceuticals,2015, 46 (12) : 1352-1365. (0)
[11] 施来顺, 万忠义. 新型甜菜碱型沥青乳化剂的合成与性能测试[J]. 山东大学学报(工学版),2008, 38 (4) : 112-126.
SHI Laishun, WAN Zhongyi. Synthesis and performance evaluation of a novel betaine-type asphalt emulsifier[J]. Journal of Shandong University (Engineering Science),2008, 38 (4) : 112-126. (0)
[12] 王丽娟, 石静, 赵方剑. 气液界面上磺基甜菜碱两性表面活性剂分子动力学模拟[J]. 山东大学学报(工学版),2014, 44 (6) : 83-89.
WANG Lijuan, SHI Jing, ZHAO Fangjian. Molecular dynamics simulation of zwitterionic surfactant sulfobetaine at the vapor/liquid interface[J]. Journal of Shandong University (Engineering Science),2014, 44 (6) : 83-89. (0)
[13] 张群, 裴梅山, 张瑾, 等. 十二烷基硫酸钠与两性表面活性剂复配体系表面性能及影响因素[J]. 日用化学工业,2006, 36 (2) : 69-72.
ZHANG Qun, PEI Meishan, ZHANG Jin, et al. Study of performance and its influencing factors of blends of sodium dodecyl sulfate and zwitterionic surfactants[J]. China Surfactant Detergent & Cosmetics,2006, 36 (2) : 69-72. (0)
[14] 周雅文, 徐宝财, 韩富, 等. 两性/阴离子表面活性剂复配体系研究[J]. 北京工商大学学报(自然科学版),2008, 26 (5) : 5-8.
ZHOU Yawen, XU Baocai, HAN Fu, et al. Study on mixed system of zwitterionic and anionic surfactants[J]. Journal of Beijing Technology and Business Uuniversity (Natural Science Edition),2008, 26 (5) : 5-8. (0)
[15] 张志庆, 徐桂英, 肖莉, 等. 十二烷基甜菜碱/十二烷基硫酸钠复配体系与原油间的界面张力[J]. 山东大学学报(理学版),2002, 37 (2) : 153-166.
ZHANG Zhiqing, XU Guiying, XIAO Li, et al. Study on interfacial tension of the mixed system of dodecyl beyaine and sodium dodecyl sulfate/crude oil[J]. Journal of Shandong University (Natural Science Edition),2002, 37 (2) : 153-166. (0)
[16] 苗领, 李玉江, 李枫. 聚合物/表面活性剂驱模拟采出水的乳化稳定性[J]. 山东大学学报(工学版),2015, 45 (2) : 89-94.
MIAO Ling, LI Yujiang, LI Feng. Emulsion stability of simulated polymer/surfactant flooding produced water[J]. Journal of Shandong University (Engineering Science),2015, 45 (2) : 89-94. (0)
[17] 马宝东, 高宝玉, 卢磊, 等. 聚合物/表面活性剂二元体系油水乳化过程[J]. 山东大学学报(工学版),2010, 40 (4) : 117-120.
MA Baodong, GAO Baoyu, LU Lei, et al. A characterization method of oil-water emulsion dynamic characteristics in surfactant/polymer flooding[J]. Journal of Shandong University (Engineering Science),2010, 40 (4) : 117-120. (0)
[18] BANERJEE C, GHOSH S, MANDAL S, et al. Exploring the photophysics of curcumin in zwitterionic micellar system: an approach to control ESIPT process in the presence of room temperature ionic liquids (RTILs) and anionic surfactant[J]. Journal of Physical Chemistry B,2014, 118 (13) : 3669-3681. (0)
[19] BANERJEE C, MANDAL S, GHOSH S, et al. Tuning the probe location on zwitterionic micellar system with variation of pH and addition of surfactants with different alkyl chains: solvent and rotational relaxation studies[J]. Journal of Physical Chemistry B,2012, 116 (36) : 11313-11322. (0)
[20] SAPELLI E, BRANDAO T A S, FIEDLER H D, et al. Fluorescence of Zn (Ⅱ) 8-hydroxyquinoline complex in the presence of aqueous micellar media: The special cetyltrimethylammonium bromide effect[J]. Journal of Colloid and Interface Science,2007, 314 (1) : 214-222. (0)
[21] MARTINEZ L, ANDRADE R, BIRGIN E G, et al. Packmol: a package for building initial configurations for molecular dynamics simulations[J]. Journal of Computational Chemistry,2009, 30 (13) : 2157-2164. (0)
[22] OOSTENBRINK C, VILLIA A, MARK A E, et al. A biomolecular force field based on the free enthalpy of hydration and solvation: the GROMOS force-field parameter sets 53A5 and 53A6[J]. Journal of Computational Chemistry,2004, 25 (13) : 1656-1676. (0)
[23] BERENDSEN H J C, POSTMA J P M, GUNSTEREN W F, et al. Molecular dynamics with coupling to an external bath[J]. Journal of Chemical Physics,1984, 81 (8) : 3684-3690. (0)
[24] HESS B, BEKKER H, BERENDSEN H J C, et al. LINCS: a linear constraint solver for molecular simulations[J]. Journal of Computational Chemistry,1997, 18 (12) : 1463-1472. (0)
[25] ESSMANN U, PERERA L, BERKOWITZ M L, et al. A smooth particle mesh Ewald method[J]. Journal of Chemical Physics,1995, 103 (19) : 8577-8593. (0)
[26] HUMPHREY W, DALKE A, SCHULTEN K. VMD-visual molecular dynamics[J]. Journal of Molecular Graphics,1996, 14 : 33-38. (0)