基于物料粉体学性质对5种制剂品种颗粒溶化性的相关性研究

中药制剂质量受多种因素影响，中间体物料粉体学性质是制剂前重要的内容，而粉末与颗粒的质量控制是保证最终产品稳定可控的关键[1-2]。颗粒溶化性是中药颗粒剂的关键质量属性，在实际大生产过程中会出现颗粒溶化性问题，最终影响产品的质量。中药生产过程中，质量控制的工艺参数及控制限设计大多源于经验和片面了解[3]。实际生产工序上存在的溶化性问题的相关影响因素主要是化学成分及提取精制工艺差异[4]。课题前期调查研究发现，通过对多品种、多批次中药物料粉体物理属性的测定，能够从数字的角度探寻中药粉体物性参数与颗粒溶化性的相关性。研究中运用SeDeM专家系统[5-8]以及基于粉体学性质对粉末与颗粒进行综合表征[9-11]。同时以粉末的休止角（α）、松装密度（Da）、振实密度（Dc）、孔隙率（Ie）、卡尔指数（IC）、豪斯纳比（IH）、含水量（HR）、吸湿性（H）、粒径（D10、D50、D60、D90）、均匀性（UN）、比表面积（SSA）、相对均齐度（Iθ）、分布范围（span）、宽度（width）、粒径＜50 μm百分比（Pf）、中间体颗粒溶化率（DR）12个物性参数构建中间体粉末的物理指纹图谱与二维相关性矩阵热图评价各品种粉末批次间及品种指标间的差异性，结合多元统计方法偏最小二乘回归（partial least squares regression method，PLS）及主成分分析（principal component analysis，PCA）[12]对影响颗粒溶化率相关物性参数进行筛选，拟预测物料在工艺单元传递中的质量变化。

中药物理指纹图谱是继中药化学指纹图谱[13]、中药生物指纹图谱[14]而发展起来的用来表征粉末物理属性的图谱，以更好评价其质量一致性，预测中间体粉末与制剂的兼容性[15-16]，达到产品质量的可控有效。PLS是1种新型的多元统计方法，可以用来分析多变量对因变量的模型分析，实现多因子的筛选[17]。PCA是1种数据简化方法，主要利用降维思维把多个指标转化为几个主要综合指标，同时保持数据集对方差贡献最大的特征。2种统计方法结合分析可以很好地预测对颗粒溶化性产生影响的关键物性属性（critical material attributes，CMA），为后期工艺改善提供指导方向。

本实验中以杏贝止咳颗粒（Xingbei Zhike Keli，XZK）、桂枝茯苓胶囊（Guizhi Fuling Jiaonang，GFJ）、腰痹通胶囊（Yaobitong Jiaonang，YJ）、天舒胶囊（Tianshu Jiaonang，TJ）及参乌益肾片（Shenwu Yishen Pian，SYP）5个品种的中间体为研究对象，分别对各品种制粒前一步物料及制粒后的颗粒（各品种批次20批）进行物料粉体学性质对颗粒溶化性的影响相关性研究。

1仪器与材料

1.1仪器

BT1001型智能粉体特性测试仪、Bettersize 2600型激光粒度分布仪，丹东百特仪器有限公司；LHS-250HC-II型恒温恒湿箱、DHG-9145A型电热鼓风干燥箱、HWS26型电热恒温水浴锅，上海一恒科学仪器有限公司；Mettler Toledo 204型电子天平，德国梅特勒–托利多公司；MYP11-2型磁力搅拌器，上海梅颖浦仪器仪表制造有限公司。

1.2材料

XZK混合粉和颗粒（XZK1～XZK20）、GFJ软材细粉和颗粒（GFJ1～GFJ20）、TJ软材细粉和颗粒（TJ1～TJ20）、YJ喷干细粉和颗粒（YJ1～YJ20）、SYP混合粉和颗粒（SYP1～SYP20），以上各品种批次中间体物料均为一一对应，所有物料均购自江苏康缘药业股份有限公司。

2方法

2.1粉末与颗粒属性测定方法

2.1.1α 使用智能粉体特性测试仪测定，采用固定底面积法，底面为直径100 mm的圆盘，通过往进料筒里进入待测粉末，启动仪器后，样品经振动过筛后流落圆盘，待形成的锥体角度稳定后，读取数据，平行测定3次后取其平均值计算。

2.1.2D a 采用智能粉体特性测试仪测定。将待测样品加入进料筒内，经振荡过筛后流落至密度容器中，待粉末流满容器后，用刮板将多余的粉末刮去，分别称定加样前与加样后容器的质量m0与m1。

Da＝( m 1 －m0)/100

2.1.3D c 采用智能粉体特性测试仪测定。测定方法为固定体积法，按提示要求安装相应组件后，测定前称定空杯质量m0，向装有延长筒的量杯中加入待测粉末，经250 Hz，1250次振动后取下延长筒，将多余样品刮去，称定装满样品的量杯质量m1，平行测定3次。

Dc＝( m 1 －m0)/100

2.1.4I е 由Da与Dc计算而得，公式如下[9,15]。

Iе＝( D c－Da)/( D a×Dc)

2.1.5IC由Da和Dc计算而得，公式如下。

IC＝( D c－Da)/ D c

2.1.6IH由Da与Dc计算得到，公式如下。

IH＝Dc/ D a

2.1.7HR精密称定2 g粉末样品，均匀平铺于恒定质量后的称量瓶（m2）中，置于105 ℃干燥箱中恒定质量，分别称定恒定质量前与恒定质量后的质量m3、m4，平行测定样品3份，计算含水量，公式如下。

HR＝( m 3－m4)/( m 3－m2)

2.1.8H 取干燥洁净的称量瓶置于恒温恒湿箱中，温度（25±2）℃，湿度（75±2）%，平衡12 h后，取出称定质量（m5），分别取平行样3份，每份约1 g均匀平铺于称量瓶中，称定质量（m6），将称量瓶敞口置于上述条件恒温恒湿箱中24 h，盖上瓶盖后取出称定质量（m7），计算H，公式如下。

H＝( m 7 －m6)/( m 6 －m5)

2.1.9Pf和Iθ使用激光粒度分布仪，以空气为介质，将粉末样品加入干法进料斗中，计算粒径小于50 μm的粉末所占比例，即为Pf。分别选取粒径355、212、100、50 μm 4个粒径节点的粉末，计算其所占比例，公式如下[18]。

Iθ＝F m /[100＋( d m －d m －1) F m －1＋( d m ＋1－d m ) F m ＋1＋( d m －

d m －2) F m －2＋( d m ＋2－d m ) F m ＋2＋…＋( d m ＋n－d m ) F m ＋n]

F m 为粒径在100～212 μm的粉末质量百分比，F m －1为粒径分布在50～100 μm的质量百分比，F m ＋1为粒径分布在212～355 μm的质量百分比，d m 、d m －1、d m ＋1为相应粒径范围的平均粒径

2.1.10粒径（D10、D50、D60、D90）、span及width 使用激光粒度分布仪测定，以空气为介质，将样品加入干法分散系统的进料斗中，分别测定累积粒度分布数达到10%、50%、60%、90%时所对应的粒径D10、D50、D60、D90及粒径比表面积（specific surface area，SSA）、均匀性（uniformity，UN）。span和width计算公式如下[9,18]。

span＝( D 90 －D10)/ D 50

width＝D90－D10

2.1.11DR 参照《中国药典》2020年版3部0104颗粒剂溶化性测定方法[19]。称取约10 g待测样品（m8），加入200 mL热水（75±2）℃搅拌5 min后，趁热滤过，取滤液置恒定质量后的蒸发皿（m9）上蒸干后，于105 ℃干燥箱中干燥至恒定质量（m10），计算DR，计算公式如下。

DR＝( m 10 －m9)/ m 8

2.2各品种批次间物理指纹图谱的构建

物理指纹图谱既能表征中间体粉末和颗粒的物理性质，还能直观反映出各品种中间体粉末批次间的差异性和相似性[2]。为方便对2级质量属性参数进行分析，通过参考文献数据[20-21]，以稳定性、均一性、流动性、堆积性及可压性作为物理指纹图谱的一级指标，12个参数为二级指标。对5个品种各20批测定数据进行范围统计，确定各指标转化公式，以便计算各参数标准化数据范围（0～10）。各二级指标转化公式见表1。

2.3 数据分析软件

采用Sigmaplot 12.5（美国Systat软件公司）软件建立各品种粉体二级属性指纹图谱，绘制相应雷达图，结合Matlab 2020a（美国MathWorks公司）软件对批次间相似度进行评价。采用GraphPad Prism 8.0（美国GraphPad Software公司）对品种批次及品种指标间相似性绘图。采用SPSS 24.0（美国IBM公司）对各品种粉体属性的物理参数进行标准化后，结合SIMCA 14.1（瑞典MKS Umetrics公司）进行PCA分析及PLS模型预测。

3结果与分析

3.1制剂粉末物理属性的数据特征

以5个品种制粒前一步粉末及制粒后颗粒为测定对象，各取20批中间体物料对相关物性参数进行测定，部分均值数据见表2。数据显示，XZK中间体混合粉α均大于46°，Da为0.24～0.31 g/mL，Dc为0.56～0.73 g/mL，HR为6.21%～7.48%，H为15.71%～25.04%，width、span及粒径小于50 μm的范围分别为36.75～63.42 μm、1.82～2.32 μm、68.32%～93.88%；GFJ软材细粉α均大于40°，TJ软材细粉、YJ喷干细粉及SYP干膏粉α均大于37°，其中YJ喷干细粉与SYP干膏粉的H分别为18.02%～27.14%、27.19%～30.59%，同时两者HR分别为3.44%～6.05%、4.84%～7.24%。

根据《中国药典》2020年版[19]中对吸湿性要求，当2%＜H＜15%时具有引湿性，XZK混合粉的粒径分布区间范围较大，且存在吸湿现象；干燥失重要求不得超过2%，说明XZK、YJ与SYP中间体粉末含水量影响较大。综上所述，各品种中间体粉体在α、HR、H及粒径分布区间范围差异较大。按照《欧洲药典》规定[20]，α＞35°、IH＞1.35且IC＞0.25时，流动性差，表明5个品种的流动性不佳，各品种物理属性及批次间差异需做进一步分析。

3.2品种批次间差异性评价

对HR、H、Pf、Iθ、span、width、IH、α、Da、Dc、Iе、IC共12个物性参数构成的物理指纹图谱二级指标进行转化，使其数值在1～10。分别对所选品种粉体批次间相似度进行评估，并最终以雷达图的形式直观可视地呈现粉末物理指纹图谱，其中阴影部分面积越大代表各品种粉末批次间差异性越小，产品生产越稳定，见图1。基于夹角余弦值方法[22]，物理指纹图谱相似度越接近于1，表明各品种批次间粉体越相近。各品种批次间相似度结果显示，XZK的混合粉相似度分布在86.52%～99.92%，YJ的喷干细粉相似度分布在87.98%～99.90%，SYP干膏粉相似度在86.94%～99.97%，GFJ相似度在90.38%～99.97%，TJ的软材细粉相似度在89.79%～99.99%。并以5个品种的12个指标转化后参数的平均值绘制雷达图作为对照图谱，其平均相似度为94.14%。结果说明各品种中间体粉末批次间粉体学性质较稳定，适合做进一步统计分析。

3.3制剂品种间物理属性相关性分析

为了对各品种物理属性相关性进行更全面分析，在原来12个物性指标基础上增加粒径指标（D10、D50、D60、D90）及UN、SSA进行分析。以5个品种各20批粉体物料构建的物性参数相关性矩阵热点图分析见图2。热图中颜色越深或色差越小代表其品种间相似性越高，反之，颜色越浅或色差越大，表明相关性弱。图2-A中反映各品种批次在18个物理属性上的相关性，结果显示，各品种在Da、Dc、Ie、IC、IH、HR、Iθ及span上存在较高相似性，在α及粒径分布范围上存在差异性，各品种指标间相关系数为?0.964～0.998。

图2-B中反映其二维矩阵相关性，结果显示，Dc、D10与α呈显著相关（P＜0.05），H、HR及SSA与Da呈显著相关（P＜0.05），Iθ及粒径D50、D60、D90、UN、span、width、Pf与吸湿性呈显著相关（P＜0.05），同时发现粉体间均匀性会对其豪斯纳比及吸湿性产生影响，表明制剂物料的吸湿性受多种物性属性影响。

3.4基于粉体学性质的PCA

各品种不同批次的18个粉体物性参数在PCA模型中的分析结果见表3。结果显示，主成分1对模型主成分变量解释率为57.06%，主成分2主变量解释率为20.76%，主成分3主变量解释率为9.23%，前2个主成分累积解释变量达77.82%，符合统计学要求[23]。故以前2个主成分参与统计学分析。

基于PCA所得到的载荷分散图见图3，图中横坐标代表主成分1，纵坐标代表主成分2，各参数指标距离原点越远，说明对相应主成分贡献率越大。可以看出，IC、IH、SSA及代表粒径性质的D10、D50、D60、D90与Iθ、width、Pf对主成分1贡献率较大；主成分2由α、Da、Dc、Ie、HR及UN、span构成。

图4的得分图中可以看出，基于主成分1与主成分2之间的关系，XZK混合粉及GFJ软材细粉分布较集中，反映批次间物性差异较小；YJ喷干细粉、TJ软材细粉及SYP干膏粉分布较离散，受物理属性影响波动较大。此外，XZK相较于其他品种相对独立，也反映干法制粒与湿法制粒过程中物料属性的差异性，为课题组后期改善颗粒溶化性提供思路。整体而言，各品种物料粉体属性在模型内，适合用于模型建立的预测。

3.5制剂中间体颗粒溶化性CMA的筛选

对18个物性参数进行标准化处理后作为自变量，以颗粒DR为因变量进行PLS模型分析，得到的模型性能指标：校正决定系数R2X＝0.953，预测决定系数R2Y＝0.852，交叉验证决定系数Q2＝0.767（R2X、R2Y和Q2越接近1，表明模型具有较好的校正或预测性能），其模型相对预测误差为8.40%小于10%，说明模型预测性能良好。图5反映的是物料粉体对颗粒溶化性响应值的相关系数。从中也可以看出，Dc、Ie、IC、IH、HR、H、SSA、Iθ及代表粒径的D50、D60、D90、width与颗粒DR呈显著正相关，粒径D10、span、Pf及UN与颗粒DR呈负相关。其粉体学性质对颗粒DR的投影变量重要响应值（variable importance for the projection，VIP）见图6。VIP值大于1，说明该变量对因变量的影响显著。可以看出，对DR影响显著的是Dc、H、SSA、HR、IH、Ie、IC及span。H、HR反映物料粉体的稳定性，Ie、IC反映粉体的可压性，Dc、span及SSA分别为粉体的堆积性、均一性、表面形态指标，因此在工艺生产上可以通过物料粉体的稳定性、压缩性、表面形态等参数指标来改善颗粒的溶化性。因此除Iθ与代表粒径的D50、D60、D90、width外，对DR呈显著影响的指标一致，即Dc、Ie、IC、IH、HR、H、SSA是影响中间体颗粒溶化性的CMA。

4讨论

本实验以5个制剂品种各20批制粒前粉末与颗粒为研究对象，以12个物性参数建立各品种中间体粉末的物理指纹图谱与相似性矩阵热图来评价各品种批次及品种间的差异性。结果表明，各品种批次间相似度波动不大，说明品种物料的属性较稳定。进而对其进行PCA及PLS分析，结果显示，Dc、Ie、IC、IH、HR、H、SSA是影响中间体颗粒溶化性的CMA，其反映物料粉体的堆积性、稳定性、流动性、可压性及表面形态会对颗粒溶化性工艺造成影响，后期将针对相关指标的关键因素进行分析与改进。

中药制剂质量问题受多重因素影响，通过对制剂物料进行粉体学性质表征，可以解析物料属性与制剂质量的关联性[24]。中药固体物料的表面性质与固态物料粒子的粒径、孔隙率、表面粗糙度和形态等因素相关[25]，因此在实际工艺生产中通过控制粒径来改变表面性质，粒子表面包覆技术作为表面改性的常用手段[26]。例如，根据中药浸膏粉末粒子的粒径、孔隙、形貌等空间构造性质与吸湿性、压缩性等面、体性质的相关性，运用辅料及工艺改性技术解决制剂成型性问题[27]。针对制剂工艺存在的析出、结块等溶化性问题，也可以考虑辅料配用及制粒方式来达到改善目的。

本研究中运用多元统计方法对关键变量进行辨识以及预测，将制剂中间体颗粒溶化性与粉体物理性质联系起来。研究中收集的品种批次样本数据远小于生产过程中产生的大批量数据，采用的PCA及PLS算法相对简单，建立的模型相对预测误差偏大。后期将考虑采用算法的更换或者次变量的剔除，来降低其相对预测误差，同时增大模型分析的样本量，提高模型预测效能，实现相关共性技术平台的搭建，以期加快中药制剂过程质量控制的规范化与现代化进程。

参考文献（略）

来 源：汪盛华，闫 明，徐芳芳，章晨峰，王团结，王振中，肖 伟．基于物料粉体学性质对5种制剂品种颗粒溶化性的相关性研究[J]. 中草药, 2022, 53(22): 7082-7090 .