摘要:传统的阿里红中齿孔酸含量测定采用高效液相色谱法(HPLC),但该方法前处理复杂,操作繁琐。为了实现对中药阿里红饮片中齿孔酸含量的快速无损监测,尝试建立基于近红外光谱(NIR)的偏最小二乘(PLS)回归模型。用来预测阿里红中齿孔酸的含量,采用传统的HPLC方法对阿里红中的齿孔酸含量进行测定,其结果作为指标值。采集近红外数据后使用五种光谱变换方法对光谱数据进行预处理,即多元散射校正(MSC)、标准正态变化(SNV)、 Savitzky-Golay平滑(7点)、一阶导数变换(FD)和二阶导数变换(SD)。通过竞争自适应重加权法(CARS)进行波长选择并对PLS模型进行优化,大大减少了光谱变量的数量,并显著提高了PLS模型的性能,尤其是SNV-CARS-PLS模型,仅占总光谱波长的5.53%,预测集的R2值为0.982 3,预测均方根误差(RMSEP)值为0.103 7%,残差预测偏差(RPD)值为5.34。通过t检验表明:该最优模型与传统HPLC法在预测阿里红中齿孔酸含量上没有显著差异。研究结果表明:近红外光谱结合竞争性自适应重加权算法对波长筛选后建立偏最小二乘回归模型对阿里红饮片中齿孔酸含量的无损检测可行。
加入收藏
阿里红是中药材, 主要来源于多孔菌科真菌药用拟层孔菌的子实体, 也被称作落叶松茸, 主要分布在新疆、 内蒙和四川等地[1]。 阿里红的性味归经为甘、 苦、 温, 具有温肺化痰、 降气平喘、 祛风除湿、 活血消肿、 利尿、 解蛇毒等功效, 常用于治疗咳嗽、 哮喘、 胃痛、 胃酸过多、 尿路结石、 肾炎、 慢性风湿性关节炎、 咽喉炎、 牙周炎以及毒蛇咬伤等症状[2]。 阿里红中的化学成分主要包括三萜类、 倍半萜类、 甾体类化合物、 多糖以及一些其他的直链烷烃类结构, 其中以三萜类化合物成分最多, 齿孔酸和去氢齿孔酸为其主要的三萜类成分。 三萜类成分是药用真菌的有效成分之一, 具有提高人体免疫功能、 抗衰老等多种保健功效。 其还含有齿孔醇、 齿孔二醇、 去氢硫色多孔菌酸、 鲨鱼烯、 去氢齿孔酮酸、 落叶松脂酸、 麦角甾醇等成分[3],主要成分化学式见表1。
为了更好地控制阿里红的质量, 实验室一般采用HPLC对阿里红中齿孔酸含量进行测定。 然而HPLC前处理繁琐,需要长时间进行分析, 设备、 操作和维护成本很高。 为了克服这些缺点, 近红外光谱技术是HPLC法的一种有利的替代策略。 近红外光通常是指在800~2 500 nm波长范围内的光, 在近红外光谱(NIR)中, 未知物质用光谱(许多波长或频率)的近红外光照射, 这些近红外光可以被样品吸收、 透射、 反射或散射, 在与样品中分子的某些键的振动频率的组合和泛音相对应的特定频率下, 光被样品吸收[4]。 偏最小二乘回归是一种线性统计算法, 在普通多元回归的基础上结合了主成分分析和正则相关分析, 可以用以解决自变量的多重共线问题[5]。 随着化学计量学的发展, 偏最小二乘法与近红外相结合已在食品、 农业、 石油和医学领域得到广泛应用[6-11]。 本研究在结合近红外光谱和偏最小二乘法建立定量模型用以预测阿里红中齿孔酸含量。
表1阿里红中主要化学成分
1、实验部分
1.1仪器
DHG-9055AD电热恒温干燥箱(上海齐欣科学仪器有限公司); 1100高效液相色谱(安捷伦); KQ-400KDE型高功率数控超声波清洗器(昆山市超声仪器有限公司); SolidSpec 3700DUV近红外光谱仪(岛津, 日本); Matlab R2021b用于数据处理、 作图和建模。
1.2样品
90%乙醇(分析纯, 天津致远化学试剂有限公司);齿孔酸对照品(NLT:≥98.0%购自四川精萃天成药物科技有限公司);乙腈(色谱纯,Merck公司, 批号:K54731230);阿里红饮片90份(新疆恩萨尔维吾尔医药饮片制药有限公司),该饮片由刘戈宇副研究员鉴定为阿里红药材。 实验前, 样品需粉碎过40目筛并置于电热恒温干燥箱中60℃干燥4 h以除去水分。
2、结果与讨论
2.1 HPLC测定阿里红中齿孔酸的含量
2.1.1色谱条件
色谱柱:Hedera ODS-2 C18柱(4.6 mm×250 mm, 5μm),如图1所示;
流动相:A相为0.1%磷酸(A), B相为乙腈; 检测波长:210 nm;流速1.0 mL·min-1;进样量10μL;
梯度洗脱时间程序:0~20 min保持A-B(40∶60), 20~28 min,线性变化至A-B(10∶90), 28~48 min,保持A-B(10∶90), 48~51 min,线性变化至A-B(40∶60), 51~53 min,保持A-B(40∶60);
柱温:30℃。
2.1.2溶液制备
齿孔酸对照品溶液制备: 精密称取齿孔酸对照品24.75 mg,加入90%乙醇溶解并定容于25 mL量瓶中, 摇匀, 即得浓度为0.97 mg·mL-1的对照品溶液, 备用。
供试品溶液制备: 称取干燥阿里红药材粗粉(过40目筛)约0.5 g,置于15 mL具塞试管中, 并用15 mL 90%乙醇处理; 样品混合物在60℃下超声处理30 min(功率200 W),过滤并收集滤液, 再加90%乙醇15 mL洗涤残渣继续超声提取2次, 合并滤液浓缩并定容至50 mL摇匀, 然后通过0.45μm有机滤膜, 取续滤液作为供试品溶液。
空白样品溶液:90%乙醇用0.45 mm滤膜过滤, 续滤液作为空白样品溶液。
图1高效液相色谱
2.1.3线性关系考察
精密量取对照品溶液1.00, 3.00, 5.00, 8.00 mL分别置于10 mL量瓶中, 加90%乙醇定容, 摇匀。 分别吸取上述各浓度对照品线性溶液和对照品溶液各10μL依次注入液相色谱仪, 按“2.1.1”项色曲线, 如图2所示得回归方程为Y=3 435.6x+440.82,r=0.999 8,结果表明在97~970μg·mL-1内, 浓度和峰面积的线性关系良好。
图2齿孔酸标准曲线
2.1.4仪器精密度试验
制备同一份供试品溶液连续进样6次, 每次精确进样10μL,记录6次峰面积并计算出RSD为0.041%,表明仪器精密度良好, 测定结果可信, 结果见表2。
2.1.5稳定性试验
精密吸取同一份供试品溶液分别于0, 4, 8, 12, 18, 24 h每次精确进样10μL,并记录每次峰面积并计算出RSD为0.057%,表明本品在室温条件下24 h内稳定性良好测定结果可信, 结果见表3。
2.1.6重复性试验
取同一份阿里红药材6个平行样, 按供试品溶液制备方法制成样品溶液, 按“2.1.1”项下色谱条件测定, 记录6次峰面积并计算出含量和RSD分别为2.41%和1.73%,表明该方法重复性良好, 测定结果可信, 结果见表4。
表2精密度试验结果
表3稳定性试验结果
表4重复性试验结果
2.1.7回收率试验
取同一份阿里红药材6个平行样约0.5 g,精密称定, 并加入一定量的对照品, 按供试溶液制备方法制成样品溶液, 精密吸取10μL,按“2.1.1”项色谱条件测定, 计算回收率。 平均回收率为100.81%, RSD为1.81% (n=6),表明回收率良好, 结果见表5。
2.1.8阿里红中齿孔酸含量计算并划分样本
取90批阿里红药材, 每一批样品平行制备3份, 按“2.1.1”项下色谱条件测定, 计算90份样品中齿孔酸的含量, 齿孔酸含量范围为1.27%~4.35%,采用K-S样本划分法, 将样品以3∶1比例进行校正集和预测集划分, 最终校正集为68个样品, 预测集为22个样品, 如表6所示。
2.2近红外光谱法建模
2.2.1光谱采集
使用近红外光谱仪(SolidSpec 3700DUV, Shimadzu, Japan)采集近红外光谱。 仪器配备了水平积分球漫反射采样元件、 近红外区的InGaAs和PbS探测器。NIRS仪器系统在测量前用白板进行校准, 收集了800~2 500 nm范围内的光谱, 分辨率为0.1 nm。 每个样品重新加载并记录3次, 计算平均光谱并用于数据处理。
表5齿孔酸加样回收率试验结果(n=6)
表6阿里红样本划分
2.2.2光谱预处理
原始光谱是药材最直接的表现形式, 光谱变换可以不同程度地削弱或消除背景噪声, 对优化特征谱带和提高谱带灵敏度至关重要。 本研究选择了五种常用的光谱变换方法进行光谱数据预处理, 包括一阶导数(FD)、 二阶导数(SD)、Savitzky-Golay平滑算法(SG)、 多元散射校正(MSC)和标准正态变化(SNV),结果如图3所示。
2.2.3近红外光谱解析
90个样品的近红外光谱特征如图3(a)所示。 在1 198、1 467、1 733、1 949、2 117和2 301 nm附近有6个主峰。1 198 nm的吸收来源于C—H伸缩和第二泛音。1 467、1 733和1 949 nm处的峰值分别来自O—H和O—H组合带的第一泛音。2 117和2 301 nm处的峰是由于脂的CH2基团中C—H和C—O伸缩振动的组合带。 图3(b)为多元散射校正(MSC)处理后的90个样品的近红外光谱, 图3(c)为标准正态变化(SNV)处理后的90个样品的近红外光谱, 图3(d)为Savitzky-Golay平滑算法(7个点)处理后的90个样品的近红外光谱图, 图3(e)为一阶导数(FD)处理后的90个样品的近红外光谱图, 图3(f)为二阶导数(SD)处理后的90个样品的近红外光谱图。 视觉观察,90个样品的NIR光谱显示出非常微小的差异, 并且所有峰的趋势和吸收强度都非常相似。 这些细微的差异可用来通过化学计量学进行量化。
2.2.4波长点选择方法
近红外光谱吸收峰重叠现象严重, 信号吸收较弱, 且背景干扰严重, 波长选择的目的就是剔除干扰和无用的变量, 并减少数据冗余和多重共线性, 提高模型的预测精度和预测效率。 波长点选择方法的特点是将每个波长点视为一个单元, 因此点选择的变量是离散的。 本研究采用竞争自适应加权采样法(CARS)进行点变量选择[12]。 竞争性自适应重加权算法(CARS)是每次通过自适应加权采样, 保留PLS模型中回归系数绝对值权重较大的点作为新的子集, 去掉权值较小的点, 再新的子集建立PLS模型, 经过多次计算, 选择PLS模型交互验证均方根误差最小子集中的波长作为特征波长。
如图4所示,CARS的波长选择基于50次迭代和十折交叉验证。 以CARS选择过程的图形表示, 图4(a)表示蒙特卡洛采样运行次数和选择率, 由于衰减指数函数的作用, 随着采样次数的增加, 在采样前期变量数快速减少, 随着采样次数的增加, 变量减少的速度减慢, 表明算法“粗选”和“精选”的过程。 图4(b)表示蒙特卡洛采样运行次数和交叉验证均方根误差, 采样次数为22次时RMSECV最小为0.123 7%。 图4(c)表示蒙特卡洛采样运行过程中各个波长的回归系数趋势。 每个波长由彩色线表示, 星号垂线表示在22次采样运行后确定的交叉验证均方误差根最小时对应的波长子集共94个波长纳入PLS模型, 占原始光谱信息的5.53%。
图3 90批次阿里红饮片的近红外光谱
图4竞争自适应重加权抽样(CARS)变量选择
22次采样运行后识别的具有最低交叉验证均方根误差(RMSE)的子集
2.2.5模型评价
使用一种普遍接受的十折交叉验证方法, 通过式(1)获得的最小预测残差平方和(PRESS)点确定模型的潜在变量数量。PLS模型评价一般由确定系数R2和均方根误差(RMSE*)确定。R2由式(2)获得,R2越接近于1,表示因变量的变化是由大部分回归线描述的。 如果直线的平方误差很大,R2越接近于0,表示因变量的变化很少由回归线描述。RMSE*从式(3)获得,RMSE*表示回归线未描述的总变化的百分比, 如果直线的平方误差很小, 意味着直线非常适合。 评价原则一般认为R2越接近1, RMSE*越接近0,模型效果越好。 性能偏差比(ratio of performance to deviation, RPD)由式(4)计算得到, 表示模型预测性能。 一般RPD越大, 模型的预测能力越强[13]。RPD的阈值选择取决于多种因素, 如背景、 目的和行业。 学者们更倾向于将RPD<1.4视为轻微或根本不具有代表性, 将1.4~3.0之间的值视为合理预测, 当RPD>3时, 表示模型预测精度高。 各项评价参数不能孤立参考, 需要综合起来评价。
表7不同前处理和CARS选择波长组合下PLS模型
图5真实含量VS预测含量
由表7所示, 建立不同前处理光谱数据和不同策略选择波长的PLS模型, 并对其建模效果进行评价。 显然在引入CARS变量选择之前, 近红外光谱吸收峰重叠现象严重, 背景干扰严重, 模型效果很差。CARS波长选择的引入剔除了干扰和无用的变量, 大大提高模型的预测精度和预测效率, 与其他模型相比,SNV+CARS模型性能最佳, 其潜变量为45个,Rp2为0.982 3, RMSEP为0.103 7%, RPD为5.34。 图5为该模型真实含量和预测含量对比。t检验表明,SNV+CARS+PLS建模方法和传统HPLC方法在预测阿里红中齿孔酸含量上没有显著差异。
2.2.6决定系数的解释
回归分析主要用于探索因变量和自变量之间的关系。 在线性回归分析中, 决定系数法(也称为R平方)用于预测和解释模型的结果。 与一元回归模型不同, 除了内部验证性能外, 多元回归还包括外部验证性能, 会导致模型选择中出现一个以上的R2。 有些模型的R
c2大于Rp2,这称为过拟合。 相反有些模型的Rp2大于Rc2,称为欠拟合。 理想情况下,Rp2的变化方向与Rc2相同, 但实际上过拟合和欠拟合是常见现象。 本文最优回归模型选择的是Rp2最大的模型。
3、结 论
采用近红外光谱和偏最小二乘相结合的方法建立阿里红饮片中齿孔酸的定量模型, 比较了CARS波长选择前后模型策略的性能。CARS波长选择后的模型有了很大的提高, 具有较高的相关系数(>0.9)。 在众多PLSR模型中,SNV+CARS模型性能最佳, 其潜变量为45个,Rp2为0.982 3, RMSEP为0.103 7%, RPD为5.34。t检验表明,SNV+CARS+PLS建模方法和传统的高效液相色谱法在预测阿里红饮片中齿孔酸含量上没有显著差异, 表明所开发的模型可作为传统方法的替代或补充。 尽管高效液相色谱法在植物药质量控制中的活性成分测定方面领先了一步, 但近红外光谱法可以在不需要任何样品前处理, 在不到5 min的时间内成功地用于测定阿里红饮片中的齿孔酸含量。 多项研究表明近红外光谱与化学计量工具相结合可以有效测定植物药物中的活性成分。 植物药的质量控制一直是一个棘手的问题, 若将红外技术与化学计量学相结合, 建立的定量模型将具备快速、 无损、 绿色和节约成本等优点, 可以为中药质量控制提供更好、 更便捷的方法。 本研究90个样品中的齿孔酸含量变化很大, 从1.27%到4.36%。 尽管R2没有达到0.999,我们开发的模型仍然可以用于筛选分析。
基金资助:国家自然科学基金联合基金项目(U1703235)资助;
文章来源:谢玉玉,陈志慧,侯雪玲,等.近红外光谱结合化学计量学快速测定阿里红饮片中齿孔酸含量[J].光谱学与光谱分析,2024,44(10):2981-2987.
分享:
慢性胃炎是一种胃黏膜的慢性炎症性疾病,发病机制复杂,涉及幽门螺杆菌感染、不良饮食习惯及遗传因素等多重诱因。该病在病理学上可表现为胃黏膜萎缩,并伴随肠上皮化生和异型增生等癌前病变特征,与胃癌的发生呈显著正相关。
2025-08-16养血安神方是由柳州市中医医院名老中医总结多年临床经验,结合广西柳州地区患者的病情和地域特点,经酸枣仁汤加减化裁而来[1],该方由炒酸枣仁、首乌藤、知母、茯神、川芎等共八味药材构成,在不改变酸枣仁汤原方主要组成及清热除烦、养血安神之功效的基础上[2],增加了首乌藤,与酸枣仁共为君药,进一步增强养血安神的功效[3]。
2025-08-07近年来,对土茯苓的需求量日益增加,人工栽培光叶菝葜尚未形成规模化、规范化的种植模式,故其表型常受种植地、栽培措施的影响,造成土茯苓药材质量的不易控制。此外,一些不良商家利用加工炮制后的土茯苓外观特征丢失、传统鉴定方法无法鉴别的特点,用混伪品冒充正品以获利[3]。
2025-07-23苦杏仁为蔷薇科植物山杏(PrunusarmeniacaL.var.ansuMaxim)、东北杏[Prunusmandshurica(Maxim.)Koehne]、西伯利亚杏(Prunus.sibiricaL.)或杏(PrunusarmeniacaL.)的干燥成熟种子,始载于《神农本草经》,具有降气止咳平喘、润肠通便的功效。苦杏仁作为药食两用的传统中药在中药制剂中被广泛使用,因此对其质量控制方面须有严格要求。
2025-06-26枳实为芸香科植物酸橙CitrusaurantiumL.及其栽培变种或甜橙C.sinensisOsbeck的干燥幼果,主产于江西、湖南、湖北、重庆等地,其中江西为枳实的道地产区。枳实味苦、辛、酸,微寒,归脾、胃经,具有破气消积、化痰散痞的功效,多用于治疗积滞内停、痞满胀痛、大便秘结、泻痢后重、胃下垂、子宫脱垂、脱肛等,其主要含有黄酮类、生物碱类等活性成分。
2025-05-20近年来,胶原蛋白特征肽鉴别方法已成功用于鹿角胶等胶类的真伪鉴别,对于鹿角胶的质量控制起到了非常重要的作用。随着研究的深入,尚有一些不足需要改进,例如《中国药典》仅收载的鹿角胶特征肽的定性鉴别方法,含量测定也仅收载了氨基酸的含量测定法,还不能真实反应鹿角胶质量。
2025-05-07现行质量标准没有含量测定项,不利于质量控制,有学者建立了宁神灵颗粒中黄芩苷或甘草酸铵的含量测定方法[3-5],但同时对上述两种成分进行质量控制的方法尚没有文献报道,本文针对两种化合物紫外最大吸收波长与极性等理化性质不同的特点,采用波长切换技术及梯度洗脱技术。
2025-04-23对20批灵芝进行DNA提取、PCR扩增、纯化、测序、输入中药材DNA条形码鉴定系统得出16批灵芝为药典品种赤芝。同时,测定灵芝中的三萜及甾醇成分含量反映其药效成分。将理化分析技术与DNA分子鉴定相结合,全面、客观地评价灵芝的质量,以期为灵芝的标准化生产、质量控制及临床应用提供科学依据。现将研究结果报道如下。
2025-03-28中药鉴定技术课程是中药专业的一门专业核心课,在学校实施中国特色学徒制人才培养的背景下,为了更紧密地贴合合作企业的人才需求,并与企业的特色产品生产实践相衔接,我们对课程进行了针对性的岗位特点开发。遵循学徒制“在岗培养、交互训教”的培养模式,本课程的目标是培养学生对企业生产中成药成分中常见中药饮片的鉴定能力。
2025-03-17秦艽的有效成分以环烯醚萜苷为主,三萜、甾体、黄酮、苯丙素类等是其有效成分[5-6],具有肝胆保护、抗氧化、抗菌、神经保护、治疗关节炎、抗肿瘤等药理活性[7-11]。黄柏的化学成分以生物碱、苯丙素、酚苷类及脂肪酸类为主,具有抗炎、抗菌、降血压、降血糖、免疫抑制、抗类风湿等功效[12-14]。
2025-03-07人气:6027
人气:3748
人气:3693
人气:3576
人气:3352
我要评论
期刊名称:光谱学与光谱分析
期刊人气:2112
主管单位:中国科学技术协会
主办单位:中国光学学会
出版地方:北京
专业分类:科学
国际刊号:1000-0593
国内刊号:11-2200/O4
邮发代号:82-68
创刊时间:1981年
发行周期:月刊
期刊开本:大16开
见刊时间:一年半以上
影响因子:0.735
影响因子:0.645
影响因子:0.874
影响因子:0.385
影响因子:0.111
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!