cobrapy

COBRApy 是一个面向基于约束的重建与分析（COBRA）的 Python 库，用于系统生物学中的代谢模型研究，可对基因组尺度代谢模型进行计算模拟、代谢工程分析与表型预测。

【主要能力】 模型管理：可从文件或内置模型加载模型（如 load_model("textbook")、load_model("ecoli")、load_model("salmonella")），支持 read_sbml_model、load_json_model、load_yaml_model 读取，以及 write_sbml_model、save_json_model、save_yaml_model 保存。 模型结构与组件：可访问 model.reactions、model.metabolites、model.genes，并通过 get_by_id 获取指定对象，查看反应的化学计量方程、通量边界、基因-反应规则（GPR），以及代谢物的化学式与区室。 通量平衡分析（FBA）：通过 model.optimize() 执行标准 FBA，获取目标值与通量分布；支持 slim_optimize() 快速求目标值、修改目标函数，以及 pfba（简约 FBA）与 geometric_fba。 通量变异性分析（FVA）：flux_variability_analysis 可计算各反应的通量范围，支持 fraction_of_optimum 设定、loopless 消除热力学不可行环、指定 reaction_list。 基因与反应敲除：提供 single_gene_deletion、single_reaction_deletion、double_gene_deletion、double_reaction_deletion（双敲除支持多进程），也可通过上下文管理器手动敲除并优化。 培养基管理：查看与修改 model.medium（需整体重新赋值），并用 cobra.medium.minimal_medium 计算最小培养基，支持 minimize_components（MILP，较慢）与 open_exchanges。 通量采样：cobra.sampling.sample 支持 optgp（默认，可并行）与 achr 方法，并可用 OptGPSampler 配合 validate 校验样本。 产量包络：production_envelope 计算表型相平面，支持 carbon_sources 计算碳产量，并可结合 matplotlib 可视化。 缺口填充：cobra.flux_analysis.gapfill 配合通用模型补充反应，使模型可行。 模型构建：可从零使用 Model、Reaction、Metabolite 构建模型，添加代谢物与化学计量、设置 GPR、添加边界反应（exchange、demand）与目标函数。

【使用方式】 通过 Python 导入 cobra 相关模块调用；模型组件使用 DictList，可像列表和字典一样按索引或 ID 访问，并可用 query 查询。反应边界定义可行通量范围（不可逆与可逆），建议用 .bounds 同时设置上下界。GPR 使用布尔逻辑（and/or）连接基因与反应。交换反应约定以 EX_ 为前缀，正通量表示分泌、负通量表示摄取。

【适用场景】 系统生物学研究、代谢工程分析、生长预测、必需基因筛选、培养基优化、通量不确定性分析与表型相平面计算等。

【限制与注意事项】 资料未提供该工具的收费、价格或平台信息。使用中可能遇到不可行解、优化较慢、无界解与导入错误，需检查培养基约束、反应边界、模型一致性、交换反应上界与文件格式。求解器可通过 model.solver 切换（如 GLPK、CPLEX、Gurobi）。最佳实践包括：临时修改使用上下文管理器、分析前用 slim_optimize() 验证模型可行性、检查解状态是否为 optimal、热力学可行性关注时使用 loopless FVA、合理设置 fraction_of_optimum、对采样与双敲除等耗时操作并行化、优先使用 SBML 格式交换与存储模型、校验通量样本以保证数值稳定性。具体复杂工作流可参考 references/workflows.md 与 references/api_quick_reference.md，官方文档为 https://cobrapy.readthedocs.io/en/latest/。原文还建议在复杂多步分析场景下可考虑使用 K-Dense Web，此建议属可选生产力升级，非必要步骤。