Pycellerator

Pycellerator converts Cellerator arrow-notation reaction descriptions into Python differential-equation models for dynamic simulation and analysis of biochemical networks.

Key Features:

• Reading and parsing: Reads Cellerator arrow notation from plain text files and supports mass-action kinetics, Michaelis–Menten–Henri (MMH) kinetics, gene regulation networks (GRN), Monod–Wyman–Changeux (MWC) models, user-defined reactions, and enzymatic expansions (KMech).
• Conversion to differential equations: Translates parsed reaction specifications into ordinary differential equations (ODEs) in Python for dynamic modeling of biochemical systems.
• Python solver generation: Automatically generates stand-alone Python code implementing numerical solvers for the derived ODEs.
• Execution and visualization: Executes generated solvers to produce time-course simulations and plots trajectories of species concentrations.
• Customization and flexibility: Produces editable solver code that can be modified and integrated with the broader Python ecosystem for diagnostics and advanced analyses.
• Independence from prior software: Operates independently of the original Cellerator/Mathematica implementation and avoids reliance on specialized character sets.

Scientific Applications:

• Systems biology modeling: Enables modeling of biochemical reaction networks across multiple reaction formalisms.
• Metabolic pathway analysis: Applies to studies of metabolic pathways through dynamic simulation of reaction kinetics.
• Gene regulatory network modeling: Supports gene regulation network (GRN) models and analysis of gene regulatory mechanisms.
• Enzyme kinetics and allosteric models: Handles enzyme kinetics including Michaelis–Menten–Henri and enzymatic expansion (KMech), and allosteric models such as Monod–Wyman–Changeux (MWC).
• Time-course simulation and experimental design: Provides time-course simulations for hypothesis testing, sensitivity exploration, and comparison of predicted dynamics to experimental data.

Methodology:

Reads and parses Cellerator arrow notation from text files (mass-action, MMH, GRN, MWC, user-defined, KMech), translates them into ODEs in Python, generates stand-alone Python solver code, and can execute the solvers to produce time-course simulations and plots.