1. 引言
壮骨关节胶囊系由壮骨关节丸改制为胶囊剂而成，壮骨关节丸收载于中国药典2020年版一部是由淫羊霍、补骨脂、桑寄生、独活、骨碎补、鸡血藤、狗脊、续断、熟地黄等12味中药组成的复方制剂 [1] ，具有补益肝肾、养血活血、舒筋活络、理气止痛等功效，用于肝肾不足、血淤气滞、经络痹阻所致的骨性关节炎、腰肌劳损，症见关节肿胀、疼痛、麻木、活动受限等疾病 [2] [3] [4] 。现代药理学研究证明：淫羊藿主要有效成分淫羊藿苷，能够显著促进人成骨细胞的增殖和分化，能够使BMP2 mRNA的表达升高 [4] ；补骨脂主要有效成分补骨脂素和异补骨脂素能提高大鼠成骨细胞内碱性磷酸酶活性，促进成骨细胞的增殖与分化 [5] [6] [7] ；桑寄生、鸡血藤、骨碎补均具有促进骨细胞生长，抗骨质疏松作用；其主要有效成分分别是槲皮素、芒柄花素、柚皮苷 [8] 。目前，对壮骨关节胶囊的研究大多为临床、药效方面 [9] [10] [11] [12] ，鲜有涉及含量测定 [13] [14] 。
近红外光谱分析(Near-Infrared Spectroscopy, NIRS)技术作为检测技术中的一朵奇葩，其具有快速、无损、绿色等特点，目前已成功应用于食品、烟和化工等诸多行业产品的分析测定中，可以用于对自动化生产进行监控、调节、反馈生产。2002年美国FDA已把该项技术作为一种标准的检测方法，2005版《中国药典》已将“近红外分光光度法指导原则”列入附录。近红外光谱技术在制药工业中的应用日趋广泛，从药物的定性、定量分析，到生产过程各个阶段(包括合成、制剂、包装等过程)的监控都体现出近红外光谱的巨大潜力。目前还没有关于近红外应用于壮骨关节胶囊的研究报道。本研究主要采用近红外光谱技术结合化学计量学方法–偏最小二乘(PLS)，对壮骨关节胶囊中流浸膏指标成分总黄酮醇苷进行检测分析，并建立定量分析模型，以期提供一种新的快速检测壮骨关节胶囊流浸膏质量的方法，能够快速、准确的监控壮骨关节胶囊流浸膏质量，为壮骨关节胶囊生产提供质量合格的壮骨关节胶囊流浸膏。
2. 实验部分
2.1. 仪器与试剂
仪器信息见表1，试剂材料见表2。
Table 1. Experimental instruments and their manufacturers and models
表1. 实验使用仪器及其厂家、型号
Table 2. Reagents used in the experiment and their manufacturers and batch numbers
表2. 实验使用试剂及其厂家、批号
2.2. 近红外漫反射光谱采集
壮骨关节胶囊流浸膏光谱采集预处理：由华润三九股份有限公司技术部负责收集，共取样65个。
光谱采集条件：采集壮骨关节胶囊流浸膏样品的漫反射光谱，光谱扫描范围4000~12000 nm，分辨率16 cm−1，每个样本扫描3次，取平均光谱如图1。实验采用内部参比进行光谱采集分析。
Figure 1. Raw near-infrared spectrum of ZhuangGuGuanJieWan fluidextract
图1. 壮骨关节胶囊流浸膏近红外原始光谱
2.3. 壮骨关节胶囊流浸膏总黄酮醇苷的测定
对照品溶液的制备取朝藿定C、淫羊藿苷对照品适量，精密称定，加稀乙醇制成每1 ml各含20 μg的混合溶液，即得。
供试品溶液的制备取本品适量，研细，混匀，取约0.45 g，精密称定，精密加入稀乙醇25 ml，称定重量，超声处理45分钟，放冷，再称定重量，用稀乙醇补足减失的重量，摇匀，滤过，取续滤液，即得。色谱检测条件见表3。
Table 3. Chromatographic detection conditions
表3. 色谱检测条件
理论板数按淫羊藿苷峰计算应不低于5000，梯度洗脱见表4。
2.4. 数据处理与模型性能评价
使用赛默飞光谱仪利用漫反射原理，对样品进行扫描，采集光谱。将光谱数据导入TQ Analyst分析软件，为保证光谱的有效信息保留，需要将光谱中的噪声、光谱偏移及背景等造成的误差进行处理，一般采用平滑、基线矫正、MSC、导数等预处理方法，进行处理。将预处理后的光谱进行分析，为保证模型的准确性，需要将对应总黄酮醇苷的波段进行筛选，缩短建模时间。
使用偏最小二乘法对建立总黄酮醇苷模型，验证方式选用运用留一交叉验证法，计算预测残差平方选择模型最佳主成分数。对建模数据的训练集和验证集的误差均方根进行考察，越接近则模型预测结果约准确。并通过模型的相关系数判定模型效果，相关系数越接近1，则模型稳定性越强。在建立模型后，还需对模型的稳定性、准确性及差值显著性进行考察。
3. 结果与讨论
3.1. 壮骨关节胶囊流浸膏总黄酮醇苷含量测定结果
按照2.3项下的方法测定65个壮骨关节胶囊流浸膏中总黄酮醇苷的含量，总黄酮醇苷含量的范围0.327~2.281 mg/g之间。
3.2. 光谱预处理
为保证采集光谱含有建模所用的有效信息，对光谱进行预处理，常见预处理有导数、SNV、MSC等处理方法，通过对比不同预处理方法，壮骨关节胶囊流浸膏光谱采用一阶导数的预处理方法，可以将光谱中的背景噪声进行剔除，并将光谱之间的差异放大，为后续利用含量和光谱的线性关系建模起到重要作用。并通过对比建立模型的误差均方根和相关系数对光谱波段的筛选，比较，最终确定建模区间为8870.95~7729.29 cm−1，此波段与总黄酮醇苷的对应性最好，结果如表5所示。
Table 5. Parameters of total flavonol glycosides content model
表5. 总黄酮醇苷含量模型参数
通过比较发现，壮骨关节胶囊流浸膏总黄酮醇苷含量模型通过进行光谱预处理处理后得到的模型有较低的RMSEC值，且相关系数R值更接近1，模型效果见图2。
Figure 2. Schematic diagram of the total flavonoid glycoside model of ZhuangGuGuanJieWan fluidextract
图2. 壮骨关节胶囊流浸膏总黄酮醇苷模型示意图
4. 结论
4.1. 准确性
Figure 3. Prediction effect of ZhuangGuGuanJieWan fluidextract
图3. 壮骨关节胶囊流浸膏预测效果图
从预测结果来看，预测值与真实值接近，符合中药分析的精度要求(如图3)。用NIR定量校正模型验证壮骨关节胶囊流浸膏总黄酮醇苷含量的预测效果优良。
4.2. 显著性
Table 6. Paired samples t-test table
表6. 成对双样本均值的t-检验表
使用SPSS分析软件对变量样本进行t检验，结果如表6所示，P > 0.05，真实值与预测值之间无显著性差异。
4.3. 差值的正态分布
Figure 4. Normal Distribution Test Histogram
图4. 正态分布检验直方图
为精确反应模型显著性，计算模型数据残差，P > 0.05 (如图4)，表明差值没有显著性。
从上述三个角度来验证壮骨关节胶囊流浸膏总黄酮醇苷模型，其结果均符合要求，可以用于快速检测壮骨关节胶囊流浸膏总黄酮醇苷含量，实现壮骨关节胶囊流浸膏质量的快速检测。
5. 结论
本研究建立了壮骨关节胶囊流浸膏总黄酮醇苷含量的近红外快速分析方法，应用PLS建立的壮骨关节胶囊流浸膏总黄酮醇苷含量的定量校正模型，模型考察的相关系数结果较好，且误差均方根接近，模型具有较高的稳定性和准确性，满足定量分析的要求。所建立的近红外模型能够在较短的时间内通过采集近红外光谱对壮骨关节胶囊流浸膏总黄酮醇苷含量进行测定，相对传统检测方法，该方法成本低、检测准确、快速。