中药色谱指纹图谱相似度评价系统

中药色谱指纹图谱相似度评价系统国家药典委员会国家药典委员会翟为民背景背景国家食品药品监督管理局于2000年先后下达“关于加国家食品药品监督管理局于2000年先后下达“关于加强中药注册管理有关事宜的通知”（国药管注［2000］157号）和“关于印发《中药注射剂色谱指纹图谱研究的技术要号）和关于印发《中药注射剂色谱指纹图谱研究的技术要求》的通知”（国药管注［2000］348号）文件。国家药典委员会负责组织实施“已批准生产的中药注射剂质量标准的提高完善及指纹图谱工作”及牵头组织“已批准多家生产的提高完善及指纹图谱工作”及牵头组织“已批准多家生产的中药注射剂，其指纹图谱联合起草工作”。国家药典委员会从2000年12月起开始组织实施，进行国家药典委员会从2000年12月起开始组织实施，进行了可行性调研、规范性技术文件起草与研究、同布进行了科技部立项。从中药注射剂的指纹图谱研究质量标准提高逐步从中药注射剂的指纹图谱研究、质量标准提高，逐步成熟、扩展，到2010版药典，部分提取物和中成药标准中都收载了指纹图谱项目，为发挥着积极的作用。收载了指纹图谱项目，为发挥着积极的作用。1论证1、论证首先进行调研，多次召开研讨会，代表涉及科研院所、中央和省级药品检验所及企业科技人员形成共识：央和省级药品检验所及企业科技人员。形成共识：①重点-色谱指纹图谱；②理论上可行②理论上可行。③通过实践，检验可行性，不断完善。谱2、修订《中药注射剂指纹图谱研究技术要求》（国药管注［2000］348号）①①优化、方法验证，细化参数、明确要求。②解决手工计算“共有峰”、“非共有峰”面积工作中大量的人为因素改变“只见树木不见森林”的局面力求量的人为因素，改变“只见树木不见森林”的局面，力求反映整体质量特性。二规范性技术文件起草与研究工作二、规范性技术文件起草与研究工作1、修订修订稿就原料药材、半成品、成品的样品收集、制备、参照物的选择，指纹图谱试验条件，对实验室的要求，对仪器的要求，试验方法和试验条件的建立及验证，色谱条件的优选，指纹图谱的建立和验证识别计算和复核标准书写格式指纹图谱的建立和验证、识别、计算和复核，标准书写格式等进行了详尽的细化规定，起草稿经广泛征求意见，先后进行了多次修改得以完成。以《中药注射剂色谱指纹图谱试验行了多次修改得以完成。以《中药注射剂色谱指纹图谱试验研究操作规程指南》（试行）的名义发布，作为局颁“技术要求”的补充。2、计算机辅助中药指纹图谱相似度计算软件的研究开发发对于成药而言，如何控制批与批之间质量的稳定，这就要求需采用某种手段能定量地表征各样品间的相似程度(或差异性)而不是真伪鉴别，从整体性考虑，采用指纹图谱，在批间稳定等方面能够发挥积极作用。也就需要建立套能整体反应各样品指纹图谱之间相似程度的计也就需要建立一套能整体反应各样品指纹图谱之间相似程度的计算方法，但原版“技术要求”采用共有峰和非共峰之比值的人工计算方法，该方法仅考虑了人为选定的几个特征峰的相似性，即将一个完整的事物割裂开进行比对，是不全面的，未能将中药内在成分的整体特性表达出来，因此，要创建一套能定量比对整体上反映中药内含成分指纹图谱相似程度的工具上反映中药内含成分指纹图谱相似程度的工具。国家药典委员会组织沈阳药科大学浙江•国家药典委员会组织沈阳药科大学、浙江大学、中南大学、清华大学、北京大学医学部药科技发中北大学中国学部医药科技开发中心、西北大学、中国药品生物制品检定所、第二军医大学八家单位，于2004年集中测试修改定型《计算机辅助中药指纹图谱相似度计算软件》。•软件的开发过程与中药注射剂指纹图谱的科研工作一直保持了同步。这种同步的应科研工作直保持了同步。这种同步的应用也给软件的测试使用创造了条件。7开发初期设计思路是将色谱数据导入•开发初期设计思路是将色谱数据导入MATLAB的专业数据处理的软件中进行相似度的计算。后考虑到知识产权保护的问题，采用自主编程，将指纹图谱数据导入，数据库管理，计程，将指纹图谱数据导入，数据库管理，计算方法选择，谱图处理等功能综合考虑，设计开发而成计开发而成。8软件基本算法与版本软件基本算法与版本一、算法选择简介二、修订改版工作二、修订改版工作三、2.0版的介绍•一、首要问题是选择算法•课题组同步进行了《指南》的制定与算法课题组同步进行了《指南》的制定与算法的评价与选择，通过大量计量化学方面的研讨与对比验证选择了夹角余弦法研讨与对比验证，选择了夹角余弦法。评价方法比较与算法选择简介中药指纹图谱是表征中药复杂成分与其质量关系的重要手段中药指纹图谱是表征中药复杂成分与其质量关系的重要手段。。在表征过在表征过程中程中，，相似度评价起着关键作用相似度评价起着关键作用。。目前常用的相似度评价方法主要有峰重叠率目前常用的相似度评价方法主要有峰重叠率评价方法比较与算法选择简介法法（（NeiNei系数法系数法））,,相关系数法相关系数法,,距离系数法距离系数法,,向量夹角余弦法与峰重叠率和共向量夹角余弦法与峰重叠率和共有峰强度结合法有峰强度结合法（（改进的改进的NeiNei系数法系数法））。。11指纹图谱波动对各种评价方法的影响指纹图谱波动对各种评价方法的影响::0.951.00ABCD0750.800.850.90095rityDE0550.600.650.700.75Simila0.000.050.100.150.200.450.500.55FlttiItitFluctuatingIntensityA.夹角余弦法夹角余弦法B.相关系数法相关系数法C.距离系数法距离系数法D.NeiNei系数法系数法E.改进的改进的NeiNei系数法系数法由于药材的产地来源储藏条件和加工工艺的不同以及仪器设备的稳定性等原结果表明：A方法受指纹图谱波动的影响最小12由于药材的产地、来源、储藏条件和加工工艺的不同，以及仪器设备的稳定性等原因使得指纹图谱产生一定的波动。举例：向痛必定粉针制剂高效液相色谱指纹图谱中随机地加入一定强度的波动，以考查波动对相似度评价方法的影响。结果如上图。小峰的缺失的影响：小峰的缺失的影响：小峰的缺失的影响：小峰的缺失的影响：091.0ABC0.70.80.9CDESimilarity040.50.6S0.30.40~5%0~3%The scope of lost peaks in fingerprint00~1%A.夹角余弦法夹角余弦法B.相关系数法相关系数法C.距离系数法距离系数法D.NeiNei系数法系数法E.改进的改进的NeiNei系数法系数法在药材加工和提取过程中，药材中的某些成分往往会有一定的损失，导致结果表明：A方法受小峰缺失的影响最小13在药材加工和提取过程中，药材中的某些成分往往会有定的损失，导致指纹图谱出现峰的缺失。假设峰的缺失主要为峰面积较小的峰，以5％作为阈值，将缺失峰分成0～1％，0～3％和0～5％三组进行考察。大峰的缺失的影响大峰的缺失的影响1.0ABC大峰的缺失的影响：大峰的缺失的影响：0.80.9larityDE0.50.60.7Simil01230.40.5NumberoflargerpeakslostinfingerprintNumberoflargerpeakslostinfingerprintA.夹角余弦法夹角余弦法B.相关系数法相关系数法C.距离系数法距离系数法D.NeiNei系数法系数法E.改进的改进的NeiNei系数法系数法结果表明A方法受大峰缺失的反映较“灵敏”结果表明：A方法受大峰缺失的反映较“灵敏”14在药材加工和提取过程中，指纹图谱中大峰的缺失并不常见。而对于同属不同种的药材，大峰的差异是有可能的。考察了大峰缺失的影响举例。数据标准化处理的影响数据标准化处理的影响：：数据标准化处理的影响数据标准化处理的影响：：数据标准化处理一般会影响各种方法的相似度评价结果。若处理前后样本的相似度评价结果的大小（等级）顺序变化较小，说明标准化处理对样本的亲疏远近关系影响较小反之影响较大采用Spearman相关系数进行检验Influenceofdataautoscalingondifferentmethods较小，反之影响较大。采用Spearman相关系数进行检验。距离系数相关系数夹角余弦gforevaluatingsimilarity距离系数相关系数夹角余弦Spearman系数0.82710.927(=0.01,2-tailed,rs=0.735)0.8420.7950.79015数据标准化对表中三种评价方法所得结果影响较小关于对照模式的计算关于对照模式的计算•平均矢量（重心）=(xx…x)/m•平均矢量（重心）=j(x1j,x2j,…,xnj)/m•中位数矢量=median((x1j,x2j,…,xnj))两种对照模式两种对照模式•批样品中不存在离群(outlier)样本时，一般推荐使用平均矢量（重心）作为对照模式矢量•样品中存在离群样本时，推荐使用具有稳健(robust)性质的中位数矢量来作为共有模式矢量，这种情况可能在计算原料药材样本时的共有模式出现出现。•采用具有稳健性质的中位数矢量来作为对照模式，就不会因为原料药材样本中出现了个离群样本就不会因为原料药材样本中出现了一个离群样本时而影响整个药材的对照模式。值得指出的是如果不存在离群样本时两种方•值得指出的是，如果不存在离群样本时，两种方法所得结果应该基本一致。一、2004年版软件与算法简介•版本:•研究版（2004A）•研究版（2004A）•检验版（2004B）研究版（2004A）检验版（2004B）应用范围科学研究指纹图谱的检验工作功能具有生成对照图谱功能不具有生成对照图谱功能结论研究版（2004A）的功能较多，且函盖了检验版（2004B）18结论研究版（2004A）的功能较多，且函盖了检验版（2004B）具有的功能•2004版以后，根据使用情况，验证并确认保持核心算法不变。保持核算法不变•根据使用中收集的反馈意见，补充必要的需求按照软件编制规范彻底梳理并重需求，按照软件编制规范，彻底梳理并重新编写了全部代码，在人机交互方面进一步完善并且为了适应管需要增加了相步完善，并且为了适应管理需要增加了相应的功能。进行了不同范围的测试，形成了2.0版。业已获得软件著作权证书。20版介绍2.0版介绍运行环境•运行环境•1、硬件：硬件•最低配置：CPU主频533MHz以上，内存64M，硬盘剩余最小空间100M屏幕分辨率硬盘剩余最小空间100M，屏幕分辨率800×600。2软件环境•2、软件环境：•支持的操作系统：windows98，windows2000、windows2003、windowsXP。windows7（32位）（64位）•下载：fpchporgcnfp.chp.org.cn•安装：•操作流程•生成对照图谱生成对照图谱•检验主界面主界面•文件列表区：显示打开的图谱文件列表。•文件详细信息区：显示选中文件信息。文件详细信息区：显示选中文件信息。•绘图区：绘制图谱，标注色谱峰。•状态栏：显示工作状态、当前用户、登录时间等信息。时间等信息。•快捷工具栏：常用按钮：•鼠标有提示鼠标有提示•工作状态切换•注意提示：注意提示：•全谱导入（单、多个）：多个*.cdf*.txt*SCP*XSCP*.SCP*.XSCP•可打开三种格式的报表：式的报表*.rdf*.qrp*.XQRP.XQRP导出文件（图谱）图谱*.txt*S*.Scp导出图像*bmp*.bmp*.wmf参照图谱信息删除选定的图谱删除所有图谱•删除选定的图谱（可多选）删除单个或多个图•删除所有图谱•多一项确认，再•删除单个或多个图谱。•也可以选定列表中多项确认，再清除。•也可以选定列表中的文件，直接按DEL键。DEL键。视图菜单视图菜单图谱间距调整图谱间距调整图谱间距•放大图谱、缩小图谱、原始大小•快捷工具栏或者菜单中选择、鼠标快捷工具栏或者菜单中选择、鼠标单谱显示多谱显示多谱显示隐藏坐标轴显示或隐藏色谱峰菜单、右键，显示或隐藏。3.3.6 隐藏坐标轴从菜单、右键，隐藏或显示。随鼠标显示峰