Pharma Mod UNIFAC预测磺胺类原料药水溶解度的热力学敏感性、参数可迁移性与可辨识性审计

《Canadian Journal of Chemical Engineering》：Thermodynamic sensitivity and transferability limitations of Pharma Mod UNIFAC for aqueous solubility prediction of sulfonamide APIs

　　

研究人员在本工作中评估了Pharma Mod UNIFAC预测磺胺类活性药物成分（active pharmaceutical ingredients, APIs）水溶解度的能力，并将新的溶解度与量热测量同可溯源文献数据相结合。在评估参数可迁移性之前，研究人员验证了每个体系的分子表示、固相参考性质以及关键水–磺胺相互作用对。直接应用时，原始Pharma Mod UNIFAC参数能够重现定性趋势，但在所评估的APIs之间不能一致迁移。对关键水–磺胺相互作用进行局部常数重参数化后，若干体系的计算溶解度得到显著改善；然而，拟合参数在结构相关化合物之间发生变化，因此应将其解释为局部有效校正，而不是对原始数据库的通用替代。温度依赖的局部回归表明，不同函数形式，如a_mk+b_mkT和a_mk+c_mkT²，在所研究温度区间内可给出几乎不可区分的有效轨迹，揭示了拟合系数的可辨识性限制。方向敏感性分析表明，对同一相互作用对进行正向和反向扰动并非对所有APIs等价。总体而言，Pharma Mod UNIFAC仍是预测水溶解度的有用且结构化的框架，但其应用于高度稀释的药物体系时需要仔细的局部验证。主要贡献是对参数可迁移性的热力学审计，而不是新的通用相互作用参数矩阵。

研究背景与意义：在高度稀释的药物水溶液中，溶解度对活度系数（activity coefficient, γ）及固相参考热力学性质高度敏感，即使活度系数模型的小偏差也可能被放大为相平衡组成的显著偏差。磺胺类活性药物成分（active pharmaceutical ingredients, APIs）具有多官能团、固相相互作用强、常温低水溶性及分子结构复杂等特征，使固液平衡（solid–liquid equilibrium, SLE）的热力学描述困难。基团贡献方法如UNIFAC及Pharma Mod UNIFAC可用于估算多组分体系活度系数，但通用基团相互作用参数能否迁移到结构复杂的API–水体系仍不清楚。本研究因此将新溶解度与量热数据同可溯源文献数据结合，评估Pharma Mod UNIFAC对磺胺类API水溶解度的预测能力、关键水–磺胺相互作用参数的可迁移性、温度依赖性与方向敏感性。研究人员发现原始参数可重现定性趋势但不能一致迁移；局部常数重参数化可显著改善若干体系，但拟合值随结构相关化合物变化，应视为局部有效校正；温度依赖形式可进一步降低偏差，但窄温区内不同函数形式近不可区分，数据主要约束局部有效a_mk(T)轨迹而非独立可迁移的a_mk、b_mk、c_mk系数；方向敏感性具有体系依赖性。其意义在于提供参数可迁移性的热力学审计框架，而非提出新的通用参数矩阵，并强调Pharma Mod UNIFAC用于高度稀释药物体系时需局部验证。

关键技术方法（不超过250字）：研究采用静态平衡法测定四种磺胺类API在水中的溶解度，并以紫外–可见光谱定量；采用差示扫描量热法（differential scanning calorimetry, DSC）测定熔点和熔化焓；以Schroeder–Van Laar固液平衡关系计算溶解度；在Aspen Plus?中调用Pharma Mod UNIFAC计算活度系数，并用其数据回归系统对关键水–磺胺相互作用对进行局部常数及温度依赖重参数化；开展结构/基团与固相参考认证、无限稀释活度系数（γ^∞）计算及正反向敏感性分析。样本来源为Sigma-Aldrich提供的磺胺类API实验样品、Milli-Q纯化水，以及经来源追溯和条件分类的文献溶解度数据；无队列样本。

研究结果：4.1 实验溶解度、量热数据与数据可溯源性：通过静态平衡实验获得sulfisoxazole、acetazolamide、ethoxzolamide和furosemide在298.15–313.15 K的溶解度，DSC获得熔点及部分熔化焓；文献数据按来源、温度、介质、不确定度和热力学量分类。结论是建立了可用于SLE分析的可溯源数据集，并强调应区分表观溶解度、总溶解度和中性物种固有溶解度。furosemide在熔化前分解，采用伪固相参考作探索性处理，三种表征方法给出的SLE摩尔分数在报告精度内相同，AARD为9.09%，但不用于支持通用迁移性结论。4.2 结构与热力学认证：通过在Aspen Plus?中核查化学式、摩尔质量、分子结构、Pharma Mod UNIFAC基团分配及固相参考性质，确认二级磺胺对应SO₂NH、一级磺胺对应SO₂NH₂，水组为H₂O^PHA，关键对为H₂O^PHA?SO₂NH或H₂O^PHA?SO₂NH₂。结论是后续偏差主要可归因于活度系数模型及参数，而非体系定义错误。4.3 原始参数直接迁移性：通过直接应用原始Pharma Mod UNIFAC参数比较计算与实验溶解度，发现部分体系可重现温度上升趋势但存在系统偏移，且预测表现因API而异。结论是仅凭磺胺基团身份不足以保证定量预测，需开展迁移性审计。4.4 常数参数表示作为基线迁移性假设：通过分析数据库关键相互作用对的常数参数表示，明确式(3)中b_mk和c_mk项未激活，常数表示是数据库参数的基线假设。结论是直接迁移偏差需区分有效相互作用水平偏移、温度依赖描述不足及更根本的局部环境迁移限制，温度回归应作为诊断扰动而非新通用参数。4.5 关键水–磺胺相互作用的局部常数重参数化：通过对各API关键对进行局部常数重估，发现计算溶解度总体改善，但响应不均一；sulfapyridine、sulfamethoxazole和acetazolamide等经常数校正后偏差大幅降低，sulfacetamide、sulfathiazole、sulfisoxazole和ethoxzolamide仍保留温度结构残差。结论是常数重参数化可识别“局部水平偏移”是否足够，但拟合常数在结构相关化合物间变化，应解释为局部有效校正而非通用替代。4.6 温度依赖局部有效轨迹：通过对常数校正后仍有温度结构残差的体系进行a_mk+b_mkT和a_mk+c_mkT²局部回归，发现若干体系残差显著降低；二级磺胺主要涉及H₂O^PHA→SO₂NH方向，一级磺胺涉及H₂O^PHA→SO₂NH₂方向，且两种函数形式在研究温区内给出近等价轨迹。结论是数据约束局部有效a_mk(T)轨迹，而非独立可迁移的a_mk、b_mk、c_mk系数。4.7 函数非唯一性与参数可辨识性：通过比较线性和二次温度依赖回归及拟合参数相关性，发现窄温区和小数据量下不同函数形式可近不可区分，拟合系数存在强相关和参数补偿。结论是局部拟合精度不等于参数可迁移性，不能将单个a_mk、b_mk、c_mk系数解释为独立可辨识或通用可迁移参数。4.8 关键相互作用方向敏感性：通过分别对正向H₂O^PHA→SO₂NH/SO₂NH₂和反向SO₂NH/SO₂NH₂→H₂O^PHA施加温度依赖扰动，发现sulfathiazole、sulfisoxazole和ethoxzolamide的正向校正显著优于反向，acetazolamide正反向近等价。结论是UNIFAC型基团相互作用非对称，方向敏感性为体系依赖，不能仅由SO₂NH或SO₂NH₂形式推断。4.9 无限稀释行为及对Pharma Mod UNIFAC迁移性的启示：通过298.15 K原始参数计算无限稀释活度系数，得到sulfisoxazole为93.52、ethoxzolamide为4425.86、acetazolamide为23.42，表明强非理想性和不利的有效溶质–水相互作用。结论是无限稀释区放大局部水–API相互作用表示差异，局部参数调整可降低特定API数值误差但不保证迁移性；主要限制在于以通用基团贡献参数表示结构特异的水–API相互作用。

讨论与结论：讨论部分指出，直接迁移、局部常数重参数化、温度依赖有效轨迹、函数非唯一性及方向敏感性共同表明，限制并非仅在于缺少可调参数，而在于高度稀释热力学区中通用基团参数难以表征结构特异的局部水–API相互作用；Pharma Mod UNIFAC仍提供识别控制溶解度之相互作用并诊断通用参数局部不足的结构化热力学基础，但用于复杂磺胺类API时须谨慎解释迁移性。结论译述：本研究将Pharma Mod UNIFAC作为预测高度稀释磺胺类API水溶解度的热力学框架进行评估，结合新溶解度与量热测量及文献数据构建可溯源SLE数据集，并关注实验溶解度数据一致性、固相参考热力学性质、分子碎片化和关键水–磺胺相互作用对。固相参考定义在SLE建模中很重要；对熔化前分解的API，实验熔化性质可能不可得，伪固相参考可用于探索性计算，但不应据此推断通用迁移性。原始参数直接应用可重现定性趋势，但标准参数化不能一致迁移到所评估磺胺类API；仅存在磺胺基团不足以保证定量预测。关键水–磺胺相互作用的局部常数重参数化显著改善若干API计算溶解度，但拟合常数在结构相关化合物间差异明显，应解释为局部有效校正而非通用替代；部分体系常数校正近足够，另一些仍存在系统残差，提示局部有效相互作用需温度依赖表示。温度依赖局部回归进一步降低选定体系偏差，但a_mk+b_mkT与a_mk+c_mkT²在窄温区内可近不可区分，数据约束局部有效a_mk(T)轨迹而非独立可迁移的a_mk、b_mk、c_mk系数，显示局部准确拟合本身不能建立参数可迁移性。方向敏感性分析表明同一水–磺胺相互作用对的正向与反向扰动对计算溶解度的影响不等价，且方向敏感性由局部分子环境控制而非仅由官能团身份决定。总体而言，本研究主要贡献不是提出新的通用参数矩阵，而是建立Pharma Mod UNIFAC参数可迁移性的热力学审计，区分直接参数迁移、局部常数校正、局部温度依赖有效相互作用轨迹、函数非唯一性和方向敏感性；Pharma Mod UNIFAC仍是有用的结构化水溶解度预测框架，但用于高度稀释药物体系时需仔细局部验证，尤其是磺胺基团处于多样分子环境中的APIs。