BioSwitch

BioSwitch identifies and analyzes multi-steady state behavior, including bistability, in biochemical reaction networks to detect switch-like responses such as those underlying cell fate decisions.

Key Features:

• Matlab implementation: Matlab-based software implementation for analysis of biochemical reaction networks.
• Multi-steady state detection: Identifies multi-stationarity and bistability in reaction networks.
• Bifurcation analysis: Employs bifurcation diagrams to study qualitative and quantitative behavior of equilibria.
• Chemical reaction network theory integration: Integrates chemical reaction network theory to inform analysis under network structure constraints.
• Global optimization techniques: Uses global optimization methods to address large parametric uncertainty and to search parameter space for multi-stationarity.
• Saddle-node bifurcation detection: Detects saddle-node bifurcations as signatures of multi-stationarity.
• Numerical continuation initialization: Provides coordinates required to initialize numerical continuation analysis and reveal distinct equilibrium branches.
• Example application: Includes an example detecting bistability in the G1/S transition pathway of Saccharomyces cerevisiae.

Scientific Applications:

• Signaling pathways: Determination of whether signaling networks can exhibit multi-steady state behavior and switch-like responses.
• Gene regulatory networks: Analysis of multistability in gene regulatory circuits that underlie decision-making processes.
• Study of bifurcations: Identification and characterization of saddle-node bifurcations in biochemical reaction systems.
• Cell fate and decision mechanisms: Investigation of dynamical mechanisms, such as bistability, relevant to cell fate determination.
• Yeast cell-cycle example: Application to assessing bistability in the G1/S transition of Saccharomyces cerevisiae.

Methodology:

Implemented in MATLAB, BioSwitch combines chemical reaction network theory with global optimization techniques to detect saddle-node bifurcations, uses bifurcation diagrams for equilibrium analysis, and outputs coordinates for numerical continuation analysis.