在中药注射剂生产过程中，由于原料或辅料批次更换，以及工艺参数波动等原因，产品质量会有不同程度的波动[1-2]。由于中药化学成分繁多、结构类型高度复杂，在治疗疾病时具有多成分、多靶标的特性，其功效取决于多种成分的共同作用[3]，况且中药注射剂给药方式的特殊性，因此，中药注射剂质量的全面表征与一致性评价是中药质量控制标准建立过程中的重点和难点，已成为制约其发展的瓶颈之一[4]。

目前，中药注射剂质量的批间一致性评价通常采用包括HPLC、GC等在内的色谱技术或色谱联用技术的定量分析及指纹图谱分析法，通过定量、半定量或定性分析[5-6]，实现化学成分的表征与量化，在此基础上进行质量一致性评价。其中中药指纹图谱是中药多成分质量综合分析及评价的重要手段[7-8]。在现阶段中药的有效成分大多数没有明确的情况下，中药指纹图谱能够有效地从整体上表征中药的整体质量，是基于中药“多成分、多靶点”的作用特性建立的一种有效的中药质量控制方法[8]。

氢核磁共振（proton nuclear magnetic resonance，1H-NMR）指纹图谱技术具有单一性、全面性、定量性和易分辨性等特点，适用于复杂混合物中各类有机物的定性定量分析[9-11]。1H-NMR指纹图谱技术的显著优势在于分析范围广，可同时对初级代谢产物（氨基酸、糖类等）与次级代谢产物（黄酮、酚酸、皂苷等药用活性成分）进行高通量分析。目前基于1H-NMR指纹图谱的中药材及其制剂的质量评价方法，一般通过对分析对象的1H-NMR全光谱进行相似度与空间距离计算[12-18]，或在1H-NMR光谱中选择部分代谢物的特征信号并计算相似度[19]，以实现一致性评价的目的。然而这些计算方法未能充分利用1H-NMR光谱中反映的样品化学信息。

本研究以丹参注射液为研究对象，考察其1H-NMR光谱中各类初生及次生代谢产物的特征信号及信号分布规律，研究基于丹参注射液1H-NMR光谱特征指纹图谱的产品质量一致性评价方法。

方法与结果

1、 样本制备与核磁分析

准确量取900 μL丹参注射液与100 μL含0.05% TSP的氘代水，混合均匀并在10 000 r/min下离心10 min，移取500 μL上清液置入5 mm标准核磁管，留待分析。

在预饱和水峰压制的Noesygppr1d的脉冲序列下进行核磁分析，以溶剂90% H2O＋10% D2O锁场，谱宽为δ 12，中心频率为δ 4.695，弛豫时间为2 s，混合时间为50 ms，90°脉冲P1为14.75 μs，空扫次数为4，扫描次数为32，采集时间为2.726 s，接收器增益为114。

在Bruker Topspin（3.6.2版，美国Bruker Biospin Corporation）软件中，对自由感应衰减（free induction decay，FID）信号进行傅里叶变换，窗函数线宽＝0.3 Hz，手动执行相位校正并以TSP的特征信号作为内部参考，将其化学位移设为δ 0.00，作为化学位移漂移矫正的基准，确保不同样品光谱之间具有可比性。

2、 光谱预处理

使用MestReNova（14.0.0版，西班牙Mestrelab Research）软件对光谱在“Peak”模式下进行分箱（Binning）宽度为δ 0.02的降维处理，处理后的光谱作为相似度计算的参考数据。分箱后，去除残留水信号影响的光谱区域（δ 4.5～5.0）及无信号的区域（δ 0.0～0.5），原始光谱由32 768个数据点降维至513个数据点。在此基础上，在Matlab中采用区间相关优化平移算法（interval correlation optimised shifting algorithm，ICOSHIFT）[20-21]工具包对齐NMR光谱信号，消除不同批次丹参注射液NMR图谱间的化学位移偏差。

3 、光谱分区

中药提取物的1H-NMR光谱按提取物中各类物质化学性质的不同，可被划分为低场区（主要为芳香族化合物信号，δ 5.50～10.00）、中场区（主要为碳水化合物区信号，δ 3.30～5.50）以及高场区（主要为小分子有机酸和氨基酸信号，δ 0.05～3.30）。为充分利用光谱中的信息，并全面表征样品的化学信息，本研究将丹参注射液的1H-NMR光谱按高、中、低场进行分区计算，标记为1、2、3区（M1、M2、M3），分别反映注射液中的氨基酸、糖类及丹参酚酸类的成分信息，如图1所示。

4、 相似度计算

4.1 相似度计算方法

本研究使用标准化欧氏距离与余弦相似度作为丹参注射液1H-NMR光谱指纹图谱的相似性计算方法[22]。相似度计算由Matlab软件（版本2019a，美国Mathworks公司）完成。为使标准化欧式距离这一指标与余弦相识度量纲一致，便于后续统一考察，本研究使用“相对欧氏距离”。样品n的1H-NMR光谱的相对欧氏距离rDE的计算公式如下。

rDEn＝DEn/DEmax

相似，相对欧氏距离越接近0说明该分区的相似度越接近于标准样品。经分析，除批次13、16，不在效期内的丹参注射液1区余弦相似度均大于0.970，效期内批次的1区余弦相似度则均大于0.985，说明余弦相似度对氨基酸区的波动敏感性不足，但仍能体现出不在效期内注射液中的氨基酸、小分子有机酸含量与近期批次有一定程度的差异；与此形成对比的是欧氏距离对样品变异具有更高的敏感性，效期外丹参注射液的1区、2区欧氏距离较大且批间波动也较大，批次20～62的3区欧式距离则明显高于标准批次。

2区余弦相似度这一指标则体现出了不在效期内的丹参注射液批次糖区光谱的无序波动，在效期内批次31～81中则均高于0.97。而欧式距离的总体趋势则随生产日期与基准批次的接近而逐渐减小：该指标在2018年及之前批次中较高，在2019年3月至6月批次中则显著降低，在0.2～0.5波动；在2019年7月至2020年1月批次中则在0.08～0.23波动，在20年3月及之后的批次中则在0.1～0.4波动。由于糖类在贮存过程中相对稳定，可以推测样品在出厂时糖类含量即不相同，而氨基酸及酚酸类物质可能在贮存过程中发生了较大变化。

3区的余弦相似度在各批次注射液中的波动更为显著，其中批次20～56的余弦相似度均低于0.75，相应的，其欧式距离则均较大，说明酚酸含量也有一定程度的批间波动。

综上，相较余弦相似度，欧式距离更能体现批次间的差异；而3区的余弦相似度对丹参注射液中酚酸的批间变异有较好的考察能力。此外，评价结果表明，丹参注射液中糖类与氨基酸类成分组成的批间差异较为显著，对注射液质量一致性的影响较大，为提高丹参注射液的用药安全性，对这2类成分应提出控制标准。

4.2 分区相似度评价方法结果的准确性验证

（1）与丹参注射液1H-NMR光谱的全谱相似度评价结果对比：计算81个样品的全1H-NMR光谱（δ 0.05～10.0）的余弦相似度与欧氏距离。全谱相似度的雷达图如图3所示。由图3可知，基于全光谱的余弦相似度与欧氏距离仅能识别部分不在效期内的注射液的异常；同时，图2与图1中2区余弦相似度及2区欧氏距离的评价结果高度雷同，这是因为糖区在丹参注射液这一品种的氢谱全谱中所占峰面积比例极大，因此基于全谱的相似度分析受糖区特征峰影响极大，对酚酸区及氨基酸区的差异识别能力因此被削弱。因此，对于丹参注射液，基于全谱的余弦相似度及欧式距离相似度评价对异常样本的识别能力不如分区相似度评价方法。

（2）基于主成分分析（principal component analysis，PCA）的一致性评价结果对比：进一步对比81个注射液的全光谱数据（XA＝81×513）及分区相似度数据（XA＝81×6）的无监督PCA模型，建模前使用中心化法对所有数据进行缩放处理。81个注射液的全光谱PCA模型包含5个主成分，其R2为0.901，Q2为0.767，前3个主成分分别解释44.6%、27.5%、10.5%的模型方差；81个注射液的分区相似度数据的PCA模型包含3个主成分，其R2为0.998，Q2为0.823，3个主成分分别解释91.5%、5.6%、2.8%的方差。建模结果显示分区相似度的PCA模型的分类性能优于全光谱PCA模型。在此基础上，为使分类结果更为明确，基于PCA结果绘制基于Ward’s最小方差法的层次聚类（hierarchical clustering analysis，HCA）树。聚类结果见图4。

由聚类图4-A可以看出，PCA-HCA无监督判别模型对基于全光谱信息的样品具有一定的聚类能力，但对各批号注射液的分离能力不高，集群间的分类界限不清晰，非相近日期生产的注射液批次（如20年3月批次与19年3月至19年6月批次）也有被聚类到同一集群（图4-A IV组）的风险。对比图4-A与4-B，显然基于分区相似度信息的模型聚类效果更好，相近批次之间的空间距离更近，类内相似度更高而类间差异更为显著。

（3）基于正交偏最小二乘判别分析（orthogonal partial least square-discriminant analysis，OPLS-DA）的批间光谱差异辨识 在“2.4.3（2）”项中对效期内样品的聚类结果，为进一步验证方法对注射液批间差异辨识的准确性，将效期内体现出差异的批次分为2类（图4中III组与IV组），对全光谱的分箱积分数据进行帕累托缩放处理后，进行OPLS- DA，以观察2类注射液中的指标差异的明确来源。OPLS-DA模型是一种有监督的判别模型，较PCA模型有更强的识别类间差异来源的能力。OPLS-DA模型相应的S线图和置换检验图如图5所示。其中，S线图用于识别不同类别注射液间的化学组成差异，置换检验图用于评估模型是否过拟合。

图5-D～F表明3个模型均未过拟合，模型性能良好。S线图的Pctr [1]反映了X变量与t1之间的协方差，而Pcorr值代表了X变量与t1之间的相关性，Pcorr的绝对值越大，表示该点对应的化合物对模型的影响越大；Pcorr的正相关系数和负相关系数分别表示2类样本中相应化合物浓度的正相关和负相关。由图5-A～C中S线图的Pcorr值可知，酚酸区的不同是2类样品的主要差异来源，这一结论与图2中2类样本的分区相似度统计结论（2类样本的3区余弦相似度与欧氏距离均有差异）完全一致。

为进一步验证这一结果的准确性，对注射液中酚酸含量进行基于HPLC-UV的定量分析[23]，并以HPLC法所得酚酸含量的定量结果为参照，考察2类注射液的酚酸含量差异。6个丹参酚酸在2类注射液中含量差异的双样本t检验结果分别为P＝0.000、0.160、0.640、0.000、0.000、0.004。6种丹参酚酸含量分布的箱型图见图6所示。由相应t检验结果可知，注射液中丹参素、紫草酸、丹酚酸B、丹酚酸A的质量浓度在III组中有显著的提升，说明2类样本在酚酸组成上存在显著差异，但是基于1H-NMR全谱的指纹图谱一致性分析（图3）并未体现出这种差异。

讨论

本研究将丹参注射液的1H-NMR光谱按代谢物分布规律进行了分区，并以余弦相似度及标准化欧氏距离为评价指标，分别计算氨基酸区、糖区及酚酸区的相似度，作为多批次丹参注射液的质量一致性评价方法，评价了同一厂家在2016年至2020年间生产的81批丹参注射液中酚酸、氨基酸及糖类成分的组成情况，结果显示，生产时间间隔较久的产品在氨基酸及糖类成分的含量与组成上存在差异；此外，随保存时间变长，丹参酚酸也会产生一定程度的降解导致产品质量波动。在此基础上，通过多角度的对比证明了所提出的1H-NMR光谱分区相似度评价方法在丹参注射液的质量一致性评价中相较传统的全光谱相似度计算法具有更强的批间化学差异识别能力。

综上所述，本研究建立的丹参注射液的1H-NMR光谱相似度评价方法可以反映出丹参注射液批间不同类别成分的波动，可用于丹参注射液的多组分质量控制及评价，同时为中药注射液的一致性评价提供了新的思路与参考方法。

利益冲突 所有作者均声明不存在利益冲突

参考文献（略）

来 源：赵 芳，李文竹，潘坚扬，杨嘉誉，瞿海斌．基于1H-NMR指纹图谱分区相似度的丹参注射剂质量评价方法 [J]. 中草药, 2022, 53(17): 5312- 5320.