吡嗪酰胺(PZA)是一种白色或类白色结晶粉末,是治疗肺结核的一线药物,也是短程化疗中最有效的抗菌药物,它与利福平、异烟肼能将结核病治疗疗程从9~12个月缩短到6个月。受水溶性低的影响,该药物的生物利用度较低,因此通过改善溶解性提高其生物利用度具有重要意义。
使用共溶剂是增大水溶性差药物溶解度的常用手段,广泛应用于制药过程和临床领域[1-2]。共溶剂现象是指在混合溶液中,当两种溶剂达到某一特定比例时,药物溶解度会达到最大值,此时的溶解度比在任一单溶剂中的溶解度都高[3-4]。共溶剂系统的原理是达到特定比例的两种溶剂与疏水性药物之间的界面张力会比任意单溶剂与疏水性药物之间的界面张力小,以此达到增大溶解度的目的[5]。
虽然吡嗪酰胺在纯溶剂中的溶解度数据已经有文献给出[6],但是其在混合溶剂中的溶解度数据尚未见报道。本文使用激光动态技术测定了吡嗪酰胺在甲醇+水和丙酮+正丙醇两种混合溶剂中的溶解度,并用热力学模型进行了关联比较,又进一步对该溶解过程的相关热力学参数进行了计算,为吡嗪酰胺药物的工业结晶分离、纯度提高提供了理论基础和数据支持。
1 实验部分 1.1 实验原料和仪器吡嗪酰胺,分析纯,北京金瑞林科技发展有限公司;甲醇、正丙醇、丙酮,分析纯,北京化工厂;实验用水为屈臣氏蒸馏水。
BT423S型电子分析天平,德国Sartorius公司;501A型超级恒温器,上海实验仪器厂;QHJ756B型磁力搅拌器,常州新析仪器厂。
1.2 吡嗪酰胺溶解度的测定使用激光动态法[7]测定吡嗪酰胺的溶解度。激光监测系统包括激光发生器、光电变压器和光强显示器;测定系统包括超级恒温器、磁力搅拌器和水银温度计。测定溶解度时,在结晶器中加入精准称量的溶剂,记录溶剂质量,结晶器敞口处用带有精密温度计的胶塞密封,记录温度,结晶器用橡胶管连接超级恒温器。打开激光检测仪器,观测并记录自动记录仪上显示的信号信息,此时溶剂内没有溶质,显示的接收信号为最强值。恒温下多次少量地加入PZA固体,由于颗粒的遮挡作用,光强会随着溶解过程不断增大,当接收的信号强度在30min内没有回升至最高点,即认为溶液饱和[8]。此时加入溶质的质量即为对应温度下PZA的溶解度。每组实验重复3次,取平均值。
吡嗪酰胺的饱和溶解度x1表示为
x1=m1M1m1M1+n∑i=2miMi | (1) |
混合溶剂中甲醇或丙酮的初始摩尔分数x20表示为
x02=m2M2m2M2+m3M3 | (2) |
式中,m1、m2、m3分别表示吡嗪酰胺、甲醇(丙酮)和水(正丙醇)的质量;M1、M2、M3分别表示吡嗪酰胺、甲醇(丙酮)和水(正丙醇)的摩尔质量[9-10]。
1.3 关联模型 1.3.1 Apelblat方程该方程是基于固液平衡关系以及溶质分子与溶液分子形成复合物的假设,由Clausius-Clapeyron方程[11]转变得到,表达式为
lnx1=A+BT/K+Cln(T/K) | (3) |
式中A、B、C均为无因次模型参数;T为系统的绝对温度。
1.3.2 λh方程1980年,Buchowski等[12]研究了溶解度、溶剂活度与温度之间的关系,得到的表达式为
ln(1+λ1−x1x1)=λh(1T−1Tm) | (4) |
式中λ、h为两个模型参数,由实验数据回归得到;T为系统的绝对温度;Tm为PZA的熔点,其值为463.02K[13]。
1.3.3 一般共存模型结合理想二元混合溶剂模型Redich-Kister (CNIBS/R-K)方程[14-15]可知,药物溶解度在恒温条件下与二元混合溶剂的组成有关,表现形式为
lnx1=x2ln(x1)2+x3ln(x1)3+x2x3n∑i=0Si(x2−x3)i | (5) |
式中x2、x3为混合溶剂的初始摩尔分数;(x1)2、(x1)3为PZA在单溶剂中的饱和溶解度;Si为模型参数,其中i的值为0~3。本文中混合溶剂为二元溶剂,此时可用(1-x2)替代x3,则有
lnx1=ln(x1)3+[ln(x1)2−ln(x1)3+S0−S1+S2]x2+(−S0+3S1−5S2)x22+(−2S1+8S2)x32+(−4S2)x42 | (6) |
方程(6)可简写为
lnx1=B0+B1x2+B2x22+B3x32+B4x42 | (7) |
式中B0、B1、B2、B3、B4均为模型参数。
1.3.4 Jouyban-Acree(J-A)模型J-A模型充分考虑了温度与混合溶剂组成对溶解度的影响[16-17],表达式为
lnx1=x2ln(x1)2+x3ln(x1)3+x2x3n∑i=0Ji(x2−x3)iT | (8) |
式中Ji为模型参数,i的值为0~3;x2、x3、(x1)2、(x1)3与上文一致。将其与Apelblat方程联立,得到新的混合模型[16]
lnx1=a2+b2T+c2lnT+(a1−a2)x2+(b1−b2+J0−J1+J2)x2T+(3J1−J0−5J2)x22T+(8J2−2J1)x32T+(−4J2)x42T+(c1−c2)x2lnT | (9) |
式中a1、b1、c1与a2、b2、c2分别为两种纯溶剂代入Apelblat方程的模型参数。方程(9)可简写为
lnx1=A1+A2T+A3lnT+A4x2+A5x2T+A6x22T+A7x32T+A8x42T+A9x2lnT | (10) |
式中A1~A9均为模型参数。
本文使用相对偏差DRD、平均偏差DMD、均方根偏差DRMSD来表示关联方程的拟合效果
DRD=x1−xcal1x1 | (11) |
DMD=1NN∑i=1|xi−xcalixi| | (12) |
DRMSD=√N∑i=1(xi−xcali)2N | (13) |
式中N为实验点数量,xi为PZA的溶解度,xical为PZA的溶解度计算值。
1.4 热力学参数的计算本文使用Van’t Hoff模型计算溶解过程的热力学参数。由该模型可知,在理想溶液中溶解过程满足
lnx1=−ΔHsolRT+ΔSsolR | (14) |
式中R为理想气体常数,8.314J/(mol·K);ΔHsol为标准摩尔溶解焓变;ΔSsol为标准摩尔溶解熵变。
可以计算出Gibbs自由能的表观值,表达式为
ΔGsol=ΔHsol−TmeanΔSsol | (15) |
式中Tmean为平均温度,计算式为
Tmean=n/n∑i=11Ti | (16) |
通过对比焓和熵对溶解过程中表观Gibbs自由能的贡献能力可以确定溶解过程的溶解原理,焓和熵的相对贡献计算方法分别为
ξH=|ΔHsol||ΔHsol|+|TmeanΔSsol|×100% | (17) |
ξS=|TmeanΔSsol||ΔHsol|+|TmeanΔSsol|×100% | (18) |
本文测量了吡嗪酰胺在水、甲醇、丙酮及正丙醇4种溶剂中的溶解度,并与文献[6]中的结果进行比较,对比结果如表 1所示,其中x1ref为文献值,x1exp为实验值。可能由于溶解平衡时间不一致以及系统误差导致两者结果有一定的偏差,但总体具有较好的一致性,证明本文实验数据真实可靠。
![]() |
下载CSV 表 1 PZA在4种溶剂中溶解度的文献值与实验值对比 Table 1 Comparison of literature values and experimental values of the solubility of PZA in four solvents |
吡嗪酰胺在两种混合溶剂中的实验值和关联模型的计算值以及关联结果的误差值均列于表 2、3,溶解度曲线如图 1、2所示。由以上结果可知,PZA的溶解度有很高的温度敏感性,会随着温度的上升而增大。所有数据的平均偏差均小于3%,均方根偏差均小于0.02%,可知修正的Apelblat方程、λh方程、一般共存模型和修正后的Jouyban-Acree方程关联溶解度数据均能较好地拟合,相对偏差的最大值分别为2.702%、9.039%、5.241%和17.058%,Apelblat方程拟合效果最好;4种模型的模型参数分别列于表 4~7。
![]() |
下载CSV 表 2 PZA在甲醇+水二元混合溶剂中的溶解度数据与关联结果 Table 2 Experimental mole fraction solubility and correlation results for PZA (x1) in methanol+water |
![]() |
下载CSV 表 3 PZA在丙酮+正丙醇二元混合溶剂中的溶解度数据与关联结果 Table 3 Experimental mole fraction solubility and correlation results for PZA (x1) in acetone+1-propanol |
![]() |
▲-0.0000; ●-0.1000; ■-0.2000; ▼-0.2997; ◆-0.3996; □-0.4999; ○-0.5999; △-0.7001; ▽-0.7998; ◇-0.8997; ☆-1.0000。 图 1 PZA在甲醇+水二元混合溶剂中的拟合Apelblat方程溶解度曲线 Fig.1 Fitting of the Apelblat solubility curve of PZA in methanol + water binary solvent |
![]() |
▲-0.0000; ●-0.1004; ■-0.2000; ▼-0.2999; ◆-0.3998; □-0.4997; ○-0.5999; △-0.7000; ▽-0.7999; ◇-0.9000; ☆-1.0000。 图 2 PZA在丙酮+正丙醇二元混合溶剂中的拟合Apelblat方程溶解度曲线 Fig.2 Fitting of the Apelblat solubility curve of PZA in acetone+1-propanol binary solvent |
![]() |
下载CSV 表 4 PZA在二元混合溶剂中的Apelblat方程溶解度数据关联参数 Table 4 Parameters of the Apelblat equation for PZA in the binary cosolvent mixtures |
![]() |
下载CSV 表 5 PZA在二元混合溶剂中的λh方程溶解度数据关联参数 Table 5 Parameters of the λh equation for PZA in the binary cosolvent mixtures |
![]() |
下载CSV 表 6 PZA在二元混合溶剂中的一般共存模型溶解度数据关联参数 Table 6 Parameters of the general cosolvency model for PZA in the binary cosolvent mixtures |
![]() |
下载CSV 表 7 PZA在二元混合溶剂中的J-A模型溶解度数据关联参数 Table 7 Parameters of the Jouyban-Acree model for PZA in the binary cosolvent mixtures |
PZA在两种混合溶剂中都会出现共溶剂现象,即当温度恒定时,在甲醇或丙酮的初始摩尔分数增大至0.6之前,PZA溶解度会随着甲醇或丙酮初始摩尔分数的增大而增大;之后,溶解度会随着甲醇或丙酮初始摩尔分数的增大而减小。出现此现象的原因为:一方面,溶质具有稳定的分子内氢键,在溶液中可形成稳定的缔合物;另一方面,混合溶剂之间会形成分子间氢键以确保它们的混溶性,两种溶剂可以互相干扰分子间氢键网络的形成,因此溶剂间分子间作用力减小,同时,纯溶剂的自缔合受到干扰,助溶剂降低了纯溶剂挤出非极性疏水化合物的能力,从而导致溶质的溶解度增加。但在文献中, 同样在甲醇+水的混合溶剂中4-氨基丁酸却没有出现共溶剂现象[18],所以我们认为共溶剂现象与溶质结构有关。
热力学参数的计算拟合结果结果如图 3、4,计算结果列于表 8。由热力学分析可知,焓变与表观Gibbs自由能的值均为正,说明该溶解过程是非自发进行的吸热过程;且在溶解过程中,对Gibbs自由能作出主要贡献的是焓变,ξH ≥ 60.82%。由图 5、6可知,斜率值是负时为熵变推动力,斜率值是正时为焓变推动力。在PZA的溶解过程中,当甲醇的初始摩尔分数在0~0.2000和0.5999~1.0000或丙酮初始摩尔分数在0.5999~1.0000范围内时,补偿图的斜率为负值,即该过程为熵变推动;甲醇的初始摩尔分数在0.2000~0.5999或丙酮初始摩尔分数在0~0.5999范围内时,补偿图的斜率为正值,即该过程为焓变推动。由图 3、4可知,PZA的溶解过程为熵变与焓变交替推动的过程。
![]() |
▲-0.0000; ●-0.1000; ■-0.2000; ▼-0.2997; ◆-0.3996; □-0.4999; ○-0.5999; △-0.7001; ▽-0.7998; ◇-0.8997; ☆-1.0000。 图 3 PZA在甲醇+水混合溶剂中溶解度(lnx1)与1/T关联数据的范特霍夫拟合直线 Fig.3 Van't Hoff plot of the mole fraction solubility (lnx1) of PZA in methanol + water binary mixed solvents against 1/T with a straight line to correlate the data |
![]() |
▲-0.0000; ●-0.1004; ■-0.2000; ▼-0.2999; ◆-0.3998; □-0.4997; ○-0.5999; △-0.7000; ▽-0.7999; ◇-0.9000; ☆-1.0000。 图 4 PZA在丙酮+正丙醇混合溶剂中溶解度(lnx1)与1/T关联数据的范特霍夫拟合直线 Fig.4 Van't Hoff plot of the mole fraction solubility (lnx1) of PZA in acetone+1-propanol binary mixed solvents against 1/T with a straight line to correlate the data |
![]() |
下载CSV 表 8 PZA恒温在二元混合溶剂中的溶解过程热力学数据 Table 8 Thermodynamic functions for the solution process of PZA in the binary solvent at constant temperature |
![]() |
图 5 302.60K下PZA在甲醇+水二元混合溶剂中的非线性焓-熵补偿图(点对应于甲醇的摩尔分数) Fig.5 ΔHsol vs. ΔGsol enthalpy-entropy compensation plot for the solubility of PZA in various methanol + water binary mixed solvents at a mean harmonic temperature of 302.60K (numbered points correspond to the mole fractions of methanol in the mixtures) |
![]() |
图 6 302.60K下PZA在丙酮+正丙醇二元混合溶剂中的非线性焓-熵补偿图(点对应于丙酮的摩尔分数) Fig.6 ΔHsol vs. ΔGsol enthalpy-entropy compensation plot for the solubility of PZA in various acetone + 1-propanol binary mixed solvents at a mean harmonic temperature of 302.60K(numbered points correspond to the mole fractions of acetone in the mixtures) |
(1) 修正的Apelblat方程、λh方程、一般共存模型和修正后的Jouyban-Acree方程均能很好地关联本文所测PZA溶解度的实验数据,其中Apelblat方程的关联效果最好。
(2) 吡嗪酰胺在两种混合溶剂中都出现了共溶剂现象,即当甲醇或丙酮的摩尔分数为0.6时,吡嗪酰胺的溶解度最大。出现共溶剂现象的原因应该与溶质分子的结构有关。
(3) 采用vant Hoff方程计算了溶解过程的热力学性质,包括Gibbs自由能、焓变和熵变,证明了该溶解过程为非自发吸热过程。通过非线性焓-熵补偿图可知,溶解过程机理与共溶剂混合物的组成相关。
[1] |
SONI K L, SOLANKI S S, MAHESHWARI R K. Solubilization of poorly water soluble drug using mixed solvency approach for aqueous injection[J]. British Journal of Pharmaceutical Research, 2014, 4(5): 549-568. DOI:10.9734/BJPR |
[2] |
AMIN K, DANNENFELSER R M, ZIELINSKI J, et al. Lyophilization of polyethylene glycol mixtures[J]. Journal of Pharmaceutical Sciences, 2004, 93(9): 2244-2249. DOI:10.1002/jps.20135 |
[3] |
YEH M K, CHANG L C, CHIOU A H J. Improving tenoxicam solubility and bioavailability by cosolvent system[J]. AAPS Pharm Sci Tech, 2009, 10(1): 166-171. DOI:10.1208/s12249-009-9189-2 |
[4] |
MIYAKO Y, KHALEF N, MATSUZAKI K, et al. Solubility enhancement of hydrophobic compounds by cosolvents:role of solute hydrophobicity on the solubilization effect[J]. International Journal of Pharmaceutics, 2010, 393: 48-54. DOI:10.1016/j.ijpharm.2010.03.059 |
[5] |
CHAUDHARI P, SHARMA P, BARHATE N, et al. Solubility enhancement of hydrophobic drugs using synergistically interacting cyclodextrins and cosolvent[J]. Current Science, 2007, 92(11): 1586-1591. |
[6] |
ZHANG K K, SHEN H, XU S J, et al. Thermodynamic study of solubility for pyrazinamide in ten solvents from T=(283.15 to 323.15) K[J]. The Journal of Chemical Thermodynamics, 2017, 112: 204-212. DOI:10.1016/j.jct.2017.04.014 |
[7] |
WANG S, WANG J K, YIN Q X. Measurement and correlation of solubility of 7-aminocephalosporanic acid in aqueous acetone mixtures[J]. Industrial and Engineering Chemistry Research, 2005, 44: 3783-3787. DOI:10.1021/ie048938e |
[8] |
WANG S, CHENG X W, LIU B L, et al. Temperature dependent solubility of sodium cyclamate in selected pure solvents and binary methanol + water mixed solvents[J]. Fluid Phase Equilibria, 2015, 390: 1-6. DOI:10.1016/j.fluid.2015.01.012 |
[9] |
ZHANG Q, YANG Y, CAO C C, et al. Thermodynamic models for determination of the solubility of dibenzothiophene in (methanol + acetonitrile) binary solvent mixtures[J]. The Journal of Chemical Thermodynamics, 2015, 80: 7-12. DOI:10.1016/j.jct.2014.08.012 |
[10] |
JI W Z, LU D Q, YANG S N, et al. Solubility and solution thermodynamics of cetilistat in water and (acetone, isopropyl alcohol, acetonitrile) binary solvent mixtures[J]. Journal of Molecular Liquids, 2016, 224: 1380-1387. DOI:10.1016/j.molliq.2016.10.127 |
[11] |
MARTHANDAN M V, ACREE W E JR. Solubility of anthracene in binary solvent mixtures containing dibutyl ether[J]. Journal of Chemical and Engineering Data, 1987, 32: 301-303. DOI:10.1021/je00049a005 |
[12] |
BUCHOWSKI H, KSIAZCZAK A, PIETRZYK S. Solvent activity along a saturation line and solubility of hydrogen-bonding solids[J]. The Journal of Physical Chemistry, 1980, 84: 975-979. DOI:10.1021/j100446a008 |
[13] |
GRADY L T, HAYS S E, KING R H, et al. Drug purity profiles[J]. Journal of Pharmaceutical Sciences, 1973, 62(3): 456-464. DOI:10.1002/jps.2600620321 |
[14] |
BARZEGAR-JALALI M, JOUYBAN-GHARAMALEKI A. A general model from theoretical cosolvency models[J]. International Journal of Pharmaceutics, 1997, 152: 247-250. DOI:10.1016/S0378-5173(97)04922-3 |
[15] |
ACREE W E JR. Mathematical representation of thermodynamic properties Part 2. Derivation of the combined nearly ideal binary solvent (NIBS) Redlich-Kister mathematical representation from a two-body and three-body interactional mixing model[J]. Thermochimica Acta, 1992, 198: 71-79. DOI:10.1016/0040-6031(92)85059-5 |
[16] |
YANG W G, FAN S M, GUO Q R, et al. Thermodynamic models for determination of the solubility of 4-(4-aminophenyl)-3-morpholinone in different pure solvents and (1, 4-dioxane+ethyl acetate) binary mixtures with temperatures from (278.15 to 333.15) K[J]. The Journal of Chemical Thermodynamics, 2016, 97: 214-220. DOI:10.1016/j.jct.2016.01.022 |
[17] |
WANG G, WANG Y L, ZHANG J, et al. Modeling and simulation of thermodynamic properties of L-alanyl-L-glutamine in different solvents[J]. Industrial and Engineering Chemistry Research, 2014, 53: 3385-3392. DOI:10.1021/ie404081c |
[18] |
ZHAO K F, YANG P, DU S C, et al. Determination and correlation of solubility and thermodynamics of mixing of 4-aminobutyric acid in mono-solvents and binary solvent mixtures[J]. The Journal of Chemical Thermodynamics, 2016, 102: 276-286. DOI:10.1016/j.jct.2016.07.014 |