Case study / scientific software

Subcellular Workflow

A modular, FAIR-compliant MATLAB framework for ODE biochemical pathway modeling, simulation, parameter estimation, global sensitivity analysis, and profile likelihood analysis. Thesis: A workflow for developing biochemical pathway models using ordinary differential equations.

Timeframe
Jan 2016-Mar 2025
Role
Architect and developer
Domain
Systems biology, neuroscience
Stack
MATLAB, Python, COPASI, Git

Problem

Make biochemical pathway modeling reproducible and adoptable.

ODE-based biological models can be hard to build, parameterize, analyze, and share when each lab relies on bespoke scripts. The goal was a workflow that supports initialization, simulation, parameter estimation, sensitivity analysis, and profile likelihood analysis while following FAIR principles.

Contribution

Designed the framework and implemented the analysis pipeline.

I architected and built the framework, implemented algorithms from the scientific literature, engineered interfaces between MATLAB solvers and COPASI, created documentation, and managed Git/GitHub collaboration across institutions.

Start
Model setup
Simulation
Parameter estimation
Sensitivity analysis
Profile likelihood
Results

Reusable ODE modeling pipeline

From setup to validated model analysis.

The workflow connects model initialization, simulation, parameter estimation, local and global sensitivity analysis, and profile likelihood analysis into a reproducible scientific software process.

Implementation details

Core capabilities

Model setup

Structured model initialization and data handling for biochemical pathway models, with attention to standard formats like SBML and SBtab.

Parameterization

Methods for parameter estimation and model calibration using MATLAB and scientific optimization tooling.

Model analysis

Local and global sensitivity analysis plus profile likelihood analysis to understand model behavior and parameter identifiability.

Validation

From simple motifs to basal ganglia signaling

LIRMO & LIREM

Ligand-receptor interactions and ligand-receptor-effector responses validated as small reference models for the workflow.

TM-M (CaMKII)

Translated complex CaMKII autophosphorylation dynamics into ODE form for medium-scale validation.

G-PIPK striatal model

Pushed the framework into larger parameter spaces: intensive parameter estimation and PLA on an Intel i9-10980XE workstation. Resulted in a peer-reviewed Neuroinformatics paper.

Workflow

FAIR-compliant scientific code with standard-tool interoperability.

The framework follows FAIR principles for findable, accessible, interoperable, and reusable research software, integrates standard formats, and uses Git workflows for cross-institutional collaboration.

Related work

See more

Browse all projects or review the CV for the broader research and publication context.

Case study · scientific software

Subcellular Workflow

A modular, FAIR-compliant MATLAB framework for ODE biochemical pathway modeling, simulation, parameter estimation, global sensitivity analysis, and profile likelihood analysis. Thesis: A workflow for developing biochemical pathway models using ordinary differential equations.

Timeframe

Jan 2016 – Mar 2025

Role

Architect and developer

Domain

Systems biology, neuroscience

Stack

MATLAB, Python, COPASI, Git

Start
Model setup
Simulation
Parameter estimation
Sensitivity analysis
Profile likelihood
Results

Setup

Model + data

Simulation

MATLAB solvers

Estimation

Optimization

Analysis

Sensitivity / PLA

Standards

SBML / SBtab

Output

Validated results

Problem

Make biochemical pathway modeling reproducible and adoptable.

ODE-based biological models can be hard to build, parameterize, analyze, and share when each lab relies on bespoke scripts. The goal was a workflow that supports initialization, simulation, parameter estimation, sensitivity analysis, and profile likelihood analysis while following FAIR principles.

Contribution

Designed the framework and implemented the analysis pipeline.

I architected and built the framework, implemented algorithms from the scientific literature, engineered interfaces between MATLAB solvers and COPASI, created documentation, and managed Git/GitHub collaboration across institutions.

Implementation details

Core capabilities

Model setup

Structured model initialization and data handling for biochemical pathway models, with attention to standard formats like SBML and SBtab.

Parameterization

Methods for parameter estimation and model calibration using MATLAB and scientific optimization tooling.

Model analysis

Local and global sensitivity analysis plus profile likelihood analysis to understand model behavior and parameter identifiability.

Validation

From simple motifs to basal ganglia signaling

LIRMO & LIREM

Ligand-receptor interactions and ligand-receptor-effector responses validated as small reference models for the workflow.

TM-M (CaMKII)

Translated complex CaMKII autophosphorylation dynamics into ODE form for medium-scale validation.

G-PIPK striatal model

Pushed the framework into larger parameter spaces: intensive parameter estimation and PLA on an Intel i9-10980XE workstation. Resulted in a peer-reviewed Neuroinformatics paper.

Workflow

FAIR-compliant scientific code with standard-tool interoperability.

The framework follows FAIR principles for findable, accessible, interoperable, and reusable research software, integrates standard formats, and uses Git workflows for cross-institutional collaboration.

Stack

Scientific software stack

MATLABSimBiologyOptimization ToolboxPythonCOPASIGit/GitHubSBMLSBtab

Related work

See more

Browse all projects or review the CV for the broader research and publication context.