StrainDesign

A COBRApy^[1]-based package for computational design of metabolic networks

The comprehensive StrainDesign package for MILP-based strain design computation with the COBRApy toolbox supports MCS, MCS with nested optimization, OptKnock ^[2], RobustKnock ^[3] and OptCouple ^[4], GPR-rule integration, gene and reaction knockouts and additions as well as regulatory interventions. The automatic lossless network and GPR compression allows strain design computations from genome-scale metabolic networks. Supported solvers are GLPK (available from COBRApy), CPLEX, Gurobi and SCIP ^[5].

Installation instructions ...

Download Jupyter notebook examples ...

Parts of the compression routine is done by efmtool’s compression function (https://csb.ethz.ch/tools/software/efmtool.html^[6]).

Installation:

The StrainDesign package is available on pip and Anaconda. To install the latest release, run:

pip install straindesign

or

conda install -c cnapy straindesign

Developer Installation:

Download the repository and run

pip install -e .

in the main folder. Through the installation with -e, updates from a ‘git pull’ are at once available in your Python envrionment without the need for a reinstallation.

JAVA_HOME path:

In some cases, installing the StrainDesign python package may fail with the error:

JVMNotFoundException: No JVM shared library file (libjli.dylib) found. Try setting up the JAVA_HOME environment variable.

In this case you need to set your JAVA_HOME environment variable

If you’re on OS X and get the error

OSError: [Errno 0] JVM DLL not found

check that your Java and the JPype library is set up correctly.

Examples:

Computation examples are provided in the different chapters of this documentation. The original Jupyer notebook files are located in the StrainDesign package at docs/source/examples.

How to cite:

Schneider P., Bekiaris P. S., von Kamp A., Klamt S. - StrainDesign: a comprehensive Python package for computational design of metabolic networks. Bioinformatics, btac632 (2022)

Contents:

• Solvers
  • 3rd party solver installation
    • CPLEX
    • Gurobi
    • SCIP
  • Solver selection
• Network Analysis
  • Flux optimization (FBA/pFBA)
    • Parsimonious FBA (pFBA)
  • Flux variability analysis (FVA)
  • Yield optimization
    • Mathematical background
• Plotting the flux space
  • Production Envelopes
  • Yield Spaces
  • Mixed Plots
  • 3D plots
    • Interactive and animated 3D plots
• Computational strain design: Growth-coupled production (GCP)
  • pGCP: potentially growth-coupled production
  • wGCP: weakly growth-coupled production
  • dGCP: directionally growth-coupled production
  • SUCP: substrate-uptake-coupled production
• Minimal Cut Sets (MCS)
  • Prerequisites
    • 1) Add and verify production pathway
  • Example 1: Strain designs with a minimum product (1,4-butanediol) yield (SUCP strain design)
  • Example 2: Enforce product (1,4-BDO) synthesis at all growth states (dGCP strain design)
  • Example 3: Suppress flux states that are optimal with respect to a pre-defined objective function (wGCP strain design)
  • Example 4: Protect flux states that are optimal with respect to a pre-defined objective function (pGCP strain design)
  • Example 5: All single gene knockouts that prohibit growth (synthetic lethals).
  • Example 6: Genome-scale strain designs with a minimum product (1,4-butanediol) yield (SUCP strain design)
  • Example 7: Suppress flux states in a toy network
  • Example 8: Suppress and protect flux states in a toy network
  • Theoretical background
• Multi-level strain optimization approaches
  • OptKnock
    • Example 9: OptKnock strain design
    • Example 10: OptKnock strain design with a tilted objective function
    • Example 11: Genome-scale OptKnock strain design
  • RobustKnock
    • Example 12: RobustKnock strain design
  • OptCouple
    • Example 13: OptCouple strain design
  • Combining nested optimization strain design with MCS
    • Example 14: Combining OptKnock with a tilted objective function and the MCS approach
• Standalone network compression
  • Standalone GPR-integraton
    • Gene perturbation studies
• CNApy interface
• StrainDesign API
  • straindesign
    • Submodules
      • straindesign.compute_strain_designs
      • straindesign.cplex_interface
      • straindesign.efmtool
      • straindesign.glpk_interface
      • straindesign.gurobi_interface
      • straindesign.indicatorConstraints
      • straindesign.lptools
      • straindesign.names
      • straindesign.networktools
      • straindesign.parse_constr
      • straindesign.pool
      • straindesign.scip_interface
      • straindesign.solver_interface
      • straindesign.strainDesignMILP
      • straindesign.strainDesignModule
      • straindesign.strainDesignProblem
      • straindesign.strainDesignSolutions
    • Package Contents
      • DisableLogger

References:

[1] Ebrahim, A., Lerman, J.A., Palsson, B.O. et al. - COBRApy: COnstraints-Based Reconstruction and Analysis for Python. BMC Syst Biol 7, 74 (2013)

[2] Burgard, A. P., Pharkya, P., & Maranas, C. D. - Optknock: a bilevel programming framework for identifying gene knockout strategies for microbial strain optimization. Biotechnology and bioengineering, 84(6), 647–657 (2003)

[3] Tepper N., Shlomi T. - Predicting metabolic engineering knockout strategies for chemical production: accounting for competing pathways, Bioinformatics. Volume 26, Issue 4, Pages 536–543 (2010)

[4] Jensen K., Broeken V., Lærke Hansen A.S., et al. - OptCouple: Joint simulation of gene knockouts, insertions and medium modifications for prediction of growth-coupled strain designs. Metabolic Engineering Communications, Volume 8 (2019)

[5] Bestuzheva K., Besançon M., Chen W.K. et al. - The SCIP Optimization Suite 8.0. Available at Optimization Online and as ZIB-Report 21-41, (2021)

[6] Marco Terzer, Jörg Stelling, Large-scale computation of elementary flux modes with bit pattern trees, Bioinformatics, Volume 24, Issue 19, (2008), Pages 2229–2235,

Indices and tables

• Index

• Module Index

• Search Page