2.1 · 网络药理学基础与 Cytoscape 网络构建

网络药理⭐⭐⭐⭐重要储备 ⭐⭐⭐|约 20 分钟13 节7 图
本章目标

掌握网络药理学"中药成分 → 靶点 → 疾病 → 通路"的完整研究流程,能独立使用 TCMSP、SwissTargetPrediction、GeneCards、STRING、Cytoscape、Metascape 等工具,产出一张"中药复方-成分-靶点-通路网络图"。

本章路线
网络药理学研究总流程
图 2.1-1网络药理学研究总流程
2.1 · 网络药理学基础与 Cytoscape 网络构建

0. 小白导读:网络药理学到底是什么?

如果你是零基础,先别被"网络药理学"这个高大上的名字吓到。用大白话说:

  • 传统药理学研究"一种药 → 一个靶点 → 一种病",像"一把钥匙开一把锁"。
  • 网络药理学研究"一个方子里的很多成分 → 很多靶点 → 很多通路 → 很多病",像"一张密密麻麻的快递物流网"。
📦生活化比喻

一个贯穿本章的比喻——快递物流

  • 中药成分 = 一批批"货物"(丹参酮、黄芩苷、川芎嗪……)

  • 靶点 = 一个个"收货地址"(AKT1、TP53、IL6……)

  • 疾病 = 一个"问题小区"(心梗、抑郁、脑缺血……)

  • 通路 = 一条条"运输干线"(PI3K-Akt 通路、TNF 信号通路……)

  • 网络图 = 一张"物流总览图",一眼看清货物从哪发出、送到哪、走哪条路

网络药理学要回答的核心问题就是:这个中药方子,到底靠哪些成分、作用在哪些靶点上、通过哪些通路,治好了这个病?

为什么学它? 你的研究方向之一是"网络药理学",这是中药现代化研究最主流的"开场白"方法——几乎所有中药复方 SCI 论文的第一张图,都是网络药理学分析图。学会它,你就掌握了中药研究的"通用入场券"。

🎯一句话总结

一句话总结:网络药理学 = 用"网络"的视角,把"药-靶-病-通路"的关系画成一张图,再用这张图解释中药为什么有效。


1. 网络药理学研究总流程

网络药理学虽然工具多、步骤多,但主线只有一条,可以概括为"五步走":

text
中药成分筛选 → 成分靶点预测 → 疾病靶点获取 → 交集/网络分析 → 富集分析
   (ADME)         (SwissTarget)     (GeneCards)      (STRING/Cytoscape)   (GO/KEGG)
📦生活化比喻

快递比喻:这五步就像"发货 → 查收货地址 → 确认问题小区 → 画物流图 → 找运输干线"。

五步对应的问题与工具速览

步骤 要回答的问题 核心工具 产出
① 成分筛选 这个中药里哪些成分"够格"当药? TCMSP / ETCM / HERB 活性成分表
② 靶点预测 这些成分会作用在哪些蛋白上? SwissTargetPrediction / SEA / PharmMapper 成分-靶点表
③ 疾病靶点 这个病和哪些基因有关? GeneCards / OMIM / DisGeNET 疾病靶点表
④ 交集与网络 药和病"撞车"在哪些靶点上?谁最关键? Venny / STRING / Cytoscape 交集靶点 + 网络图
⑤ 富集分析 这些靶点主要干哪些活、走哪条路? Metascape / DAVID / 微生信 GO/KEGG 富集图
💡提示

记住这条主线,后面每一节都是在给这条主线"填细节"。


2. 中药活性成分筛选(ADME)

截图
TCMSP 中药成分搜索页
图 2.1-2TCMSP 中药成分搜索页截图

2.1 为什么要筛选?

一味中药里有几十上百种化学成分,但不是所有成分都能成为"药"。一个成分要口服起效,得先过"吸收、分布、代谢、排泄"这几关——这就是 ADME

  • Absorption(吸收):吃进去能不能被肠道吸收?
  • Distribution(分布):吸收后能不能到达病灶?
  • Metabolism(代谢):会不会被肝脏迅速代谢掉?
  • Excretion(排泄):会不会太快被排出去?
📦生活化比喻

快递比喻:ADME 就像"货物能不能顺利通关"。有些货物(成分)过不了海关(吸收差)、半路被扣(代谢快),就送不到收货地址(靶点)。我们只挑"能通关的货"来研究。

2.2 两个核心筛选指标:OB 与 DL

网络药理学论文中最常用的两个"通关标准"是:

指标 全称 含义 常用阈值
OB Oral Bioavailability(口服生物利用度) 药物口服后进入血液循环的比例 ≥ 30%
DL Drug-Likeness(类药性) 化合物与已知药物结构的相似程度 ≥ 0.18
  • OB ≥ 30%:说明这个成分口服后"能进到血里",有机会起效。
  • DL ≥ 0.18:说明这个成分"长得像药",有成为药物的潜质。
注意

️ 这两个阈值是经验值,来自 TCMSP 数据库的官方推荐,也是绝大多数网络药理学论文的通行标准。部分研究还会叠加 Lipinski 五规则(分子量 ≤ 500、脂水分配系数 AlogP ≤ 5、氢键供体 ≤ 5、氢键受体 ≤ 10 等)进一步筛选。

2.3 工具:TCMSP 及其替代数据库

TCMSP(中药系统药理学数据库)是网络药理学的"第一站",提供中药成分、ADME 参数、靶点等信息。

  • 新版地址:https://tcmsp-e.com/ (若打不开,可尝试旧版 https://old.tcmsp-e.com/tcmsp.php 或更换浏览器)
  • 基本用法:搜索中药名(如"丹参 Salvia miltiorrhiza")→ 查看成分列表 → 按 OB/DL 阈值筛选 → 导出成分及其靶点。
注意

TCMSP 偶有打不开的情况(服务器在国外、访问不稳定),这是正常现象。遇到这种情况,用下面这些替代数据库顶上:

数据库 地址 特点
ETCM(中医药整合药理学平台) http://www.tcmip.cn/ 含药材、方剂、成分、靶标、疾病五大库,ETCM 2.0 收录 2000+ 药材
HERB(本草组鉴) http://herb.ac.cn/ 高通量中药实验参考数据库,含转录组数据,适合找"有实验证据"的成分
SymMap(症状映射库) http://www.symmap.org/ 把中医症状与西医症状、成分、靶点关联起来
BATMAN-TCM 2.0 http://bionet.ncpsb.org.cn/batman-tcm 基于化学相似性预测成分-靶点关系
💡提示

小白建议:入门阶段主攻 TCMSP,把它用熟;ETCM 作为"备用粮仓",TCMSP 打不开时顶上。


3. 成分靶点预测

截图
SwissTargetPrediction 靶点预测
图 2.1-3SwissTargetPrediction 靶点预测截图

筛选出活性成分后,下一步要知道:这些成分会作用在哪些蛋白(靶点)上?

📦生活化比喻

快递比喻:成分是"货物",靶点预测就是"查这批货的收货地址"。一个成分往往有多个收货地址(多靶点),这正是中药"多靶点"特性的来源。

3.1 三大靶点预测工具

工具 地址 原理 特点
SwissTargetPrediction http://www.swisstargetprediction.ch/ 基于配体 2D/3D 结构相似性"反向找靶" 最常用;输出靶点概率排名,可限定物种为 Homo sapiens
SEA http://sea.bkslab.org/ 基于化学相似性集合方法 输入一组配体,预测其共同靶点
PharmMapper http://www.lilab-ecust.cn/pharmmapper/ 基于药效团映射 输入小分子结构,输出潜在靶点

SwissTargetPrediction 使用要点: 1. 输入成分的 SMILES 结构式(可从 TCMSP 或 PubChem 获取)。 2. 物种选择 Homo sapiens(人类),因为我们要研究的是人体疾病。 3. 提交后得到靶点列表,按概率(Probability)从高到低排序。 4. 一般取 Probability > 0 或排名靠前的靶点;更严格可设阈值(如 Probability > 0.5)。

💡提示

SMILES 是什么? 一种用一行字符表示分子结构的"分子身份证",例如阿司匹林是 CC(=O)Oc1ccccc1C(=O)O。TCMSP 和 PubChem 都能直接复制。

3.2 靶点标准化:UniProt 基因名

不同数据库返回的靶点名称格式不统一(有的是蛋白名、有的是基因名、有的是 ID),必须统一转换成"官方基因名",否则后面合并分析会乱套。

  • 工具:UniProt(https://www.uniprot.org/)的 Retrieve/ID mapping 功能
  • 操作:粘贴靶点列表 → 选择"从蛋白名/ID 转成基因名" → 下载结果
  • 产出:一张"成分 → 基因名"的对照表,去重后得到成分靶点基因列表
注意

这一步最容易出错:靶点名称不统一会导致后续交集分析"对不上号"。养成"拿到靶点先转 UniProt 基因名"的习惯。


4. 疾病靶点获取

知道了"药"的靶点,还要知道"病"的靶点——这个疾病和哪些基因有关?

📦生活化比喻

快递比喻:疾病靶点 = "问题小区里的住户名单"。我们要找的是"既在发货单上、又住在问题小区"的那批人(交集靶点)。

4.1 三大疾病靶点数据库

数据库 地址 特点 使用场景
GeneCards https://www.genecards.org/ 人类基因综合数据库,收录 35 万+ 基因 网络药理学最常用;输入疾病名(如 "myocardial infarction")
OMIM https://www.omim.org/ 人类孟德尔遗传在线数据库 查遗传病相关基因
DisGeNET https://www.disgenet.org/ 基因-疾病关联综合数据库 整合多来源的基因-疾病关联

GeneCards 使用要点: 1. 搜索疾病英文名(如 "myocardial infarction")。 2. 得到基因列表,按 Relevance score(相关度)排序。 3. 一般取 Relevance score > 中位数(或按论文惯例取前若干)的基因作为疾病靶点。 4. 导出基因名列表,同样用 UniProt 标准化。

💡提示

疾病名要用英文。中文疾病名先翻译成英文(可用 MeSH 词表查标准译名,见模块四 4.1)。


5. 药物-疾病交集分析

现在手里有两份名单: - 成分靶点(药能作用的蛋白) - 疾病靶点(病相关的基因)

两者"撞车"的部分,就是药物-疾病交集靶点——这些靶点既被中药成分作用,又和疾病相关,是网络药理学的"核心研究对象"。

📦生活化比喻

快递比喻:交集靶点 = "既在发货单上、又住在问题小区"的那批收货人。网络药理学后面所有分析,都围绕这批人展开。

5.1 工具:Venny / Jvenn

  • Venny:https://bioinfogp.cnb.csic.es/tools/venny/ (经典在线工具)
  • Jvenn:https://jvenn.toulouse.inra.fr/app/example.html (交互式)

用法:把两份基因列表分别粘贴进两个圈 → 点击提交 → 得到韦恩图 → 点击中间交集区域,导出交集基因列表。

text
成分靶点                    疾病靶点
    ┌──────────────┐          ┌──────────────┐
    │              │  交集靶点  │              │
    │    A         │◄────────►│      B       │
    │              │   (核心)   │              │
    └──────────────┘          └──────────────┘
💡提示

交集靶点通常有几十个。后面 PPI 网络、富集分析都用这批交集靶点,而不是全部靶点。


6. PPI 蛋白互作网络构建

截图
STRING 数据库 PPI 网络构建
图 2.1-4STRING 数据库 PPI 网络构建截图

交集靶点不是孤立的——蛋白之间会"互相认识、互相配合"。PPI(Protein-Protein Interaction,蛋白-蛋白互作)网络就是把这些"认识关系"画成一张网。

📦生活化比喻

快递比喻:交集靶点是"收货人",PPI 网络是"收货人之间的熟人关系网"。谁和谁关系铁(互作强),谁就是网络里的"关键人物"。

6.1 工具:STRING 数据库

  • 地址:https://string-db.org/
  • 当前版本:STRING 12.5(2025 年更新,收录 5900 万+ 蛋白质、275 亿+ 条互作)

STRING 使用要点: 1. 粘贴交集靶点基因列表(物种选 Homo sapiens)。 2. 点击 Search → 得到 PPI 网络图。 3. 设置置信度阈值(confidence score,一般取 0.4 或 0.9,0.9 更严格)。 4. 导出数据: - TSV 格式(互作关系表,用于 Cytoscape 导入) - 网络图截图(用于论文)

注意

导出 TSV 很重要:STRING 网页上的图可以直接截图,但要做拓扑分析(找核心靶点),必须把互作关系表(TSV)下载下来,导入 Cytoscape。


7. Cytoscape 网络可视化与核心靶点筛选

截图
Cytoscape 网络可视化界面
图 2.1-5Cytoscape 网络可视化界面截图

7.1 Cytoscape 是什么?

Cytoscape 是网络可视化的"行业标准"软件,把"谁和谁有关系"画成漂亮的网络图,还能做拓扑分析。

  • 官网:https://www.cytoscape.org/
  • 当前稳定版:Cytoscape 3.10.x(3.10.0 已正式发布)
📦生活化比喻

快递比喻:STRING 给你"关系清单",Cytoscape 把清单画成"物流总览图",还能标出谁是"关键枢纽"。

7.2 基本操作流程

  1. 导入数据:File → Import → Network from File,选择从 STRING 下载的 TSV 文件。
  2. 美化网络: - 节点大小 = 度值(Degree,连接数越多越大) - 节点颜色 = 按拓扑指标渐变 - 布局:Layout → 选择 yFiles Organic / Circular 等
  3. 拓扑分析(找核心靶点):安装 CytoNCA 插件(Apps → App Manager → 搜索 CytoNCA),计算以下指标:
指标 英文 含义 通俗解释
Degree 一个节点连接了多少条边 "认识多少人"
介数 Betweenness 有多少最短路径经过该节点 "多少条路必经此地"
紧密度 Closeness 到其他所有节点的平均距离 "离大家有多近"
  • 综合 Degree、Betweenness、Closeness 都高的节点 = 核心靶点(Hub 基因)。 4. 模块聚类:安装 MCODE 插件,对网络聚类,找出"抱团的功能模块"(往往对应蛋白复合物或功能通路)。
💡提示

论文常见写法:取 Degree 排名前 10-20 的靶点作为核心靶点,用柱状图展示,并说明"这些靶点可能是复方发挥药效的关键"。

7.3 网络图的类型

网络药理学论文通常画三层网络图

text
复方 → 中药 → 成分 → 靶点 → 通路
  • 成分-靶点网络:中药成分与交集靶点的连接关系
  • 靶点-通路网络:交集靶点与富集通路的连接关系
  • 复方-成分-靶点-通路综合网络:把上面所有层叠在一起(本章交付物)

8. GO 与 KEGG 富集分析

复方-成分-靶点-通路网络图
图 2.1-6复方-成分-靶点-通路网络图

拿到交集靶点后,最后一个问题:这些靶点主要在细胞里干什么活、走哪条信号通路?

📦生活化比喻

快递比喻:富集分析 = "统计这批收货人都是干什么职业的、主要走哪几条运输干线"。如果发现很多靶点都跟"炎症"有关,说明这个药可能通过抗炎起效。

8.1 两个核心概念

GO(Gene Ontology,基因本体):把基因功能分成三大类——

GO 类别 英文 回答的问题
BP Biological Process(生物过程) 参与什么生理过程?(如炎症反应、凋亡)
CC Cellular Component(细胞组分) 在细胞的哪里工作?(如细胞膜、线粒体)
MF Molecular Function(分子功能) 有什么分子功能?(如蛋白结合、酶活性)

KEGG(京都基因与基因组百科全书):核心是通路数据库(Pathway),回答"这些靶点主要富集在哪些信号通路"(如 PI3K-Akt 信号通路、TNF 信号通路)。

8.2 四大富集分析工具

工具 地址 特点 适合谁
Metascape https://metascape.org/ 一站式,整合 40+ 数据库,粘贴基因列表一键分析 新手首选
DAVID https://david.ncifcrf.gov/ 经典工具,历史最悠久、引用最广 需要详细参数时
微生信 http://bioinformatics.com.cn/ 中文在线平台,支持气泡图、圈图等可视化 中文界面友好
ChiPlot https://www.chiplot.online/ 免费在线绘图,支持富集气泡图、热图 出图漂亮

Metascape 使用要点: 1. 粘贴交集靶点基因列表 → 物种选 Homo sapiens → 提交。 2. 自动完成 GO BP/CC/MF 与 KEGG 富集,输出富集表(含 P 值、富集倍数)。 3. 下载结果表,用微生信或 ChiPlot 画气泡图(横轴富集倍数、纵轴通路名、气泡大小代表基因数、颜色代表 P 值)。

注意

富集分析要关注 P 值:一般取 P < 0.05(或校正后 q 值 < 0.05)的通路作为显著富集通路。论文里通常展示排名前 10-20 的通路


9. 完整流程案例:丹参治疗心肌梗死

AI 绘图
PPI 网络艺术图
图 2.1-7PPI 网络艺术图AI 绘图

把前面所有步骤串起来,用"丹参治疗心肌梗死"走一遍完整流程(这也是你交付物的模板):

步骤 操作 示例结果
① 成分筛选 TCMSP 搜"丹参",OB≥30%、DL≥0.18 得到 60+ 活性成分(如丹参酮IIA、隐丹参酮、丹酚酸B)
② 靶点预测 SwissTargetPrediction 预测成分靶点,UniProt 标准化 得到 200+ 成分靶点基因
③ 疾病靶点 GeneCards 搜 "myocardial infarction",取 Relevance 前 50% 得到 1000+ 疾病靶点
④ 交集分析 Venny 取交集 得到 80+ 交集靶点
⑤ PPI 网络 STRING 导入交集靶点,置信度 0.9,导出 TSV 得到 PPI 网络
⑥ 核心靶点 Cytoscape + CytoNCA 算 Degree/Betweenness/Closeness 得到 AKT1、TP53、IL6、TNF 等核心靶点
⑦ 富集分析 Metascape 做 GO/KEGG 富集于 PI3K-Akt、TNF、IL-17 等通路
⑧ 出图 Cytoscape 画复方-成分-靶点-通路网络图 一张综合网络图

论文里的"一句话解读"模板

🎯一句话总结

"网络药理学分析显示,丹参的 60 个活性成分可能通过作用于 AKT1、TP53、IL6 等核心靶点,调控 PI3K-Akt 与 TNF 信号通路,从而发挥抗炎、抗凋亡作用,治疗心肌梗死。"

这就是网络药理学的"标准答案"——学会这个模板,你就能读懂并复现绝大多数网络药理学论文。


10. 常见坑

表现 解决办法
TCMSP 打不开 页面加载失败 换浏览器 / 用旧版链接 / 改用 ETCM、HERB
靶点名称不统一 交集分析对不上号 所有靶点先过 UniProt 转成官方基因名
疾病靶点太多 上万个基因没法分析 按 Relevance score 取前 50%,或加阈值
STRING 图导不出数据 只有截图没有 TSV 在 STRING 下载 TSV 格式互作表,再导入 Cytoscape
Cytoscape 没装插件 找不到 CytoNCA/MCODE Apps → App Manager 搜索安装
富集结果全是"无关"通路 富集到泛泛的通路 检查基因列表是否标准化、是否用了交集靶点而非全部靶点
网络图太乱 一团乱麻看不清 用 yFiles Organic 布局、按 Degree 调节点大小、隐藏孤立节点

11. 练习与交付物

练习(必须亲手做)

  1. 成分筛选:用 TCMSP 检索一味你感兴趣的中药(如"黄芪"),按 OB≥30%、DL≥0.18 筛选,导出活性成分表。
  2. 靶点预测:取 3 个活性成分,用 SwissTargetPrediction 预测靶点,用 UniProt 标准化。
  3. 疾病靶点:用 GeneCards 检索一个疾病(如"脑缺血 cerebral ischemia"),导出疾病靶点。
  4. 交集分析:用 Venny 取成分靶点与疾病靶点的交集。
  5. PPI 网络:用 STRING 构建交集靶点的 PPI 网络,导出 TSV。
  6. 核心靶点:用 Cytoscape + CytoNCA 计算拓扑指标,找出 Top 10 核心靶点。
  7. 富集分析:用 Metascape 做 GO/KEGG 富集,用微生信画气泡图。

交付物(验收标准)

笔记

用 Cytoscape 完成一张"中药复方-成分-靶点-通路网络图",输出:

  1. 核心靶点表(Top 10,含 Degree/Betweenness/Closeness 数值)

  2. GO/KEGG 富集图(Top 10 通路气泡图)

  3. 简要文字解读(150-300 字,说明核心靶点与通路如何解释药效)


12. 速查表

数据库与工具速查

用途 工具 地址
中药成分与 ADME TCMSP https://tcmsp-e.com/
成分库替代 ETCM http://www.tcmip.cn/
成分库替代 HERB http://herb.ac.cn/
靶点预测 SwissTargetPrediction http://www.swisstargetprediction.ch/
靶点预测 SEA http://sea.bkslab.org/
靶点预测 PharmMapper http://www.lilab-ecust.cn/pharmmapper/
基因名标准化 UniProt https://www.uniprot.org/
疾病靶点 GeneCards https://www.genecards.org/
疾病靶点 OMIM https://www.omim.org/
疾病靶点 DisGeNET https://www.disgenet.org/
交集分析 Venny https://bioinfogp.cnb.csic.es/tools/venny/
PPI 网络 STRING https://string-db.org/
网络可视化 Cytoscape https://www.cytoscape.org/
富集分析 Metascape https://metascape.org/
富集分析 DAVID https://david.ncifcrf.gov/
富集分析+绘图 微生信 http://bioinformatics.com.cn/
科研绘图 ChiPlot https://www.chiplot.online/

关键参数速查

参数 常用阈值 说明
OB ≥ 30% 口服生物利用度
DL ≥ 0.18 类药性
STRING 置信度 0.4 / 0.9 0.9 更严格
富集 P 值 < 0.05 显著富集
核心靶点 Top 10-20 按 Degree 等指标

下一章:2.2 · 分子对接基础与 PyMOL 可视化

配套学习资源

共 4 条 · 点击跳转