CN114566292B - 应用内源性物质预测肾移植药物剂量的模型的构建方法 - Google Patents
应用内源性物质预测肾移植药物剂量的模型的构建方法 Download PDFInfo
- Publication number
- CN114566292B CN114566292B CN202210111071.3A CN202210111071A CN114566292B CN 114566292 B CN114566292 B CN 114566292B CN 202210111071 A CN202210111071 A CN 202210111071A CN 114566292 B CN114566292 B CN 114566292B
- Authority
- CN
- China
- Prior art keywords
- model
- sample
- metabolic
- acetonitrile
- detected
- 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
- 239000003814 drug Substances 0.000 title claims abstract description 25
- 229940079593 drug Drugs 0.000 title claims abstract description 19
- 239000000126 substance Substances 0.000 title claims abstract description 12
- 210000003734 kidney Anatomy 0.000 title abstract description 8
- 238000002054 transplantation Methods 0.000 title abstract description 8
- 238000010276 construction Methods 0.000 title description 3
- 230000002503 metabolic effect Effects 0.000 claims abstract description 64
- 238000000034 method Methods 0.000 claims abstract description 50
- 230000000694 effects Effects 0.000 claims abstract description 45
- 238000012216 screening Methods 0.000 claims abstract description 33
- 239000003018 immunosuppressive agent Substances 0.000 claims abstract description 30
- 229960003444 immunosuppressant agent Drugs 0.000 claims abstract description 23
- 230000001861 immunosuppressant effect Effects 0.000 claims abstract description 22
- 210000004369 blood Anatomy 0.000 claims abstract description 10
- 239000008280 blood Substances 0.000 claims abstract description 10
- WEVYAHXRMPXWCK-UHFFFAOYSA-N Acetonitrile Chemical group CC#N WEVYAHXRMPXWCK-UHFFFAOYSA-N 0.000 claims description 69
- 150000002500 ions Chemical class 0.000 claims description 43
- 239000000523 sample Substances 0.000 claims description 41
- 230000036961 partial effect Effects 0.000 claims description 26
- 210000001744 T-lymphocyte Anatomy 0.000 claims description 21
- 239000012634 fragment Substances 0.000 claims description 18
- DDRJAANPRJIHGJ-UHFFFAOYSA-N creatinine Chemical compound CN1CC(=O)NC1=N DDRJAANPRJIHGJ-UHFFFAOYSA-N 0.000 claims description 16
- 238000007477 logistic regression Methods 0.000 claims description 14
- 239000002207 metabolite Substances 0.000 claims description 12
- 239000002904 solvent Substances 0.000 claims description 12
- 230000004044 response Effects 0.000 claims description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 10
- FERIUCNNQQJTOY-UHFFFAOYSA-N Butyric acid Chemical compound CCCC(O)=O FERIUCNNQQJTOY-UHFFFAOYSA-N 0.000 claims description 10
- 239000013062 quality control Sample Substances 0.000 claims description 10
- 229940109239 creatinine Drugs 0.000 claims description 8
- 210000002966 serum Anatomy 0.000 claims description 8
- 238000010219 correlation analysis Methods 0.000 claims description 6
- 239000004475 Arginine Substances 0.000 claims description 5
- AGPKZVBTJJNPAG-WHFBIAKZSA-N L-isoleucine Chemical compound CC[C@H](C)[C@H](N)C(O)=O AGPKZVBTJJNPAG-WHFBIAKZSA-N 0.000 claims description 5
- 229930182844 L-isoleucine Natural products 0.000 claims description 5
- QIVBCDIJIAJPQS-VIFPVBQESA-N L-tryptophane Chemical compound C1=CC=C2C(C[C@H](N)C(O)=O)=CNC2=C1 QIVBCDIJIAJPQS-VIFPVBQESA-N 0.000 claims description 5
- UMOXFSXIFQOWTD-LURJTMIESA-N N(2)-succinyl-L-arginine Chemical compound NC(=N)NCCC[C@@H](C(O)=O)NC(=O)CCC(O)=O UMOXFSXIFQOWTD-LURJTMIESA-N 0.000 claims description 5
- QIVBCDIJIAJPQS-UHFFFAOYSA-N Tryptophan Natural products C1=CC=C2C(CC(N)C(O)=O)=CNC2=C1 QIVBCDIJIAJPQS-UHFFFAOYSA-N 0.000 claims description 5
- JAZBEHYOTPTENJ-JLNKQSITSA-N all-cis-5,8,11,14,17-icosapentaenoic acid Chemical compound CC\C=C/C\C=C/C\C=C/C\C=C/C\C=C/CCCC(O)=O JAZBEHYOTPTENJ-JLNKQSITSA-N 0.000 claims description 5
- ODKSFYDXXFIFQN-UHFFFAOYSA-N arginine Natural products OC(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-N 0.000 claims description 5
- JAZBEHYOTPTENJ-UHFFFAOYSA-N eicosapentaenoic acid Natural products CCC=CCC=CCC=CCC=CCC=CCCCC(O)=O JAZBEHYOTPTENJ-UHFFFAOYSA-N 0.000 claims description 5
- 229960005135 eicosapentaenoic acid Drugs 0.000 claims description 5
- 235000020673 eicosapentaenoic acid Nutrition 0.000 claims description 5
- 229940125721 immunosuppressive agent Drugs 0.000 claims description 5
- 229960000310 isoleucine Drugs 0.000 claims description 5
- 230000003907 kidney function Effects 0.000 claims description 5
- 229910052757 nitrogen Inorganic materials 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- BPYKTIZUTYGOLE-IFADSCNNSA-N Bilirubin Chemical compound N1C(=O)C(C)=C(C=C)\C1=C\C1=C(C)C(CCC(O)=O)=C(CC2=C(C(C)=C(\C=C/3C(=C(C=C)C(=O)N\3)C)N2)CCC(O)=O)N1 BPYKTIZUTYGOLE-IFADSCNNSA-N 0.000 claims description 4
- QJJXYPPXXYFBGM-LFZNUXCKSA-N Tacrolimus Chemical compound C1C[C@@H](O)[C@H](OC)C[C@@H]1\C=C(/C)[C@@H]1[C@H](C)[C@@H](O)CC(=O)[C@H](CC=C)/C=C(C)/C[C@H](C)C[C@H](OC)[C@H]([C@H](C[C@H]2C)OC)O[C@@]2(O)C(=O)C(=O)N2CCCC[C@H]2C(=O)O1 QJJXYPPXXYFBGM-LFZNUXCKSA-N 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- BDAGIHXWWSANSR-UHFFFAOYSA-N methanoic acid Natural products OC=O BDAGIHXWWSANSR-UHFFFAOYSA-N 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 108090000623 proteins and genes Proteins 0.000 claims description 4
- 238000003908 quality control method Methods 0.000 claims description 4
- 229960001967 tacrolimus Drugs 0.000 claims description 4
- QJJXYPPXXYFBGM-SHYZHZOCSA-N tacrolimus Natural products CO[C@H]1C[C@H](CC[C@@H]1O)C=C(C)[C@H]2OC(=O)[C@H]3CCCCN3C(=O)C(=O)[C@@]4(O)O[C@@H]([C@H](C[C@H]4C)OC)[C@@H](C[C@H](C)CC(=C[C@@H](CC=C)C(=O)C[C@H](O)[C@H]2C)C)OC QJJXYPPXXYFBGM-SHYZHZOCSA-N 0.000 claims description 4
- 102000009027 Albumins Human genes 0.000 claims description 3
- 108010088751 Albumins Proteins 0.000 claims description 3
- JWODVNNWDACELT-UHFFFAOYSA-N C(C)(=O)O.COC=1C=C2C=CNC2=CC1 Chemical compound C(C)(=O)O.COC=1C=C2C=CNC2=CC1 JWODVNNWDACELT-UHFFFAOYSA-N 0.000 claims description 3
- ODKSFYDXXFIFQN-BYPYZUCNSA-P L-argininium(2+) Chemical compound NC(=[NH2+])NCCC[C@H]([NH3+])C(O)=O ODKSFYDXXFIFQN-BYPYZUCNSA-P 0.000 claims description 3
- 238000003657 Likelihood-ratio test Methods 0.000 claims description 3
- 238000007664 blowing Methods 0.000 claims description 3
- 238000010828 elution Methods 0.000 claims description 3
- 238000005534 hematocrit Methods 0.000 claims description 3
- 230000003908 liver function Effects 0.000 claims description 3
- 102000004169 proteins and genes Human genes 0.000 claims description 3
- 239000000243 solution Substances 0.000 claims description 3
- OSWFIVFLDKOXQC-UHFFFAOYSA-N 4-(3-methoxyphenyl)aniline Chemical compound COC1=CC=CC(C=2C=CC(N)=CC=2)=C1 OSWFIVFLDKOXQC-UHFFFAOYSA-N 0.000 claims description 2
- 239000007864 aqueous solution Substances 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 235000019253 formic acid Nutrition 0.000 claims description 2
- 238000012417 linear regression Methods 0.000 claims description 2
- 238000004949 mass spectrometry Methods 0.000 claims description 2
- 238000001819 mass spectrum Methods 0.000 claims description 2
- 230000001376 precipitating effect Effects 0.000 claims description 2
- 238000004885 tandem mass spectrometry Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 208000030453 Drug-Related Side Effects and Adverse reaction Diseases 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 8
- 230000002829 reductive effect Effects 0.000 description 4
- 229940046731 calcineurin inhibitors Drugs 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000001506 immunosuppresive effect Effects 0.000 description 3
- 238000001558 permutation test Methods 0.000 description 3
- 230000001225 therapeutic effect Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- COCNDHOPIHDTHK-UHFFFAOYSA-N 5-methoxyindole-3-acetic acid Chemical compound COC1=CC=C2NC=C(CC(O)=O)C2=C1 COCNDHOPIHDTHK-UHFFFAOYSA-N 0.000 description 2
- 102000012192 Cystatin C Human genes 0.000 description 2
- 108010061642 Cystatin C Proteins 0.000 description 2
- 206010062016 Immunosuppression Diseases 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000024924 glomerular filtration Effects 0.000 description 2
- -1 i.e. Substances 0.000 description 2
- 230000004060 metabolic process Effects 0.000 description 2
- 239000012071 phase Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 150000003384 small molecules Chemical class 0.000 description 2
- 206010067484 Adverse reaction Diseases 0.000 description 1
- PMATZTZNYRCHOR-CGLBZJNRSA-N Cyclosporin A Chemical compound CC[C@@H]1NC(=O)[C@H]([C@H](O)[C@H](C)C\C=C\C)N(C)C(=O)[C@H](C(C)C)N(C)C(=O)[C@H](CC(C)C)N(C)C(=O)[C@H](CC(C)C)N(C)C(=O)[C@@H](C)NC(=O)[C@H](C)NC(=O)[C@H](CC(C)C)N(C)C(=O)[C@H](C(C)C)NC(=O)[C@H](CC(C)C)N(C)C(=O)CN(C)C1=O PMATZTZNYRCHOR-CGLBZJNRSA-N 0.000 description 1
- 108010036949 Cyclosporine Proteins 0.000 description 1
- 206010071602 Genetic polymorphism Diseases 0.000 description 1
- 108010022337 Leucine Enkephalin Proteins 0.000 description 1
- 238000010220 Pearson correlation analysis Methods 0.000 description 1
- 206010052779 Transplant rejections Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000006838 adverse reaction Effects 0.000 description 1
- 150000001412 amines Chemical class 0.000 description 1
- 229940024606 amino acid Drugs 0.000 description 1
- 150000001413 amino acids Chemical class 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000035 biogenic effect Effects 0.000 description 1
- 239000012496 blank sample Substances 0.000 description 1
- 239000012482 calibration solution Substances 0.000 description 1
- 150000001735 carboxylic acids Chemical class 0.000 description 1
- 229960001265 ciclosporin Drugs 0.000 description 1
- 229930182912 cyclosporin Natural products 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000037213 diet Effects 0.000 description 1
- 235000005911 diet Nutrition 0.000 description 1
- 235000014113 dietary fatty acids Nutrition 0.000 description 1
- 150000002016 disaccharides Chemical class 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000857 drug effect Effects 0.000 description 1
- 238000000132 electrospray ionisation Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011067 equilibration Methods 0.000 description 1
- 229930195729 fatty acid Natural products 0.000 description 1
- 239000000194 fatty acid Substances 0.000 description 1
- 150000004665 fatty acids Chemical class 0.000 description 1
- HQVFCQRVQFYGRJ-UHFFFAOYSA-N formic acid;hydrate Chemical compound O.OC=O HQVFCQRVQFYGRJ-UHFFFAOYSA-N 0.000 description 1
- 239000005556 hormone Substances 0.000 description 1
- 229940088597 hormone Drugs 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000543 intermediate Substances 0.000 description 1
- URLZCHNOLZSCCA-UHFFFAOYSA-N leu-enkephalin Chemical compound C=1C=C(O)C=CC=1CC(N)C(=O)NCC(=O)NCC(=O)NC(C(=O)NC(CC(C)C)C(O)=O)CC1=CC=CC=C1 URLZCHNOLZSCCA-UHFFFAOYSA-N 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 238000002705 metabolomic analysis Methods 0.000 description 1
- 230000001431 metabolomic effect Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 150000002772 monosaccharides Chemical class 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000003285 pharmacodynamic effect Effects 0.000 description 1
- 230000002974 pharmacogenomic effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- ZAHRKKWIAAJSAO-UHFFFAOYSA-N rapamycin Natural products COCC(O)C(=C/C(C)C(=O)CC(OC(=O)C1CCCCN1C(=O)C(=O)C2(O)OC(CC(OC)C(=CC=CC=CC(C)CC(C)C(=O)C)C)CCC2C)C(C)CC3CCC(O)C(C3)OC)C ZAHRKKWIAAJSAO-UHFFFAOYSA-N 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 229960002930 sirolimus Drugs 0.000 description 1
- QFJCIRLUMZQUOT-HPLJOQBZSA-N sirolimus Chemical compound C1C[C@@H](O)[C@H](OC)C[C@@H]1C[C@@H](C)[C@H]1OC(=O)[C@@H]2CCCCN2C(=O)C(=O)[C@](O)(O2)[C@H](C)CC[C@H]2C[C@H](OC)/C(C)=C/C=C/C=C/[C@@H](C)C[C@@H](C)C(=O)[C@H](OC)[C@H](O)/C(C)=C/[C@@H](C)C(=O)C1 QFJCIRLUMZQUOT-HPLJOQBZSA-N 0.000 description 1
- 230000000391 smoking effect Effects 0.000 description 1
- 150000003431 steroids Chemical class 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 210000002700 urine Anatomy 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- 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
- G16H70/00—ICT specially adapted for the handling or processing of medical references
- G16H70/40—ICT specially adapted for the handling or processing of medical references relating to drugs, e.g. their side effects or intended usage
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
- G16B15/30—Drug targeting using structural data; Docking or binding prediction
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Physics & Mathematics (AREA)
- Medicinal Chemistry (AREA)
- Pharmacology & Pharmacy (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Toxicology (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Biophysics (AREA)
- Crystallography & Structural Chemistry (AREA)
- Epidemiology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明涉及一种应用内源性物质预测肾移植药物剂量的模型的构建方法，方法包括：收集个体使用免疫抑制剂后的血液样本和使用效果，根据血液样本中的代谢特征和临床特征结合使用效果筛选出与使用效果强相关的代谢特征和临床特征，然后使用筛选后的代谢特征和临床特征建立免疫抑制剂的药物剂量预测模型。本发明的方法所构建的模型能够预测个体所适用的给药量，避免潜在的药物不良反应，有助于个体化给药。
Description
技术领域
本发明属于药物预测技术领域，具体涉及一种应用内源性物质预测肾移植药物剂量的模型的构建方法。
背景技术
钙调神经磷酸酶抑制剂 (CNIs) 自20世纪80年代问世以来，已经成为目前临床常用的治疗干预手段，改善了移植患者的结局，提高了移植患者的生存率。目前，钙调神经磷酸酶抑制剂主要有环孢素、他克莫司、西罗莫司等。
然而，目前免疫抑制剂的治疗窗窄，不同个体间药动学和药效学差异大。剂量过高容易免疫抑制过度，增加严重感染的风险；剂量过低，容易免疫抑制不足，从而增加移植器官发生排斥反应的风险。因此给临床应用带来一定的困扰。
随着药物代谢组学技术的发展，通过分析给药前血浆或尿液代谢概况来预测药物不良反应或药物疗效已经成为研究的热点。药物吸收，分布，代谢和消除是体内的基本过程，可能受遗传 (如遗传多态性)，生理状态和环境因素 (如饮食，年龄和吸烟) 的影响。内源性小分子作为基因与环境相互作用的最终产物，其代谢特征能够兼容更多的信息，而另一方面，部分免疫抑制剂主要通过抑制T淋巴细胞发挥免疫抑制作用，因此T淋巴细胞含量可以表征免疫抑制剂的药物生效情况。
现有技术中，若想检测药物生效情况需要监测免疫抑制剂的血药浓度，尤其是全血谷浓度，及时调整给药剂量，增加疗效，降低不良反应。但是，监测谷浓度需要严格采集定时血样，这对于临床繁忙的工作环境难以操作和实现。因此，需要一种方法能够建立一个预测免疫抑制剂优选剂量的模型。
发明内容
基于现有技术中存在的上述缺点和不足，本发明的目的之一是至少解决现有技术中存在的上述问题之一或多个，换言之，本发明的目的之一是提供满足前述需求之一或多个的一种应用内源性物质预测肾移植药物剂量的模型的构建方法。
为了达到上述发明目的，本发明采用以下技术方案：
一种应用内源性物质预测肾移植药物剂量的模型的构建方法，包括以下步骤：
S1、收集待测样本，将所述待测样本中的蛋白质沉淀；所述待测样本为个体使用免疫抑制剂后的血液样本；
S2、将所述个体使用所述免疫抑制剂的效果分为高反应组、正常组和低反应组，所述效果根据T淋巴细胞含量判断；
S3、使用质谱仪分析所述待测样本的代谢特征，所述代谢特征包括所述待测样本中的若干种内源性代谢物质含量；
S4、通过偏最小二乘法分析所述代谢特征与所述效果分组的相关性，筛选出与所述效果分组具有强相关性的代谢特征；
S5、获取所述个体使用所述免疫抑制剂后的临床特征，所述临床特征包括所述个体的年龄、BMI、肾功能指标、肝功能指标、红细胞压积、白蛋白含量参数；
S6、通过分析所述临床特征与所述效果分组的相关性，筛选出与所述效果分组具有强相关性的临床特征；
S7、根据相关性分析和多项逻辑回归模型，以筛选后的所述代谢特征及所述临床特征与所述效果分组的相关性建立所述免疫抑制剂的药物剂量预测模型。
作为优选方案，所述使用质谱仪分析所述待测样本的代谢特征，具体包括如下步骤：
S31、以随机方式分析所述待测样本，首先连续检测3-5个质控样本，然后每10个所述待测样本插入1个质控样本，且同时检测正离子和负离子模式；
S32、在所述同时检测正离子和负离子模式下，通过预设数据库获得待测样本的碎片离子峰作为所述代谢特征。
作为优选方案，所述获得待测样本的碎片离子峰作为所述代谢特征之后还包括数据预处理过程S221：
S321、对所述代谢特征进行峰提取、峰对齐、缺失值筛选、缺失值填补，并基于质控样本进行局部线性回归校正。
作为优选方案，所述步骤S4具体为通过两次独立的偏最小二乘法分析所述待测样本的代谢特征数据与所述效果的相关性进行筛选，包括以下步骤：
S41、以所有检测鉴定到的碎片离子峰为基础变量，目标效应值为结局变量，建立第一个独立的偏最小二乘法模型;
S42、通过第一个独立的偏最小二乘法模型，获得所述碎片离子峰对第一个独立的偏最小二乘法模型的首轮贡献指数；
S43、筛选首轮贡献指数大于1.5的碎片离子峰作为候选变量，建立第二个独立的偏最小二乘法模型；
S44、通过第二个独立的偏最小二乘法模型，获得所述候选变量对第二个独立的偏最小二乘法模型的二轮贡献指数；
S45、筛选二轮贡献指数大于1的碎片离子峰作为目标变量，以目标变量离子峰对应的内源性代谢物质作为筛选后的所述代谢特征。
作为优选方案，所述步骤S45之后还包括步骤S46：
S46、通过相关性分析确定所述代谢特征与效果分组的相关性，以相关性的显著水平小于0.05为标准，缩小并优化代谢特征的范围。
作为优选方案，所述步骤S7具体包括
S71、根据多项逻辑回归模型，从所述初选临床特征中进一步筛选得到筛选后的所述临床特征；通过多项逻辑回归，选用似然比检验实现变量筛选，将显著水平小于0.05的所述初选临床特征纳入最终回归方程。
S72、根据所述初选临床特征建立多项逻辑回归方程。
作为优选方案，所述方法还包括步骤S8、使用正交最小二乘法判别模型验证模型的预测能力；具体包括：
S81、使用筛选后的所述代谢特征及所述临床特征形成正交最小二乘法判别模型的协变量；
S82、根据所述协变量构建正交最小二乘法判别模型，验证所述模型的预测能力。
作为优选方案，所述使用质谱仪分析所述待测样本的代谢特征的参数具体为：
运用UPLC-Q-TOF/MSE质谱仪分析，飞行质谱采用ESI源，毛细管电压2.5 kv，锥孔电压40 v，离子源温度100 ℃，脱溶剂温度450 ℃，锥孔反吹氮气流速50 L/h，脱溶剂气(氮气) 流速800 L/h；扫描时间0.2 s，MS/MS二级质谱分析时碰撞能量设定为20-30 ev，质量扫描范围为50-1200 m/z；
色谱柱选用100 mm×2.1 mm、1.8 μm的 HSS T3柱；流动相为溶剂A和溶剂B，所述溶剂A为0.1 %甲酸水溶液，所述溶剂B为乙腈溶液，所述溶剂B中乙腈的梯度洗脱设置如下：0-0.5 min，5 % 乙腈；0.5-2 min，5-20 % 乙腈；2-3.5 min，20-27.5 % 乙腈；3.5-4 min，27.5-70 % 乙腈；4-7.5 min，70-75 %乙腈；7.5-8.5,75-95 % 乙腈，8.5-13.5 min，95 %乙腈，13.5-16 min，95-5 %乙腈，16-18 min，5 % 乙腈；流速为0.3 mL/min，柱温为35 ℃，进样量2 μL。
作为优选方案，所述免疫抑制剂为他克莫司；所述若干种内源性代谢物质包括血清肌酐、胆红素原、L-异亮氨酸、5-甲氧基吲哚乙酸、二十碳五烯酸、N2-琥珀酰精氨酸、色氨酸精氨酸和丁酸。
另一方面，本发明提供一种免疫抑制剂剂量预测的方法，使用如上述任一项的方法所构建的模型，包括如下步骤：
S1、将目标样本的参数代入所述模型，获得疗效预测结果；
S2、根据所述疗效预测结果迭代药物剂量，获得不同疗效预测对应的药物剂量，实现剂量预测
本发明与现有技术相比，有益效果是：
（1）本发明的方法所构建的模型能够预测个体所适用的给药量，避免潜在的药物不良反应，有助于个体化给药；
（2）本发明的方法所构建的模型，通过筛选代谢特征和临床特征，提高了变量的准确程度，具有优秀的预测能力，减少了过拟合风险。
附图说明
图1是本发明实施例的第一次偏最小二乘法模型的相关图；
图2是本发明实施例的第一次偏最小二乘法模型的载荷图；
图3是本发明实施例的第二次偏最小二乘法模型的相关图；
图4是本发明实施例的置换检验的结果分析图；
图5是本发明的正交最小二乘法判别模型的分析结果。
具体实施方式
为了更清楚地说明本发明实施例，下面将对照附图说明本发明的具体实施方式。显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图，并获得其他的实施方式。
实施例：本实施例的一种应用内源性物质预测肾移植药物剂量的模型的构建方法，首先收集个体使用免疫抑制剂后采集的血液样本，然后根据血液样本中代谢特征、临床特征与使用效果的相关性筛选出具有强相关的代谢特征和临床特征，并将二者结合进行预测模型的构建。
具体的，本实施例的方法首先进行步骤S1、收集待测样本，将待测样本中的蛋白质沉淀；待测样本为个体使用免疫抑制剂后采集的血液样本。在本实施例中，免疫抑制剂可以具体为他克莫司。
S2、将个体使用免疫抑制剂的效果分为高反应组、正常组和低反应组，效果根据T淋巴细胞含量判断。
T淋巴细胞的含量是免疫抑制剂疗效的体现，在移植排斥反应中起着重要的作用。血清肌酐、胱抑素C和肾小球率过滤是评价肾移植后肾脏功能的重要指标，经过临床研究后发现T淋巴细胞与肾脏功能指标中的血清肌酐、胱抑素C和肾小球滤过率有显著相关性，表明T淋巴细胞含量可表示肾移植患者在免疫抑制剂发挥免疫抑制作用后的疗效结果。
根据免疫抑制剂疗效的情况可以将用药效果分为三组，低反应组即所测T淋巴细胞百分比高于参考范围值上限 (n=16)、正常反应组即所测T淋巴细胞百分比在参考范围值内 (n=73) 和高反应组即所测T淋巴细胞百分比低于参考范围值下限 (n=20)。其中在高反应组更容易发生延迟复功 (P<0.05)，低反应组与正常组在急性排斥反应发生率方面无显著差异，但是在低反应组急性排斥反应发生率远高于正常组。
S3、使用质谱仪分析待测样本的代谢特征，代谢特征主要包括待测样本中的若干种内源性代谢物质含量。内源性代谢物质即机体代谢过程中所形成的产物或中间产物，具体包括氨基酸类、脂类、脂肪酸、类固醇、激素、生物胺、单糖、双糖、羧酸等
具体来说，药物可以影响代谢特征，而临床特征可以影响药物效果，因此结合临床特征和代谢特征可以更好的反应药物疗效。
进一步的，使用质谱仪分析待测样本可以包括以下步骤：
S31、以随机方式分析待测样本，首先连续检测3-5个质控样本，然后每10个待测样本插入1个质控样本，且同时检测正离子和负离子模式；
详细的，待测样本的分析可以使用以下方案：非靶向代谢组学分析在UPLC-Q-TOF质谱仪上进行，质谱采用MSE continuum模式，离子源：电喷雾电离模式 (ESI)，正离子和负离子模式均检测，毛细管电压2.5 kv，锥孔电压40 v，离子源温度100 ℃，脱溶剂温度450℃，锥孔反吹氮气流速50 L/h，脱溶剂气 (氮气) 流速800 L/h。扫描时间0.2 s；扫描间隔0.015 s，低通道碰撞能量6 eV；高通道碰撞能量20-30 ev。质量扫描范围为50-1200 m/z。
液相条件如下：样品在HSS T3柱 (100 mm×2.1 mm，1.8 μm) 上分离。流动相为0.1 %甲酸水溶液 (溶剂A) 和乙腈 (溶剂B)，梯度洗脱0-0.5 min，5 % B；0.5-2 min，5-20% B；2-3.5 min，20-27.5 % B；3.5-4 min，27.5-70 % B；4-7.5 min，70-75 % B；7.5-8.5,75-95 % B，8.5-13.5 min，95 % B，13.5-16 min，95-5 % B，16-18 min，5 % B。流速为0.3mL/min，柱温为35 ℃，进样量2 μL，自动进样器温度为10 ℃。
亮氨酸-脑啡肽 (5 ng/mL) 作为校正液随行于在整个分析中，用于判断仪器质量轴是否有偏差。此外，3-5个空白样品在正式进样前运行，用于平衡色谱柱并分析样品引入的背景离子。待测样本以随机方式于UPLC-QTOF-MS分析。待测样品分析前连续检测3-5个质控样品，然后每10个样本后插入1个质控样本。
S32、在同时检测正离子和负离子模式下，通过预设数据库获得待测样本的碎片离子峰作为代谢特征，在每个样品中，测量的离子峰面积表示该代谢物的相对强度。
为了确保数据的可靠性，本实施例中还进一步的对数据了进行以下预处理：S321、对代谢特征进行峰提取(选择合适的加合离子提取峰)、删除样本中缺失值超过80 %的代谢物，用k-最近邻算法 (KNN) 填充缺失值，QC的相对标准偏差（RSD）小于30 %，并根据质控样品进行局部多项式拟合。经过数据预处理后，质控样品(QC) 峰的变异系数（CV）得到降低，更加可靠。
得到代谢特征的数据后，进行步骤S4、通过偏最小二乘法分析代谢特征与效果分组的相关性，筛选出与效果分组具有强相关性的代谢特征；
具体的，为了选择与T淋巴细胞显著相关的代谢物，本实施例选择连续两次偏最小二乘法分析，从而步骤S4可以包括以下过程：
S41、以所有检测鉴定到的碎片离子峰为基础变量，目标效应值为结局变量，建立第一个独立的偏最小二乘法模型：以所有检测鉴定到的碎片离子峰为基础，建立初始PLS模型，PLS模型的相关图如图1所示。以每个点代表一个个体，X变量是所有检测到的离子峰，Y变量是T淋巴细胞百分比。
S42、通过第一个独立的偏最小二乘法模型，获得碎片离子峰对第一个独立的偏最小二乘法模型的首轮贡献指数；
模型的载荷图如图2所示，载荷图显示了每个预测变量和响应变量之间的相关性。图2右上角的X变量与T淋巴细胞百分比呈正相关，而左下角的变量与T淋巴细胞百分比呈负相关，VIP是代表X变量对PLS模型贡献的指数。
S43、筛选首轮贡献指数大于1.5的碎片离子峰作为候选变量，建立第二个独立的偏最小二乘法模型，来筛选具有较强相关性的代谢特征 (VIP>1.5) 以预测个体T淋巴细胞百分比。
S44、通过第二个独立的偏最小二乘法模型，获得候选变量对第二个独立的偏最小二乘法模型的二轮贡献指数；如图3所示，以上一轮筛选的代谢特征为X变量再次建立PLS模型，其模型置换检验的结果分析图如图4所示。通过置换检验 (n=200) 进行模型内部验证，以检验模型是否过度拟合，若没有则模型可靠。
进一步的，一般非靶向代谢组学采集的数据比较庞大，可以获得成千上万的代谢特征，如果仅使用一种筛选方法，例如仅使用采用贡献指数大于1的一次筛选，得到筛选后的数据仍然比较多，而两次独立筛选可以产生叠加效应，通过数据迭代分析，获得两次独立筛选的筛选条件，可以获得更加精确的目标代谢物。因而，通过步骤S42进行第一轮筛选以后，从成千上万个（非靶向代谢组学的）代谢特征中筛选得到数百个具有相关性的代谢特征，接下来需要进一步缩减范围，得到更具相关性的代谢特征。重复上一轮的筛选过程，但VIP更加严格，进行步骤S45、筛选二轮贡献指数大于1的碎片离子峰作为目标变量，以目标变量离子峰所对应的内源性代谢物质作为筛选后的代谢特征，经过这一步筛选的代谢特征具有更强的相关性，两次筛选后得到的代谢特征可以缩减到一百种左右，进一步缩减后续物质鉴定的工作量，提高预测效率及精确度。
进一步的，为了加强所筛选出代谢特征与T淋巴细胞百分比的相关性，还可以进行步骤S46、通过相关性分析确定代谢特征与效果分组的相关性，以相关性的显著水平小于0.05为标准，缩小并优化代谢特征的范围。
筛选完具有强相关的代谢特征，还需要筛选具有相关性的临床特征，因此进行步骤S5、采集个体使用免疫抑制剂后的临床特征，临床特征包括年龄、BMI、肾功能指标、肝功能指标、红细胞压积、白蛋白含量参数；临床信息也能够解释一部分个体疗效差异，因此也需要考虑被纳入到剂量预测模型中。
S6、通过相关性分析筛选出与效果分组具有强相关性的临床特征；在本实施例中，可以使用pearson相关分析进行筛选，筛选得到的临床特征为年龄和血清肌酐，这两个临床特征与T淋巴细胞百分比显著相关(P<0.05)。
筛选完代谢特征和临床特征后，进行步骤S7、根据相关性分析和多项逻辑回归模型，以筛选后的代谢特征及临床特征与效果分组的相关性建立免疫抑制剂的药物剂量预测模型
进一步的，步骤S7包括以下步骤：
S71、根据多项逻辑回归模型，从初选临床特征中进一步筛选得到筛选后的临床特征；通过多项逻辑回归，选用似然比检验实现变量筛选，将显著水平小于0.05的初选临床特征纳入最终回归方程。
在本实施例中，初选的代谢特征和临床特征包括血清肌酐、胆红素原、L-异亮氨酸、5-甲氧基吲哚乙酸、二十碳五烯酸、N2-琥珀酰精氨酸、色氨酸精氨酸和丁酸。
S72、根据初选临床特征建立多项逻辑回归方程。
拟合情况如下表所示：
多项逻辑回归模型拟合情况
在本实施例中根据上述多项逻辑回归模型，可以得到以下方程：
G1=LOG[P(正常反应组)/P(低反应组)]=-3.875-0.113(日剂量)-0.753(Log10 (血清肌酐))+0.126(胆红素原)-0.432(L-异亮氨酸)-0.068(5-甲氧基吲哚乙酸)+0.079(二十碳五烯酸)+0.013 (N2-琥珀酰精氨酸)+0.147(色氨酸精氨酸)-0.11(丁酸)
G2=LOG[P(高反应组)/P(低反应组)]=-13.151-0.588(日剂量)+4.461(Log10 (血清肌酐))+0.301(胆红素原)-1.362(L-异亮氨酸)+1.156(5-甲氧基吲哚乙酸)+0.522(二十碳五烯酸)+0.014 (N2-琥珀酰精氨酸)+0.305(色氨酸精氨酸)-0.117(丁酸)
G3=0 (对照组)
根据三组G1、G2、G3的值，带入以下公式，最终可得到三个疗效组的概率。
P1=exp(G1)/[exp(G1)+exp(G2)+exp(G3)]
P2=exp(G2)/[exp(G1)+exp(G2)+exp(G3)]
P3=exp(G3)/[exp(G1)+exp(G2)+exp(G3)]
上述方程能够较好地解释解释个体服用免疫抑制剂服用后T淋巴细胞含量的差异，进而根据疗效可以递归反馈进行剂量调整。
作为上述模型预测准确度的验证，本实施例的方法还包括步骤S8、使用正交最小二乘法判别模型验证模型的预测能力。
步骤S8可以具体包括：S71、使用筛选后的所述代谢特征及所述临床特征形成正交最小二乘法判别模型的协变量；
S72、根据所述协变量构建正交最小二乘法判别模型。
根据所选特异性内源物质和临床特征建立正交最小二乘法判别模型，用正交最小二乘法判别模型能够上述筛选的代谢特征和临床特征的预测能力。模型的分析结果用图5表示，每个圆圈代表一个个体，得分图显示了聚类的概况。所选择的变量将高低疗效组清楚地分开，R2Y为0.574，Q2为0.415。同时，执行了200次迭代的置换检验，以防模型过度拟合结果表明，模型是稳定的，没有过度拟合的风险。由正交最小二乘法判别模型可以发现高反应组和低反应组中的这几个变量有显著差异，即该模型可以用于免疫抑制剂给药剂量的预测。
本实施例还提供一种免疫抑制剂药物剂量预测的方法，应用上述方法所构建的模型进行免疫抑制剂给药剂量的预测，具体包括：将给药后目标样本中与模型相对应的特征参数代入该模型，获得T细胞含量的疗效预测结果。然后根据该T细胞含量偏高/低迭代给药剂量，从而获得不同疗效所对应的给药剂量，从而实现给药剂量的精准预测。
应当说明的是，上述实施例仅是对本发明的优选实施例及原理进行了详细说明，对本领域的普通技术人员而言，依据本发明提供的思想，在具体实施方式上会有改变之处，而这些改变也应视为本发明的保护范围。
Claims (8)
1.一种应用内源性物质预测肾移植药物剂量的模型的构建方法，其特征在于，包括以下步骤：
S1、收集待测样本，将所述待测样本中的蛋白质沉淀；所述待测样本为个体使用免疫抑制剂后的血液样本；
S2、将所述个体使用所述免疫抑制剂的效果分为高反应组、正常组和低反应组，所述效果根据T淋巴细胞含量判断；
S3、使用质谱仪分析所述待测样本的代谢特征，所述代谢特征包括所述待测样本中的若干种内源性代谢物质含量；
S4、通过两次独立的偏最小二乘法分析所述待测样本的代谢特征数据与所述效果的相关性进行筛选，筛选出与所述效果分组具有强相关性的代谢特征；
S5、获取所述个体使用所述免疫抑制剂后的临床特征，所述临床特征包括所述个体的年龄、BMI、肾功能指标、肝功能指标、红细胞压积、白蛋白含量参数；
S6、通过分析所述临床特征与所述效果分组的相关性，筛选出与所述效果分组具有强相关性的临床特征；
S71、根据多项逻辑回归模型，从初选临床特征中进一步筛选得到筛选后的所述临床特征；通过多项逻辑回归，选用似然比检验实现变量筛选，将显著水平小于0.05的所述初选临床特征纳入最终回归方程；
S72、根据所述初选临床特征建立多项逻辑回归方程；
S73、根据所述多项逻辑回归方程建立药物剂量预测模型；
所述步骤S4具体包括以下步骤：
S41、以所有检测鉴定到的碎片离子峰为基础变量，目标效应值为结局变量，建立第一个独立的偏最小二乘法模型;
S42、通过第一个独立的偏最小二乘法模型，获得所述碎片离子峰对第一个独立的偏最小二乘法模型的首轮贡献指数；
S43、筛选首轮贡献指数大于1.5的碎片离子峰作为候选变量，建立第二个独立的偏最小二乘法模型；
S44、通过第二个独立的偏最小二乘法模型，获得所述候选变量对第二个独立的偏最小二乘法模型的二轮贡献指数；
S45、筛选二轮贡献指数大于1的碎片离子峰作为目标变量，以目标变量离子峰对应的内源性代谢物质作为筛选后的所述代谢特征。
2.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法，其特征在于，所述使用质谱仪分析所述待测样本的代谢特征，具体包括如下步骤：
S31、以随机方式分析所述待测样本，首先连续检测3-5个质控样本，然后每10个所述待测样本插入1个质控样本，且同时检测正离子和负离子模式；
S32、在所述同时检测正离子和负离子模式下，通过预设数据库获得待测样本的碎片离子峰作为所述代谢特征。
3.如权利要求2所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法，其特征在于，所述获得待测样本的碎片离子峰作为所述代谢特征之后还包括数据预处理过程S221：
S321、对所述代谢特征进行峰提取、峰对齐、缺失值筛选、缺失值填补，并基于质控样本进行局部线性回归校正。
4.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法，其特征在于，所述步骤S45之后还包括步骤S46：
S46、通过相关性分析确定所述代谢特征与效果分组的相关性，以相关性的显著水平小于0.05为标准，缩小并优化代谢特征的范围。
5.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法，其特征在于，所述方法还包括步骤S8、使用正交最小二乘法判别模型验证模型的预测能力；具体包括：
S81、使用筛选后的所述代谢特征及所述临床特征形成正交最小二乘法判别模型的协变量；
S82、根据所述协变量构建正交最小二乘法判别模型，验证所述模型的预测能力。
6.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法，其特征在于，所述使用质谱仪分析所述待测样本的代谢特征的参数具体为：
运用UPLC-Q-TOF/MSE质谱仪分析，飞行质谱采用ESI源，毛细管电压2.5 kv，锥孔电压40 v，离子源温度100 ℃，脱溶剂温度450 ℃，锥孔反吹氮气流速50 L/h，脱溶剂气流速800L/h；扫描时间0.2 s，MS/MS二级质谱分析时碰撞能量设定为20-30 ev，质量扫描范围为50-1200 m/z；
色谱柱选用100 mm×2.1 mm、1.8 μm的 HSS T3柱；流动相为溶剂A和溶剂B，所述溶剂A为0.1 %甲酸水溶液，所述溶剂B为乙腈溶液，所述溶剂B中乙腈的梯度洗脱设置如下：0-0.5min，5 % 乙腈；0.5-2 min，5-20 % 乙腈；2-3.5 min，20-27.5 % 乙腈；3.5-4 min，27.5-70% 乙腈；4-7.5 min，70-75 %乙腈；7.5-8.5,75-95 % 乙腈，8.5-13.5 min，95 % 乙腈，13.5-16 min，95-5 %乙腈，16-18 min，5 % 乙腈；流速为0.3 mL/min，柱温为35 ℃，进样量2 μL。
7.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法，其特征在于，所述免疫抑制剂为他克莫司；所述若干种内源性代谢物质包括血清肌酐、胆红素原、L-异亮氨酸、5-甲氧基吲哚乙酸、二十碳五烯酸、N2-琥珀酰精氨酸、色氨酸精氨酸和丁酸。
8.一种药物剂量预测的方法，其特征在于，使用如权利要求1-7任一项所述的方法所构建的模型，包括如下步骤：
S1、将目标样本的参数代入所述模型，获得疗效预测结果；
S2、根据所述疗效预测结果迭代药物剂量，获得不同疗效预测对应的药物剂量，实现剂量预测。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210111071.3A CN114566292B (zh) | 2022-01-29 | 2022-01-29 | 应用内源性物质预测肾移植药物剂量的模型的构建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210111071.3A CN114566292B (zh) | 2022-01-29 | 2022-01-29 | 应用内源性物质预测肾移植药物剂量的模型的构建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114566292A CN114566292A (zh) | 2022-05-31 |
| CN114566292B true CN114566292B (zh) | 2022-11-11 |
Family
ID=81714166
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210111071.3A Active CN114566292B (zh) | 2022-01-29 | 2022-01-29 | 应用内源性物质预测肾移植药物剂量的模型的构建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114566292B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010016110A1 (ja) * | 2008-08-05 | 2010-02-11 | 富士通株式会社 | 重回帰分析による予測モデルの作成方法、作成システムおよび作成プログラム |
| CN107247095A (zh) * | 2017-05-25 | 2017-10-13 | 武汉大学 | 一种基于代谢组学分析的用于抗氧化衰老药物筛选的大鼠模型的构建方法 |
| CN108280321A (zh) * | 2018-01-29 | 2018-07-13 | 沈阳药科大学 | 基于细胞代谢轮廓分析构建抗肿瘤候选化合物作用机制预测模型的方法 |
| CN112786145A (zh) * | 2021-03-04 | 2021-05-11 | 华中科技大学同济医学院附属协和医院 | 一种器官移植患者他克莫司用药剂量精准预测方法 |
| CN113035369A (zh) * | 2021-03-10 | 2021-06-25 | 浙江大学 | 一种肾移植抗感染药物剂量预测模型的构建方法 |
| CN113889221A (zh) * | 2021-12-08 | 2022-01-04 | 首都医科大学宣武医院 | 重症肌无力患者他克莫司治疗剂量精准选择的系统及其应用 |
| CN113930495A (zh) * | 2021-10-18 | 2022-01-14 | 上海市第一人民医院 | 一种肝移植术后他克莫司初始剂量的预测模型及其个体化应用 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101361038B1 (ko) * | 2013-02-21 | 2014-02-12 | 경북대학교 산학협력단 | 신장 이식 후 환자의 뇨 시료로부터 사구체여과율을 예측하는 방법 |
| CN109668984A (zh) * | 2019-01-30 | 2019-04-23 | 山西医科大学第医院 | 一种基于代谢组学的喉癌血清判别模型的构建方法 |
| CN113671078B (zh) * | 2021-08-18 | 2023-03-03 | 郑州大学第一附属医院 | 一种基于代谢组学的脓毒症预后模型建立方法 |
-
2022
- 2022-01-29 CN CN202210111071.3A patent/CN114566292B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010016110A1 (ja) * | 2008-08-05 | 2010-02-11 | 富士通株式会社 | 重回帰分析による予測モデルの作成方法、作成システムおよび作成プログラム |
| CN107247095A (zh) * | 2017-05-25 | 2017-10-13 | 武汉大学 | 一种基于代谢组学分析的用于抗氧化衰老药物筛选的大鼠模型的构建方法 |
| CN108280321A (zh) * | 2018-01-29 | 2018-07-13 | 沈阳药科大学 | 基于细胞代谢轮廓分析构建抗肿瘤候选化合物作用机制预测模型的方法 |
| CN112786145A (zh) * | 2021-03-04 | 2021-05-11 | 华中科技大学同济医学院附属协和医院 | 一种器官移植患者他克莫司用药剂量精准预测方法 |
| CN113035369A (zh) * | 2021-03-10 | 2021-06-25 | 浙江大学 | 一种肾移植抗感染药物剂量预测模型的构建方法 |
| CN113930495A (zh) * | 2021-10-18 | 2022-01-14 | 上海市第一人民医院 | 一种肝移植术后他克莫司初始剂量的预测模型及其个体化应用 |
| CN113889221A (zh) * | 2021-12-08 | 2022-01-04 | 首都医科大学宣武医院 | 重症肌无力患者他克莫司治疗剂量精准选择的系统及其应用 |
Non-Patent Citations (5)
| Title |
|---|
| A Population Pharmacokinetic Model to Predict the Individual Starting Dose of Tacrolimus Following Pediatric Renal Transplantation;Louise M. Andrews等;《CLINICAL PHARMACOKINETICS》;20180430;第57卷(第4期);全文 * |
| an Integrative approach for Identifying a Metabolic Phenotype Predictive of Individualized Pharmacokinetics of Tacrolimus;PB Phapale等;《Clinical pharmacology & therapeutics》;20100430;第87卷;全文 * |
| 肾移植受者中他克莫司群体药动学分析及应用;彭凯;《中国优秀硕士学位论文全文数据库（医药卫生科技辑）》;20210215;全文 * |
| 肾移植患者西罗莫司的群体药代动力学模型研究;张钰等;《南京医科大学学报(自然科学版)》;20170915(第09期);全文 * |
| 药物代谢组学研究进展;葛纯等;《药学进展》;20170425(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114566292A (zh) | 2022-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106970228B (zh) | 用于蛋白质或多肽的混合物的从上到下多路复用质谱分析的方法 | |
| Gika et al. | Sample preparation prior to the LC–MS-based metabolomics/metabonomics of blood-derived samples | |
| CN112881547B (zh) | 一种肝硬化及肝炎人群早期肝癌诊断标志物的筛选方法 | |
| US20070176088A1 (en) | Feature selection in mass spectral data | |
| IL150840A (en) | Unintentional analysis method of complex samples | |
| Matsuda et al. | Assessment of metabolome annotation quality: a method for evaluating the false discovery rate of elemental composition searches | |
| van der Greef et al. | The role of metabolomics in systems biology: a new vision for drug discovery and development | |
| Van Riper et al. | Mass spectrometry-based proteomics: basic principles and emerging technologies and directions | |
| Luo et al. | Spatial metabolomics reveals skeletal myofiber subtypes | |
| CN114566292B (zh) | 应用内源性物质预测肾移植药物剂量的模型的构建方法 | |
| CN113406226B (zh) | 一种基于非靶向代谢组学检测gist患者血浆中伊马替尼代谢物的方法 | |
| CN112305118B (zh) | L-辛酰基肉碱作为疾病诊断的生物标志物 | |
| CN112669958A (zh) | 代谢物作为疾病诊断的生物标志物 | |
| Mahmud et al. | Comparison of global metabolite extraction strategies for soybeans using UHPLC-HRMS | |
| EP3911951A1 (en) | Urine metabolomics based method of detecting renal allograft injury | |
| CN115206420B (zh) | 一种精神分裂症异常基因-代谢调控网络的构建方法及应用 | |
| CN111896728A (zh) | 一种检测柴胡石油醚部位抗抑郁量-效与量-毒关系的组合标志物及其应用 | |
| CN116482380A (zh) | 抑郁症的生物标志物的用途 | |
| CN112630344B (zh) | 代谢标志物在脑梗死中的用途 | |
| CN114354827A (zh) | 代谢标志物及其在制备2型糖尿病的风险预测试剂盒中的应用和试剂盒 | |
| CN111383708B (zh) | 基于化学基因组学的小分子靶标预测算法及其应用 | |
| Duncan et al. | State-of-the-art capillary electrophoresis mass spectrometry methods for analyzing the polar metabolome | |
| CN112083110A (zh) | 一种何首乌导致药物肝损伤标志物的筛选和应用 | |
| CN112834652B (zh) | 急性主动脉夹层患者特异性生物标志组合物及其用途 | |
| AU2021106589A4 (en) | Method for Analyzing Differential Metabolites of Traditional Fermented Food Based on Non-targeted Metabolomics |
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 |