COBRApy 约束代谢分析
将以下提示词粘贴到你的 AI 对话框中:
请根据 https://skillhub.cn/install/skillhub.md,安装 @org-02qudk26/cobrapy。
技能介绍
要解决什么问题
在系统生物学和代谢工程里,你往往手里是一份基因组规模的代谢网络(反应、代谢物、基因–反应规则),真正需要的是:在给定培养基和边界条件下,细胞能否生长、各反应通量是多少、敲掉哪个基因会停长、怎样改培养基或网络才更利于产物。手工在电子表格里做 FBA 不可扩展,也不方便做 FVA、敲除筛选和通量采样。
技能如何工作
COBRApy 把 约束重建与分析(COBRA) 落成可复用的 Python 工作流:
- 模型管理:从文件或 BiGG 等仓库加载模型,优先用 SBML 交换;用
DictList统一访问反应、代谢物与基因。 - FBA 与变体:标准 FBA、
slim_optimize()(只要目标值时更快)、简约 FBA、几何 FBA;优化后看解状态是否为optimal。 - FVA 与不确定性:
fraction_of_optimum控制次优空间;需要热力学可行时考虑 loopless FVA。 - 扰动与表型:单/双基因或反应敲除、交换反应(
EX_前缀,经model.medium管摄取/分泌)、生长与最小培养基、生产包络与 通量采样。 - 模型维护:间隙填充、从头构建;临时改边界或培养基时用 上下文管理器,避免污染全局状态。
典型路径:加载模型 → model.slim_optimize() 验可行 → FBA/FVA 或敲除筛选 → 按需采样或培养基优化。
适用边界与注意点
- 不可行/无界 多与培养基、反应上下界或交换反应上界有关;慢时可换 GLPK、CPLEX、Gurobi 等求解器(
model.solver)。 - 大规模双敲除、采样等适合 并行;采样结果要做数值稳定性检查。
- 复杂多步管线可对照技能内
references/workflows.md与 官方文档;本介绍不替代完整 API 说明。
使用场景
- 从 BiGG 或 SBML 加载基因组规模模型后,用 FBA 判断某培养基组合下理论生长率是否大于零
- 对候选基因做单敲除或双敲除筛选,找出使目标产物通量上升且仍满足生长约束的突变组合
- 在固定生物量目标下做 FVA,并设置 fraction_of_optimum 查看次优通量范围内某途径是否仍可行
- 比较两种碳源摄取上限下的生产包络或通量采样结果,评估代谢工程改造的可行区间
适合人员
- 做底盘菌株改造、需要在计算模型里先验评估敲除与培养基方案的代谢工程师
- 维护实验室或合作方 SBML 模型、负责 FBA/FVA 与可行性检查的系统生物学研究生
- 用 COBRA 方法写论文或报告、要把「加载—优化—敲除/培养基」流程脚本化的生物信息学研究者
- 在课程或课题中教授通量分析、希望用标准 Python API 演示交换反应与 model.medium 设置的教师