CN107480467A - 一种判别或比较药物作用模块的方法 - Google Patents
一种判别或比较药物作用模块的方法 Download PDFInfo
- Publication number
- CN107480467A CN107480467A CN201610399510.XA CN201610399510A CN107480467A CN 107480467 A CN107480467 A CN 107480467A CN 201610399510 A CN201610399510 A CN 201610399510A CN 107480467 A CN107480467 A CN 107480467A
- Authority
- CN
- China
- Prior art keywords
- module
- networks
- bio
- network
- topological
- 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
Links
Classifications
-
- 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
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
-
- 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
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
-
- 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
-
- 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
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
- G16B40/30—Unsupervised data analysis
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Biophysics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Biotechnology (AREA)
- Physiology (AREA)
- Molecular Biology (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioethics (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明提供一种判别或比较药物作用模块的方法,所述方法包括以下步骤:(1)对两个生物网络进行模块的划分,分别得到模块集合;(2)计算两个模块集合中各个模块的一个或多个拓扑结构特征,得到单一或综合的拓扑指标,并由此计算两个模块集合中对应模块的拓扑参数变化;(3)将对应模块的拓扑参数变化与显著性界值进行比较,当两个生物网络分别为药物干预前和干预后的生物网络时,如果拓扑参数变化≥显著性界值,则所述模块为药物的作用模块;当两个生物网络分别为第一药物干预后和第二药物干预后的生物网络时,如果拓扑参数变化≥显著性界值,则所述模块为第一药物和第二药物的作用差异模块。
Description
技术领域
本发明属于生物信息技术领域。具体而言,本发明涉及基于生物网络判别和比较药物靶点模块的方法。
背景技术
系统生物学和网络药理学的兴起,使药物机制研究及新药研发从传统的“单组分、单靶点、单疾病”研究模式向“多组分、多靶点、多途径”的方向发展。整合基因组、转录组、蛋白质组和代谢组等高通量数据,结合各种数学模型和算法,在网络的背景下研究疾病、药物、靶点数据及其相互作用关系也成为疾病和药物作用机制研究的热点。
许多复杂疾病的发生、发展都与一系列相互作用的基因或蛋白相关,疾病表型是不同的生物学过程在一个复杂网络中相互作用的表现。而药物则是通过作用于疾病网络中的多个靶点,对各靶点的作用产生协同效应,从而对疾病的发生、发展进行干预,最终达到治疗效果。因此,在疾病研究和药物研究领域,网络分析方法越来越多地被应用于识别复杂疾病的生物标志物及药物的作用靶点等。生物网络的分析也为更好的理解复杂疾病的病理机制、系统揭示药物的药理机制及新药的研发提供了新的思路。
模块性(modularity)是系统生物网络的一个重要特性。在系统生物网络的背景下,通过各种计算模型和算法对网络进行模块化分析,在模块的层次下对网络进行解构和分析成为网络药理学和系统生物学新的研究方向。模块指生物网络中紧密连接的功能集团,相关文献中亦表示为聚类(cluster)、亚网络(subnetwork)、社团结构(community),子集(subset)、模体(motif)等。模块在生物网络中具有紧密连接的拓扑结构特征,同时也具有稳定性和功能性的特点。研究者基于网络聚类、启发式搜索、种子扩展、网络拓扑结构、矩阵分解等原理发展了大量的模块划分方法,模块化研究成为生物网络研究的热点。
疾病病理和药物作用机制同样具有模块化的特征。“单组分、单靶点”的药物研究模式忽视了靶点基因或蛋白之间的相互作用关系,因此,一些研究尝试从模块的角度来筛选具有特异的靶点和较低副作用的药物。模块药理学(Modular pharmacology,MP)的概念框架也被提出,从模块化角度揭示药物与疾病间关系,提出通过药物的多靶点、模块化设计来治疗复杂疾病,这种以模块而不是单个基因或蛋白为靶点的方法为研究药物作用机制及新药研发提供了新的策略,而如何在复杂的生物网络中识别药物的靶点模块就成为一种挑战。因此,结合生物网络模块的拓扑结构、稳定性和功能性特点,探寻药物靶点模块识别方法具有重要的意义。
发明内容
本发明提供了一种基于生物网络判别和比较药物作用模块的方法,本方法整合了模块的拓扑结构分析、统计分析来综合判别药物干预后生物网络的作用模块(又称“应答模块”),组合生物功能分析则可以识别其是否为药物的靶点模块;此外,本发明的方法还可以通过比较不同药物干预后的生物网络的差异模块,来对比分析不同药物的应答差异。
定义
本文所用术语“生物网络”是指由节点与边组成的网络,由G(V,E)表示,V表示网络中的节点的集合,节点可以是基因、蛋白、化合物、代谢物等药物可能作用的靶点;E表示网络中的边的集合,边是节点之间的相互作用关系,例如可以是基因的共表达或相互作用关系、蛋白的相互作用关系、转录调节关系等靶点间可能的相互作用关系。常见的生物网络包括基因共表达网络、蛋白质相互作用网络、代谢网络、基因转录调控网络、非编码基因调控网络及各种跨组学的调控及相互作用网络等。以基因共表达网络为例,基因共表达网络是以基因之间表达谱数据的相关性为基础而构建的,常使用图模型来描述基因之间的关系,图中的节点代表基因,边表示两个基因之间的共表达相互作用关系。
本文所用术语“模块”是指生物网络中由节点和边组成的紧密连接的功能集团,其涵义等同于生物网络中的“子集”、“聚类”、“社团”、“子网”、“模体”等。
本文所用术语“药物的作用模块”是指与药物作用之前的生物网络中对应模块相比较,拓扑结构参数发生显著性变化的模块,其可与术语“药物的应答模块”互换使用。
本文所用术语“药物的靶点模块”是指将药物的作用模块经功能富集而确定确实具有相关功能的模块。
本文所用术语“药物的差异模块”是指将两种药物作用之后获得的生物网络中对应模块相比较,拓扑结构参数发生显著性变化的模块。
本文所用术语“药物的保守模块”是指将两种药物作用之后获得的生物网络中对应模块相比较,拓扑结构参数未发生显著性变化的模块。
本发明的具体技术方案如下:
本发明提供一种判别或比较药物作用模块的方法,所述方法包括以下步骤:
(1)对第一生物网络Gi(Vi,Ei)和第二生物网络Gj(Vj,Ej)进行模块的划分,分别得到模块集合Gi(Mi)和Gj(Mj);其中第一生物网络Gi(Vi,Ei)为药物干预前的生物网络,而第二生物网络Gj(Vj,Ej)为药物干预后的生物网络;或者第一生物网络Gi(Vi,Ei)为第一药物干预后的生物网络,而第二生物网络Gj(Vj,Ej)为第二药物干预后的生物网络;
(2)计算模块集合Gi(Mi)和Gj(Mj)中各个模块的一个或多个拓扑结构特征M(t),得到单一或综合的拓扑指标T=M(t1,t2...tm),m为拓扑结构特征数目,并由此计算两个模块集合中对应模块的拓扑参数变化△t/T;
(3)将对应模块的拓扑参数变化△t/T与显著性界值W进行比较,当第一生物网络Gi(Vi,Ei)为药物干预前的生物网络,而第二生物网络Gj(Vj,Ej)为药物干预后的生物网络时,如果△t/T≥显著性界值W,则所述模块为所述药物的作用模块;当第一生物网络Gi(Vi,Ei)为第一药物干预后的生物网络,而第二生物网络Gj(Vj,Ej)为第二药物干预后的生物网络时,如果△t/T≥显著性界值W,则所述模块为第一药物和第二药物的作用差异模块。
优选地,本发明的方法可以用于判别所述药物的作用模块是否是药物的靶标模块,即该方法还包括:
(4)将所述药物的作用模块进行生物功能富集分析,根据功能富集结果判断所述模块是否是所述药物的靶点模块。
在本发明方法的步骤(1)中,所述药物、第一药物、第二药物为同类型的用以预防和治疗疾病的中药或化学药。优选地,所述第一药物、第二药物为同类型的中药或化学药的药效组分。
所述第一生物网络和第二生物网络中的Vi或Vj为网络中节点的集合,Ei或Ej为网络中边的集合。优选地,所述第一生物网络和第二生物网络为生物分子网络,例如选自基因共表达网络、蛋白质相互作用网络、基因转录调控网络、非编码基因调控网络及代谢网络。更优选地,所述第一生物网络和第二生物网络为相同类型的基因共表达网络或蛋白质相互作用网络;更优选地,所述生物分子网络的节点数目≥200、≥500、≥1000或≥5000。
优选地,采用选自聚类、启发式搜索、种子扩展、矩阵分解、拓扑网络等网络分析算法进行模块划分;更优选地,采用选自分层聚类算法、WGCNA、MCODE、MCL、CFinder、CPM、SPC、G-N algorithm、ModuLand、DME、MINE、SVD的方法进行模块划分。
在本发明方法的步骤(2)中,所述拓扑结构特征M(t)选自模块性(modularity)、度(degree)、连接度(connectivity)、密度(density)、聚类系数(clustering coefficient)、介数中心性(betweenness)、网络直径(diameter)、最短路径(shortest path)、熵(entropy)和特征路径长度(characteristic path length)。
优选地,拓扑结构特征数目m为1-10,例如可以是1-10中的任何整数个。
优选地,拓扑结构特征数目m为2,拓扑结构特征M(t)包括密度和连接度;或者优选地,拓扑结构特征数目m为3,拓扑结构特征M(t)包括密度、连接度和熵。
优选地,计算拓扑参数变化△t/T包括:
计算第一生物网络中每个模块的拓扑参数与第二生物网络中相应模块的拓扑参数相比的变化量,即各个拓扑参数的△t/T值,采用离差标准化方法将各个拓扑参数的△t/T进行标准化处理,经过离差标准化后,各个拓扑参数的△t/T应在(0,1)之间。
在本发明方法的步骤(3)中,利用置换检验(permutation test)或自助重采样(bootstrap resampling)程序计算模块拓扑参数变化的显著性界值W;
优选地,显著性界值W通过以下方式获得:
利用置换检验或自助重采样方法得到大于1000个与第一生物网络中每个模块对应的随机模块,计算每个随机模块与第二网络相应模块的标准化△t/T值,进而计算所有随机模块的标准化△t/T值的正态分布概率,显著性界值W等于所述正态分布概率中p≤0.05时对应的标准化△t/T值。
在本发明方法的步骤(4)中,对药物的作用模块进行生物功能富集分析,以判断其是否是药物的靶点模块。模块的功能富集分析工具较多,可根据模块具体情况选择适合的工具。根据模块的功能富集结果可以判断药物的作用模块是否具有相应的生物学功能。例如,可以通过基因本体(GO)功能和/或通路(pathway)分析进行生物学功能富集分析。一般根据超几何分布计算功能富集的p值。GO的功能富集是计算模块组成基因在GO分类中某一(几)个特定的分支的超几何分布关系,返回模块的GO功能富集p值,小的p值则表示这个模块在该GO中出现了富集。同样,通路(pathway)富集分析是计算模块组成基因在同一条Pathway中的超几何分布关系,返回模块在该通路的富集p值,小的p值则表示这个模块在该通路中出现了富集。
本发明的方法的示意图见图1。
在生物网络中,对于判别与比较药物的作用模块一直缺少有效的方法,如何判断药物干预前后生物网络中模块的变化及其显著性一直存在挑战。本发明提出了一种利用网络的拓扑参数、统计和功能分析方法来判别与比较药物的作用靶点模块的方法,该方法能够有效地发现药物的作用模块,比较不同药物的作用模块的异同,为以模块为药物靶点的药物研发模式提供了新的策略。
具体地,药物干预后必然会对生物网络产生影响,而怎么判断哪些模块对药物产生应答,也就是怎样选择药物的靶点模块就是其中的关键问题。模块本身也是由节点和边组成的小网络,所以模块本身具有网络的拓扑结构特征,如模块会有特定的密度、连接度、最短路径、特征路径长度、模块性等。因此,一方面,本发明的方法利用了在特定药物干预前和干预后的网络中模块的拓扑结构变化,通过对比分析来判断该模块是否是药物的作用模块,进而通过功能富集分析,可以从中选择药物的靶点模块。另一方面,本发明的方法还可以利用在不同药物干预后的网络中模块的拓扑结构异同,通过对比分析来判断不同药物的不同作用模块,从而判断药物的作用异同。由此,本发明的方法是通过观察药物干预前后模块的拓扑结构变化(而不是单个基因、蛋白表达水平的变化)来反应药物的作用机制。
附图说明
以下,结合附图来详细说明本发明的实施方案,其中:
图1示出本发明所述药物靶点模块识别方法的流程图;
图2示出实施例1中BA组数据的基因共表达模式;
图3示出实施例1中BA组的基因共表达模块;
图4示出实施例1中BA组与模型组相比得到的应答模块(W>0.2);
图5示出实施例1中BA组与JA组相比得到的差异模块(W>0.2);
图6示出实施例1中BA组的16号模块示例及富集到的生物学功能;
图7示出实施例1中BA组的靶点模块示例;
图8示出实施例1中BA组的特征性模块示例;
图9示出实施例2中判别出的BA药物靶点模块中各节点在已知数据库中的相互作用关系。
图10示出实施例2中判别出的药物靶点模块内基因的表达情况。
图11示出实施例2中判别出的药物靶点模块内基因的相互作用关系情况。
图12示出实施例3中效应组的模块划分结果。
图13示出实施例3中效应组与对照组相比的抗PD-1治疗黑色素瘤的应答模块。
图14示出实施例4中效应组的模块划分结果。
图15示出实施例4中效应组与对照组相比的英夫利昔单抗治疗RA的应答模块。
具体实施方式
以下参照具体的实施例来说明本发明。本领域技术人员能够理解,这些实施例仅用于说明本发明,其不以任何方式限制本发明的范围。
下述实施例中的实验方法,如无特殊说明,均为常规方法。下述实施例中所用的药材原料、试剂材料等,如无特殊说明,均为市售购买产品。
实施例1清开灵有效组分干预小鼠脑缺血模型的靶点模块识别
本实例以精制清开灵有效组分黄芩苷(BA)、栀子苷(JA)、胆酸(UA)干预小鼠脑缺血模型的基因表达谱数据为例,构建不同药物干预前和干预后的基因共表达网络,并且以BA组的数据来具体说明本方法的实施。
整个实施过程的示意图见图1。
1.获得基因共表达数据
应用精制清开灵有效组分黄芩苷(BA)、栀子苷(JA)、胆酸(UA)分别对小鼠脑缺血模型(MCAO)进行药物干预。
小鼠脑缺血模型的构建:小鼠MCAO模型制备采用Zea-Longa线栓法,稍加改进,致备小鼠脑缺血1.5h再灌注24h的模型。小鼠术前禁食12h,自由饮水,经以10%水合氯醛溶液为麻醉剂,采用腹腔注射的方式,以0.4ml/100g体重的剂量,麻醉小鼠,仰卧固定,颈部正中切口,分离并暴露左侧颈总动脉,颈内外动脉,用镊子小心剥离出颈总动脉、颈外动脉,分离迷走神经。用动脉夹夹住颈总动脉,暂时阻断颈总动脉与颈内动脉的血流。用消毒后的缝合线在颈外动脉远心端结扎,并剪断颈外动脉,确保颈外动脉残端长度长于0.5cm。在颈内动脉近端备线,远端放置动脉夹,颈总动脉分叉膨大处切口,松开动脉夹,向颈内动脉插入栓线18~20㎜,在大脑中动脉起始端堵塞大脑中动脉,然后将颈总动脉连同鱼线一起结扎,外留10㎜线头,缝合皮肤。脑缺血1.5h时退出线栓,使血流再通。24h后观察小鼠的神经功能评分判断造模是否成功。
将构建得到的脑缺血模型小鼠随机分成5组,包括sham组、vehicle组、BA组、JA组和UA组,缺血后2h,用药组按2ml药液/kg体重分别经尾静脉注射BA(5mg/m1)、UA(7mg/m1)和JA(25mg/m1)。Sham组和vehicle组小鼠接经尾静脉注射0.9%NaCl(2ml/kg体重)。然后按一步法分别抽提小鼠海马RNA,利用Bioanalyzer基因芯片对RNA质量进行评估。利用小鼠Oligo片段的芯片对干预后小鼠的基因表达谱进行分析。所用芯片由与脑缺血相关的374个c-DNA组成。本实例所用基因表达谱数据共有4组,即模型组(Vehicle)、黄芩苷治疗组(BA)、栀子苷治疗组(JA)及胆酸治疗组(UA)。
将模型组基因表达谱数据作为药物干预前的数据,BA组、JA组和UA组基因表达谱数据作为药物干预后的数据。每组基因表达谱数据都由12个样本的374个基因组成(Tbp、Zeb1、Pou2f1、Foxb1、Creb1、Camk2g、Csf1、F5、Hspd1、Matn2、Mt1、Adamts1、Klf6、Dffa、Rgs18、Rhoa、Kcnmb1、Pdcd11、Pdpk1、Casp8ap2、Mogat1、Rps26、Ak1、Csnk2a2、Dkk2、Ppm1e、Tnfrsf22、Trp53i11、Smpd3、Grin1、Cdk5、Jund、E2f1、Apoe、Il1b、Prkar1b、Il7r、Ngfb、Rela、Ifnar1、Adcy6、Bak1、Fzd6、Prkch、Rgs4、Actg1、Gck、Rgs9、Sox9、Rgs1、Dgke、Rgs20、Map2k2、Pin1、Prkcn、Dgkz、Csnk1g1、Dusp4、Il11、Grb2、Shc1、Syk、Sim2、Ywhah、Fgf13、Bid、Gstm2、Rarg、Pou3f1、Camk2b、Mapkapk2、Tcf4、Sos1、Stat5a、Vegfb、Bad、Etv3、Id1、Lcat、Nf1、Gsn、Bbc3、Clu、Capn9、Ercc5、Comt、Ctsl、Amph、Vegfc、Bax、Cyp51、Sox10、Nfyc、Gata2、Id3、Lef1、Pou6f1、6330503C03Rik、Ech1、Ccl4、Itm2a、Hspa1a、Cbx3、Klf10、Idh3g、Gpx2、Map2k5、Daxx、E2f3、Fgf12、Ikbkg、Btrc、Ikbkap、Ifnar2、Cdk5、Psmb1、Sufu、Gab1、Sox30、Pxn、Pygo2、Ctnnb1、Grin2a、Il5ra、Cdk4、Bcl2l1、Actb、Myb、Prkca、Csf2rb2、Gnaq、B-raf、Wnt6、Adcy7、Cacna1b、Fzd7、Prkcm、Rock1、Adcy8、Prkcc、Sub1、Tuba1b、Rgs6、Plcb1、Mknk1、Diablo、Mef2c、Lrp1b、Dgkg、Rgs12、Serpina5、Hspb1、Ppm1b、Dlk1、Cdc42、Fadd、Mdfi、Fgf11、Map3k4、Klk1b3、Il6ra、Tgfb2、Wnt11、Ccna1、Map2k6、Htr1f、Zmat3、Bnip3、Tsg101、Vim、Srf、D14Abb1e、Cdh11、Vdac2、Tfdp1、Gak、Ccna2、Vegfa、Vegfa、Hdac1、Srebf1、Stch、E2f1、Nfatc1、Gna12、Gna13、Cacnb3、Zic1、Pou4f3、Tcf12、Ldb1、Capns1、Fxyd2、Gcgr、LOC100304588、Syt11、Gadd45a、Pbx2、Ier3、Mapk9、Ctnnbip1、Fgf15、Smad3、Nlk、Mecp2、Sigirr、Rgs18、Ptk2b、Sap30bp、Pcmt1、Tcf3、Braf、Ankrd6、Rgs5、Rap1gap、Adcy1、Grin2b、Gap43、Map2k1、Mapk10、Tgfb1、Lta、Rps6ka1、Wnt3、Rara、Prkcd、Atf4、Adcyap1r1、Cycs、Hint1、Rdx、Src、Adcy9、Prkce、Shcbp1、Elk3、Rgs14、Rgs17、Dusp10、Tubb3、Cyc1、Dusp16、Plcg2、Fzd10、Dgkd、Stat3、Mapk14、Map2k4、Htr1a、Map3k2、Frat1、Casp7、Eef2k、Thbd、Rarb、Camk4、Htr2c、E2f5、Met、Htr7、Camk2b、Stat6、Sod1、Efna4、Vdac3、Adora1、Bmp1、Vdac1、Grb2、Igfbp2、Top2b、Rpl35、Bdnf、Ppp3cb、Raf1、Cpe、Cacnb3、0610007C21Rik、Gna14、Gna11、Tuba1a、Zic3、Mlx、Id4、Ldb2、Sepp1、Prodh、S100a9、Pgam2、Rcan1、Abcc5、Ccr5、Ap1m1、Map3k5、Csnk1e、Axin1、Freq、Sh2b1、Rps6ka4、Wif1、Nkd1、Pam、Crem、Tgm2、Barhl1、Tradd、Plcd4、Ppp2r4、Otud7b、Rgs7、Casp2、Junb、Il2rg、Bad、Il1a、Egr1、Pdgfa、Gapdh、Eif4e、Apc、Prkcz、Parp1、Egfr、Prkcb1、Rgs2、Traf2、Ccr3、Rgs16、Smpd1、Tbp、Dgka、Mos、B230120H23Rik、Eif4e2、Rgs19、Adcy3、Creb5、Taf7、Pik3ca、Stat1、Il15、Atf3、Dvl3、Map3k3、Casp4、Kcnq1、Ptp4a3、fosB、Wnt3a、Calm1、Htr3a、Crkl、Casp3、Lhx1、Camk4、Selenbp2、Tcfe2a、Scg5、Pold3、Mmp2、Farp2、Pold2、Pold1、Gpx4、App、Mlh3、Rbl2、Tpp2、Cdh3、Fmo2、Pold4、Arf1、Sox1、Arhgef1),数据构成了374×12的基因表达矩阵,部分数据元素的格式如下:
Vehicle组:
BA组:
Array ID | Sam 1 | Sam 2 | Sam 3 | Sam 4 | Sam 5 | Sam 6 |
Sc02R02C10 | 1.019 | 0.966 | 1.019 | 1.121 | 1.146 | 0.967 |
Sc03R02C10 | 0.79 | 0.948 | 1.106 | 0.758 | 1.016 | 0.995 |
Sc04R02C10 | 0.931 | 1.003 | 0.873 | 0.712 | 1.02 | 0.833 |
Sc01R03C01 | 1.089 | 1.137 | 1.368 | 1.326 | 1.282 | 0.991 |
Sc02R03C01 | 1.269 | 1.482 | 1.451 | 1.522 | 1.044 | 1.321 |
Sc03R03C01 | 1.21 | 0.625 | 1.082 | 0.903 | 0.83 | 0.876 |
Sc04R03C01 | 1.822 | 2.376 | 0.804 | 1.226 | 2.714 | 1.944 |
Sc01R03C02 | 1.782 | 2.199 | 1.929 | 2.252 | 4.277 | 4.575 |
每一个单元格内的数值为该基因(行)在该样本(列)中的表达值,根据网络构建的方法,在这个表达谱中,具有同表达模式的基因可以被聚类在一起,形成共表达基因,这些共表达的基因聚类成为“基因共表达模块”。
2.网络构建和模块划分
用加权共表达网络分析(WGCNA)工具分别构建3个药物组数据的基因共表达网络(节点数目均为374个)并划分模块。经过模块识别BA组一共得到23个模块,模块的大小在3-149个节点(基因)之间;JA组得到42个模块,模块大小在3-46个节点之间;UA组得到15个模块,模块大小在3-29个节点之间。
以BA组的模块为例,其整个网络所有基因的表达情况如图2所示,可以看到明显的共表达模块的存在。基于分层聚类算法得到BA组的模块,每一种颜色代表一个模块,结果见图3。
3.模块拓扑参数计算
计算药物干预前(模型组)、BA干预后(BA组)的不同状态下模块的拓扑结构特征变化。
选择密度、连接度作为模块的拓扑指标,分别计算BA干预组和模型组模块的密度、连接度。将BA干预组与模型组进行对照,观察两组模块在密度、连接度指标中的变化。
密度:density=mean(vectorize Matrix(A))
连接度:
4.计算拓扑结构变化的显著性。
用置换检验(Permutation test)对两组模块的拓扑结构变化进行显著性检验,对BA组23个模块随机置换1000次,与模型组相比,计算每个模块的标准化△t/T,进而得到标准化△t/T在1000次随机置换中的正态分布概率,当统计p值≤0.05时对应的标准化△t/T值最小值,即药物干预后模块拓扑参数显著变化的界值。当p≤0.05时,W=0.2为模块显著变化的界值。与模型组相比,BA干预组中△t/T>0.2的模块即可认为是药物的应答模块。
将BA组与模型组相比,共发现11个BA的应答模块,结果见图4。
还发现不同药物相对于模型组的应答模块存在差别,这说明不同药物干预后产生的生物网络中存在的模块具有拓扑结构差异。因此,本发明的方法还可用于不同药物组之间的相比,△t/T>0.2的模块为两组之间的差异模块。将BA组与JA组相比,可发现9个差异模块,结果见图5。
5.通过功能分析判别与比较药物的靶点模块
利用功能富集分析工具(DAVIDE等工具)对三个药物组的应答模块进行GO功能和KEGG通路的功能富集分析,以p<0.05为富集分析显著性阈值。以BA组为例,BA组的11个应答模块中,BA-16模块富集到最多功能,包括41个GO功能和9条KEGG通路,BA-6模块没有富集到任何功能。如BA-16模块富集到的生物学功能如下:
可以看到BA-16模块与细胞信号、胎盘发育、磷代谢过程等生物学功能相关,而且BA-16模块与结肠癌和MAPK信号通路相关,因此可以推断BA干预BA-16应答模块的机制与如上的功能有关。图6为BA-16模块及其富集到的功能。
对得到的BA、JA、UA三个药物组的应答模块都进行生物学功能的富集。发现,与模型组相比,BA组得到11个应答模块,其中3个模块未能富集到任何的生物学功能,其余的8个模块为BA的靶点模块;与模型组相比,JA组得到22个应答模块,其中3个模块未能富集到任何的生物学功能,其余的19个模块为JA的靶点模块;与模型组相比,UA组得到8个应答模块,其中1个模块未能富集到任何的生物学功能,其余的7个模块为UA的靶点模块。BA组的靶点模块(BA-5、BA-15)示例见图7。
BA、JA、UA三组间两两比较,如果一组的模块与其他两组相比较均不同,则该模块为该组的特征性模块。经比较发现,BA有5个特征性模块,特征性模块(BA-9、BA-21)示例见图8;JA有8个特征性模块,UA有2个特征性模块。
通过模块的功能富集分析可以得到模块相关的生物学功能,即可能是药物的作用机制。
实施例2对判别出的药物的靶点模块进行验证
1.基于数据库已知知识的验证
通过String(functional protein association networks)数据库可以查找已知的模块内基因(蛋白)的各种相互作用关系。如BA-16模块,将模块内的节点输入到数据库中,可以得到节点的相互作用关系,参见图9。
由图9可以看到,BA-16模块内的多个节点有多条边相连,每一条边代表一种相互作用关系的证据类型,证明已知这些节点确实存在着相互作用关系。
2.基于实验的验证
采用实验方法,通过基因的表达量来证实靶点模块的组成意义。
用Western Blot方法对各治疗组的靶点模块基因进行验证。BA-21模块中,VEGF和B230120H23Rik的表达水平与Vehicle组相比具有显著性差异,分别参见图10A和图10B。
利用免疫共沉淀(Co-Immunoprecipitation)证实模块内基因编码蛋白的相互作用关系。BA-16模块内Met和Egfr蛋白的免疫共沉淀结果表明二者有直接的相互作用关系(图11B),BA-9模块内的Fmo2和Jund的免疫共沉淀结果表明二者有相互作用关系(图11A)。
实施例3在抗PD-1治疗黑色素瘤基因表达网络中的应用
本实例将本发明的方法应用于抗PD-1治疗黑色素瘤的mRNA表达数据,用以判别抗PD-1治疗黑色素瘤的靶点模块。将抗PD-1治疗产生反应的病例分为一组(效应组),将未产生反应的病例分为一组(对照组),利用WGCNA方法分别构建两组的基因共表达网络。
本实例中所用数据下载自GEO(Gene Expression Omnibus)数据库(https://www.ncbi.nlm.nih.gov/geo/,数据ID为GSE78220),该数据集为抗PD-1治疗黑色素瘤的mRNA表达数据,样本来源为人,共有28例样本,其中15例样本为效应组,13例样本为对照组。数据中共包含mRNA数量为25268。具体实施过程如下:
1.mRNA表达数据的整理
将下载的数据整理成WGCNA分析所需的数据表格格式,表格的列为mRNA,行为治疗样本,相对应的单元格内数值为该mRNA(行)在该样本(列)中的表达值,基于mRNA的共表达情况,进行网络构建和模块划分。
2.网络构建和模块划分
同样用WGCNA工具构建效应组的mRNA共表达网络并划分模块。在网络构建过程中,有3080个mRNA由于过多的缺失值而被移除,剩余的22188个mRNA被用来构建网络(所构建网络的节点数为21976个)和模块划分。在效应组中共得到30个模块,模块的大小在38-2900个节点(mRNA)之间。经WGCNA识别得到效应组的模块划分结果见图12。
3.模块拓扑参数计算
计算效应组和对照组模块的拓扑结构特征变化。选择密度、连接度作为模块的拓扑指标,分别计算效应组和对照组模块的密度、连接度。将效应组与对照组进行比较,观察两组模块在密度、连接度指标中的变化。
4.计算拓扑结构变化的显著性。
用自助重采样程序(bootstrap resampling)对两组模块的拓扑结构变化进行显著性检验,对效应组的模块组成基因进行10000次重采样,与对照组相比,计算每个模块的标准化△t/T,进而得到标准化△t/T在10000次重采样过程中的正态分布概率,当统计p值≤0.05时对应的标准化△t/T值最小值,即模块拓扑参数显著变化界值。当p≤0.05时,W=0.26为模块显著变化的界值。与对照组相比,抗PD-1效应组中△t/T>0.26的模块即可认为是药物的应答模块。
将效应组与对照组相比,共发现8个抗PD-1治疗黑色素瘤的应答模块,结果见图13。
5.通过功能分析判别抗PD-1治疗黑色素瘤的靶点模块
利用功能富集分析工具(DAVIDE 6.7工具)对效应组的应答模块进行GO功能和KEGG通路的功能富集分析,以p<0.05为富集分析显著性阈值。在抗PD-1治疗效应组中的8个应答模块中都能够富集到相应的功能,抗PD-1治疗黑色素瘤的靶点模块。
实施例4在英夫利昔单抗(infliximab)治疗类风湿性关节炎(RA)的基因表达网络中的应用
本实例将本发明应用于英夫利昔单抗治疗RA的mRNA表达数据,用3个拓扑参数(密度、连接度、熵)来判别英夫利昔单抗治疗RA的靶点模块。将病例按照英夫利昔单抗治疗前、后分为对照组和治疗组,利用WGCNA方法构建治疗组和对照组的基因共表达网络。
本实例中所用数据下载自GEO(Gene Expression Omnibus)数据库(https://www.ncbi.nlm.nih.gov/geo/,数据ID为GSE8350),该数据集为英夫利昔单抗治疗RA的mRNA表达数据,样本来源为人,本实例选用其中36例样本,其中18例为对照组,另18例为治疗组。具体实施过程如下:
1.mRNA表达数据的整理
数据整理过程参见实施例3。
2.网络构建和模块划分
同样用WGCNA工具构建效应组的mRNA共表达网络并划分模块。共有776个mRNA被用来构建网络(所构建网络的节点数为776个)和模块划分,在治疗组中共得到46个模块,模块的大小在3-73个节点(mRNA)之间。经WGCNA识别得到效应组的模块划分结果见图14。
4.模块拓扑参数计算
计算效应组和对照组模块的拓扑结构特征变化。选择密度、连接度、熵作为模块的拓扑指标,分别计算效应组和对照组模块的密度、连接度、熵。将治疗组与对照组进行比较,观察两组模块在密度、连接度、熵指标中的变化。
4.计算拓扑结构变化的显著性。
用置换检验(Permutation test)对两组模块的拓扑结构变化进行显著性检验,对治疗组的模块进行1000次置换,与对照组相比,计算每个模块的标准化△t/T,进而得到标准化△t/T在1000次置换过程中的正态分布概率,当统计p值≤0.05时对应的标准化△t/T值最小值,即模块拓扑参数显著变化界值。当p≤0.05时,W=0.15为模块显著变化的界值。与对照组相比,英夫利昔单抗治疗组中△t/T>0.15的模块即可认为是药物的应答模块。
将治疗组与对照组相比,共发现14个英夫利昔单抗治疗RA的应答模块,结果见图15。
5.通过功能分析判别英夫利昔单抗治疗RA的靶点模块
利用功能富集分析工具(DAVIDE 6.7工具)对效应组的应答模块进行GO功能和KEGG通路的功能富集分析,以p<0.05为富集分析显著性阈值。在英夫利昔单抗治疗RA的14个应答模块中,有2个模块未富集到任何生物学功能,其余12个模块为英夫利昔单抗治疗RA的靶点模块。
实施例5在黄芩苷治疗小鼠脑缺血模型蛋白相互作用网络中的应用
本实例以精制清开灵有效组分黄芩苷(BA)干预小鼠脑缺血模型的蛋白相互作用网络为例进行发明的应用。与模型组对照,得到BA干预小鼠脑缺血模型的差异表达基因,基于String数据库进行差异表达基因的网络映射,从而构建BA干预小鼠脑缺血模型的蛋白相互作用网络,用本发明的方法判别BA干预小鼠脑缺血蛋白相互作用网络的靶点模块。
本实施例与实施例1的差别在于所用数据来源数量不同,且网络类型为蛋白相互作用网络,因此与实施例1中的模块划分结果和应答模块判别结果均存在一定差异。下面仅对与实施例1中差异部分进行描述。
1.蛋白相互作用网络数据。
BA干预小鼠脑缺血模型的差异表达基因获取参见实施例1。利用String数据库,通过差异表达基因的映射构建蛋白相互作用网络,所构建的网络共包含2229个蛋白(节点)。
2.模块划分
BA干预小鼠脑缺血模型的蛋白相互作用网络中,用MCODE方法进行模块划分,共得到49个模块。
3.模块拓扑参数计算
计算BA干预后模块的拓扑结构特征变化。
选择密度、连接度作为模块的拓扑指标,分别计算BA干预组和模型组模块的密度、连接度。将BA干预组组成模块与模型组中相应节点进行对照,观察两组模块在密度、连接度指标中的变化。
4.计算拓扑结构变化的显著性。
用置换检验(Permutation test)对两组模块的拓扑结构变化进行显著性检验,对BA组49个模块随机置换1000次,与模型组相比,计算每个模块的标准化△t/T,进而得到标准化△t/T在1000次随机置换中的正态分布概率,当统计p值≤0.05时对应的标准化△t/T值最小值,即药物干预后模块拓扑参数显著变化的界值。当p≤0.05时,W=0.11为模块显著变化的界值。与模型组相比,BA干预组中△t/T>0.11的模块即可认为是药物的应答模块。
将BA组与模型组相比,共发现13个应答模块。
5.通过功能分析有个6应答模块没有富集到相应功能,其余7个模块为BA在蛋白相互作用网络是的靶点模块。
通过以上实例分析,证明了本发明能够有效地在生物网络中识别药物的靶点模块,还可以通过不同药物组间的比较得到不同药物的保守模块和差异模块。
以上对本发明具体实施方式的描述并不限制本发明,本领域技术人员可以根据本发明作出各种改变或变形,只要不脱离本发明的精神,均应属于本发明所附权利要求的范围。
Claims (10)
1.一种判别或比较药物作用模块的方法,所述方法包括以下步骤:
(1)对第一生物网络Gi(Vi,Ei)和第二生物网络Gj(Vj,Ej)进行模块的划分,分别得到模块集合Gi(Mi)和Gj(Mj);其中第一生物网络Gi(Vi,Ei)为药物干预前的生物网络,而第二生物网络Gj(Vj,Ej)为药物干预后的生物网络;或者第一生物网络Gi(Vi,Ei)为第一药物干预后的生物网络,而第二生物网络Gj(Vj,Ej)为第二药物干预后的生物网络;
(2)计算模块集合Gi(Mi)和Gj(Mj)中各个模块的一个或多个拓扑结构特征M(t),得到单一或综合的拓扑指标T=M(t1,t2…tm),m为拓扑结构特征数目,并由此计算两个模块集合中对应模块的拓扑参数变化Δt/T;
(3)将对应模块的拓扑参数变化Δt/T与显著性界值W进行比较,当第一生物网络Gi(Vi,Ei)为药物干预前的生物网络,而第二生物网络Gj(Vj,Ej)为药物干预后的生物网络时,如果Δt/T≥显著性界值W,则所述模块为所述药物的作用模块;当第一生物网络Gi(Vi,Ei)为第一药物干预后的生物网络,而第二生物网络Gj(Vj,Ej)为第二药物干预后的生物网络时,如果Δt/T≥显著性界值W,则所述模块为第一药物和第二药物的作用差异模块。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
(4)将所述药物的作用模块进行生物功能富集分析,根据功能富集结果判断所述模块是否是所述药物的靶点模块。
3.根据权利要求1或2所述的方法,其特征在于,步骤(1)中,所述药物、第一药物、第二药物为同类型的用以预防和治疗疾病的中药或化学药;优选地,所述第一药物、第二药物为同类型的中药或化学药的药效组分。
4.根据权利要求1至3中任一项所述的方法,其特征在于,步骤(1)中,所述第一生物网络和第二生物网络为生物分子网络,例如选自基因共表达网络、蛋白质相互作用网络、基因转录调控网络、非编码基因调控网络及代谢网络;更优选地,所述第一生物网络和第二生物网络为相同类型的基因共表达网络或蛋白质相互作用网络;
更优选地,所述生物分子网络的节点数目≥200、≥500、≥1000或≥5000。
5.根据权利要求1至4中任一项所述的方法,其特征在于,步骤(1)中,采用选自聚类、启发式搜索、种子扩展、矩阵分解、拓扑网络等网络分析算法进行模块划分;更优选地,采用选自分层聚类算法、WGCNA、MCODE、MCL、CFinder、CPM、SPC、G-N algorithm、ModuLand、DME、MINE、SVD的方法进行模块划分。
6.根据权利要求1至5中任一项所述的方法,其特征在于,步骤(2)中,所述拓扑结构特征M(t)选自模块性(modularity)、度(degree)、连接度(connectivity)、密度(density)、聚类系数(clustering coefficient)、介数中心性(betweenness)、网络直径(diameter)、最短路径(shortest path)、熵(entropy)和特征路径长度(characteristic path length);
优选地,拓扑结构特征数目m为1-10;
更优选地,拓扑结构特征数目m为2,拓扑结构特征M(t)包括密度和连接度;或者优选地,拓扑结构特征数目m为3,拓扑结构特征M(t)包括密度、连接度和熵。
7.根据权利要求1至6中任一项所述的方法,其特征在于,步骤(2)中,计算拓扑参数变化Δt/T包括:
计算第一生物网络中每个模块的拓扑参数与第二生物网络中相应模块的拓扑参数相比的变化量,即各个拓扑参数的Δt/T值,采用离差标准化方法将各个拓扑参数的Δt/T进行标准化处理,经过离差标准化后,各个拓扑参数的Δt/T应在(0,1)之间。
8.根据权利要求1至7中任一项所述的方法,其特征在于,步骤(3)中,利用置换检验(permutation test)或自助重采样(bootstrap resampling)程序计算模块拓扑参数变化的显著性界值W。
9.根据权利要求1至8中任一项所述的方法,其特征在于,步骤(3)中,显著性界值W通过以下方式获得:
利用置换检验或自助重采样方法得到大于1000个与第一生物网络中每个模块对应的随机模块,计算每个随机模块与第二网络相应模块的标准化Δt/T值,进而计算所有随机模块的标准化Δt/T值的正态分布概率,显著性界值W等于所述正态分布概率中p≤0.05时对应的标准化Δt/T值。
10.根据权利要求1至9中任一项所述的方法,其特征在于,步骤(4)中,通过基因本体(GO)功能和/或通路(pathway)分析进行生物学功能富集分析。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610399510.XA CN107480467B (zh) | 2016-06-07 | 2016-06-07 | 一种判别或比较药物作用模块的方法 |
EP16904507.7A EP3467689A4 (en) | 2016-06-07 | 2016-12-14 | METHOD FOR DIFFERENTIATING OR COMPARING AN ACTIVE SUBSTANCE ACTIVITY MODULE |
PCT/CN2016/109895 WO2017211059A1 (zh) | 2016-06-07 | 2016-12-14 | 一种判别或比较药物作用模块的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610399510.XA CN107480467B (zh) | 2016-06-07 | 2016-06-07 | 一种判别或比较药物作用模块的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107480467A true CN107480467A (zh) | 2017-12-15 |
CN107480467B CN107480467B (zh) | 2020-11-03 |
Family
ID=60578315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610399510.XA Active CN107480467B (zh) | 2016-06-07 | 2016-06-07 | 一种判别或比较药物作用模块的方法 |
Country Status (3)
Country | Link |
---|---|
EP (1) | EP3467689A4 (zh) |
CN (1) | CN107480467B (zh) |
WO (1) | WO2017211059A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109192252A (zh) * | 2018-08-23 | 2019-01-11 | 南开大学 | 共表达周期昼夜节律的转录组学在药物作用机制发现中的用途 |
CN110211634A (zh) * | 2018-02-05 | 2019-09-06 | 深圳华大基因科技服务有限公司 | 一种多组学数据联合分析的方法 |
CN110349625A (zh) * | 2019-07-23 | 2019-10-18 | 中国科学院心理研究所 | 一种人类大脑基因表达时空常模的建立方法 |
CN110957002A (zh) * | 2019-12-17 | 2020-04-03 | 电子科技大学 | 一种基于协同矩阵分解的药物靶点相互作用关系预测方法 |
CN112053742A (zh) * | 2020-07-23 | 2020-12-08 | 中南大学湘雅医院 | 分子靶标蛋白的筛选方法、装置、计算机设备和存储介质 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113567639B (zh) * | 2021-07-13 | 2023-05-16 | 中国食品药品检定研究院 | 中药材质量综合评价方法 |
CN114067906B (zh) * | 2021-11-15 | 2024-04-05 | 扬州大学 | 一种融合多源生物信息的关键蛋白质识别方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049678A (zh) * | 2012-11-23 | 2013-04-17 | 中国科学院自动化研究所 | 基于蛋白质交互作用网络的异病同治分子机理分析方法 |
CN103065066A (zh) * | 2013-01-22 | 2013-04-24 | 四川大学 | 基于药物组合网络的药物联合作用预测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013034300A2 (en) * | 2011-09-09 | 2013-03-14 | Philip Morris Products S.A | Systems and methods for network-based biological activity assessment |
CN103514381B (zh) * | 2013-07-22 | 2016-05-18 | 湖南大学 | 整合拓扑属性和功能的蛋白质生物网络模体识别方法 |
WO2015054266A1 (en) * | 2013-10-08 | 2015-04-16 | The Regents Of The University Of California | Predictive optimization of network system response |
CN103778349B (zh) * | 2014-01-29 | 2017-02-15 | 思博奥科生物信息科技(北京)有限公司 | 一种基于功能模块的生物分子网络分析的方法 |
-
2016
- 2016-06-07 CN CN201610399510.XA patent/CN107480467B/zh active Active
- 2016-12-14 WO PCT/CN2016/109895 patent/WO2017211059A1/zh unknown
- 2016-12-14 EP EP16904507.7A patent/EP3467689A4/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049678A (zh) * | 2012-11-23 | 2013-04-17 | 中国科学院自动化研究所 | 基于蛋白质交互作用网络的异病同治分子机理分析方法 |
CN103065066A (zh) * | 2013-01-22 | 2013-04-24 | 四川大学 | 基于药物组合网络的药物联合作用预测方法 |
Non-Patent Citations (4)
Title |
---|
PETER LANGFELDER ET AL: "Is My Network Module Preserved and Reproducible", 《PLOS COMPUTATIONAL BIOLOGY》 * |
张莹莹 等: "药物--疾病复杂网络的模块化解构", 《中国药理学通报》 * |
陈寅萤 等: "网络分析方法在疾病和药物研究中的应用", 《中国中药杂志》 * |
陈寅萤: "多维度多组分干预脑缺血药理机制的比较研究", 《中国博士学位论文全文数据库医药卫生科技辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110211634A (zh) * | 2018-02-05 | 2019-09-06 | 深圳华大基因科技服务有限公司 | 一种多组学数据联合分析的方法 |
CN110211634B (zh) * | 2018-02-05 | 2022-04-05 | 深圳华大基因科技服务有限公司 | 一种多组学数据联合分析的方法 |
CN109192252A (zh) * | 2018-08-23 | 2019-01-11 | 南开大学 | 共表达周期昼夜节律的转录组学在药物作用机制发现中的用途 |
CN110349625A (zh) * | 2019-07-23 | 2019-10-18 | 中国科学院心理研究所 | 一种人类大脑基因表达时空常模的建立方法 |
CN110957002A (zh) * | 2019-12-17 | 2020-04-03 | 电子科技大学 | 一种基于协同矩阵分解的药物靶点相互作用关系预测方法 |
CN112053742A (zh) * | 2020-07-23 | 2020-12-08 | 中南大学湘雅医院 | 分子靶标蛋白的筛选方法、装置、计算机设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107480467B (zh) | 2020-11-03 |
EP3467689A4 (en) | 2020-02-26 |
WO2017211059A1 (zh) | 2017-12-14 |
EP3467689A1 (en) | 2019-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107480467A (zh) | 一种判别或比较药物作用模块的方法 | |
Zhang et al. | A fast linear neighborhood similarity-based network link inference method to predict microRNA-disease associations | |
EP3191605B1 (en) | A droplet-based method and apparatus for composite single-cell nucleic acid analysis | |
CN106503482B (zh) | 一种用于定量分析药物干预前后生物分子网络中模块变化的方法 | |
EP2387758B1 (en) | Evolutionary clustering algorithm | |
Thareja et al. | A review of data mining optimization techniques for bioinformatics applications | |
Mallik et al. | WeCoMXP: Weighted connectivity measure integrating Co-methylation, Co-expression and protein-protein interactions for gene-module detection | |
Marchetti et al. | A methodology based on MP theory for gene expression analysis | |
Griffin et al. | Exploiting microarrays to reveal differential gene expression in the nervous system | |
CN113921084A (zh) | 疾病相关非编码rna调控轴多维靶向预测方法及系统 | |
Manganaro et al. | Adding semantics to gene expression profiles: new tools for drug discovery | |
Madhamshettiwar et al. | RMaNI: regulatory module network inference framework | |
Chen et al. | Identification of Important Biological Pathways for Ischemic Stroke Prediction through a Mathematical Programming Optimisation Model-DIGS | |
CN107038350B (zh) | 一种药物的长非编码rna靶点预测方法和系统 | |
Parvez et al. | The Role of AI, Fuzzy Logic System in Computational Biology and Bioinformatics | |
Kourou et al. | Identification of differentially expressed genes through a meta-analysis approach for oral cancer classification | |
WO2012123374A2 (en) | Method for robust comparison of data | |
Wimalarasan | CO-expression Analysis of RNA-sequence Data from Parkinson's Disease Patients | |
Sikora et al. | Quality improvement of rule-based gene group descriptions using information about GO terms importance occurring in premises of determined rules | |
Tripathi et al. | Genetic algorithm based clustering for gene-gene interaction in episodic memory | |
Anik et al. | Prediction of a Gene Regulatory Network in Cancer Cells | |
Benabderrahmane | Biomedical knowledge extraction using fuzzy differential profiles and semantic ranking | |
Stanfield | Comprehensive Characterization of the Transcriptional Signaling of Human Parturition through Integrative Analysis of Myometrial Tissues and Cell Lines | |
Cavallaro | Neuronal apoptosis revealed by genomic analysis: Integrating gene expression profiles with functional information | |
Nathoo | Senior Honors Thesis in Computational Biology April 2021 |
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 |