cobrapy

COBRApy - Metabolic Modeling

Models the "metabolism" of a cell as a linear optimization problem. Used to predict bacterial growth under different conditions or design GMO strains.

When to Use

• Predicting microbial growth rates under different nutrient conditions.
• Designing metabolic engineering strategies (knockouts, additions).
• Understanding metabolic flux distributions.
• Comparing metabolic capabilities across organisms.
• Identifying essential genes and reactions.

Core Principles

Flux Balance Analysis (FBA)

Optimizes metabolic fluxes to maximize biomass production (or other objectives) subject to stoichiometric constraints.

Gene-Protein-Reaction (GPR)

Genes encode proteins (enzymes) that catalyze reactions. Knockouts affect reaction availability.

Constraints

Reaction bounds (lower/upper limits) represent enzyme capacity or nutrient availability.

Quick Reference

Standard Imports

import cobra
from cobra.io import load_model, save_model

Basic Patterns

# 1. Load model (e.g., E. coli)
model = cobra.io.load_model("iJO1366")
# Or: model = cobra.io.read_sbml_model("model.xml")

# 2. Run Flux Balance Analysis (FBA)
solution = model.optimize()
print(f"Growth rate: {solution.objective_value:.4f}")
print(f"Status: {solution.status}")

# 3. Knockout simulation (Gene essentiality)
with model:
    model.genes.get_by_id("b0002").knock_out()
    print(f"Growth after knockout: {model.optimize().objective_value:.4f}")

# 4. Change medium (nutrient availability)
model.medium = {
    'EX_glc__D_e': 10.0,  # Glucose uptake
    'EX_o2_e': 1000.0     # Oxygen
}
solution = model.optimize()

Critical Rules

✅ DO

• Check solution status - Ensure status is 'optimal' before using results.
• Use context managers - Wrap modifications in with model: to avoid permanent changes.
• Set appropriate bounds - Reaction bounds should reflect biological reality.
• Validate model - Use model.validate() to check for common issues.