CN112748200B - 一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法 - Google Patents
一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法 Download PDFInfo
- Publication number
- CN112748200B CN112748200B CN202011560084.6A CN202011560084A CN112748200B CN 112748200 B CN112748200 B CN 112748200B CN 202011560084 A CN202011560084 A CN 202011560084A CN 112748200 B CN112748200 B CN 112748200B
- Authority
- CN
- China
- Prior art keywords
- sample
- spectrum
- near infrared
- conjugated estrogens
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 135
- 229940035811 conjugated estrogen Drugs 0.000 title claims abstract description 65
- 210000002700 urine Anatomy 0.000 title claims abstract description 61
- 230000008569 process Effects 0.000 title claims abstract description 49
- 238000004440 column chromatography Methods 0.000 title claims abstract description 36
- 238000012544 monitoring process Methods 0.000 title claims abstract description 22
- 230000002159 abnormal effect Effects 0.000 claims abstract description 55
- 239000003480 eluent Substances 0.000 claims abstract description 45
- 238000001514 detection method Methods 0.000 claims abstract description 30
- 238000010828 elution Methods 0.000 claims abstract description 21
- 230000003595 spectral effect Effects 0.000 claims abstract description 17
- 238000003908 quality control method Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 107
- 238000012937 correction Methods 0.000 claims description 58
- WEVYAHXRMPXWCK-UHFFFAOYSA-N Acetonitrile Chemical compound CC#N WEVYAHXRMPXWCK-UHFFFAOYSA-N 0.000 claims description 51
- QTTMOCOWZLSYSV-QWAPEVOJSA-M equilin sodium sulfate Chemical compound [Na+].[O-]S(=O)(=O)OC1=CC=C2[C@H]3CC[C@](C)(C(CC4)=O)[C@@H]4C3=CCC2=C1 QTTMOCOWZLSYSV-QWAPEVOJSA-M 0.000 claims description 51
- PMZURENOXWZQFD-UHFFFAOYSA-L Sodium Sulfate Chemical compound [Na+].[Na+].[O-]S([O-])(=O)=O PMZURENOXWZQFD-UHFFFAOYSA-L 0.000 claims description 37
- 229910052938 sodium sulfate Inorganic materials 0.000 claims description 37
- 235000011152 sodium sulphate Nutrition 0.000 claims description 37
- 239000011159 matrix material Substances 0.000 claims description 32
- VUCAHVBMSFIGAI-ZFINNJDLSA-M estrone sodium sulfate Chemical compound [Na+].[O-]S(=O)(=O)OC1=CC=C2[C@H]3CC[C@](C)(C(CC4)=O)[C@@H]4[C@@H]3CCC2=C1 VUCAHVBMSFIGAI-ZFINNJDLSA-M 0.000 claims description 28
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 claims description 27
- 238000002329 infrared spectrum Methods 0.000 claims description 25
- 239000007864 aqueous solution Substances 0.000 claims description 24
- NLLMJANWPUQQTA-SPUZQDLCSA-N estra-1,3,5(10),7-tetraene-3,17alpha-diol Chemical compound OC1=CC=C2[C@H]3CC[C@](C)([C@@H](CC4)O)[C@@H]4C3=CCC2=C1 NLLMJANWPUQQTA-SPUZQDLCSA-N 0.000 claims description 24
- DNXHEGUUPJUMQT-CBZIJGRNSA-N Estrone Chemical compound OC1=CC=C2[C@H]3CC[C@](C)(C(CC4)=O)[C@@H]4[C@@H]3CCC2=C1 DNXHEGUUPJUMQT-CBZIJGRNSA-N 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000005457 optimization Methods 0.000 claims description 16
- 238000004811 liquid chromatography Methods 0.000 claims description 14
- BNIILDVGGAEEIG-UHFFFAOYSA-L disodium hydrogen phosphate Chemical compound [Na+].[Na+].OP([O-])([O-])=O BNIILDVGGAEEIG-UHFFFAOYSA-L 0.000 claims description 12
- 229910000403 monosodium phosphate Inorganic materials 0.000 claims description 12
- 235000019799 monosodium phosphate Nutrition 0.000 claims description 12
- AJPJDKMHJJGVTQ-UHFFFAOYSA-M sodium dihydrogen phosphate Chemical compound [Na+].OP(O)([O-])=O AJPJDKMHJJGVTQ-UHFFFAOYSA-M 0.000 claims description 12
- 238000010187 selection method Methods 0.000 claims description 11
- 239000011259 mixed solution Substances 0.000 claims description 8
- 230000035935 pregnancy Effects 0.000 claims description 8
- 230000003287 optical effect Effects 0.000 claims description 7
- 239000011347 resin Substances 0.000 claims description 7
- 229920005989 resin Polymers 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 5
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical group CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 239000000839 emulsion Substances 0.000 claims description 3
- 229910052736 halogen Inorganic materials 0.000 claims description 3
- 150000002367 halogens Chemical class 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000002203 pretreatment Methods 0.000 claims description 3
- 239000011550 stock solution Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000011426 transformation method Methods 0.000 claims description 3
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 claims description 3
- 229910052721 tungsten Inorganic materials 0.000 claims description 3
- 239000010937 tungsten Substances 0.000 claims description 3
- 230000005526 G1 to G0 transition Effects 0.000 claims description 2
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 claims 1
- KXNLCSXBJCPWGL-UHFFFAOYSA-N [Ga].[As].[In] Chemical compound [Ga].[As].[In] KXNLCSXBJCPWGL-UHFFFAOYSA-N 0.000 claims 1
- 238000004128 high performance liquid chromatography Methods 0.000 abstract description 3
- 239000000463 material Substances 0.000 abstract description 3
- 238000005070 sampling Methods 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 89
- 229960000747 sodium estrone sulfate Drugs 0.000 description 11
- 238000012795 verification Methods 0.000 description 8
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 description 7
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 description 7
- 229910052708 sodium Inorganic materials 0.000 description 7
- 239000011734 sodium Substances 0.000 description 7
- 239000000243 solution Substances 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- NBIIXXVUZAFLBC-UHFFFAOYSA-N Phosphoric acid Chemical compound OP(O)(O)=O NBIIXXVUZAFLBC-UHFFFAOYSA-N 0.000 description 2
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 229960000074 biopharmaceutical Drugs 0.000 description 2
- 239000003153 chemical reaction reagent Substances 0.000 description 2
- 238000004587 chromatography analysis Methods 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 125000000524 functional group Chemical group 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 239000000741 silica gel Substances 0.000 description 2
- 229910002027 silica gel Inorganic materials 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 208000011580 syndromic disease Diseases 0.000 description 2
- NWUYHJFMYQTDRP-UHFFFAOYSA-N 1,2-bis(ethenyl)benzene;1-ethenyl-2-ethylbenzene;styrene Chemical compound C=CC1=CC=CC=C1.CCC1=CC=CC=C1C=C.C=CC1=CC=CC=C1C=C NWUYHJFMYQTDRP-UHFFFAOYSA-N 0.000 description 1
- 239000004925 Acrylic resin Substances 0.000 description 1
- 206010003439 Artificial menopause Diseases 0.000 description 1
- 208000001132 Osteoporosis Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 229910000147 aluminium phosphate Inorganic materials 0.000 description 1
- 239000003957 anion exchange resin Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000001076 estrogenic effect Effects 0.000 description 1
- 239000003456 ion exchange resin Substances 0.000 description 1
- 229920003303 ion-exchange polymer Polymers 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003960 organic solvent Substances 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 125000001453 quaternary ammonium group Chemical group 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000012488 sample solution Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 150000003431 steroids Chemical class 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3577—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/26—Conditioning of the fluid carrier; Flow patterns
- G01N30/28—Control of physical parameters of the fluid carrier
- G01N30/32—Control of physical parameters of the fluid carrier of pressure or speed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/26—Conditioning of the fluid carrier; Flow patterns
- G01N30/28—Control of physical parameters of the fluid carrier
- G01N30/34—Control of physical parameters of the fluid carrier of fluid composition, e.g. gradient
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/50—Conditioning of the sorbent material or stationary liquid
- G01N30/52—Physical parameters
- G01N30/54—Temperature
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/88—Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/483—Physical analysis of biological material
- G01N33/487—Physical analysis of biological material of liquid biological material
- G01N33/493—Physical analysis of biological material of liquid biological material urine
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/26—Conditioning of the fluid carrier; Flow patterns
- G01N30/28—Control of physical parameters of the fluid carrier
- G01N30/32—Control of physical parameters of the fluid carrier of pressure or speed
- G01N2030/324—Control of physical parameters of the fluid carrier of pressure or speed speed, flow rate
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/88—Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86
- G01N2030/8809—Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86 analysis specially adapted for the sample
- G01N2030/8813—Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86 analysis specially adapted for the sample biological materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/62—Detectors specially adapted therefor
- G01N30/74—Optical detectors
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Urology & Nephrology (AREA)
- Molecular Biology (AREA)
- Hematology (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Biophysics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Endocrinology (AREA)
- Biotechnology (AREA)
- Cell Biology (AREA)
- Microbiology (AREA)
Abstract
本发明提供了一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法,涉及质量监控技术领域。本发明提供的近红外质量监控方法可以快速检测从孕马尿中获取结合雌激素成分的柱层析工艺环节中孕马尿洗脱液的质量,比起常规取样后HPLC检测,时间上更省时,且无污染,节约了大量的人力物力。从另一个角度讲,针对从孕马尿中获取结合雌激素成分的柱层析工艺的质量监控,一方面,可以对柱层析洗脱的起点和终点进行判断,另一方面,可以对柱层析过程进行主要质控指标成分的监控。本发明采用基于L1‑PCA的马氏距离法剔除异常光谱值,能够显著提高检测结果的准确性。
Description
技术领域
本发明涉及质量监控技术领域,具体涉及一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法。
背景技术
众所周知,结合雌激素是治疗更年期综合症的有效药物,尤其是天然的结合雌激素混合物具有确切的功效和可靠的安全性,在临床上不仅可以用于治疗和预防女性生理或人工绝经期后出现的更年期综合症,而且能够预防和治疗骨质疏松症,得到人们长期的应用和认可。
早期报道的结合雌激素混合物提取方法的专利有美国专利:US2429398、US2519516、US2696265、US2711988、US2834712等,其多为采用有机溶剂萃取的方法。上世纪60年代开始出现采用活性炭、离子交换树脂、反向硅胶等方法从孕马尿中提取结合雌激素混合物的方法。目前的提取方法主要是利用结合雌激素为甾体结构物质的疏水性来分离制备;其中,国内外陆续有多个结合雌激素分离制备相关的发明专利被公开,包括反相硅胶、含有各种功能团的大孔树脂、苯乙烯-二乙烯聚合物非极性树脂、聚丙烯酸酯树脂、含季铵官能团的强碱性阴离子交换树脂等吸附树脂等。
采用大孔吸附树脂富集、提取是从孕马尿中提取结合雌激素的关键工艺步骤,常规的工艺监测方法是取样后,回到实验室采用HPLC检测雌酮硫酸钠和马烯雌酮硫酸钠等的含量,通常需要几小时甚至1天后才能获得结果,滞后于柱层析过程,无法实现柱层析工艺的过程控制。
在近红外光谱分析领域,异常光谱的剔除方法常采用马氏距离法,但马氏距离法要求总样本数要大于样本的维数,处理繁琐。
发明内容
本发明的目的在于提供一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法,本发明提供的方法能够快速检测从孕马尿中获取结合雌激素成分的柱层析工艺中孕马尿洗脱液的质量,相比于常规取样后再进行液相色谱检测,更加省时,且无污染,节约了大量人力物力。本发明采用基于L1-PCA的马氏距离法剔除异常光谱值,能够显著提高检测结果的准确性。
为了实现上述发明目的,本发明提供以下技术方案:
本发明提供了一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法,包括以下步骤:
采集孕马尿原液柱层析得到的洗脱液,作为待测样品;
将所述待测样品进行近红外光谱扫描,得到原始光谱数据,采用基于L1-PCA的马氏距离法剔除所述原始光谱数据中的异常光谱值,将剔除异常光谱值后的光谱数据导入校正模型,得到所述待测样品中结合雌激素的含量;
所述校正模型为真实值和测定值之间的线性方程,所述测定值为剔除异常光谱值后的近红外光谱数据;
所述结合雌激素包括17α-二氢马烯雌酮硫酸钠、马烯雌酮硫酸钠和雌酮硫酸钠中的一种或多种。
优选地,所述校正模型的建立方法包括以下步骤:
(1)将孕马尿原液进行柱层析,得到孕马尿洗脱液样本;
(2)将所述孕马尿洗脱液样本进行液相色谱检测,得到孕马尿洗脱液样本中结合雌激素的真实含量值;
(3)将步骤(1)所述的孕马尿洗脱液样本进行近红外光谱扫描,得到原始样本光谱数据,采用基于L1-PCA的马氏距离法剔除异常样本光谱值,采集孕马尿洗脱液样本的光谱数据;
(4)对步骤(3)采集的光谱数据进行预处理,然后对所述预处理后的光谱数据进行波段选择,得到特征波段;采用偏最小二乘法对特征波段的光谱数据和孕马尿洗脱液样本中对应的结合雌激素真实含量值做回归拟合,建立校正模型;
所述步骤(2)和步骤(3)没有时间顺序的限制。
优选地,不同结合雌激素的校正模型为:
17α-二氢马烯雌酮硫酸钠的校正模型为:y=0.9173x+0.0128;
马烯雌酮硫酸钠的校正模型为:y=0.9079x+0.0258;
雌酮硫酸钠的校正模型为:y=0.9151x+0.0396;
马烯雌酮硫酸钠+雌酮硫酸钠的校正模型为:y=0.9148x+0.0636;
上述校正模型中,x代表真实值,y代表预测值。
优选地,当所述待测样品中结合雌激素的总含量大于0.001mg/mL时,判定为孕马尿柱层析洗脱的起点;
当所述待测样品中结合雌激素的总含量小于0.001mg/mL时,判定为孕马尿柱层析洗脱的终点。
优选地,所述基于L1-PCA的马氏距离法剔除异常光谱值包括:
将所述原始光谱数据构造光谱矩阵;
按照式I所示计算式,采用L1-PCA算法求解所述光谱矩阵,得到光谱主成分;
按照式II所示计算式将所述主成分构造协方差矩阵;
将所述方差矩阵按照式III所示计算式计算得到马氏距离;
设置阈值和剔除异常光谱值;
式I中,X'为n×m的光谱样本矩阵,n为样本数,m为每条光谱所采集的数据点数;U为投影矩阵;V为系数矩阵;L1为矩阵的1范数;
S=T'T/n,式II;
式II中,T'为T的转置,n为样本数,T的计算方法为:获得光谱数据的信号子空间P后,然后P根据计算均值光谱向量μ,对于P矩阵的每个样本减去均值向量μ,得到T;
式III中,P为光谱数据的信号子空间;μ为均值光谱向量;S为由T构造的样本信号子空间协方差矩阵;
所述阈值为2~3。
优选地,步骤(2)所述液相色谱检测的参数包括:
色谱柱:C18色谱柱;
色谱柱规格:250mm×4.6mm,5μm,100A;
流动相:包括A相和B相,其中A相为磷酸二氢钠水溶液、乙腈和甲醇的混合溶液,所述磷酸二氢钠水溶液的浓度为20mmol/L,pH值为3.5;所述A相中磷酸二氢钠水溶液、乙腈和甲醇的体积比为17:2:1;B相为磷酸氢二钠水溶液和乙腈的混合溶液,所述磷酸氢二钠水溶液的浓度为10mmol/L,pH值为3.5;所述B相中磷酸氢二钠水溶液和乙腈的体积比为3:7;
流动相洗脱程序为:0~18min,A相体积分数由70%降低为67%;18~23min,A相体积分数由67%降低为20%;23~28min,A相体积分数由20%升高为70%;28~35min,A相体积分数稳定在70%;
流速:1.0mL/min;
柱温:40℃;
检测波长:205nm;
进样量:1μL。
优选地,所述近红外光谱扫描包括:
采用在线或离线检测,以空气为背景,采用透射测样方式,波长检测范围为10000~4000cm-1,扫描次数为32次,分辨率为8cm-1,光程为2mm,每个待测样品重复扫描3~5次,取平均值作为原始光谱数据;
或者,依据光栅扫描分光原理,以卤钨灯为光源;光谱范围1000~1800nm,检测器:InGaAs检测器,分辨率:8cm-1;扫描次数:32次;光程1mm。
优选地,所述步骤(4)中预处理的方法包括:卷积平滑、一阶卷积求导、二阶卷积求导、多元散射校正、标准正态变量变换和归一化法中的一种或几种相结合的方法。
优选地,所述步骤(4)中波段选择的方法包括全波长、相关系数法选择波长区间、相关成分法选择波长区间、迭代优化波长选择方法1或迭代优化波长选择方法2。
本发明提供了一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法,本发明建立校正模型,该校正模型为真实值和测定值之间的线性方程,所述测定值为剔除异常光谱值后的近红外光谱数据;本发明采用基于L1-PCA的马氏距离法剔除异常光谱值,能够显著提高检测结果的准确性。本发明采用基于L1-PCA的马氏距离法,通过数据降维,消除众多信息共存中相互重叠的信息部分,更方便用于小样本数量的处理,并且抑制了拖尾噪声,提高信号的识别度。另外,当提取的特征数量较少时,基于L1-PCA的马氏距离法更适用于异常光谱的剔除过程。
本发明的方法可以快速检测从孕马尿中获取结合雌激素成分的柱层析工艺环节中孕马尿洗脱液的质量,比起常规取样后HPLC检测,时间上更省时,且无污染,节约了大量的人力物力。从另一个角度讲,针对从孕马尿中获取结合雌激素成分的柱层析工艺的质量监控,一方面,可以对柱层析洗脱的起点和终点进行判断,另一方面,可以对柱层析过程进行主要质控指标成分(17α-二氢马烯雌酮硫酸钠、马烯雌酮硫酸钠、雌酮硫酸钠、马烯雌酮硫酸钠+雌酮硫酸钠)的监控。
附图说明
图1为液相色谱检测中17α-二氢马烯雌酮硫酸钠的含量数据图;
图2为液相色谱检测中马烯雌酮硫酸钠的含量数据图;
图3为液相色谱检测中雌酮硫酸钠的含量数据图;
图4为液相色谱检测中马烯雌酮硫酸钠+雌酮硫酸钠的含量数据图;
图5为结合雌激素近红外光谱图;
图6为基于L1-PCA的马氏距离法计算异常光谱的结果图;
图7为基于马氏距离法计算异常光谱的结果图;
图8为17α-二氢马烯雌酮硫酸钠建模样本集的含量趋势图;
图9为20181211-2批次样本17α-二氢马烯雌酮硫酸钠的预测趋势图;
图10为剔除图9中的异常样本后得到的预测趋势图;
图11为马烯雌酮硫酸钠建模样本集的含量趋势图;
图12为20181211-2批次样本马烯雌酮硫酸钠的预测趋势图;
图13为剔除图12中的异常样本后得到的预测趋势图;
图14为雌酮硫酸钠建模样本集的含量趋势图;
图15为20181211-2批次样本雌酮硫酸钠的预测趋势图;
图16为剔除图15中的异常样本后得到的预测趋势图;
图17为马烯雌酮硫酸钠+雌酮硫酸钠建模样本集的含量趋势图;
图18为20181211-2批次样本马烯雌酮硫酸钠+雌酮硫酸钠的预测趋势图;
图19为剔除图17中的异常样本后得到的预测趋势图。
具体实施方式
本发明提供了一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法,包括以下步骤:
采集孕马尿原液柱层析得到的洗脱液,作为待测样品;
将所述待测样品进行近红外光谱扫描,得到原始光谱数据,采用基于L1-PCA的马氏距离法剔除所述原始光谱数据中的异常光谱值,将剔除异常光谱值后的光谱数据导入校正模型,得到所述待测样品中结合雌激素的含量;
所述校正模型为真实值和测定值之间的线性方程,所述测定值为剔除异常光谱值后的近红外光谱数据;
所述结合雌激素包括17α-二氢马烯雌酮硫酸钠、马烯雌酮硫酸钠和雌酮硫酸钠中的一种或多种。
本发明采集孕马尿原液柱层析得到的洗脱液,作为待测样品。在本发明中,所述柱层析的固定相优选为大孔树脂,流动相优选为乙醇。本发明对所述柱层析的具体工艺参数没有特殊要求,采用本领域技术人员所熟知的工艺即可。
得到待测样品后,本发明将所述待测样品进行近红外光谱扫描,得到原始光谱数据,采用基于L1-PCA的马氏距离法剔除所述原始光谱数据中的异常光谱值,将剔除异常光谱值后的光谱数据导入校正模型,得到所述待测样品中结合雌激素的含量;所述校正模型为真实值和测定值之间的线性方程,所述测定值为剔除异常光谱值后的近红外光谱数据;所述结合雌激素包括17α-二氢马烯雌酮硫酸钠、马烯雌酮硫酸钠和雌酮硫酸钠中的一种或多种。
在本发明中,所述近红外光谱扫描优选包括:
采用在线或离线检测,以空气为背景,采用透射测样方式,波长检测范围为10000~4000cm-1,扫描次数为32次,分辨率为8cm-1,光程为2mm,每个供试品溶液重复扫描3~5次,取平均值作为光谱数据;
或者,依据光栅扫描分光原理,以卤钨灯为光源;光谱范围1000~1800nm,检测器:InGaAs检测器,分辨率:8cm-1;扫描次数:32次;光程1mm。
在本发明中,每次扫描平均需要3~5秒。
在本发明中,所述基于L1-PCA的马氏距离法剔除异常光谱值优选包括:
将所述原始光谱数据构造光谱矩阵;
按照式I所示计算式,采用L1-PCA算法求解所述光谱矩阵,得到光谱主成分;
按照式II所示计算式将所述主成分构造协方差矩阵;
将所述方差矩阵按照式III所示计算式计算得到马氏距离;
设置阈值和剔除异常光谱值;
式I中,X'为n×m的光谱样本矩阵,n为样本数,m为每条光谱所采集的数据点数;U为投影矩阵;V为系数矩阵;L1为矩阵的1范数;
S=T'T/n,式II;
式II中,T'为T的转置,n为样本数,T的计算方法为:获得光谱数据的信号子空间P后,然后P根据计算均值光谱向量μ,对于P矩阵的每个样本减去均值向量μ,得到T;
式III中,P为光谱数据的信号子空间;μ为均值光谱向量;S为由T构造的样本信号子空间协方差矩阵;
所述阈值为2~3。
在本发明的具体实施例中,所述采用L1-PCA算法求解所述光谱矩阵,得到光谱主成分,即是求解优化问题,在求解式I的的最优化问题时,由于L1范数构成的目标函数并不是一个凸函数,因此不直接通过凸优化算法来求解,令U与V交替已知,该代价函数就变成一个凸函数,然后通过凸优化算法进行求解。
在本发明中,所述按照式II所示计算式将所述主成分构造协方差矩阵中所选主成分相应特征值占所有特征值之和的95%以上。
在本发明的具体实施例中,计算所得马氏距离即为L1范数约束主成分分析后的马氏距离;
在本发明的具体实施例中,所述阈值为2.5。本发明根据阈值的范围剔除异常样本光谱值。
在本发明中,所述校正模型的建立方法优选包括以下步骤:
(1)将孕马尿原液进行柱层析,得到孕马尿洗脱液样本;
(2)将所述孕马尿洗脱液样本进行液相色谱检测,得到孕马尿洗脱液样本中结合雌激素的真实含量值;
(3)将步骤(1)所述的孕马尿洗脱液样本进行近红外光谱扫描,得到原始样本光谱数据,采用基于L1-PCA的马氏距离法剔除异常样本光谱值,采集孕马尿洗脱液样本的光谱数据;
(4)对步骤(3)采集的光谱数据进行预处理,然后对所述预处理后的光谱数据进行波段选择,得到特征波段;采用偏最小二乘法对特征波段的光谱数据和孕马尿洗脱液样本中对应的结合雌激素真实含量值做回归拟合,建立校正模型;
所述步骤(2)和步骤(3)没有时间顺序的限制。
在本发明中,建立校正模型时柱层析的工艺与采集待测样品时柱层析的工艺相同,这里就不再赘述。
在本发明中,所述液相色谱检测的参数优选包括:
色谱柱:C18色谱柱;
色谱柱规格:250mm×4.6mm,5μm,100A;
流动相:包括A相和B相,其中A相为磷酸二氢钠水溶液、乙腈和甲醇的混合溶液,所述磷酸二氢钠水溶液的浓度为20mmol/L,pH值为3.5;所述A相中磷酸二氢钠水溶液、乙腈和甲醇的体积比为17:2:1;B相为磷酸氢二钠水溶液和乙腈的混合溶液,所述磷酸氢二钠水溶液的浓度为10mmol/L,pH值为3.5;所述B相中磷酸氢二钠水溶液和乙腈的体积比为3:7;
流动相洗脱程序为:0~18min,A相体积分数由70%降低为67%;18~23min,A相体积分数由67%降低为20%;23~28min,A相体积分数由20%升高为70%;28~35min,A相体积分数稳定在70%;
流速:1.0mL/min;
柱温:40℃;
检测波长:205nm;
进样量:1μL;
不同结合雌激素在同一色谱条件下的不同保留时间处出峰。
本发明优选在对所述孕马尿洗脱液样本进行液相色谱检测后,剔除异常数据值,得到孕马尿洗脱液样本中结合雌激素的真实含量值。本发明对所述剔除异常数据值的方法没有特殊要求,采用本领域技术人员所熟知的方法即可。在本发明的具体实施例中,所述液相色谱检测所得异常数据可以直观观察到,直接剔除即可。
本发明将所述的孕马尿洗脱液样本进行近红外光谱扫描,得到原始样本光谱数据,采用基于L1-PCA的马氏距离法剔除异常样本光谱值,采集孕马尿洗脱液样本的光谱数据。在本发明中,所述近红外光谱扫描的参数,以及基于L1-PCA的马氏距离法剔除异常样本光谱值的方法与前文检测待测样品时的参数一致,这里不再赘述。
得到孕马尿洗脱液样本的光谱数据后,本发明对采集的光谱数据进行预处理,然后对所述预处理后的光谱数据进行波段选择,得到特征波段;采用偏最小二乘法对特征波段的光谱数据和孕马尿洗脱液样本中对应的结合雌激素真实含量值做回归拟合计算,建立校正模型。
在本发明中,所述预处理的方法优选包括:卷积平滑、一阶卷积求导、二阶卷积求导、多元散射校正、标准正态变量变换和归一化法中的一种或几种相结合的方法,更优选为卷积平滑。
在本发明中,所述波段选择的方法优选包括全波长、相关系数法选择波长区间、相关成分法选择波长区间、迭代优化波长选择方法1或迭代优化波长选择方法2,更优选为迭代优化波长选择方法1。在本发明中,所述迭代优化波长选择方法1是指对N个波长区间段进行全排列组合,使用每一组合进行建模,最后选择SECV最小的为本次优化的最佳模型;所述迭代优化波长选择方法2是指从N个波长区间段中选择M段组合成一段光谱,进行建模,即N选M,对所有的可能组合进行建模,最后选择SECV最小的为本次优化的最佳模型,本部分选择的N为10,M为1、2或3。
在本发明中,所述结合雌激素包括17α-二氢马烯雌酮硫酸钠、马烯雌酮硫酸钠和雌酮硫酸钠中的一种或多种,优选为17α-二氢马烯雌酮硫酸钠、马烯雌酮硫酸钠、雌酮硫酸钠以及马烯雌酮硫酸钠+雌酮硫酸钠;其中,马烯雌酮硫酸钠+雌酮硫酸钠是指以二者的含量和为指标建立校正模型。
在本发明中,不同结合雌激素的校正模型优选为:
17α-二氢马烯雌酮硫酸钠的校正模型为:y=0.9173x+0.0128;
马烯雌酮硫酸钠的校正模型为:y=0.9079x+0.0258;
雌酮硫酸钠的校正模型为:y=0.9151x+0.0396;
马烯雌酮硫酸钠+雌酮硫酸钠的校正模型为:y=0.9148x+0.0636;
上述校正模型中,x代表真实值,y代表预测值。
在本发明的具体实施例中,不同结合雌激素的校正模型如表1所示:
表1不同结合雌激素的校正模型
表1中预测-真实拟合方程中x代表真实值,y代表预测值。
在本发明的具体实施例中,将柱层析过程中的孕马尿洗脱液进行近红外光谱扫描后得到的光谱数据作为预测值,导入所述校正模型,得到柱层析过程中孕马尿洗脱液中结合雌激素的真实值,据此,实现孕马尿柱层析过程的质量监控。
在本发明中,当孕马尿洗脱液中结合雌激素的含量大于0.001mg/mL时,判定为孕马尿柱层析洗脱的起点;当孕马尿洗脱液中结合雌激素的含量小于0.001mg/mL时,判定为孕马尿柱层析洗脱的终点。
采用本发明提供的方法能够及时判断柱层析洗脱的起点和终点,精确控制柱层析过程。
下面将结合本发明中的实施例,对本发明中的技术方案进行清楚、完整地描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例中采用的实验仪器:
美国Waters 2996高效液相色谱仪(包括梯度泵G1311A、自动进样器G1329A、柱恒温系统G1316A、二极管阵列检测器DAD-G1315B、色谱工作站);
采用的实验试剂:
磷酸(分析级,广州化学试剂厂),甲醇、乙腈(色谱级,德国默克公司),水(屈臣氏有限公司)。
实验材料:
孕马尿洗脱液样本(180个,新疆新姿源生物制药有限责任公司提供),17α-二氢马烯雌酮硫酸钠、马烯雌酮硫酸钠、雌酮硫酸钠的混合标准品,由新疆新姿源生物制药有限责任公司提供。
实施例1
(1)将孕马尿原液经大孔树脂洗脱,得到的不同批次的孕马尿洗脱液样本;
(2)将所述孕马尿洗脱液样本进行液相色谱检测,得到孕马尿洗脱液样本中结合雌激素的真实含量值;所述液相色谱检测的参数为:
色谱柱:Sharpsil-UC18;
色谱柱规格:250mm×4.6mm,5μm,100A;
流动相:A相和B相,其中A相为磷酸二氢钠水溶液、乙腈和甲醇的混合溶液,所述A相中磷酸二氢钠水溶液的浓度为20mmol/L,pH值为3.5;所述A相中磷酸二氢钠水溶液、乙腈和甲醇的体积比为17:2:1;B相为磷酸氢二钠水溶液和乙腈的混合溶液,所述B相中磷酸氢二钠水溶液的浓度为10mmol/L,pH值为3.5;所述B相中磷酸氢二钠水溶液和乙腈的体积比为3:7;
流动相洗脱程序为:0~18min,A相体积分数由70%降低为67%;18~23min,A相体积分数由67%降低为20%;23~28min,A相体积分数由20%升高为70%;28~35min,A相体积分数稳定在70%;
流速:1.0mL/min;
柱温:40℃;
检测波长:205nm;
进样量:1μL;
不同结合雌激素在同一色谱条件下的不同保留时间处出峰。
所得孕马尿洗脱液样本中结合雌激素的真实含量值如表2所示;
表2结合雌激素中各质控指标成分的含量值结果表
上述得到171个样品的测定结果,从各批次的趋势图进行分析,存在有测定的异常值,如图1~4所示,需要剔除圈出的异常含量数据,方可进行校正模型的建立。
本发明采用聚光NIR1500对所述孕马尿洗脱液样本进行近红外光谱扫描,采用基于L1-PCA的马氏距离法剔除异常样本光谱值,采集孕马尿洗脱液样本的光谱数据;采用离线检测,以空气为背景,采用透射测样方式,其中波长检测范围为10000~4000cm-1,扫描次数为64次,分辨率为8cm-1,光程为2mm;每个孕马尿样品重复扫描4次,每测一次平均需要4秒,取平均值作为光谱数据;采用卷积平滑对采集的光谱数据进行预处理,然后采用迭代优化波长选择方法1的波段选择方法,采用偏最小二乘法对特征波段的光谱数据和孕马尿洗脱液样本中对应的结合雌激素真实含量值做回归拟合计算,建立校正模型,具体如下:
经聚光NIR1500对孕马尿洗脱液样本进行NIR光谱的采集,得出结果如图5所示,从图5中看出,存在有异常光谱,采用基于L1-PCA的马氏距离法,设置阈值2~3,可以找出MTC-20181209-1-1、MTC-20181210-2-1为异常光谱,剔除这两个异常光谱,再进行校正模型的建立,如图6所示。以采用马氏距离法计算异常光谱的结果为对比例,如图7所示。由图6~7对比可以看出,马氏距离法无法计算出异常光谱数据,而采用采用基于L1-PCA的马氏距离法可以准确地识别出异常光谱数据,有利于提高检测结果的准确性。
(1)17α-二氢马烯雌酮硫酸钠建模过程中,采用采集的5批结合雌激素样本,即MTC20181209-1、MTC20181209-2、MTC20181210-1、MTC20181210-2、MTC20181211-1作为校正集,将1批的结合雌激素样本MTC20181211-2作为验证集,同时,采用L1-PCA的马氏距离法,设置阈值2~3,剔除异常光谱,再进行校正模型的建立和未知样本的预测,如表3所示;
表3 17α-二氢马烯雌酮硫酸钠校正集样本
17α-二氢马烯雌酮硫酸钠建模样本集的含量趋势图如图8所示。
17α-二氢马烯雌酮硫酸钠的校正模型如表4所示;
表4 17α-二氢马烯雌酮硫酸钠的校正模型
表4中预测-真实拟合方程中x代表真实值,y代表预测值。
17α-二氢马烯雌酮硫酸钠样本的预测结果:
采用已建立的校正模型对20181211-2批次的样本进行含量预测,结果如表5所示:
表5 20181211-2批次样本17α-二氢马烯雌酮硫酸钠的预测结果表
20181211-2批次样本17α-二氢马烯雌酮硫酸钠的预测趋势图如图9所示。通过上述对17α-二氢马烯雌酮硫酸钠的校正模型的验证,结果表明20181211-2批次中预测值与真实值的含量趋势一致。但在预测过程中,存在有异常点,即洗脱过程中的第6个点,出现预测值产生较大偏差的原因可能是因为近红外光谱采集过程中出现的误差、或者样品放置时间太久影响最终测定含量导致的。除去洗脱过程中的第6个点,得到的预测趋势图如图10所示。由图10可以看出,20181211-2批次中预测值与真实值的含量趋势一致。
(2)马烯雌酮硫酸钠建模过程中,采用采集的5批的结合雌激素样本作为校正集,将1批的结合雌激素样本作为验证集,采用L1-PCA的马氏距离法,设置阈值2~3,剔除异常光谱,再进行校正模型的建立和未知样本的预测,如表6所示;
马烯雌酮硫酸钠建模样本集的含量趋势图如图11所示。
马烯雌酮硫酸钠的校正模型如表6所示;
表6马烯雌酮硫酸钠的校正模型
表6中预测-真实拟合方程中x代表真实值,y代表预测值。
马烯雌酮硫酸钠样本的预测结果:
采用已建立的校正模型对20181211-2批次的样本进行含量预测,结果如表7所示:
表7 20181211-2批次样本马烯雌酮硫酸钠的预测结果表
20181211-2批次样本马烯雌酮硫酸钠的预测趋势图如图12所示。通过上述对马烯雌酮硫酸钠校正模型的验证,结果表明20181211-2批次中预测值与真实值的含量趋势一致。但在预测过程中,存在有异常点,即洗脱过程中的第6个点和第25个点,出现预测值产生较大偏差的原因可能是因为近红外光谱采集过程中出现的误差、或者样品放置时间太久影响最终测定含量导致的。除去洗脱过程中的第6个点和第25个点,得到的预测趋势图如图13所示。由图13可以看出,20181211-2批次中预测值与真实值的含量趋势一致。
(3)雌酮硫酸钠建模过程中,采用采集的5批的结合雌激素样本作为校正集,将1批的结合雌激素样本作为验证集,采用L1-PCA的马氏距离法,设置阈值2~3,剔除异常光谱,再进行模型的建立和未知样本的预测,如表8所示;
雌酮硫酸钠建模样本集的含量趋势图如图14所示。
雌酮硫酸钠的校正模型如表8所示;
表8雌酮硫酸钠的校正模型
表8中预测-真实拟合方程中x代表真实值,y代表预测值。
雌酮硫酸钠样本的预测结果:
采用已建立的校正模型对20181211-2批次的样本进行含量预测,结果如表9所示:
表9 20181211-2批次样本雌酮硫酸钠的预测结果表
20181211-2批次样本雌酮硫酸钠的预测趋势图如图15所示。通过上述对雌酮硫酸钠校正模型的验证,结果表明20181211-2批次中预测值与真实值的含量趋势一致。但在预测过程中,存在有异常点,即洗脱过程中的第6个点和第25个点,出现预测值产生较大偏差的原因可能是因为近红外光谱采集过程中出现的误差、或者样品放置时间太久影响最终测定含量导致的。除去洗脱过程中的第6个点和第25个点,得到的预测趋势图如图16所示。由图16可以看出,20181211-2批次中预测值与真实值的含量趋势一致。
(4)马烯雌酮硫酸钠+雌酮硫酸钠建模过程中,采用采集的5批的结合雌激素样本作为校正集,将1批的结合雌激素样本作为验证集,采用基于L1-PCA的马氏距离法,设置阈值2~3,剔除异常光谱,再进行模型的建立和未知样本的预测,如表10所示;
马烯雌酮硫酸钠+雌酮硫酸钠建模样本集的含量趋势图如图17所示。
马烯雌酮硫酸钠+雌酮硫酸钠的校正模型如表10所示;
表10马烯雌酮硫酸钠+雌酮硫酸钠的校正模型
表10中预测-真实拟合方程中x代表真实值,y代表预测值。
马烯雌酮硫酸钠+雌酮硫酸钠样本的预测结果:
采用已建立的校正模型对20181211-2批次的样本进行含量预测,结果如表11所示:
表11 20181211-2批次样本马烯雌酮硫酸钠+雌酮硫酸钠的预测结果表
20181211-2批次样本马烯雌酮硫酸钠+雌酮硫酸钠的预测趋势图如图18所示。通过上述对马烯雌酮硫酸钠+雌酮硫酸钠校正模型的验证,结果表明20181211-2批次中预测值与真实值的含量趋势一致。但在预测过程中,存在有异常点,即洗脱过程中的第6个点和第25个点,出现预测值产生较大偏差的原因可能是因为近红外光谱采集过程中出现的误差、或者样品放置时间太久影响最终测定含量导致的。除去洗脱过程中的第5个点、第6个点和第25个点,得到的预测趋势图如图19所示。由图19可以看出,20181211-2批次中预测值与真实值的含量趋势一致。
对比例
与实施例1基本相同,不同之处仅在于,未剔除异常光谱数据。包含异常光谱数据进行模型建立,得出如下:
17α-二氢马烯雌酮硫酸钠的校正模型如表12所示:
表12 17α-二氢马烯雌酮硫酸钠的校正模型
表12中预测-真实拟合方程中x代表真实值,y代表预测值。
17α-二氢马烯雌酮硫酸钠样本的预测结果:
采用已建立的校正模型对20181211-2批次的样本进行含量预测,结果如表13所示:
表13 20181211-2批次样本17α-二氢马烯雌酮硫酸钠的预测结果表
马烯雌酮硫酸钠的校正模型表14所示:
表14马烯雌酮硫酸钠的校正模型
表14中预测-真实拟合方程中x代表真实值,y代表预测值。
马烯雌酮硫酸钠样本的预测结果:
采用已建立的校正模型对20181211-2批次的样本进行含量预测,结果如表15所示:
表15 20181211-2批次样本马烯雌酮硫酸钠的预测结果表
雌酮硫酸钠的校正模型如表16所示:
表16雌酮硫酸钠的校正模型
表16中预测-真实拟合方程中x代表真实值,y代表预测值。
雌酮硫酸钠样本的预测结果:
采用已建立的校正模型对20181211-2批次的样本进行含量预测,结果如表17所示:
表17 20181211-2批次样本雌酮硫酸钠的预测结果表
马烯雌酮硫酸钠+雌酮硫酸钠的校正模型如表18所示:
表18马烯雌酮硫酸钠+雌酮硫酸钠的校正模型
表18中预测-真实拟合方程中x代表真实值,y代表预测值。
马烯雌酮硫酸钠+雌酮硫酸钠样本的预测结果:
采用已建立的校正模型对20181211-2批次的样本进行含量预测,结果如表19所示:
表19 20181211-2批次样本马烯雌酮硫酸钠+雌酮硫酸钠的预测结果表
上述包含异常光谱的预测值与实施例1中不包含异常光谱的预测值比较起来,绝对偏差更大,预测不够准确。根据实施例的记载内容可以看出,本发明提供的方法准确度较高,可以快速检测孕马尿柱层析工艺环节中孕马尿洗脱液的质量。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法,包括以下步骤:
采集孕马尿原液柱层析得到的洗脱液,作为待测样品;所述柱层析的固定相为大孔树脂,流动相为乙醇;
将所述待测样品进行近红外光谱扫描,得到原始光谱数据,采用基于L1-PCA的马氏距离法剔除所述原始光谱数据中的异常光谱值,将剔除异常光谱值后的光谱数据导入校正模型,得到所述待测样品中结合雌激素的含量;
所述校正模型为真实值和测定值之间的线性方程,所述测定值为剔除异常光谱值后的近红外光谱数据;
所述结合雌激素包括17α-二氢马烯雌酮硫酸钠、马烯雌酮硫酸钠和雌酮硫酸钠中的一种或多种;
所述基于L1-PCA的马氏距离法剔除异常光谱值包括:
将所述原始光谱数据构造光谱矩阵;
按照式I所示计算式,采用L1-PCA算法求解所述光谱矩阵,得到光谱主成分;
按照式II所示计算式将所述主成分构造协方差矩阵;
将所述方差矩阵按照式III所示计算式计算得到马氏距离;
设置阈值和剔除异常光谱值;
式I中,X'为n×m的光谱样本矩阵,n为样本数,m为每条光谱所采集的数据点数;U为投影矩阵;V为系数矩阵;L1为矩阵的1范数;
S=T'T/n,式II;
式II中,T'为T的转置,n为样本数,T的计算方法为:获得光谱数据的信号子空间P后,然后P根据计算均值光谱向量μ,对于P矩阵的每个样本减去均值向量μ,得到T;
式III中,P为光谱数据的信号子空间;μ为均值光谱向量;S为由T构造的样本信号子空间协方差矩阵;
所述阈值为2~3;
不同结合雌激素的校正模型为:
17α-二氢马烯雌酮硫酸钠的校正模型为:y=0.9173x+0.0128;
马烯雌酮硫酸钠的校正模型为:y=0.9079x+0.0258;
雌酮硫酸钠的校正模型为:y=0.9151x+0.0396;
马烯雌酮硫酸钠+雌酮硫酸钠的校正模型为:y=0.9148x+0.0636;
上述校正模型中,x代表真实值,y代表预测值;
当所述待测样品中结合雌激素的总含量大于0.001mg/mL时,判定为孕马尿柱层析洗脱的起点;
当所述待测样品中结合雌激素的总含量小于0.001mg/mL时,判定为孕马尿柱层析洗脱的终点。
2.根据权利要求1所述的近红外质量监控方法,其特征在于,所述校正模型的建立方法包括以下步骤:
(1)将孕马尿原液进行柱层析,得到孕马尿洗脱液样本;
(2)将所述孕马尿洗脱液样本进行液相色谱检测,得到孕马尿洗脱液样本中结合雌激素的真实含量值;
(3)将步骤(1)所述的孕马尿洗脱液样本进行近红外光谱扫描,得到原始样本光谱数据,采用基于L1-PCA的马氏距离法剔除异常样本光谱值,采集孕马尿洗脱液样本的光谱数据;
(4)对步骤(3)采集的光谱数据进行预处理,然后对所述预处理后的光谱数据进行波段选择,得到特征波段;采用偏最小二乘法对特征波段的光谱数据和孕马尿洗脱液样本中对应的结合雌激素真实含量值做回归拟合,建立校正模型;
所述步骤(2)和步骤(3)没有时间顺序的限制。
3.根据权利要求2所述的近红外质量监控方法,其特征在于,步骤(2)所述液相色谱检测的参数包括:
色谱柱:C18色谱柱;
色谱柱规格:250mm×4.6mm,5μm,100A;
流动相:包括A相和B相,其中A相为磷酸二氢钠水溶液、乙腈和甲醇的混合溶液,所述磷酸二氢钠水溶液的浓度为20mmol/L,pH值为3.5;所述A相中磷酸二氢钠水溶液、乙腈和甲醇的体积比为17:2:1;B相为磷酸氢二钠水溶液和乙腈的混合溶液,所述磷酸氢二钠水溶液的浓度为10mmol/L,pH值为3.5;所述B相中磷酸氢二钠水溶液和乙腈的体积比为3:7;
流动相洗脱程序为:0~18min,A相体积分数由70%降低为67%;18~23min,A相体积分数由67%降低为20%;23~28min,A相体积分数由20%升高为70%;28~35min,A相体积分数稳定在70%;
流速:1.0mL/min;
柱温:40℃;
检测波长:205nm;
进样量:1μL。
4.根据权利要求1或2所述的近红外质量监控方法,其特征在于,所述近红外光谱扫描包括:
采用在线或离线检测,以空气为背景,采用透射测样方式,波长检测范围为10000~4000cm-1,扫描次数为32次,分辨率为8cm-1,光程为2mm,每个待测样品重复扫描3~5次,取平均值作为原始光谱数据;
或者,依据光栅扫描分光原理,以卤钨灯为光源;光谱范围1000~1800nm,检测器:InGaAs检测器,分辨率:8cm-1;扫描次数:32次;光程1mm。
5.根据权利要求2所述的近红外质量监控方法,其特征在于,所述步骤(4)中预处理的方法包括:卷积平滑、一阶卷积求导、二阶卷积求导、多元散射校正、标准正态变量变换和归一化法中的一种或几种相结合的方法。
6.根据权利要求2所述的近红外质量监控方法,其特征在于,所述步骤(4)中波段选择的方法包括全波长、相关系数法选择波长区间、相关成分法选择波长区间、迭代优化波长选择方法1或迭代优化波长选择方法2。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011560084.6A CN112748200B (zh) | 2020-12-25 | 2020-12-25 | 一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法 |
US17/510,667 US20220206019A1 (en) | 2020-12-25 | 2021-10-26 | NEAR-INFRARED (NIR) QUALITY MONITORING METHOD USED IN COLUMN CHROMATOGRAPHY FOR EXTRACTING CONJUGATED ESTROGENS (CEs) FROM PREGNANT MARE URINE (PMU) |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011560084.6A CN112748200B (zh) | 2020-12-25 | 2020-12-25 | 一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112748200A CN112748200A (zh) | 2021-05-04 |
CN112748200B true CN112748200B (zh) | 2023-06-23 |
Family
ID=75645883
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011560084.6A Active CN112748200B (zh) | 2020-12-25 | 2020-12-25 | 一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20220206019A1 (zh) |
CN (1) | CN112748200B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114184712B (zh) * | 2021-12-27 | 2023-09-26 | 南通联亚药业股份有限公司 | 一种结合雌激素片中游离甾体杂质含量的检测方法 |
CN115274001B (zh) * | 2022-07-25 | 2023-04-25 | 江南大学 | 基于近红外光谱的有机物在线检测结果的可靠性评价方法 |
CN117171566B (zh) * | 2022-08-17 | 2024-05-10 | 无锡迅杰光远科技有限公司 | 样品含量识别方法、装置和存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070222100A1 (en) * | 2006-03-21 | 2007-09-27 | Huber Engineered Woods L.L.C. | Method and system using NIR spectroscopy for in-line monitoring and controlling content in continuous production of engineered wood products |
CN105334185B (zh) * | 2015-09-14 | 2018-04-03 | 上海创和亿电子科技发展有限公司 | 基于光谱投影判别的近红外模型维护方法 |
CN111579523B (zh) * | 2020-04-17 | 2024-03-08 | 新疆新姿源生物制药有限责任公司 | 一种利用近红外光谱法快速检测孕马尿中结合雌激素含量的方法 |
CN111458442B (zh) * | 2020-06-04 | 2022-08-12 | 新疆新姿源生物制药有限责任公司 | 一种孕马尿中结合雌激素含量的测定方法 |
-
2020
- 2020-12-25 CN CN202011560084.6A patent/CN112748200B/zh active Active
-
2021
- 2021-10-26 US US17/510,667 patent/US20220206019A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
CN112748200A (zh) | 2021-05-04 |
US20220206019A1 (en) | 2022-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112748200B (zh) | 一种用于从孕马尿中获取结合雌激素成分的柱层析工艺的近红外质量监控方法 | |
CN111579523B (zh) | 一种利用近红外光谱法快速检测孕马尿中结合雌激素含量的方法 | |
CN101984343B (zh) | 一种中药大孔树脂分离纯化过程关键点的判别方法 | |
US20020063208A1 (en) | Median filter for liquid chromatography-mass spectrometry data | |
CN105651892B (zh) | 一种基于gc-msqqq的卷烟指纹图谱分析方法 | |
CN102106939A (zh) | 一种六味地黄丸浓缩丸提取浓缩液质量控制方法 | |
CN111024643B (zh) | 一种麻花艽药材品质评价的近红外光谱检测方法 | |
CN111488926A (zh) | 一种基于优化模型的土壤有机质测定方法 | |
CN106353275A (zh) | 一种基于紫外光谱快速测定麻黄药材生物碱的方法 | |
Van Veen et al. | Kalman filtering for data reduction in inductively coupled plasma atomic emission spectrometry | |
CN110806456B (zh) | 一种UPLC-HRMS Profile模式非靶向代谢轮廓数据自动解析的方法 | |
CN109001182B (zh) | 封闭容器中酒精含量的拉曼光谱无损测定方法 | |
CN201449370U (zh) | 多波长核酸蛋白层析分离和检测装置 | |
Seifi et al. | Performance assessment of chemometric resolution methods utilized for extraction of pure components from overlapped signals in gas chromatography–mass spectrometry | |
CN102393378B (zh) | 一种用近红外光谱技术快速检测聚乙烯醇生产中醇解度和聚合度的方法 | |
CN107655849A (zh) | 一种凉茶近红外在线检测方法 | |
Kuligowski et al. | Application of point-to-point matching algorithms for background correction in on-line liquid chromatography–Fourier transform infrared spectrometry (LC–FTIR) | |
Lu et al. | Identification of flurochloridone metabolites in rat urine using liquid chromatography/high resolution mass spectrometry | |
CN115047099A (zh) | 一种银杏叶提取物洗脱过程中总萜内酯含量的快速测定方法 | |
Wiberg et al. | Rapid determination of lidocaine solutions with non-column chromatographic diode array UV spectroscopy and multivariate calibration | |
CN111650142A (zh) | 含铀有机相中亚硝酸浓度分析方法 | |
CN105044019A (zh) | 海产品处理过程中牛磺酸含量的快速检测方法 | |
CN108318447B (zh) | 一种基于近红外光谱技术快速测定吡喹酮对映体组成的方法 | |
Asprilla et al. | Near-infrared spectroscopy: a rapid alternative technique to reducing sugars determination in juice of sugarcane (Saccharum officinarum L.) | |
Cerdá et al. | Multicomponent analysis of environmental matrices |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |