CN110167659B - 利用基质辅助激光解吸/离子化飞行时间质谱仪进行时间对强度分布分析 - Google Patents
利用基质辅助激光解吸/离子化飞行时间质谱仪进行时间对强度分布分析 Download PDFInfo
- Publication number
- CN110167659B CN110167659B CN201780062773.4A CN201780062773A CN110167659B CN 110167659 B CN110167659 B CN 110167659B CN 201780062773 A CN201780062773 A CN 201780062773A CN 110167659 B CN110167659 B CN 110167659B
- Authority
- CN
- China
- Prior art keywords
- time
- flight
- data
- intensity
- spectrometer
- 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
- 238000004458 analytical method Methods 0.000 title claims description 28
- 239000011159 matrix material Substances 0.000 title claims description 13
- 238000003795 desorption Methods 0.000 title claims description 10
- 239000002245 particle Substances 0.000 claims abstract description 147
- 239000000523 sample Substances 0.000 claims abstract description 103
- 238000000034 method Methods 0.000 claims abstract description 60
- 238000012360 testing method Methods 0.000 claims abstract description 44
- 239000013074 reference sample Substances 0.000 claims abstract description 7
- 238000001840 matrix-assisted laser desorption--ionisation time-of-flight mass spectrometry Methods 0.000 claims description 44
- 238000003860 storage Methods 0.000 claims description 31
- 230000008569 process Effects 0.000 claims description 26
- 238000003745 diagnosis Methods 0.000 claims description 22
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 238000013135 deep learning Methods 0.000 claims description 11
- 238000013473 artificial intelligence Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000002955 isolation Methods 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 4
- 241001465754 Metazoa Species 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 11
- 238000004590 computer program Methods 0.000 abstract description 3
- 150000002500 ions Chemical class 0.000 description 27
- 238000009826 distribution Methods 0.000 description 21
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 17
- 201000010099 disease Diseases 0.000 description 16
- 230000000694 effects Effects 0.000 description 15
- 238000005259 measurement Methods 0.000 description 13
- 230000005684 electric field Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 210000004369 blood Anatomy 0.000 description 6
- 239000008280 blood Substances 0.000 description 6
- 150000004676 glycans Chemical class 0.000 description 6
- 108090000623 proteins and genes Proteins 0.000 description 6
- 102000004169 proteins and genes Human genes 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 239000000090 biomarker Substances 0.000 description 4
- 238000013500 data storage Methods 0.000 description 4
- 238000002360 preparation method Methods 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 206010033128 Ovarian cancer Diseases 0.000 description 3
- 206010061535 Ovarian neoplasm Diseases 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 238000000816 matrix-assisted laser desorption--ionisation Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000005381 potential energy Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000005119 centrifugation Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000022811 deglycosylation Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000001819 mass spectrum Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000012887 quadratic function Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000001712 DNA sequencing Methods 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000001124 body fluid Anatomy 0.000 description 1
- 239000010839 body fluid Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000004925 denaturation Methods 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 208000035475 disorder Diseases 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000000752 ionisation method Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000002207 metabolite Substances 0.000 description 1
- 230000000813 microbial effect Effects 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000002414 normal-phase solid-phase extraction Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000013610 patient sample Substances 0.000 description 1
- 238000000059 patterning Methods 0.000 description 1
- 108090000765 processed proteins & peptides Proteins 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 210000002966 serum Anatomy 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000001269 time-of-flight mass spectrometry Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/26—Mass spectrometers or separator tubes
- H01J49/34—Dynamic spectrometers
- H01J49/40—Time-of-flight spectrometers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2462—Approximate or statistical queries
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/40—Information retrieval; Database structures therefor; File system structures therefor of multimedia data, e.g. slideshows comprising image and additional audio data
- G06F16/43—Querying
- G06F16/432—Query formulation
- G06F16/434—Query formulation using image data, e.g. images, photos, pictures taken by a user
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/0027—Methods for using particle spectrometers
- H01J49/0036—Step by step routines describing the handling of the data generated during a measurement
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/02—Details
- H01J49/025—Detectors specially adapted to particle spectrometers
- H01J49/027—Detectors specially adapted to particle spectrometers detecting image current induced by the movement of charged particles
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/02—Details
- H01J49/04—Arrangements for introducing or extracting samples to be analysed, e.g. vacuum locks; Arrangements for external adjustment of electron- or ion-optical components
- H01J49/0409—Sample holders or containers
- H01J49/0418—Sample holders or containers for laser desorption, e.g. matrix-assisted laser desorption/ionisation [MALDI] plates or surface enhanced laser desorption/ionisation [SELDI] plates
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/02—Details
- H01J49/10—Ion sources; Ion guns
- H01J49/16—Ion sources; Ion guns using surface ionisation, e.g. field-, thermionic- or photo-emission
- H01J49/161—Ion sources; Ion guns using surface ionisation, e.g. field-, thermionic- or photo-emission using photoionisation, e.g. by laser
- H01J49/164—Laser desorption/ionisation, e.g. matrix-assisted laser desorption/ionisation [MALDI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
Abstract
公开了一种设备、方法或计算机程序。可以接收样本的光谱仪测试数据以进行处理。该光谱仪测试数据可以包括以穿过飞行管的离子化粒子的时间和强度为单位的飞行时间数据。可将光谱仪测试数据与参考库匹配,以确定样本的特征信息。参考库可包含以光谱仪检测过去检测到的预存参考样本的离子化粒子的时间和强度为单位的光谱仪样本数据。光谱仪参考数据具有匹配与接收到的光谱仪测试数据关联的已知特征。
Description
相关申请的交叉引用
本申请要求2016年8月22日提交的第62/377,768号美国临时专利申请的优先权,其全部内容兹通过引用并入本案。
背景技术
基质辅助激光解吸/离子化飞行时间质谱仪(MALDI-TOF MS) 是光谱仪的一个例子,用于检测样本的属性。例如,样本可以是人血,并且检测到的属性可以是被抽血的人的疾病或其他医学病症。除了疾病诊断之外,光谱仪还有其他应用。
然而,在临床环境中使用MALDI-TOF MS进行疾病诊断存在挑战。一个挑战是质量分析数据的再现性差。诸如,在癌症诊断或疾病状态检查的情况下,当诊断需要除质荷比信息之外的强度信息时,每个质荷比的强度的再现性可能是诊断准确性的关键因素。
再现性差的另一个原因可能来自MALDI-TOF在离子光学中的固有特性。MALDI-TOFMS可以以高灵敏度相对迅速地分析样本。在某些情况下,MALDI-TOF MS可能是相对较差的定量分析工具,因为所有离子化粒子在室入口处不具有相同的初始速度并且/或由于粒子的不同入射角而不在室中行进相同的距离。因此,随着时间的推移可能存在不可避免的数据散布,从而降低了数据再现性。
在某些情况下,在样本发生MALDI-TOF MS离子化后,数据被校准并调整为质荷比(m/z),其可以表示为具有统计变化的质量峰值。当光谱仪测量离子化粒子的飞行时间时,可利用二次等式将测量的飞行时间转换为质荷比。不幸的是,由于光谱仪的瑕疵,离子化粒子的飞行时间存在某些偏差。因此,相同类型的离子化粒子可具有统计学分布的飞行时间测量结果。对于针对来自参考库的匹配布型而优化的光谱仪数据,飞行时间测量结果的统计变化应具有最小化的标准误差。当利用二次等式将飞行时间测量结果转换为质荷比时,飞行时间测量结果的标准偏差受到影响,这可能使得难以使属性与样本关联 (例如,在MALDI-TOF MS中难以用人血液样本产生诊断)。
发明内容
实施例涉及设备、方法或计算机程序。在实施例中,可以接收样本的光谱仪测试数据来进行处理。光谱仪测试数据可包含以穿过光谱仪的飞行管的离子化粒子的时间和强度为单位的飞行时间数据。光谱仪测试数据可以与参考库匹配,以确定样本的特征信息。参考库可包含以光谱仪过去检测到的预存参考样本的离子化粒子的时间和强度为单位的光谱仪样本数据。光谱仪参考数据具有匹配运算与接收到的光谱仪测试数据关联的已知特征。
附图说明
示例性图1示出根据实施例的疾病诊断实验室。
示例性图2是根据实施例的MALDI-TOF MS系统的示意图。
示例性图3是根据实施例将样本处理单元、MALDI-TOF MS单元和诊断单元包含在一个系统中的集成系统的系统图。
示例性图4是根据实施例包含集成于一个系统中的样本处理单元和MALDI-TOF MS单元,而诊断单元作为独立单元提供的集成诊断系统的系统图。
示例性图5示出根据实施例的MALDI板,在该MALDI板上,激光脉冲照射样本板上的样点(spot)。
示例性图6是根据实施例的MALDI-TOF MS的硬件图。
示例性图7示出根据实施例对MALDI室中的样本施照的激光脉冲。
示例性图8示出根据实施例通过激光辐照导致样本粒子离子化。
示例性图9、10A-B和11A-B示出根据实施例的假设匹配分析。
具体实施方式
生物标志物是存在于血液、其他体液或组织中的生物分子,其是正常过程或异常过程的或者病症或疾病的标志。例如,糖蛋白CA-125 是指示存在癌症的生物标志物。因此,通常测量并且评估生物标志物,以认定特定疾病的存在或进展或者观察身体对疾病或病症的治疗的反应好不好。生物标志物在蛋白质、肽、脂质、聚糖或代谢物中的存在或其数量水平的变化能够由质谱仪测量。
在众多类型的质谱仪中,基质辅助激光解吸/离子化飞行时间质谱(MALDI-TOFMS)是采用软离子化技术的分析工具(例如示例性图1-8)。MALDI-TOF MS可在实验室环境下用于快速并且准确地分析生物分子并将其应用扩展到临床领域,诸如,微生物检测和诸如癌症的疾病诊断。例如,如图6所示,样本可埋置于基质607中,激光器601 可以对样本发射激光脉冲。基质607可以吸收来自激光器601的能量,并且样本的分子可以是离子化分子615a-c。然后,靠近基质607的电极605产生的电场可使离子化分子615a-c加速进入飞行管603内。离子化粒子615a-c可以飞过飞行管603并且可以检测器613检测到的离子化粒子615a-c的飞行时间测量结果的分布方式,检测离子化粒子 615a-c。在实施例中,在与参考库匹配(例如,由图3和4所示的诊断单元303和403)之前,不将飞行时间测量结果(例如,图11B中所示的分布布型1109和1111)转换为质荷比对强度(m/z)。
影响数据再现性的因素可以是MALDI-TOF MS系统的测量灵敏度或测量过程。虽然MALDI-TOF MS能够以高灵敏度快速分析样本,但可能存在定量分析复杂化问题,因为由于离子化过程中的瑕疵,检测到的分布布型的相对标准偏差(RSD)可能较高。在实施例中,可以使数据更可再现的方式校准、标准化、归一化并且/或者以其他方式操控光谱仪数据。
示例性图1示出根据实施例的疾病诊断实验室,其中样本处理设施101包含彼此分开的多个样本处理工具、MALDI-TOF MS系统102 和诊断软件系统103。例如,为了提取用于卵巢癌诊断的聚糖,将患者的血清加入到多孔板111内,以进行样本接收过程和蛋白质变性过程112,然后是采用酶的去糖基化过程113。然后,进行蛋白质去除过程114、干燥及离心过程、聚糖提取过程115和点样(spotting)过程116。 MALDI-TOF MS系统102分析所点样本,以产生至少一个聚糖布型。诊断软件103将该样本的聚糖布型与一个或多个预存聚糖布型进行比较,以认定存在卵巢癌和卵巢癌的进展。示例性图2是根据实施例的 MALDI-TOF MS系统的示意图。
示例性图3是根据实施例将样本处理单元、MALDI-TOF MS单元和诊断单元包含在一个系统中的集成系统的系统图。样本可以由所选模块进行组合处理。在样本制备系统301中,在自动样本制备单元311 中,根据诊断或筛选目的,样本经过预定的预设工序。在实施例中,对于聚糖提取,可选择多个处理模块,多个处理模块用于样本接收、蛋白质变性、去糖基化、蛋白质去除、干燥、离心、固相提取和/或点样。在样本制备后,样本装载器312将样本装载到板306上,并在样本干燥器307中干燥。
然后,根据实施例,可将样本提供给具有离子飞行室321和/或高电压真空发生器322的MALDI-TOF MS单元302。MALDI-TOF MS中的处理单元323可以认定离子化粒子(例如图6所示离子化粒子 615a-c)的飞行时间和检测器检测到的相应强度分布。出于疾病诊断的目的,根据实施例,可以重新组织这样获取的飞行时间和强度数据,以建立标准飞行时间列表,在标准飞行时间列表中引入飞行时间分布中心,该中心是强度平衡并且均衡之处。标准飞行时间列表可基于机器准确性和其他相关考虑因素。对每次激光照射存储的光谱数据也可用于建立标准飞行时间列表。然后,诊断单元303可将来自患者样本的光谱与预存的光谱进行比较,并分析两个光谱的布型差异。然后,诊断单元可以认定存在疾病及其进展。
示例性图4是根据实施例包含集成于一个系统中的样本处理单元和MALDI-TOF MS单元,而诊断单元403作为独立单元提供的集成诊断系统的系统图。示例性图4示出根据实施例样本制备单元401和 MALDI-TOF 402集成在一起,而诊断单元403作为独立单元分开布置的集成疾病诊断系统。
在实施例中,诊断单元可以使用参考库。参考库可与诊断单元安置在一起,也可与诊断单元分开。诊断单元可与光谱仪安置在一起,也可与光谱仪分开。在实施例中,参考库可存储于储存装置、基质辅助激光解吸/离子化飞行时间质谱仪(MALDI-TOF MS)、位于光谱仪中的数据储存装置、与光谱仪分开的数据储存装置、通过网络与光谱仪通信的数据储存装置、云储存系统和/或通过互联网连接与光谱仪通信的数据储存装置。
实施例涉及设备、方法或计算机程序。在实施例中,可以接收样本的光谱仪测试数据进行处理(例如,在诊断单元103、303和/或403 处)。光谱仪测试数据可以包含以离子化粒子穿过飞行管的时间和强度为单位的飞行时间数据。可将光谱仪测试数据与参考库匹配来确定样本的特征信息。参考库可包含以光谱仪过去检测到的预存参考样本的离子化粒子的时间和强度为单位的光谱仪样本数据。光谱仪参考数据具有匹配运算与收到的光谱仪测试数据关联的已知特征。
在实施例中,光谱仪测试数据是质谱仪测试数据并且/或者该光谱仪是质谱仪。在实施例中,光谱仪是基质辅助激光解吸/离子化飞行时间质谱仪(MALDI-TOF MS)。
在实施例中,样本包含生物分子并且/或者样本的特征信息包含样本的生物分析信息。生物分析信息可以是人、动物、植物和/或生物体的医学诊断。
示例性图5示出根据实施例的MALDI板,在该MALDI板上,激光脉冲照射样本板上的样点。
示例性图6-7示出根据实施例的光谱仪硬件。激光器601(例如,图6)可对靶707发射激光脉冲701(例如,图7)。靶707可包含样本材料709(例如,人血液和/或其他物质)。样本材料709可吸收来自激光脉冲701的能量并产生具有代表样本材料709的特征和/或属性的离子化粒子715a-e。在实施例中,激光701可以以大约30度的角度投射到样本709上,但是也可以采用所有其他角度。
如图6所示,电极605可产生较接近靶607的电场609。响应于激光器601在靶607处产生的任何离子化粒子615a-c都可被磁性地吸入飞行管603。实际上,离子化粒子615a-c被加速进入飞行管603中并沿着飞行管603的长度611投射并且被检测器613接收。离子化粒子615a可具有比离子化粒子615b-c大的质量。因此,较小的离子化粒子615b-c 可在飞行管603中比较大的离子化粒子615a更快地行进。如在检测器 613所做的测量,在来自激光器601的一个或多个激光脉冲的定时,检测器613都测量每个离子化粒子615a-c的飞行时间。本领域普通技术人员明白可以使用检测器613之外的不同类型的检测器。
由于来自激光器601的每个激光脉冲可以产生较大量的离子化粒子,因此检测器613可以有效测量由一个或多个时间同步激光脉冲产生的所有离子化粒子的飞行时间测量结果的统计分布(例如,图9、 10A-B和11A-B)。离子化粒子615a-c可具有相同或基本相同的动能,因为,与它们在进入飞行管603之前被电极605加速。穿过电极605的范围后,离子化粒子615a-c可以进入基本上没有电场的飞行管603。飞行管603可具有预定长度611。一旦离子化粒子615a-c进入没有电场的飞行管603,离子化粒子615a-c可取决于其质量具有不同速度。与较小的离子化粒子615b-c相比,较大的离子化粒子615a穿过飞行管 603花费较长的时间。
示例性图8示出根据实施例通过辐照激光束801导致样本807离子化。激光束801可以照射样本807。样本807可汽化为离子化粒子(例如,离子化粒子815或817)。当样本807吸收激光801并且部分样本807 离子化时,一些能量传递到离子化粒子815或817。电压施加到电极803 和805,将离子化分子815或817吸入819飞行管中(例如,图6所示的飞行管603)。
在实施例中,匹配过程包括在不将飞行时间数据转换为质荷比对强度布型的情况下,分析飞行时间数据。在实施例中,匹配运算可以处理飞行时间数据,以最小化峰值散布,而当将飞行时间数据转换为质荷比对强度布型时,二次运算不必要地散布峰值。在实施例中,峰值散布的最小化可以最小化飞行时间数据的相对标准偏差。根据实施例,最小化的相对标准偏差可以最大化光谱仪的再现性或准确性中的至少一个。在实施例中,匹配过程利用飞行时间数据来最小化相邻峰值之间不希望的重叠。在实施例中,匹配过程利用飞行时间数据来优化光谱仪的分箱(binning)和/或校准准确性。在实施例中,匹配可以确定样本的特征信息。
在实施例中,可以由人工智能和/或深度学习算法确定匹配。在实施例中,匹配可以优化诊断、筛选或认定医学病症中的至少一个。在实施例中,匹配可以包括如果存在非常相似的飞行时间布型,则判定将光谱仪测试数据与光谱仪参考数据进行匹配。在实施例中,可根据预定阈值执行对非常相似的飞行时间布型进行匹配的判定。在实施例中,可根据动态阈值执行对非常相似的飞行时间布型进行匹配的判定。根据实施例,可利用人工智能和/或深度学习算法确定动态阈值。
在实施例中,离子化粒子可以由激光器产生,该激光器配置成照射靶区,以使放置于靶区中的样本离子化。飞行管的第一端可靠近至少一个电极,该电极配置成使离子化粒子加速进入飞行管中。飞行管的第二对置端可靠近检测器,该检测器测量离子化粒子通过飞行管的飞行时间和离子化粒子的强度。
在实施例中,以时间和强度为单位的飞行时间数据优化对离子化粒子的属性变化的隔离。根据实施例,每个离子化粒子的属性可包含通过至少一个电极的每个离子化粒子的加速效率。根据实施例,每个离子化粒子的属性可包含进入飞行管的至少一个离子化粒子的延迟。根据实施例,每个离子化粒子的属性可包含飞行管内的至少一个离子化粒子的飞行路径的变化。
在实施例中,以时间和强度为单位的飞行时间数据可补偿样本中的物理变化。在实施例中,以时间和强度为单位的飞行时间数据可以优化数据再现性。在实施例中,以时间和强度为单位的飞行时间数据可以将诊断准确性提高到最高。
实施例涉及在匹配分析之前不将飞行时间测量结果转换为质荷比对强度布型(m/z),而使飞行时间测试数据预定质量范围和/或分箱(bin)(例如,与已知特征关联的)最佳地一致。因此,在实施例中,利用飞行时间对强度布型,而不利用转换的质荷比对强度布型,测试数据可与参考库中的准确匹配项最佳地一致。因此,在实施例中,光谱仪能够根据血液样本更准确地诊断疾病(或其他应用)。根据实施例,飞行时间数据及其相应的散布可以直线地而非二次方地增加,并且因此可以相对于任何比质量,产生更对称的时间范围。根据实施例,因为这些范围可更加对称或者至少与指数速率相比,其至少是增大的,所以与高质量和低质量的固定质量分箱校准质量值相比,固定时间分箱可更准确地校准长时间和短时间的飞行时间值。在实施例中,与比质荷比对强度布型相比,飞行时间对强度布型可具有更理想的相对标准偏差(RSD)。
示例性图9、10A-B和11A-B示出根据实施例的假设匹配分析。
例如,根据实施例,当它们所代表的峰值被隔开和/或不同时,例如图9中的匹配分析#1所示,质量分箱可能更有效。在该示例中,测试样本布型峰值901与库样本的质量分箱903匹配。
然而,当质量峰值更靠近在一起和/或它们的相应分箱重叠时,可能会出现问题。该问题示于图10A和10B中的匹配分析#2和匹配分析#3。图10A中的匹配分析#2示出两个相近峰值1005和1007之间重叠的质量分箱1001和1003。在该示例中,峰值1005可以准确地与分箱 1001关联,因为质量对强度峰值1005与分箱1001相交。然而,在该示例中,峰值1007可以与分箱1001和1003两者进行匹配,因此可能产生不准确的诊断分析。同样地,图10B中的测试样本#3可能产生不准确的诊断分析,因为峰值1009与分箱1001和1003两者相交。
在实施例中,如果将飞行时间测试数据转换为质荷比数据,则在不进一步采用数据截断估计技术的情况下,有效地区分用哪个质荷 (m/z)比峰值进行校准可能具有挑战性,这反过来可能具有降低诊断准确性的效果。
示例性图11A-B示出根据实施例在示例匹配分析#4中与参考分箱(例如,分箱1101和1103)关联的飞行时间测试数据(例如,峰值 1105和1107)。与其布型被表达为质荷比对强度的匹配分析#1、#2 和#3不同,将匹配分析#4表达为飞行时间对强度布型。根据实施例,由于匹配分析#4是飞行时间对强度布型,所以用于确定峰值1105和 1107的实际分布1109和111分别具有最小相对标准偏差,因为飞行时间数据不受二次函数的影响,二次函数具有不必要的散布分布1109 和1111。在实施例中,由于分布1109和1111未不必要地散布,所以明显地,关联峰值1105和1107可分别与分箱1101匹配,而不与分箱1103 匹配。
实施例涉及使用时间对强度分析而不是质量对强度分析。实施例可以使用均匀间隔的时间特征。在实施例中,量纲效应(例如,室的直径和长度)和/或优化技术可使用在组间和组内(例如,聚类、归类和分类)具有交叉验证约束的深度学习技术。
实施例可以涉及数据操控,用于补偿MALDI-TOF中的瓶颈效应和非直线或非均匀的行进路径。在实施例中,可以以足以使数据可再现的方式采集并且处理数据,这样其能够用于准确诊断。
在MALDI-TOF系统中,进入飞行时间(TOF)室的离子不能均匀地进入室。进入TOF室的离子不会以相同的距离行进到检测器,因为在到达检测器之前,离子可能以某个角度进入TOF室或以其他方式沿着室壁反弹。这些因素的组合可能导致不一致的质量布型。为了补偿这些限制,实施例可以将TOF室的长度视为检测器处测量的室中的任意粒子的最小行进距离。可基于粒子在内室的弹性平面上未经历反弹、经历一次反弹、二次反弹和/或更多次反弹的假设估计室的有效长度(例如,粒子的预期平均行进距离)。
最大有效长度可用作时间范围或被视为相同质量的质量范围。如果我们假设大多数粒子在室中没有经历反弹,则任何离子粒子的最大可能行进长度可以是L*[1+(D/L)^2]^(1/2),而如果在室内只经历一次反弹并且经历一次以上反弹的离子粒子的数量小到足以忽略不计,则粒子的最大可能行进长度可以是L*[1+4*(D/L)^2]^(1/2)。如果(例如) 室的直径仅比室长度小2.5%,则在该时间范围内测量的任何强度信号被认为是离子粒子中具有相同质量的离子粒子的时间散布Δt如果没有反弹则是1.00031,而如果有一次反弹则是1.00125。
实施例可以将一起移动的离子看作分组,使得以(m/z)c=∑ (m/z)j*Ij计算质量中心(例如,代表性质量),其中代表性质量(m/z)c是质量的加权平均值,每个质量(m/z)j具有其自己的强度Ij。所有质量被看作相同质量的质量的容许质量范围Δm可以是2*t*Δt,其中(Δ t)^2小到足以忽略不计。这可相当于说Δm/m=2*Δt/t。由于在没有电场的室中(1/2)*m*v^2=z*U,且v=L/t,因此可以假设L/t [=SQRT(2*U/(m/z))]对于相同的m/z彼此相等。对于假设所有粒子来自相同质量粒子的任何质量范围,代表性质量mc[(m/z)c=∑(m/z)j*Ij/∑Ij] 可以是由满足L/t=Le/(t+Δt)的量纲等价范围界定的粒子组的加权平均质量,其中Le是有效室长度并且等于粒子在漂移室中的实际行进距离。
室长度L是任意粒子在室中的最短行进距离,而有效管长度Le是 L*SQRT(1+β(D/L)^2),其中D是室的内径,β是行进最短长度L的粒子的量纲常数为0,对于无反弹和穿过中心的粒子(例如,通过管入口的中心中的缝隙进入并到达室环的检测平面的边缘的粒子),量纲常数为1/4,无反弹并且穿透边缘的粒子(例如,粒子进入管边缘的缝隙的粒子)的量纲常数为1,一次反弹的粒子的量纲常数为4,等等。可能希望将室设计成粒子漂移,而不在室的壁上反弹,可以假量纲常数β在1/4和1之间。因为Le/L=(t+Δt)/t=SQRT(1+β(D/L)^2),所以t+ Δt=t*SQRT(1+β(D/L)^2)。与t=L/SQRT(2*U/(m/z))结合可得到t+Δ t=L*SQRT[(1+β(D/L)^2)(m/z)/(2*U)]。
实施例涉及数据操控,以补偿质量确定和校准中的固有不一致性。校准质量数据中可能存在复杂性。因为质荷比取决于时间的二次函数(m/z=k*t^2),所以相对质量散布Δ(m/z)/(m/z)可以是相对时间散布的两倍(Δt/t),并且质量散布Δ(m/z)随着质荷比(m/z)在给定时间间隔或范围Δt内变大而增大,[(m/z)=k*t^2其中k=2*U/(L^2),Δ (m/z)=2*k*t*Δt=2*SQRT[k*(m/z)]*Δt。Δ(m/z)/(m/z)=2*k*t*Δ t/(k*t^2)=2*Δt/t]。
为了消除质量数据的这种不均匀间隔和不对称性质,实施例可以利用飞行时间对强度数据而不利用质荷比对强度数据进行施照。在某些情况下,质量分箱技术或质量校准可产生不充分并且/或不准确的结果,因为其是任意选择的并且分箱有时于相邻峰之间重叠。实施例可利用基于深度学习过程的算法确定最佳质量分箱数或时间分箱数,以用于质量的确定和校准。
在实施例中,可以以足以使数据可再现的方式采集并且/或处理数据,这样其能够用于准确诊断。例如,质谱可表示为曲线图,其中 x轴表示样本的质荷比,y轴表示强度或相对强度与在给定时间或给定时间分箱(转换为质量)到达传感器的粒子数之比。假定通过质谱仪的大多数离子具有+1的电荷,可以认为质荷比(m/z)代表穿过飞行时间(TOF)室的离子的质量。当离子化分子一路穿过室并击中传感器时,机器可以计算分子到达传感器花费多长时间。该时间信息可用于计算传感器检测到的离子质量。该质量可根据势能等于动能的能量守恒定律确定。
例如,下面是势能和动能的公式:
Ep=zU
通过将这些等式合并,我们可以得到下面的关系:
在实施例中,从MALDI-TOF MS获得的m/z对强度数据本质上可不一致。这是因为在由激光器或其他光源离子化后,进入飞行时间 (MALDI-TOF MS中的“TOF”)的粒子在进入时实际上可能不具有均匀的初始速度。另外,它们可不同时进入TOF室。TOF室的入口可由缝隙构成,粒子通过该缝隙进入。因为有如此多的粒子进入一系列窄缝隙,所以粒子不沿缝隙的轴线均匀地排列直线。
在实施例中,代替相似直线竞走的同等开始位置,可能存在粒子在缝隙周围聚集的瓶颈效应(例如,相似于大量跑步者如何沿着狭窄的道路不均匀地开始马拉松比赛)。进入速度和进入定位的这种差异影响MALDI-TOF结果,因为如前所述,m/z可以由下面的等式确定:
L表示TOF管的长度,t表示粒子到达检测器花费的时间。不一致的粒子速度和/或不一致的粒子开始位置(瓶颈效应引起的)可能影响t。粒子在入口处的角度影响粒子在室中的行进长度。离子化粒子 (即使是同一样本)也不可能以均匀时间到达检测器,因为理论上如此。在实施例中,MALDI-TOF MS数据可不呈现为不同峰值或隔开峰值,而是呈现为m/z值范围内的峰值沿x轴的散布。上述瓶颈效应可能是为什么与一致可再现方式能够获得的相异峰值不同,m/z对强度和时间对强度产生质量或时间的散布的核心原因之一。通过将穿过TOF 室行进的离子看作离子分组移动而不是单独离子移动,实施例消除这种固有数据差异。
在实施例中,支持该概念的潜在假设可以是每个离子分组维持由如下表示的平衡状态:
1)对于质量m,(m/z)c=∑(m/z)j*Ij/∑Ij,其中代表性质量(m/z)c是每个质量的质量加权平均值,(m/z)j具有自己的强度Ij。
2)对于一般情况,∑[m(i,j)*I(j)]=m(i,c)*∑I(j),其中m(i,c)表示离子分组的质量中心i,m(i,j)表示质量i周围的第j个质量粒子,并且I (j)表示质量i周围的第j个粒子的强度。
在实施例中,通过产生如下,可对m(i,c)求解:
m(i,c)=∑[m(i,j)*I(j)]/∑I(j)
利用该质量中心信息,实施例可对每个质量获得标准形式的质量分布。该数据可以由MALDI-TOF MS操作获取,在该MALDI-TOF MS 操作中,激光器逐次地、逐点地和/或逐样本地照射样本。通过将各种操作类型的每个质量的标准质量分布进行组合,可将标准质量分布布型构建为离散或连续分布布型。该标准(离散或连续)质量布型可产生较高的诊断准确性,因为根据实施例,在不截断数据的情况下,消除或调整MALDI-TOF MS的许多固有不准确性。
仅TOF室的长度不足以表示上述等式中的L:
m/z=2*U/(L2)*t2
该等式可依赖于TOF室的长度准确地反映粒子穿过TOF室到达检测器的距离的假设。换句话说,在这种假设下,进入飞行时间管的粒子直接向检测器行进。
由于粒子的初始速度不均匀、瓶颈效应和粒子的不同入射角,粒子通常不从入口缝隙直线地穿过TOF室到达检测器。相反,粒子可在到达检测器之前倾斜地进入缝隙并且甚至通常从TOF室壁反弹多次。
因此,TOF室的长度不是离子化粒子行进距离的综合测量值,而是其能够行进的最小距离。在实施例中,表示行进距离的等式是:
L*√1+β*(D/L)2
其中β是常数,表示每个粒子前往检测器的途中经历的反弹(β=0,行进距离为L,β=1/4、1、4等等)。此假设/等式可提高数据准确性,因为与假设的L=飞行时间管长度的假设相比,其可更准确地描述粒子在质谱中经过的一个或多个可能路径。该等式仅表示指出消除离子或粒子在TOF室中前往检测器的非直线路径的方式的一些实施例,因为从本公开能够设想其他实施例。
改进的(或有效)离子路径等式可用于消除离子的非直线行进路径,通过更准确地确定确切质量或代表性质量,消除瓶颈效应的离子分组假设和/或质量强度效应的考虑改善数据分析,因此其改善数据再现性和诊断准确性。与仅直线路径的假设相比,与该离子路径等式和相似概念有关的实施例可以更准确地描述粒子在TOF室中可经历的实际路径。与该离子分组假设有关的实施例可消除当离子进入TOF 室时离子不相同的开始位置。实施例可包含质量产量密度(样本的深度)效应/强度效应分布布型,其还提高数据再现性和诊断准确性。根据实施例,这些数据操控技术既可应用于测试样本,又可应用于数据库中用于参考(比较)的样本,以将两者之间的数据标准化和归一性提高到最高。根据实施例,这样最大化的数据归一性可减少外部因素并可提高数据再现性和诊断准确性。
在MALDI-TOF中,任何比质量都可能不会以相同的m/z值重复。测量的质量值随测量而不同。为了消除这些变化的m/z峰值,出现于特定范围内的质量峰值(围绕每个质量的几个分箱或每个质量水平的 ppm范围)可聚合成代表性质量(有时是质量的中心),并且然后利用一个或多个已知质量的一种或多种校准物质进行校准或调整。在实施例中,为了更加准确,除了标准参考校准物质之外,还可利用DNA 测序数据库对样本计算理论上可获得的蛋白质质量,并且还可用于标准校准基础。
例如,每个质量等级的百万分率(ppm)范围意味着,例如,10,000 道尔顿质量的1000ppm指将10,000*(1-1000*10-6)道尔顿与10,000* (1+1000*10-6)道尔顿范围内的质量计为相同质量。
然而,ppm质量校准技术或其他任意质量间隔存在固有缺陷,并且因此可能由于几个主要原因而产生不准确的数据再现性。例如,(1) 较高质量可产生较宽质量散布,(2)附近峰值的分箱可重叠,并且/ 或(3)用于校准m/z的分箱本征上具有任意性质。
较高的质量可产生较宽的质量散布。MALDI-TOF MS检测器最初不可将数据作为m/z的函数记录可加剧这种现象。相反,MALDI-TOF MS检测器记录粒子的飞行时间,并且然后根据下面的公式将此时间转换为m/z:
通过这种关系,m/z与t2有关,这是二次关系。例如,如果粒子花费1秒钟通过TOF室,则其相应的m/z为作为比较示例,如果另一个粒子花费2秒钟通过TOF室,则其相应的m/z为同样,对于3秒通过TOF,m/z是因为m/z是平方地增大(作为时间的函数)而非线性地增大,所以质量的散布也可呈指数增大(而不是保持相对一致或线性地增大)。
结果,峰值周围的质量分布通常可以是不对称的,因为对于相同的时间散布,较大的质量通常可比较小的质量具有较大的散布。因此,在校准不同质量的粒子的质量时,固定的ppm范围方面不具有相同的效果。
无论是利用ppm分箱还是利用时间分箱来校准质量,无论使用何种范围或判据都可能存在缺陷,因为其可被任意选择。除了经验信息,指出100ppm或86ppm是否是较好的校准分箱范围的信息有限。实施例涉及利用优化的间隔和/或散布计算来减轻任意校准中固有的不确定性,其中利用最高独立性约束优化质量间隔。实施例可采用基于深度学习的算法(例如,人工智能),从而利用算法来确定理想质量或时间分箱,以实现最佳数据再现性和诊断准确性。根据实施例,最佳质量和/或时间分箱可因样本与样本、疾病与疾病以及其他参数而不同。根据实施例,可优化系统以针对每个单独样本/情况选择理想的校准分箱,并且可优化系统以选择一般的理想校准分箱以及其他迭代。根据实施例,可以在测试样本和与其进行比较的数据库样本之间归一化选择的校准分箱。
由于瓶颈效应和引入数据不稳定的其他因素,自然可以看到质量在成分的质量中心周围的散布。在实施例中,当将样本的质量布型数据与用于诊断分析的其他样本中的质量布型数据进行比较时,要做第一件事情是对于其他参考样本数据的相应质量,找到测试样本的标准质量。根据实施例,如果质量匹配不准确,则分析就可能不准确。在实施例中,我们能够将测试样本的每个质量与参考标准化参考质量匹配,但是通常希望测试样本的所有质量与其他样本的参考数据库的相应质量准确匹配是不可能的,也是不切实际的,因为数据布型中质量 (或物质)的校准或调整能够同时改变其他质量值,或者用于匹配的某些质量不能用于理论质量信息。在实施例中,此时我们可以尽可能多地开发样本的标准质量列表。在实施例中,第i个质量m(i)的质量中心mc(i)可以与理论质量一起用于确认和进一步校准。
实施例涉及考虑到质量的密度(强度)效应,而非忽略它们。数据匹配问题的另一个来源是仪器本身。由于MALDI-TOF MS的固有特性,每当激光照射在样本上时,数据布型都可不同。由于仪器的几何形状和进入没有电场的室内的粒子的非均匀动能等,相同质量(物质) 的每个粒子不会一致地到达检测器。MALDI-TOF MS可以将每个质量的每次激光照射的强度值(将代表性质量周围的质量修改为代表性质量)累加或均化为代表性质量的方式显示数据。在该过程中,所看到的质量值可以是忽略质量的强度效应(密度效应)的质量值,其中强度与该物质或质量的粒子数成比例。同时,可以忽略与激光照射关联的质量和强度分布效应,这可能导致关于一个或多个质量的数量或模式的不正确信息。
即使质量校准期间有小的改变也可对诊断的准确性有不利影响。为了解决该问题,实施例建议使用由质量和强度的原始离散值创建的连续分布概念。强度可视为出现次数。可对每个质量创建概率密度函数,并且连续质量分布布型可用于与其他样本的一个或多个数据布型进行比较。这仅是标准质量分布的可能迭代之一。每个质量可具有其自己的标准质量分布,并且可将样本的质量分布数据布型存储为库 DB的一部分。总之,每个标准质量的质量密度函数(分布)、每个标准质量的相应强度密度函数(分布)和强度对标准质量分布布型可以构建为每个样本的库。
根据实施例,这些连续分布布型可用于较大数据的再现性和更准确诊断。
在实施例中,假定通过质谱仪的离子的大部分电荷是+1,则质荷比可表示通过室的离子的质量。当离子化分子通过室并击中传感器时,可以测量该质量。机器计算分子到达传感器花费多长时间。
在飞行时间质谱中,势能等于动能(根据能量守恒定律)。
质荷比与时间的平方相关。根据实施例,假定MALDI-TOF测量离子化分子行进的时间,使用质荷比可能是误导性的。例如,如果 m/z=1,m/z=4,则行进的时间量等于1秒和2秒。如果m/z=81,m/z=85,则这两个单元的行进时间是9秒和因此,在质量-强度布型中,每个峰值都可以是不对称的。根据实施例,看到时间对称征的机会比看到质量对称征的机会更多。实施例涉及能够解释质量分布的集中趋势的布型的对称征,这反过使分析更简单并且/或更快捷。
实施例涉及布型中的部分信号噪声可能来自相邻质量粒子之间的重叠效应。根据实施例,重叠可能是不可避免的,因为相同质量的粒子不能以相同的动能(相同的速度)进入没有电场的漂移区,并且 /或者它们的轨迹可能不同。在实施例中,即使对粒子施加的势能相同(=z*U),粒子的飞行距离仍将介于没有电场的漂移室的长度L与 L*√1+β*(D/L)2之间,其中是与每个粒子的反弹次数和粒子的进入角关联的常数{β=0(最小行程长度L,1/4(中心进入,边缘到达,无反弹),1(边缘进入,边缘到达,无反弹)和4(最大长度,有一次反弹)。在实施例中,如果室内表面可以认为是完全弹性的,并且因此可以认为最多仅一次反弹(换句话说,一次以上反弹到内室上的粒子的数量可忽略不计,并且粒子在检测器中的轨迹可忽略不计)。因此,根据实施例,相对散布Δt/t可以是β*(D/L)2,这又可以是TOF 系统的最大质量间隔。
对所公开的实施例进行各种修改和变更对于本领域技术人员明显且显而易见。所公开的实施例旨在覆盖明显且显而易见的修改和变更,只要这些修改和变更在所附权利要求及其等同的范围内。
Claims (73)
1.一种利用基质辅助激光解吸/离子化飞行时间质谱仪进行时间对强度分布分析的方法,包括:
接收样本的光谱仪测试数据,包括以离子化粒子通过飞行管行进的时间和强度为单位的飞行时间数据;
将所述光谱仪测试数据与参考库匹配,以确定所述样本的特征信息,其中所述参考库包括以至少一个光谱仪检测到的预存参考样本的离子化粒子的飞行时间和强度为单位的光谱仪样本数据,并且其中所述光谱仪参考数据具有所述匹配与接收到的光谱仪测试数据关联的已知特征;
其中在不将所述飞行时间数据转换为质荷比和强度的单位的情况下,所述匹配处理所述飞行时间数据;
其中所述匹配处理所述飞行时间数据,以最小化由于在将飞行时间数据转换为质荷比和强度的单位时的二次运算而与所述飞行时间数据关联的峰值的散布。
2.根据权利要求1所述的方法,其中峰值的所述散布的所述最小化是最小化以飞行时间和强度为单位的所述飞行时间数据的相对标准偏差。
3.根据权利要求2所述的方法,其中所述最小化的相对标准偏差是最大化以飞行时间和强度为单位的飞行时间数据的再现性或准确性中的至少一个。
4.根据权利要求1所述的方法,其中所述匹配处理所述飞行时间数据,以最小化与以飞行时间和强度为单位的所述飞行时间数据关联的相邻峰值之间的重叠。
5.根据权利要求4所述的方法,其中所述匹配处理所述飞行时间数据,以优化分箱或校准准确性中的至少一个。
6.根据权利要求1所述的方法,其中:
所述光谱仪测试数据是质谱仪测试数据;
所述至少一个光谱仪是至少一个质谱仪。
7.根据权利要求6所述的方法,其中所述至少一个光谱仪包括基质辅助激光解吸/离子化飞行时间质谱仪(MALDI-TOF MS)。
8.根据权利要求1所述的方法,其中所述匹配确定所述样本的特征信息。
9.根据权利要求1所述的方法,其中:
所述样本包含生物分子;
所述样本的所述特征信息包括所述样本的生物分析信息。
10.根据权利要求1所述的方法,其中所述匹配由人工智能或深度学习算法中的至少一个确定。
11.根据权利要求1所述的方法,其中所述匹配包括如果存在基本相似的飞行时间对强度布型,则判定将所述光谱仪测试数据与所述光谱仪参考数据进行匹配。
12.根据权利要求11所述的方法,其中根据预定阈值执行对匹配所述基本相似的飞行时间对强度布型的判定。
13.根据权利要求11所述的方法,其中根据动态阈值执行对匹配所述基本相似的飞行时间对强度布型的判定。
14.根据权利要求13所述的方法,其中所述动态阈值由人工智能或深度学习算法中的至少一个确定。
15.根据权利要求1所述的方法,其中:
所述离子化粒子由激光器产生,所述激光器配置成照射靶区以离子化放置于所述靶区中的样本;
所述飞行管的第一端靠近至少一个电极,所述电极配置成使离子化粒子加速进入所述飞行管内;以及
所述飞行管的第二对置端靠近检测器,所述检测器测量所述离子化粒子通过所述飞行管的飞行时间和所述离子化粒子的强度。
16.根据权利要求1所述的方法,其中以飞行时间和强度为单位的所述飞行时间数据优化对所述离子化粒子的属性的变化的隔离。
17.根据权利要求16所述的方法,其中每个所述离子化粒子的所述属性包括每个离子化粒子通过至少一个电极的加速效率。
18.根据权利要求16所述的方法,其中每个所述离子化粒子的所述属性包括所述离子化粒子中的至少一个进入所述飞行管的延迟。
19.根据权利要求16所述的方法,其中每个所述离子化粒子的所述属性包括所述离子化粒子中的至少一个在所述飞行管内的飞行路径的变化。
20.根据权利要求1所述的方法,其中以飞行时间和强度为单位的所述飞行时间数据补偿所述样本中的物理变化。
21.根据权利要求1所述的方法,其中以飞行时间和强度为单位的所述飞行时间数据优化数据再现性。
22.根据权利要求1所述的方法,其中以飞行时间和强度为单位的所述飞行时间数据使诊断准确性最大化。
23.根据权利要求1所述的方法,其中所述参考库存储于储存装置中。
24.一种利用基质辅助激光解吸/离子化飞行时间质谱仪进行时间对强度分布分析的设备,包括:
至少一个处理器;
接收单元,所述接收单元配置成利用所述至少一个处理器接收样本的光谱仪测试数据,所述样本的光谱仪测试数据包括以离子化粒子通过飞行管行进的时间和强度为单位的飞行时间数据;
匹配单元,所述匹配单元配置成将所述光谱仪测试数据与参考库匹配,以利用所述至少一个处理器确定所述样本的特征信息,其中所述参考库包括以至少一个光谱仪检测到的预存参考样本的离子化粒子的时间和强度为单位的光谱仪样本数据,并且其中所述光谱仪参考数据具有所述匹配单元与接收到的光谱仪测试数据关联的已知特征,
其中在不将所述飞行时间数据转换为质荷比和强度的单位的情况下,所述匹配处理所述飞行时间数据;
其中所述匹配处理所述飞行时间数据,以最小化由于在将飞行时间数据转换为质荷比和强度的单位时的二次运算而与所述飞行时间数据关联的峰值的散布。
25.根据权利要求24所述的设备,其中峰值的所述散布的所述最小化是最小化以飞行时间和强度为单位的所述飞行时间数据的相对标准偏差。
26.根据权利要求25所述的设备,其中所述最小化的相对标准偏差是最大化以飞行时间和强度为单位的飞行时间数据的再现性或准确性中的至少一个。
27.根据权利要求24所述的设备,其中所述匹配处理所述飞行时间数据,以最小化与以飞行时间和强度为单位的所述飞行时间数据关联的相邻峰值之间的重叠。
28.根据权利要求27所述的设备,其中所述匹配处理所述飞行时间数据,以优化分箱或校准准确性中的至少一个。
29.根据权利要求24所述的设备,其中:
所述光谱仪测试数据是质谱仪测试数据;
所述至少一个光谱仪是至少一个质谱仪。
30.根据权利要求29所述的设备,其中所述至少一个光谱仪包括基质辅助激光解吸/离子化飞行时间质谱仪(MALDI-TOF MS)。
31.根据权利要求24所述的设备,其中所述匹配确定所述样本的特征信息。
32.根据权利要求24所述的设备,其中:
所述样本包含生物分子;
所述样本的所述特征信息包括所述样本的生物分析信息。
33.根据权利要求32所述的设备,其中所述生物分析信息是人、动物、植物或生物体中的至少一个的医学诊断。
34.根据权利要求24所述的设备,其中所述匹配由人工智能或深度学习算法中的至少一个确定。
35.根据权利要求24所述的设备,所述匹配优化诊断、筛选或认定医学病症中的至少一个。
36.根据权利要求24所述的设备,其中所述匹配包括如果存在基本相似的飞行时间对强度布型,则判定将所述光谱仪测试数据与所述光谱仪参考数据进行匹配。
37.根据权利要求36所述的设备,其中根据预定阈值执行对匹配所述基本相似的飞行时间对强度布型的判定。
38.根据权利要求36所述的设备,其中根据动态阈值执行对匹配所述基本相似的飞行时间对强度布型的判定。
39.根据权利要求38所述的设备,其中所述动态阈值由人工智能或深度学习算法中的至少一个确定。
40.根据权利要求24所述的设备,其中:
所述离子化粒子由激光器产生,所述激光器配置成照射靶区以离子化放置于所述靶区中的样本;
所述飞行管的第一端靠近至少一个电极,所述电极配置成使离子化粒子加速进入所述飞行管内;以及
所述飞行管的第二对置端靠近检测器,所述检测器测量所述离子化粒子通过所述飞行管的飞行时间和所述离子化粒子的强度。
41.根据权利要求24所述的设备,其中以飞行时间和强度为单位的所述飞行时间数据优化对所述离子化粒子的属性的变化的隔离。
42.根据权利要求41所述的设备,其中每个所述离子化粒子的所述属性包括每个离子化粒子通过至少一个电极的加速效率。
43.根据权利要求41所述的设备,其中每个所述离子化粒子的所述属性包括所述离子化粒子中的至少一个进入所述飞行管的延迟。
44.根据权利要求41所述的设备,其中每个所述离子化粒子的所述属性包括所述离子化粒子中的至少一个在所述飞行管内的飞行路径的变化。
45.根据权利要求24所述的设备,其中以飞行时间和强度为单位的所述飞行时间数据补偿所述样本中的物理变化。
46.根据权利要求24所述的设备,其中以飞行时间和强度为单位的所述飞行时间数据优化数据再现性。
47.根据权利要求24所述的设备,其中以飞行时间和强度为单位的所述飞行时间数据使诊断准确性最大化。
48.根据权利要求24所述的设备,其中所述参考库存储于储存装置中。
49.一种利用基质辅助激光解吸/离子化飞行时间质谱仪进行时间对强度分布分析的计算机可读存储介质,包括存储有计算机可读程序代码的计算机可读硬件储存装置,所述程序代码含有计算机系统的一个或多个处理器可执行的指令,以实现对目标评估损坏的方法,所述方法包括:
接收样本的光谱仪测试数据,其包括以离子化粒子通过飞行管行进的时间和强度为单位的飞行时间数据;
将所述光谱仪测试数据与参考库进行匹配,以确定所述样本的特征信息,其中所述参考库包括以至少一个光谱仪检测到的预存参考样本的离子化粒子的时间和强度为单位的光谱仪样本数据,并且其中所述光谱仪参考数据具有所述匹配与接收到的光谱仪测试数据关联的已知特征;
其中在不将所述飞行时间数据转换为质荷比和强度的单位的情况下,所述匹配处理所述飞行时间数据;
其中所述匹配处理所述飞行时间数据,以最小化由于在将飞行时间数据转换为质荷比和强度的单位时的二次运算而与所述飞行时间数据关联的峰值的散布。
50.根据权利要求49所述的计算机可读存储介质,其中峰值的所述散布的所述最小化是最小化以飞行时间和强度为单位的所述飞行时间数据的相对标准偏差。
51.根据权利要求50所述的计算机可读存储介质,其中所述最小化的相对标准偏差是最大化以飞行时间和强度为单位的飞行时间数据的再现性或准确性中的至少一个。
52.根据权利要求49所述的计算机可读存储介质,其中所述匹配处理所述飞行时间数据,以最小化与以飞行时间和强度为单位的所述飞行时间数据关联的相邻峰值之间的重叠。
53.根据权利要求52所述的计算机可读存储介质,其中所述匹配处理所述飞行时间数据,以优化分箱或校准准确性中的至少一个。
54.根据权利要求49所述的计算机可读存储介质,其中:
所述光谱仪测试数据是质谱仪测试数据;
所述至少一个光谱仪是至少一个质谱仪。
55.根据权利要求54所述的计算机可读存储介质,其中所述至少一个光谱仪包括基质辅助激光解吸/离子化飞行时间质谱仪(MALDI-TOF MS)。
56.根据权利要求49所述的计算机可读存储介质,其中所述匹配确定所述样本的特征信息。
57.根据权利要求49所述的计算机可读存储介质,其中:
所述样本包含生物分子;
所述样本的所述特征信息包括所述样本的生物分析信息。
58.根据权利要求57所述的计算机可读存储介质,其中所述生物分析信息是人、动物、植物或生物体中的至少一个的医学诊断。
59.根据权利要求49所述的计算机可读存储介质,其中所述匹配由人工智能或深度学习算法中的至少一个确定。
60.根据权利要求49所述的计算机可读存储介质,所述匹配优化诊断、筛选或认定医学病症中的至少一个。
61.根据权利要求49所述的计算机可读存储介质,其中所述匹配包括如果存在基本相似的飞行时间对强度布型,则判定将所述光谱仪测试数据与所述光谱仪参考数据进行匹配。
62.根据权利要求61所述的计算机可读存储介质,其中根据预定阈值执行对匹配所述基本相似的飞行时间对强度布型的判定。
63.根据权利要求61所述的计算机可读存储介质,其中根据动态阈值执行对匹配所述基本相似的飞行时间对强度布型的判定。
64.根据权利要求63所述的计算机可读存储介质,其中所述动态阈值由人工智能或深度学习算法中的至少一个确定。
65.根据权利要求49所述的计算机可读存储介质,其中:
所述离子化粒子由激光器产生,所述激光器配置成照射靶区以离子化放置于所述靶区中的样本;
所述飞行管的第一端靠近至少一个电极,所述电极配置成使离子化粒子加速进入所述飞行管内;以及
所述飞行管的第二对置端靠近检测器,所述检测器测量所述离子化粒子通过所述飞行管的飞行时间和所述离子化粒子的强度。
66.根据权利要求49所述的计算机可读存储介质,其中以飞行时间和强度为单位的所述飞行时间数据优化对所述离子化粒子的属性的变化的隔离。
67.根据权利要求66所述的计算机可读存储介质,其中每个所述离子化粒子的所述属性包括每个离子化粒子通过至少一个电极的加速效率。
68.根据权利要求66所述的计算机可读存储介质,其中每个所述离子化粒子的所述属性包括所述离子化粒子中的至少一个进入所述飞行管的延迟。
69.根据权利要求66所述的计算机可读存储介质,其中每个所述离子化粒子的所述属性包括所述离子化粒子中的至少一个在所述飞行管内的飞行路径的变化。
70.根据权利要求49所述的计算机可读存储介质,其中以飞行时间和强度为单位的所述飞行时间数据补偿所述样本中的物理变化。
71.根据权利要求49所述的计算机可读存储介质,其中以飞行时间和强度为单位的所述飞行时间数据优化数据再现性。
72.根据权利要求49所述的计算机可读存储介质,其中以飞行时间和强度为单位的所述飞行时间数据使诊断准确性最大化。
73.根据权利要求49所述的计算机可读存储介质,其中所述参考库存储于储存装置中。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201662377768P | 2016-08-22 | 2016-08-22 | |
US62/377,768 | 2016-08-22 | ||
PCT/US2017/047833 WO2018039136A1 (en) | 2016-08-22 | 2017-08-21 | Time versus intensity distribution analysis using a matrix-assisted laser desorption/ionization time-of-flight mass spectrometer |
US15/682,166 US10497553B2 (en) | 2016-08-22 | 2017-08-21 | Time versus intensity distribution analysis using a matrix-assisted laser desorption/ionization time-of-flight mass spectrometer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110167659A CN110167659A (zh) | 2019-08-23 |
CN110167659B true CN110167659B (zh) | 2023-01-31 |
Family
ID=61763766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780062773.4A Active CN110167659B (zh) | 2016-08-22 | 2017-08-21 | 利用基质辅助激光解吸/离子化飞行时间质谱仪进行时间对强度分布分析 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10497553B2 (zh) |
EP (1) | EP3493895A4 (zh) |
KR (1) | KR102256075B1 (zh) |
CN (1) | CN110167659B (zh) |
WO (1) | WO2018039136A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111665291A (zh) * | 2020-04-27 | 2020-09-15 | 浙江迪谱诊断技术有限公司 | 一种飞行时间核酸质谱激光时序算法及控制方法 |
CN113219042A (zh) * | 2020-12-03 | 2021-08-06 | 深圳市步锐生物科技有限公司 | 一种用于人体呼出气体中各成分分析检测的装置及其方法 |
CN116660358B (zh) * | 2023-08-01 | 2023-11-24 | 浙江迪谱诊断技术有限公司 | 一种高分辨飞行时间质谱检测方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7091479B2 (en) | 2000-05-30 | 2006-08-15 | The Johns Hopkins University | Threat identification in time of flight mass spectrometry using maximum likelihood |
US6683301B2 (en) * | 2001-01-29 | 2004-01-27 | Analytica Of Branford, Inc. | Charged particle trapping in near-surface potential wells |
JP2005536714A (ja) * | 2001-11-13 | 2005-12-02 | カプリオン ファーマシューティカルズ インコーポレーティッド | 質量強度プロファイリングシステムおよびその使用法 |
US20030027231A1 (en) * | 2002-01-15 | 2003-02-06 | Bryden Wayne A. | Methods for using mass spectrometry to identify and classify filamentous fungi, yeasts, molds and pollen |
JP4620446B2 (ja) * | 2004-12-24 | 2011-01-26 | 株式会社日立ハイテクノロジーズ | 質量分析方法、質量分析システム、診断システム、検査システム及び質量分析プログラム |
CA2787504A1 (en) * | 2010-10-07 | 2012-04-12 | The Government Of United States Of America As Represented By The Secreta Ry Of The Department Of Health And Human Services Center For Disease Con | Use of detector response curves to optimize settings for mass spectrometry |
DE102012102875B4 (de) * | 2011-04-04 | 2024-04-18 | Wisconsin Alumni Research Foundation | Vorläuferauswahl mit einem Artificial-Intelligence-Algorithmus erhöht Abdeckung und Reproduzierbarkeit von proteomischen Proben |
NL2009015C2 (en) * | 2012-04-10 | 2013-10-15 | Biosparq B V | Method for classification of a sample on the basis of spectral data, method for creating a database and method for using this database, and corresponding computer program, data storage medium and system. |
US10197551B2 (en) * | 2014-05-30 | 2019-02-05 | ASTA, Co. Ltd. | Integrated disease diagnostic system using matrix-assisted laser desorption/ionization time-of-flight mass spectrometer |
CN106663589B (zh) * | 2014-08-29 | 2019-06-14 | 生物梅里埃有限公司 | 具有延迟时间变化的maldi-tof质谱仪及相关方法 |
-
2017
- 2017-08-21 US US15/682,166 patent/US10497553B2/en active Active
- 2017-08-21 EP EP17844222.4A patent/EP3493895A4/en active Pending
- 2017-08-21 KR KR1020197008144A patent/KR102256075B1/ko active IP Right Grant
- 2017-08-21 WO PCT/US2017/047833 patent/WO2018039136A1/en unknown
- 2017-08-21 CN CN201780062773.4A patent/CN110167659B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
KR20190075048A (ko) | 2019-06-28 |
US10497553B2 (en) | 2019-12-03 |
EP3493895A4 (en) | 2020-07-15 |
KR102256075B1 (ko) | 2021-05-25 |
EP3493895A1 (en) | 2019-06-12 |
WO2018039136A1 (en) | 2018-03-01 |
CN110167659A (zh) | 2019-08-23 |
US20180053645A1 (en) | 2018-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210151306A1 (en) | Shot-to-shot sampling using a matrix-assisted laser desorption/ionization time-of-flight mass spectrometer | |
CN105321795B (zh) | 改进的成像质谱分析方法和装置 | |
RU2633797C2 (ru) | Способ классификации образца на основании спектральных данных, способ создания базы данных, способ использования этой базы данных и соответсвующие компьютерная программа, носитель данных и система | |
CN110167659B (zh) | 利用基质辅助激光解吸/离子化飞行时间质谱仪进行时间对强度分布分析 | |
US9043164B2 (en) | Method of generating a mass spectrum having improved resolving power | |
US10153145B2 (en) | Method of mass spectrometry and a mass spectrometer | |
US11423331B2 (en) | Analytical data analysis method and analytical data analyzer | |
US20220059332A1 (en) | Multiplex charge detection mass spectrometry | |
US10197551B2 (en) | Integrated disease diagnostic system using matrix-assisted laser desorption/ionization time-of-flight mass spectrometer | |
US20240087864A1 (en) | Analysis system using matrix-assisted laser desorption/ionization time-of-flight mass spectrometer | |
US6365893B1 (en) | Internal calibration of time to mass conversion in time-of-flight mass spectrometry | |
KR20190054994A (ko) | 말디토프 질량 분석에 의한 항생제 내성 판별 장치 및 방법 | |
CN110178200B (zh) | 采用基质辅助激光解吸/离子化飞行时间质谱仪逐次发射抽样 | |
CN109155230B (zh) | 对敞开式电离斑的测量与验证 | |
CN116635975A (zh) | 用于将样品定时引入质谱仪中的方法和系统 | |
US20210217526A1 (en) | Cancer diagnosis using optimal clustering with successive deconvolution | |
US7653493B1 (en) | Proteomic sample analysis and systems therefor | |
Dexter | Developing computational methods for fundamentals and metrology of mass spectrometry imaging | |
US10359440B2 (en) | Mass spectrometry method and apparatus for clinical diagnostic applications | |
WO2015063821A1 (ja) | 分析システム | |
JPWO2023037536A5 (zh) | ||
Buddi | Statistical signal processing of ESI-TOF-MS for biomarker discovery |
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 |