CRNT4SBML

CRNT4SBML detects bistability in biochemical reaction networks by applying Chemical Reaction Network Theory (CRNT) and computational bifurcation analysis.

Key Features:

• SBML parsing: Parses Systems Biology Markup Language (SBML) models to extract reactions and species.
• C-graph construction: Constructs C-graph representations of biochemical networks under the assumption of mass action kinetics.
• Mass action kinetics assumption: Assumes mass action kinetics for reaction rate laws used in analysis.
• CRNT-based analysis: Applies Chemical Reaction Network Theory (CRNT) to assess whether networks can exhibit bistability independent of specific kinetic constants and concentrations.
• Optimization-based saddle-node search: Uses an optimization technique to search for saddle-node bifurcation points in network parameter space.
• Continuation analysis and bifurcation diagram generation: Performs continuation with respect to a specified parameter to extend steady-state branches and generate bifurcation diagrams.
• S-shaped bifurcation detection: Identifies S-shaped bifurcation diagrams indicative of bistability for specified parameters.
• Python implementation: Implemented as a Python package for programmatic analysis.

Scientific Applications:

• Bistability detection: Determining whether biochemical reaction networks can exhibit multiple stable steady states (bistability).
• Signaling pathway analysis: Investigating switching behavior in signaling pathways related to cellular decision-making.
• Identification of switching points: Locating saddle-node bifurcations that serve as switching points between steady states.
• Parameter-independent assessment: Assessing potential for bistability without exhaustive experimental testing of kinetic constants and molecule concentrations.

Methodology:

Parses SBML models, constructs C-graph representations assuming mass action kinetics, performs an optimization-based search for saddle-node bifurcations, and conducts continuation analysis with respect to a specified parameter to produce bifurcation diagrams.