用途:COBRApy 是用于代谢模型基于约束的重建与分析(COBRA)的 Python 库,服务于系统生物学研究,可用于基因组规模代谢模型的计算模拟、代谢工程分析与表型预测。
主要能力:模型管理支持从内置或文件加载与保存模型(SBML、JSON、YAML);提供通量平衡分析(FBA)、简约 FBA、几何 FBA;通量变异性分析(FVA,含 loopless 模式);单/双基因与反应敲除分析;最小培养基计算;通量采样(OptGP、ACHR);产率包络(表型相平面);缺口填补;以及从零构建模型。模型组件通过 DictList 访问反应、代谢物和基因,可通过 get_by_id 或 query 检索,并能查看反应方程、边界、GPR 逻辑等属性。
使用方式:典型流程包括加载模型、运行优化、查看通量与状态。示例:import 相关函数后使用 model = load_model("textbook") 加载内置模型,solution = model.optimize() 执行 FBA,solution.objective_value 查看目标值;FVA 使用 flux_variability_analysis(model),可设置 fraction_of_optimum 和 loopless;敲除可用 single_gene_deletion(model)、double_gene_deletion(model, processes=4),或用 with model 上下文临时敲除并自动还原;采样使用 from cobra.sampling import sample 后 sample(model, n=1000, method="optgp", processes=4);培养基通过 model.medium 字典修改并整体重新赋值,最小培养基使用 minimal_medium。
适用场景:代谢工程靶点筛选、必需基因预测、培养基与生长条件优化、通量分布不确定性分析及代谢表型相平面研究。
限制:需保证模型可行且解状态为 optimal;可能遇到不可行、无界或求解缓慢问题,可检查培养基、边界与求解器设置;双敲除与采样计算量大,建议并行处理;loopless FVA 用于热力学可行性要求高的场景。模型交换与长期存储建议使用 SBML 格式;仅需目标值时用 slim_optimize() 提升性能。具体收费、价格与平台信息暂未提供。