评价中药色谱指纹谱相似性的新方法向量夹角法(1)

大连化物所评价中药色谱指纹谱相似性的新方法：向量夹角法中国科学院大连化学物理研究所:王龙星大连化物所前言中医药现代化的重要基础是中药材的质量保障各种色谱指纹谱在中药质量控制中发挥着越来越重要的作用。在不知道中药材的主要有效成分或对其缺乏全面了解的情况下用指纹谱整体对药材进行评价是一个有效的、有前途的方法。大连化物所前言目前根据整体谱图进行质量控制的数据处理方法主要有模式识别与人工神经网络的方法，它们给出的答案是一种分类及定性的信息，不直接定量比较两张指纹图谱间的相似度。但在很多情况下我们需要利用这种相似度的计算提供定量的评价信息，为中药质量控制提供依据。基于这种情况我们提出了一种定量直接比较两张指纹谱相似度的新方法：向量夹角法大连化物所向量夹角法的基本原理n个数值组成的行(x1,x2,x3,…xn)称为n维向量简记为大写字母X定义为向量X的模而是向量X与Y的内积向量X及Y的向量夹角余弦按照下式计算向量夹角余弦约接近1说明两向量相似度越高2232221...nxxxxX＝nnyxyxyxYX...2211＝YXYX?=qcos大连化物所向量夹角法的基本原理每个色谱指纹图谱都可以看作一组对应保留时间下的峰高(或峰面积)的数值，可以把这组数值看作多维空间中的向量两个指纹图谱间相似性的问题转化为多维空间的两个向量的相似性问题，利用上式能够定量表征指纹图谱间的相似性大连化物所向量夹角法计算相似度的基本步骤一在统一的样品处理与分析测试条件下获得色谱指纹图谱；二进行谱图积分获得保留时间、峰高等信息三根据分析时的基线噪声水平，将信噪比5的色谱峰除去不作考虑；四用计算机自动进行两个指纹图谱之间的峰匹配，并利用上式计算相似度大连化物所色谱指纹谱图自动峰匹配问题的解决要利用向量夹角公式计算相似度就必须进行两个指纹谱的峰匹配，判断共有峰与特有峰并把对应的数值填入向量的对应位置，但实际的指纹图谱往往有几十、上百个色谱峰，手工进行峰匹配是非常繁琐的。我们使用VB6.0编制了一个小程序，能利用两个色谱图的积分报告自动进行峰匹配，判断共有峰及特有峰，并把对应的数值填入向量的对应位置、计算相似度。大连化物所色谱指纹谱图自动峰匹配问题的解决随着分析仪器技术的进步，分析仪器上保留时间的重复性越来越好，只要严格控制色谱条件就能够保证色谱峰的出峰次序不变并且保留时间差异较小我们的程序是在下列三条原则的基础上进行自动峰匹配：(自动峰匹配的流程图见下图)一.同一分析条件下两张指纹谱对应色谱峰的保留时间差异在合理的范围以内。(对液相色谱保留时间差异小于1分是容易做到的)二.同一分析条件下两张指纹谱出峰次序不变三.在满足一、二的同时寻找使相似度最大(峰高最接近)的峰匹配结果。大连化物所色谱指纹谱图自动峰匹配问题的解决大连化物所应用实例:吴茱萸液相谱图相似度比较我们对11种不同产地与炮制方法的吴茱萸样品(沈阳药科大学提供)进行了液相色谱分析样品处理：称取吴茱萸4.00克加100ml40％乙醇水溶液，回流煎煮2小时，冷却后放入冰箱冷冻4小时,过滤定容至100ml。过0.45μm滤膜待分析。液相分析条件:INTERSILODS3柱(4.6mm×250mm)Alliance2690液相色谱仪，检测波长为280nm。进样量2μL。流动相A为甲醇，流动相B为含0.3％醋酸氨的水溶液。二元梯度分离，流速1.0ml/min，梯度条件如下：0min-5minA:5%,5min-50minA:5％－100%，50min－60minA:100%。大连化物所实例:不同吴茱萸样品相似度比较061218243036424854广西上林汤洗7编湖南怀化湖南新晃甘草炮制广西灌阳四川广元广西百色广西文市镇广西灌阳茶用广西上林甘草炮制广西凌云仰村贵州遵义大连化物所实例:不同吴茱萸样品相似度比较表1：相似度计算结果Table1:ResultsofsimilarityofChromatogramNo.Origin(产地与炮制方法)Similarity1Shanglin,GuangXi(广西上林汤洗7遍)42.30%2Huaihua,HuNan(湖南怀化)78.60%3Xinhuang,HuNan(湖南新晃甘草炮制)84.90%4Guanyang,GuangXi(广西灌阳县)85.20%5Guangyuan,SiChuan(四川广元)85.20%6Baise,GuangXi(广西百色)85.80%7Wenshizhen,GuangXi(广西文市镇)87.50%8Guanyang,GuangXi(广西灌阳茶用)87.90%9Shanglin,GuangXi(广西上林甘草炮制)89.70%10Lingyun,GuangXi(广西凌云仰村)92.10%11Zuni,GuiZhou(贵州遵义)100%大连化物所实例:不同吴茱萸样品相似度比较据文献报道贵州遵义产吴茱萸是一种道地吴茱萸品种，以此为基准比较、计算了其它样品和它的相似度结果如表1。所有吴茱萸样品除汤洗7遍吴茱萸外都比较相似，相似度都在78％以上，而汤洗7遍吴茱萸的相似性只有42.3%，这说明该方法能够很好的区分汤洗7遍这种炮制方法。从图1的直观比较可以看出，自上而下谱图与贵州遵义样品越来越相似，而相似度的计算结果正好越来越高，从42.3%上升到92.1%。本程序的相似度计算结果与直观的比较是一致的，这说明该计算方法是合理的大连化物所实例:不同吴茱萸样品相似度比较人工比较时通常只比较大一点的色谱峰而忽略小峰，而计算机比较则可以把小色谱峰也考虑在内，这样的计算结果能反应谱图的细微差别以广西灌阳样品与四川广元样品的相似度计算为例，说明小峰的考虑与否对相似度计算的影响。大连化物所实例:不同吴茱萸样品相似度比较061218243036424854如果只考虑信噪比30的大色谱峰相似度＝95％如果考虑信噪比5的所有色谱峰相似度＝86％计算机自动峰匹配比人工比较更能反应指纹谱的细微差别广西灌阳四川广元大连化物所总结与展望向量夹角余弦法数学概念简单、数据处理过程简便快速，在定量计算两张色谱指纹谱相似度方面有其优势。一般来说不同品种药材其指纹图谱差异较大而相同品种的药材指纹图谱差异较小因此本方法也可以用作区分不同的药材。如果利用多种道地药材的标准提取物建立了标准谱图，那么通过待测样品与标准谱图的相似度比较有可能判断待测样品究竟是哪一种道地药材，其相似度如何，这也是我们今后进一步研究的目标