CN114566292A - 应用内源性物质预测肾移植药物剂量的模型的构建方法 - Google Patents

应用内源性物质预测肾移植药物剂量的模型的构建方法 Download PDF

Info

Publication number
CN114566292A
CN114566292A CN202210111071.3A CN202210111071A CN114566292A CN 114566292 A CN114566292 A CN 114566292A CN 202210111071 A CN202210111071 A CN 202210111071A CN 114566292 A CN114566292 A CN 114566292A
Authority
CN
China
Prior art keywords
model
metabolic
sample
correlation
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.)
Granted
Application number
CN202210111071.3A
Other languages
English (en)
Other versions
CN114566292B (zh
Inventor
楼燕
洪东升
何晓颖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
First Affiliated Hospital of Zhejiang University School of Medicine
Original Assignee
First Affiliated Hospital of Zhejiang University School of Medicine
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by First Affiliated Hospital of Zhejiang University School of Medicine filed Critical First Affiliated Hospital of Zhejiang University School of Medicine
Priority to CN202210111071.3A priority Critical patent/CN114566292B/zh
Publication of CN114566292A publication Critical patent/CN114566292A/zh
Application granted granted Critical
Publication of CN114566292B publication Critical patent/CN114566292B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H70/00ICT specially adapted for the handling or processing of medical references
    • G16H70/40ICT specially adapted for the handling or processing of medical references relating to drugs, e.g. their side effects or intended usage
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B15/00ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
    • G16B15/30Drug 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;
色谱柱选用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.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、根据初选临床特征建立多项逻辑回归方程。
拟合情况如下表所示:
多项逻辑回归模型拟合情况
Figure DEST_PATH_IMAGE002
在本实施例中根据上述多项逻辑回归模型,可以得到以下方程:
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 (10)

1.一种应用内源性物质预测肾移植药物剂量的模型的构建方法,其特征在于,包括以下步骤:
S1、收集待测样本,将所述待测样本中的蛋白质沉淀;所述待测样本为个体使用免疫抑制剂后的血液样本;
S2、将所述个体使用所述免疫抑制剂的效果分为高反应组、正常组和低反应组,所述效果根据T淋巴细胞含量判断;
S3、使用质谱仪分析所述待测样本的代谢特征,所述代谢特征包括所述待测样本中的若干种内源性代谢物质含量;
S4、通过偏最小二乘法分析所述代谢特征与所述效果分组的相关性,筛选出与所述效果分组具有强相关性的代谢特征;
S5、获取所述个体使用所述免疫抑制剂后的临床特征,所述临床特征包括所述个体的年龄、BMI、肾功能指标、肝功能指标、红细胞压积、白蛋白含量参数;
S6、通过分析所述临床特征与所述效果分组的相关性,筛选出与所述效果分组具有强相关性的临床特征;
S7、根据相关性分析和多项逻辑回归模型,以筛选后的所述代谢特征及所述临床特征与所述效果分组的相关性建立所述免疫抑制剂的药物剂量预测模型。
2.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法,其特征在于,所述使用质谱仪分析所述待测样本的代谢特征,具体包括如下步骤:
S31、以随机方式分析所述待测样本,首先连续检测3-5个质控样本,然后每10个所述待测样本插入1个质控样本,且同时检测正离子和负离子模式;
S32、在所述同时检测正离子和负离子模式下,通过预设数据库获得待测样本的碎片离子峰作为所述代谢特征。
3.如权利要求2所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法,其特征在于,所述获得待测样本的碎片离子峰作为所述代谢特征之后还包括数据预处理过程S221:
S321、对所述代谢特征进行峰提取、峰对齐、缺失值筛选、缺失值填补,并基于质控样本进行局部线性回归校正。
4.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法,其特征在于,所述步骤S4具体为通过两次独立的偏最小二乘法分析所述待测样本的代谢特征数据与所述效果的相关性进行筛选,包括以下步骤:
S41、以所有检测鉴定到的碎片离子峰为基础变量,目标效应值为结局变量,建立第一个独立的偏最小二乘法模型;
S42、通过第一个独立的偏最小二乘法模型,获得所述碎片离子峰对第一个独立的偏最小二乘法模型的首轮贡献指数;
S43、筛选首轮贡献指数大于1.5的碎片离子峰作为候选变量,建立第二个独立的偏最小二乘法模型;
S44、通过第二个独立的偏最小二乘法模型,获得所述候选变量对第二个独立的偏最小二乘法模型的二轮贡献指数;
S45、筛选二轮贡献指数大于1的碎片离子峰作为目标变量,以目目标变量离子峰对应的内源性代谢物质作为筛选后的所述代谢特征。
5.如权利要求4所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法,其特征在于,所述步骤S45之后还包括步骤S46:
S46、通过相关性分析确定所述代谢特征与效果分组的相关性,以相关性的显著水平小于0.05为标准,缩小并优化代谢特征的范围。
6.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法,其特征在于,所述步骤S7具体包括:
S71、根据多项逻辑回归模型,从所述初选临床特征中进一步筛选得到筛选后的所述临床特征;通过多项逻辑回归,选用似然比检验实现变量筛选,将显著水平小于0.05的所述初选临床特征纳入最终回归方程;
S72、根据所述初选临床特征建立多项逻辑回归方程。
7.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法,其特征在于,所述方法还包括步骤S8、使用正交最小二乘法判别模型验证模型的预测能力;具体包括:
S81、使用筛选后的所述代谢特征及所述临床特征形成正交最小二乘法判别模型的协变量;
S82、根据所述协变量构建正交最小二乘法判别模型,验证所述模型的预测能力。
8.如权利要求1所述的一种基于代谢组学构建他克莫司剂量预测模型的方法,其特征在于,所述使用质谱仪分析所述待测样本的代谢特征的参数具体为:
运用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;
色谱柱选用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.5min,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.3 mL/min,柱温为35 ℃,进样量2 μL。
9.如权利要求1所述的一种应用内源性物质预测肾移植药物剂量的模型的构建方法,其特征在于,所述免疫抑制剂为他克莫司;所述若干种内源性代谢物质包括血清肌酐、胆红素原、L-异亮氨酸、5-甲氧基吲哚乙酸、二十碳五烯酸、N2-琥珀酰精氨酸、色氨酸精氨酸和丁酸。
10.一种药物剂量预测的方法,其特征在于,使用如权利要求1-9任一项所述的方法所构建的模型,包括如下步骤:
S1、将目标样本的参数代入所述模型,获得疗效预测结果;
S2、根据所述疗效预测结果迭代药物剂量,获得不同疗效预测对应的药物剂量,实现剂量预测。
CN202210111071.3A 2022-01-29 2022-01-29 应用内源性物质预测肾移植药物剂量的模型的构建方法 Active CN114566292B (zh)

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 true CN114566292A (zh) 2022-05-31
CN114566292B 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 (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010016110A1 (ja) * 2008-08-05 2010-02-11 富士通株式会社 重回帰分析による予測モデルの作成方法、作成システムおよび作成プログラム
US20140235503A1 (en) * 2013-02-21 2014-08-21 Kyungpook National University Industry-Academic Cooperation Foundation Prediction method of glomerular filtration rate from urine samples after kidney transplantation
CN107247095A (zh) * 2017-05-25 2017-10-13 武汉大学 一种基于代谢组学分析的用于抗氧化衰老药物筛选的大鼠模型的构建方法
CN108280321A (zh) * 2018-01-29 2018-07-13 沈阳药科大学 基于细胞代谢轮廓分析构建抗肿瘤候选化合物作用机制预测模型的方法
CN109668984A (zh) * 2019-01-30 2019-04-23 山西医科大学第医院 一种基于代谢组学的喉癌血清判别模型的构建方法
CN112786145A (zh) * 2021-03-04 2021-05-11 华中科技大学同济医学院附属协和医院 一种器官移植患者他克莫司用药剂量精准预测方法
CN113035369A (zh) * 2021-03-10 2021-06-25 浙江大学 一种肾移植抗感染药物剂量预测模型的构建方法
CN113671078A (zh) * 2021-08-18 2021-11-19 郑州大学第一附属医院 一种基于代谢组学的脓毒症预后模型建立方法
CN113889221A (zh) * 2021-12-08 2022-01-04 首都医科大学宣武医院 重症肌无力患者他克莫司治疗剂量精准选择的系统及其应用
CN113930495A (zh) * 2021-10-18 2022-01-14 上海市第一人民医院 一种肝移植术后他克莫司初始剂量的预测模型及其个体化应用

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010016110A1 (ja) * 2008-08-05 2010-02-11 富士通株式会社 重回帰分析による予測モデルの作成方法、作成システムおよび作成プログラム
US20140235503A1 (en) * 2013-02-21 2014-08-21 Kyungpook National University Industry-Academic Cooperation Foundation Prediction method of glomerular filtration rate from urine samples after kidney transplantation
CN107247095A (zh) * 2017-05-25 2017-10-13 武汉大学 一种基于代谢组学分析的用于抗氧化衰老药物筛选的大鼠模型的构建方法
CN108280321A (zh) * 2018-01-29 2018-07-13 沈阳药科大学 基于细胞代谢轮廓分析构建抗肿瘤候选化合物作用机制预测模型的方法
CN109668984A (zh) * 2019-01-30 2019-04-23 山西医科大学第医院 一种基于代谢组学的喉癌血清判别模型的构建方法
CN112786145A (zh) * 2021-03-04 2021-05-11 华中科技大学同济医学院附属协和医院 一种器官移植患者他克莫司用药剂量精准预测方法
CN113035369A (zh) * 2021-03-10 2021-06-25 浙江大学 一种肾移植抗感染药物剂量预测模型的构建方法
CN113671078A (zh) * 2021-08-18 2021-11-19 郑州大学第一附属医院 一种基于代谢组学的脓毒症预后模型建立方法
CN113930495A (zh) * 2021-10-18 2022-01-14 上海市第一人民医院 一种肝移植术后他克莫司初始剂量的预测模型及其个体化应用
CN113889221A (zh) * 2021-12-08 2022-01-04 首都医科大学宣武医院 重症肌无力患者他克莫司治疗剂量精准选择的系统及其应用

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
LOUISE M. ANDREWS等: "A Population Pharmacokinetic Model to Predict the Individual Starting Dose of Tacrolimus Following Pediatric Renal Transplantation", 《CLINICAL PHARMACOKINETICS》 *
PB PHAPALE等: "an Integrative approach for Identifying a Metabolic Phenotype Predictive of Individualized Pharmacokinetics of Tacrolimus", 《CLINICAL PHARMACOLOGY & THERAPEUTICS》 *
张钰等: "肾移植患者西罗莫司的群体药代动力学模型研究", 《南京医科大学学报(自然科学版)》 *
彭凯: "肾移植受者中他克莫司群体药动学分析及应用", 《中国优秀硕士学位论文全文数据库(医药卫生科技辑)》 *
荆宁宁等: "免疫抑制剂在肾移植患者中的群体药动学研究进展", 《中国药学杂志》 *
葛纯等: "药物代谢组学研究进展", 《药学进展》 *
辛华雯等: "相关基因多态性与肾移植术后他克莫司个体化用药的关系", 《医药导报》 *

Also Published As

Publication number Publication date
CN114566292B (zh) 2022-11-11

Similar Documents

Publication Publication Date Title
Zimmer et al. Advances in proteomics data analysis and display using an accurate mass and time tag approach
Gika et al. Sample preparation prior to the LC–MS-based metabolomics/metabonomics of blood-derived samples
CN106970228B (zh) 用于蛋白质或多肽的混合物的从上到下多路复用质谱分析的方法
Paša-Tolić et al. Proteomic analyses using an accurate mass and time tag strategy
US6940065B2 (en) Method for characterizing biomolecules utilizing a result driven strategy
Taylor Therapeutic drug monitoring of immunosuppressant drugs by high-performance liquid chromatography–mass spectrometry
CN112881547B (zh) 一种肝硬化及肝炎人群早期肝癌诊断标志物的筛选方法
Van Riper et al. Mass spectrometry-based proteomics: basic principles and emerging technologies and directions
van der Greef et al. The role of metabolomics in systems biology: a new vision for drug discovery and development
Matsuda et al. Assessment of metabolome annotation quality: a method for evaluating the false discovery rate of elemental composition searches
Luo et al. Spatial metabolomics reveals skeletal myofiber subtypes
US20090269796A1 (en) Methods of detecting and treating myocardial ischemia and myocardial infarction
Michelhaugh et al. Using artificial intelligence to better predict and develop biomarkers
Valentine et al. Developing IMS–IMS–MS for rapid characterization of abundant proteins in human plasma
Higgins et al. Principles of phosphoproteomics and applications in cancer research
CN114566292B (zh) 应用内源性物质预测肾移植药物剂量的模型的构建方法
CN112305118B (zh) L-辛酰基肉碱作为疾病诊断的生物标志物
CN112669958A (zh) 代谢物作为疾病诊断的生物标志物
Merchant Mass spectrometry in chronic kidney disease research
Lee et al. Precursor mass prediction by clustering ionization products in LC-MS-based metabolomics
CN112630344B (zh) 代谢标志物在脑梗死中的用途
CN110068624B (zh) 生物标志物用于预测血管再狭窄以及检测试剂盒
Duncan et al. State-of-the-art capillary electrophoresis mass spectrometry methods for analyzing the polar metabolome
Furlani et al. Liquid chromatography-mass spectrometry for clinical metabolomics: an overview
CN112083110A (zh) 一种何首乌导致药物肝损伤标志物的筛选和应用

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