仿制药冒充原研药的一种检测方法

摘要

本发明提供了仿制药冒充原研药的一种检测方法，包括五个步骤：1、获得待测样品的拉曼光谱数据；2、依据原研药的标准拉曼光谱数据及待测样品拉曼光谱数据计算待测样品的HQI值，根据HQI值对待测样品进行首次判定；3、选出待测样品特征峰和原研药特征峰；4、计算待测样品特征峰的个数与原研药特征峰的个数的比值t，根据t对待测样品进行二次判定；5、计算匹配特征峰的个数和原研药特征峰的个数的比值s，根据s对待测样品进行三次判定。本发明的技术方案准确性高，合理性强，显著降低判别的假阳性率，并且能够满足现场快检的需要。

说明书

技术领域

本发明涉及仿制药冒充原研药的一种检测方法，属于医药分析技术领域。

背景技术

仿制药与原研药通常含有相同的有效成分(API)，有的甚至所采用的辅料 也相同。为了谋取非法暴利，一些不法商家采用价格低廉的仿制药来冒充原研 药进行销售，此举严重侵害了研发者的知识产权。并且，由于API和辅料的来 源不同等原因，原研药和仿制药之间也会存在一定的质量差异，消费者使用该 类假药还可能产生严重的用药安全问题。

由于仿制药与原研药含有相同API和辅料，且规格相同，当采用药品法定 标准对该类假药进行检测时常得到假阳性结果，即将仿制药错判为原研药，为 药检工作人员对该类假药的检测和判定带来极大的挑战。

由于拉曼光谱特征性强，几乎所有的物质都有其特征性的拉曼光谱峰，因 此一些分析工作者尝试采用拉曼光谱技术对是否采用仿制药冒充原研药进行分 析检测。

目前采用拉曼光谱对仿制药冒充原研药进行快检的方法主要是：计算待测 样品的拉曼光谱与其API对照品的拉曼光谱(事先存储于本地数据库中)的相 似度(夹角余弦、欧氏距离等)，然后将相似度与事先设定的阈值(通常设定为 0.95)进行比较，从而完成对待测样品的判别。但该方法灵敏度较低，当检测一 些和原研药相似度较高的仿制药(如相似度大于0.95)时，由于仿制药与原研 药光谱的差异仅体现在局部的一些位置或者吸收峰上，这些差异对于基于全谱 的相似度数值上的改变较小，该方法难以适用。

CN201410081628.9公开了仿制药冒充原研药的一种检测方法，该方法是将 待检药品的特征峰与原研药的特征峰进行峰匹配，根据原研药特征峰的个数和 匹配特征峰的个数的差值，判定待检药品时原研药还是仿制药。但由于药品的 组分复杂，API、辅料及杂质的光谱峰可能会重叠性，待检药品光谱中能与原研 药品特征峰匹配的光谱峰也可能是待检药品中其他组分产生的光谱峰，因此会 导致判断结果出现假阳性。并且，该方法是根据峰位筛选出匹配特征峰，当待 检药品的特征峰因为其他物质的影响变宽或出现肩峰时，就对匹配特征峰的筛 选产生影响，从而导致检测结果不准确。

发明内容

本发明的目的是提供一种仿制药冒充原研药的检测方法。

为了实现上述目的，本发明所采用的技术方案是：

仿制药冒充原研药的一种检测方法，待测样品的拉曼光谱数据由拉曼光谱 仪检测得到，所得的拉曼光谱数据采用计算机进行处理和计算，该计算机中安 装有数据预处理模块和HQI值计算模块，其特征在于，包括以下五个步骤：

步骤一：待测样品的拉曼光谱数据的获得

采用拉曼光谱仪测得待测样品的原始拉曼光谱数据，采用计算机对原始拉 曼光谱数据进行预处理，得到待测样品的拉曼光谱数据，该预处理包括特征波 段的选取、对特征波段进行9点平滑、对拉曼光谱数据进行基线校正、和对校 正后的拉曼光谱数据进行矢量归一化处理；

步骤二：对待测样品进行首次判定

依据原研药的标准拉曼光谱数据及待测样品拉曼光谱数据，采用计算机计 算待测样品的HQI值，当HQI值小于k，则判定该待测样品为仿制药，当HQI 值大于等于k，则转入步骤三，k为0.90-0.98之间的任意数值；

步骤三：选峰处理

选出待测样品的拉曼光谱数据中强度大于最强峰强度的p％的散射峰作为 待测样品特征峰，p为1～10之间的任意数值，选出原研药的标准拉曼光谱数据 中强度大于最强峰强度的p％的散射峰作为原研药特征峰；

步骤四：对待测样品进行二次判定

计算待测样品特征峰的个数与原研药特征峰的个数的比值t，当t不为 80％-120％之间的数值，则判定该待测样品为仿制药，当t为80％-120％之间的 数值，转入步骤五进行三次判定；

步骤五：对待测样品进行三次判定

筛选出待测样品特征峰中，峰位和原研药特征峰的峰位波数相差-w～w的待 测样品特征峰作为匹配特征峰，w为2cm^-1～4cm^-1之间的任意数值，

计算匹配特征峰的个数和原研药特征峰的个数的比值s，当s为90％-100％ 之间的任意数值，则判定为该药品为原研药，反之，则判定为该药品为仿制药。

本发明的技术方案的进一步特征在于：k为0.95。

本发明的技术方案的进一步特征在于：p为5。

本发明的技术方案的进一步特征在于：步骤四中，t不为95％-100％之间的 数值时，判定该待测样品为仿制药，当t为95％-100％之间的数值时，转入步骤 五进行三次判定。

本发明的技术方案的进一步特征在于：w为3cm^-1。

本发明的技术方案的进一步特征在于：当s为98％-100％之间的任意数值时， 判定为该药品为原研药，反之，则判定为该药品为仿制药。

与背景技术相比，本发明所提供的技术方案将标准拉曼光谱数据作为建立 判别分析方法的依据，从原研药的标准拉曼光谱数据的整体轮廓和特征峰特点 出发，将待测样品与原研药的拉曼光谱进行整体相似度和特征峰个数比较，其 优点和积极效果如下：

1.准确性高

本发明的技术方案将待测样品与原研药的拉曼光谱进行整体相似度比较， 找出不符合首次判定条件的待测样品，再通过二次判定和三次判定，完成对待 测样品的检测，与背景技术相比，对仿制药冒充原研药的检测准确率高。

2.显著降低假阳性率

本发明的技术方案将待测样品的拉曼光谱图与原研药的标准拉曼光谱数据 进行比较，采用的相似度分析方法可以对光谱整体轮廓与原研药光谱差异较大 的样品进行快速、准确地判别。对于相似度极大的样品采用特征峰的分析方法， 能获得较准确的分析结果。两种分析方法相结合既从API的角度出发的同时， 也将其他辅料光谱信息也包括在内，同时整体光谱还考虑到光谱峰强度、及形 状对判别结果的影响，使判定结果的假阳性率显著降低。

3.能够满足现场快检的需要

本发明的技术方案采用拉曼光谱仪即可实现，方法简单，检测快速，能够 在基层一线检验中实现对仿制药冒充原研药的现场快速检测。

4.更具有合理性

本发明的技术方案采用待检药品的特征峰的个数与原研药品的特征峰个数 的比值进行三次判定，由于不同品种药品其原研药特征峰的个数也不尽相同， 从几个至几十个不等，采用百分比的结果进行比较，与背景技术中采用峰个数 差异的绝对值进行比较更具有合理性。

附图说明

图1为本发明的技术方案在实施例一中的流程图；

图2为实施例一中十一个待测样品的拉曼光谱数据图；

图3为实施例一中步骤三的选峰处理的说明图；

图4为实施例二中的选峰处理说明图；以及

图5为实施例三中的选峰处理说明图。

具体实施方式

以下结合附图，对本发明所涉及的仿制药冒充原研药的一种检测方法做进 一步说明。

<实施例一>

一、样品和分析系统

1、样品选取：选取杭州默沙东制药有限公司生产的辛伐他汀片做为本实施 例的原研药；选取杭州默沙东制药有限公司生产的其他批次的辛伐他汀片和其 他十个厂家生产的辛伐他汀片仿制药作为待测样品，十一个待测样品标记为依 次标记为待测样品一(来源：厂家1的仿制药)、待测样品二(来源：厂家2的 仿制药)、待测样品三(来源：默沙东)、待测样品四(来源：厂家4的仿制药)、 待测样品五(来源：厂家5的仿制药)、待测样品六(来源：厂家6的仿制药)、 待测样品七(来源：厂家7的仿制药)、待测样品八(厂家8的仿制药)、待测 样品九(来源：厂家9的仿制药)、待测样品十(来源：厂家10的仿制药)、待 测样品十一(来源：厂家11的仿制药)。

2、分析系统：在本实施例中，采用B&WTek i-Raman-785型拉曼光谱仪获 得相关的拉曼光谱数据，所得的拉曼光谱数据采用计算机进行处理和计算。该 计算机的硬件系统为Dual E2180＠2.00 GHz芯片，该计算机 中安装有数据预处理模块、平均值计算模块和HQI值计算模块，HQI值为快速 匹配系数。

二、待测样品的检测方法

在本实施例中，采用的HQI值计算模块为matlab7.0软件。

在本实施例中，原研药的标准拉曼光谱数据的获得方法为：取原研药，刮 除表面的包膜，采用拉曼光谱仪分别测得三组原研药的原始拉曼光谱数据，采 用数据预处理模块分别对三组原研药的原始拉曼光谱数据进行预处理，包括： 选取拉曼光谱数据中200cm^-1～1800cm^-1的波段作为特征波段，对特征波段采用 Savitzky-Golay算法9点平滑，并用Rubberband算法对拉曼光谱数据进行基线 校正，将校正后的拉曼光谱数据进行矢量归一化处理，得到预处理后的三组拉 曼光谱数据，采用计算机计算预处理后的三组拉曼光谱数据的平均值，得到原 研药的标准拉曼光谱数据。

原研药的标准拉曼光谱数据还可以从数据库中查询得到。

图1为本发明的技术方案在实施例一中的流程图。

图2为实施例一中十一个待测样品的拉曼光谱数据图。

如图1所示，本实施例所提供的仿制药冒充原研药的一种检测方法，包括 以下五个步骤：

步骤一(S1)：待测样品的拉曼光谱数据的获得

取待测样品一，刮除其表面的包膜，采用拉曼光谱仪测得待测样品一的原 始拉曼光谱数据，采用计算机对待测样品一的原始拉曼光谱数据进行预处理， 预处理方法与原研药的原始拉曼光谱数据进行预处理方法相同，得到如图2所 示的待测样品一的拉曼光谱数据。

参照上述方法，分别得到如图2所示的其他十个待测样品的拉曼光谱数据。

步骤二(S2)：对待测样品进行首次判定

依据标准拉曼光谱数据及待测样品拉曼光谱数据，采用计算机计算待测样 品的HQI值。当HQI值小于k，则判定该待测样品为仿制药，当HQI值大于等 于k，则转入步骤三，k为0.90-0.98之间的任意数值。

在本实施例中，k为0.95。采用matlab7.0软件计算得到十一个待测样品的 HQI值如表1所示：

表1十一个待测样品的HQI值

根据表1中每个待测样品的HQI值，对十一个待测样品的进行首次判定， 其中，待测样品四、待测样品六、待测样品七、待测样品八、待测样品十及待 测样品十一的HQI＜0.95，判定这六个待测样品均为仿制药；待测样品一、待测 样品二、待测样品三、待测样品五及待测样品九的HQI≥0.95，转入步骤三。

步骤三(S3)：选峰处理

选出待测样品的拉曼光谱数据中强度大于最强峰强度的p％的散射峰作为 待测样品特征峰，p为1～10之间的任意数值，选出原研药的标准拉曼光谱数据 中强度大于最强峰强度的p％的散射峰作为原研药特征峰。在本实施例中，p为 5。

图3为实施例一中步骤三的选峰处理的说明图。

取待测样品一的拉曼光谱数据，如图3所示，选出强度大于最强峰强度的 5％的散射峰作为待测样品一的特征峰，共得到39个待测样品一的特征峰。参照 上述方法，如图3所示，分别选出了33个待测样品二的特征峰、36个待测样品 三的特征峰、34的待测样品五的特征峰及32个待测样品九的特征峰。

取原研药的标准拉曼光谱数据，选出强度大于最强峰强度的5％的散射峰作 为原研药特征峰，共得到36个原研药特征峰。

步骤四(S4)：对待测样品进行再次判定

计算待测样品特征峰的个数与原研药特征峰的个数的比值t，在本实施例中， 当t不为95％-105％之间的数值，则判定该待测样品为仿制药，当t为95％-105％ 之间的数值，转入步骤五进行三次判定。

在本实施例中，首先计算得到待测样品一的特征峰的个数与原研药特征峰 的个数的比值t，为108％，因此判定该待测样品即待测样品一为仿制药。依据 该方法，分别计算待测样品二、待测样品三、待测样品五及待测样品九的特征 峰的个数与原研药特征峰的个数的比值t，分别为92％、100％、96％及89％，因 此判定待测样品二和待测样品九就为仿制药，将待测样品三和待测样品五转入 步骤五进行三次判定。

步骤五(S5)：对待测样品进行三次判定

筛选出待测样品特征峰中，峰位和原研药特征峰的峰位波数相差-w～w的待 测样品特征峰作为匹配特征峰，w为2cm^-1～4cm^-1之间的任意数值。在本实施例 中，w为3cm^-1，筛选出的待测样品三的匹配特征峰的个数为36个，筛选出的 待测样品五的匹配特征峰的个数为34个。

计算匹配特征峰的个数和原研药特征峰的个数的比值s，当s为98％～100％ 之间的任意数值，则判定为该待测样品为原研药，反之，则判定为该待测样品 为仿制药。

在本实施例中，待测样品三的匹配特征峰的个数与原研药品特征峰的个数 的比值s为1，因此，判定待测样品三为原研药；待测样品五的匹配特征峰的个 数与原研药品特征峰的个数的比值s为94.44％，因此，判定待测样品五为仿制 药。

以上判定结果均与实际相符。

在本实施例中，k为0.95。当k为0.98时，首次判定的判定结果与k为0.95 的首次判定的判定结果相同。当k为0.90时，首次判定的判定结果为：待测样 品四、八、十、和十一为仿制药，待测样品一、二、三、五、六、七、和九为 待测样品，然后，参照步骤三、步骤四及步骤五的操作对七个待测样品进行二 次判定和三次判定即可，最终的判定结果与上述结果相符。

<实施例二>

本实施例使用的分析仪器与实施例一相同。选用的待测样品有五个，分别 为实施例一中的：待测样品一(来源：厂家1的仿制药)、待测样品二(来源： 厂家2的仿制药)、待测样品三(来源：默沙东)、待测样品五(来源：厂家5 的仿制药)、和待测样品九(来源：厂家9的仿制药)。

图4为实施例二中的选峰处理的说明图。

参照实施例一所提供的步骤一至步骤五的检测方法，在本实施例，p为1， 其他与实施例一的方法相同，对五个待测样品的进行步骤三的选峰处理，处理 结果如图4所示，待测样品一共筛选到47个的特征峰，待测样品二共筛选到47 个的特征峰，待测样品三共筛选到48个的特征峰，待测样品五共筛选到47个 的特征峰，待测样品九共筛选到44个的特征峰。

当p为1时，筛选出的原研药的标准拉曼光谱数据中的原研药特征峰共有 48个。

分别计算五个待测样品的特征峰的个数与原研药特征峰的个数的比值t，结 果如表2所示：

表2五个待测样品的特征峰的个数与原研药特征峰的个数的比值

待测样品 待测样品一 待测样品二 待测样品三 待测样品五 待测样品九 t 97.9％ 97.9％ 100％ 97.9％ 91.7％

依据实施例一中所提供的对t的判定方法进行判定，待测样品九为仿制药， 将其他三个待测样品转入步骤五中进行三次判定。

在本实施例中，w为3cm^-1，筛选出的待测样品三的匹配特征峰的个数为48 个，筛选出的待测样品一、二、五的匹配特征峰的个数分别为42、45、46个。

计算匹配特征峰的个数和原研药特征峰的个数的比值s，当s为98％～100％ 之间的任意数值，则判定为该待测样品为原研药，反之，则判定为该待测样品 为仿制药。

在本实施例中，待测样品三的匹配特征峰的个数与原研药品特征峰的个数 的比值s为1，因此，判定待测样品三为原研药；而其他三个待测样品的匹配特 征峰的个数与原研药品特征峰的个数的比值s皆不在98％～100％之间，因此，判 定这三个待测样品都为仿制药。

以上判定结果均与实际相符。

<实施例三>

本实施例使用的分析仪器与实施例一相同。选用的待测样品有五个，分别 为实施例一中的：待测样品一(来源：厂家1的仿制药)、待测样品二(来源： 厂家2的仿制药)、待测样品三(来源：默沙东)、待测样品五(来源：厂家5 的仿制药)、和待测样品九(来源：厂家9的仿制药)。

图5为实施例三中的选峰处理的说明图。

参照实施例一所提供的骤一至步骤五的检测方法，在本实施例，p为10， 其他与实施例一的方法相同，对五个待测样品的进行步骤三的选峰处理，处理 结果如图4所示，待测样品一共筛选到36个的特征峰，待测样品二共筛选到28 个的特征峰，待测样品三共筛选到29个的特征峰，待测样品五共筛选到28个 的特征峰，待测样品九共筛选到24个的特征峰。

当p为10时，筛选出的原研药的标准拉曼光谱数据中的原研药特征峰共有 29个。

分别计算五个待测样品的特征峰的个数与原研药特征峰的个数的比值t，结 果如表3所示：

表3五个待测样品的特征峰的个数与原研药特征峰的个数的比值

待测样品 待测样品一 待测样品二 待测样品三 待测样品五 待测样品九 t 124.1％ 96.6％ 100％ 96.6％ 82.8％

依据实施例一中所提供的对t的判定方法进行判定，待测样品一和待测样品 九为仿制药，将其他三个待测样品转入步骤五中进行三次判定。

在本实施例中，w为3cm^-1，筛选出的待测样品三的匹配特征峰的个数为29 个，筛选出的待测样品二、五的匹配特征峰的个数分别为27、27个。

计算匹配特征峰的个数和原研药特征峰的个数的比值s，当s为98％～100％ 之间的任意数值，则判定为该待测样品为原研药，反之，则判定为该待测样品 为仿制药。

在本实施例中，待测样品三的匹配特征峰的个数与原研药品特征峰的个数 的比值s为1，因此，判定待测样品三为原研药；而其他两个待测样品的匹配特 征峰的个数与原研药品特征峰的个数的比值s皆不在98％～100％之间，因此，判 定这两个待测样品都为仿制药。

以上判定结果均与实际相符。

<实施例四>

本实施例使用的分析仪器与实施例一相同。选用的待测样品有五个，分别 为实施例一中的：待测样品一(来源：厂家1的仿制药)、待测样品二(来源： 厂家2的仿制药)、待测样品三(来源：默沙东)、待测样品五(来源：厂家5 的仿制药)、和待测样品九(来源：厂家9的仿制药)。

参照实施例一所提供的步骤一至步骤五的检测方法，在本实施例，步骤四 中当t不为80％-120％之间的数值，则判定该待测样品为仿制药，当t为80％-120％ 之间的数值，转入步骤五进行三次判定，其他与实施例一的方法相同，对五个 待测样品依次进行首次判定、二次判定和三次判定。

当p＝5时，五个待测样品筛选出的待测样品特征峰的个数如4所示：

表4五个待测样品的特征峰的个数

待测样品 待测样品一 待测样品二 待测样品三 待测样品五 待测样品九 特征峰个数 39 33 36 34 32

分别计算每个待测样品的特征峰个数与原研药特征峰个数的比值t，结果如 表5所示：

表5五个待测样品的特征峰的个数与原研药特征峰的个数的比值

待测样品 待测样品一 待测样品二 待测样品三 待测样品五 待测样品九 t 108.3％ 91.8％ 100％ 94.4％ 88.9％

根据表5的结果，五个待测样品的t值均在80％-120％之间，则转入步骤五 进行三次判定。在本实施例中，w为3cm^-1，筛选出的待测样品三的匹配特征峰 的个数为36个，筛选出的待测样品一、二、五、九的匹配特征峰的个数分别为 33、31、32、28个。

计算匹配特征峰的个数和原研药特征峰的个数的比值s，当s为98％～100％ 之间的任意数值，则判定为该待测样品为原研药，反之，则判定为该待测样品 为仿制药。

在本实施例中，待测样品三的匹配特征峰的个数与原研药品特征峰的个数 的比值s为1，因此，判定待测样品三为原研药；而其他四个待测样品的匹配特 征峰的个数与原研药品特征峰的个数的比值s皆不在98％～100％之间，因此，判 定这两个待测样品都为仿制药。

以上判定结果均与实际相符。

<实施例五>

本实施例使用的分析仪器与实施例一相同。选用的待测样品有五个，分别 为实施例一中的：待测样品一(来源：厂家1的仿制药)、待测样品二(来源： 厂家2的仿制药)、待测样品三(来源：默沙东)、待测样品五(来源：厂家5 的仿制药)、和待测样品九(来源：厂家9的仿制药)。

参照实施例一所提供的步骤一至步骤五的检测方法，在本实施例，步骤五 中当s为90％～100％之间的任意数值，则判定为该待测样品为原研药，反之，则 判定为该待测样品为仿制药。对五个待测样品依次进行首次判定、二次判定和 三次判定。

在本实施例中，待测样品一和待测样品三的匹配特征峰的个数与原研药品 特征峰的个数的比值s在90％～100％之间，因此，判定待测样品一和待测样品三 为原研药。其他三个待测样品的匹配特征峰的个数与原研药品特征峰的个数的 比值s皆不在90％～100％之间，因此，判定其他三个待测样品都为仿制药。

判定结果除一个与实际不符，其他均与实际相同。

与背景技术相比，上述实施例一至五所提供的技术方案将标准拉曼光谱数 据作为建立判别分析方法的依据，从原研药的标准拉曼光谱数据的整体轮廓和 特征峰特点出发，将待测样品与原研药的拉曼光谱进行整体相似度和特征峰个 数比较，其优点和积极效果如下：

1.准确性高

上述实施例一至五所涉及的技术方案首先将待测样品与原研药的拉曼光谱 进行整体相似度比较，找出不符合首次判定条件的待测样品，再将该待测样品 的特征峰个数与原研药的特征峰个数进行比较，通过再次判定和三次判定，完 成对待测样品的判定，与背景技术相比，本方法对仿制药冒充原研药的鉴别准 确率高。

2.显著降低假阳性率

上述实施例一至五的技术方案将待测样品的拉曼光谱图与原研药的标准拉 曼光谱数据进行比较，采用的相似度分析方法可以对光谱整体轮廓与原研药光 谱差异较大的样品进行快速、准确地判别。对于相似度极大的样品采用特征峰 的分析方法，能获得较准确的分析结果。两种分析方法相结合既从API的角度 出发的同时，也将其他辅料光谱信息也包括在内，同时整体光谱还考虑到光谱 峰强度、及形状对判别结果的影响，使判定结果的假阳性率显著降低。

3.能够满足现场快检的需要

上述实施例一至五的技术方案采用拉曼光谱仪即可实现，方法简单，检测 快速，能够在基层一线检验中实现对仿制药冒充原研药的现场快速检测。

4.更具有合理性

上述实施例一至五所提供的技术方案采用待检药品的特征峰的个数与原研 药品的特征峰个数的比值进行三次判定，由于不同品种药品其原研药特征峰的 个数也不尽相同，从几个至几十个不等，采用百分比的结果进行比较，与背景 技术中采用峰个数差异的绝对值进行比较更具有合理性。

当然，本发明所涉及的仿制药冒充原研药的一种检测方法并不仅仅限定于 上述实施例一至五中的内容。以上内容只是按示例的方式而不是限制的方式来 描述，这些设定和数值可以有变化，而依据本发明的技术方案所作的任何等效 变换，均属于本发明的保护范围。