COBRApy - Python 代谢建模与通量分析

技能概述

COBRApy 是用于约束基础重建和分析（COBRA）基因组规模代谢模型的 Python 库，支持加载 SBML 模型、执行 FBA/FVA 分析、基因敲除模拟和代谢工程优化，是系统生物学和代谢工程研究的核心计算工具。

适用场景

1. 代谢模型构建与分析 加载和操作基因组规模代谢模型（SBML/JSON/YAML 格式），从本地文件或 BiGG/BioModels 远程仓库获取模型，检查模型结构和代谢物、反应、基因的关联关系。

2. 通量平衡分析与表型预测 执行标准 FBA、parsimonious FBA 和 geometric FBA 优化，预测细菌生长速率；进行通量变异分析（FVA）确定反应通量范围；使用通量采样探索可行解空间。

3. 代谢工程与基因敲除研究 单/双基因或反应敲除筛选，评估必需基因和代谢瓶颈；计算最小培养基组分；生成生产包络线，优化代谢途径设计；执行 gap-filling 修复不可行模型。

核心功能

1. 模型管理与格式转换 从文件或远程仓库（BiGG、BioModels）加载 SBML/JSON/YAML 格式代谢模型；将模型导出为标准格式进行长期存储或交换；内置教科书模型（textbook/e_coli_core）和大肠杆菌 iJO1366/iML1515 等基因级模型。

2. 代谢模拟与通量分析 通量平衡分析（FBA）优化目标函数；通量变异分析（FVA）在指定最优性范围内探索通量边界；单/双基因或反应删除研究；几何 FBA 和 parsimonious FBA 寻找中心解；使用 OptGP 或 ACHR 方法进行通量采样。

3. 模型构建与优化工具 从零创建 Model、Reaction、Metabolite 对象构建自定义代谢网络；设置边界反应和交换反应；配置基因-反应规则（GPR）布尔逻辑；使用上下文管理器临时修改模型并自动回滚；通过 GLPK/CPLEX/Gurobi 求解器执行优化。

常见问题

COBRApy 需要什么 Python 版本和求解器？

COBRApy 0.31.1 要求 Python 3.9 或更高版本（已放弃 3.8）。默认使用 GLPK 求解器（通过 swiglpk 自动安装），可选用 CPLEX 或 Gurobi 处理大规模问题。安装命令：uv pip install "cobra==0.31.1"。

如何加载和运行 FBA 分析？

使用 from cobra.io import load_model 加载模型（内置 textbook 或 e_coli_core 教程模型；远程需网络），然后调用 solution = model.optimize() 执行 FBA，通过 solution.objective_value 获取最优目标值（如生长速率），使用 solution.fluxes 访问各反应通量。

可以从哪些来源获取代谢模型？

COBRApy 支持三种来源：本地 SBML/JSON/YAML 文件（read_sbml_model() 等）；内置教科书模型（textbook、iJO1366、salmonella）；远程模型仓库 BiGG Models 和 BioModels（load_model("iML1515") 首次获取需网络，后续使用缓存）。SBML 是推荐交换格式。