A Probabilistic Approach to Explore Signal Execution Mechanisms With Limited Experimental Data

Abstract

Mathematical models of biochemical reaction networks are central to the study of dynamic cellular processes and hypothesis generation that informs experimentation and validation. Unfortunately, model parameters are often not available and sparse experimental data leads to challenges in model calibration and parameter estimation. This can in turn lead to unreliable mechanistic interpretations of experimental data and the generation of poorly conceived hypotheses for experimental validation. To address this challenge, we evaluate whether a Bayesian-inspired probability-based approach, that relies on expected values for quantities of interest calculated from available information regarding the reaction network topology and parameters can be used to qualitatively explore hypothetical biochemical network execution mechanisms in the context of limited available data. We test our approach on a model of extrinsic apoptosis execution to identify preferred signal execution modes across varying conditions. Apoptosis signal processing can take place either through a mitochondria independent (Type I) mode or a mitochondria dependent (Type II) mode. We first show that in silico knockouts, represented by model subnetworks, successfully identify the most likely execution mode for specific concentrations of key molecular regulators. We then show that changes in molecular regulator concentrations alter the overall reaction flux through the network by shifting the primary route of signal flow between the direct caspase and mitochondrial pathways. Our work thus demonstrates that probabilistic approaches can be used to explore the qualitative dynamic behavior of model biochemical systems even with missing or sparse data.

Unified Tumor Growth Mechanisms from Multimodel Inference and Dataset Integration

“…(3) Consensus across the field led to labeling SCLC phenotypic subtypes by the dominant transcription factor expressed in that subtype. (4) Subtype with transcriptional signature intermediate between NE and Non-NE, named SCLC-A2. (5) Phenotypic transitions occur in a hierarchical manner from SCLC-A to SCLC-N to SCLC-Y cells.…”

Distinct execution modes of a biochemical necroptosis model explain cell type-specific responses and variability to cell-death cues

“…We estimated parameter values using PyDREAM, 48 a Python implementation of the DiffeRential E volution A daptive M etropolis (DREAM) method. 80 We utilized pMLKL Western blot data at the two highest TNF doses (100 and 10 ng/ml) and defined a multi-objective cost function, where Θ is the parameter set, x m ( t,d ) and x e ( t,d ) are model-predicted and experimentally measured pMLKL concentrations, respectively, at time t and TNF dose d , and σ ( t ) = 0.1· x e ( t,d ) (following previous studies 49,81,82 ). Parameter sampling was performed using five Monte Carlo chains, each run for 50,000 iterations, the first 25,000 of which were considered burn-in and discarded, resulting in 125,000 parameter sets.…”

“…where Q is the parameter set, xm(t,d) and xe(t,d) are model-predicted and experimentally measured pMLKL concentrations, respectively, at time t and TNF dose d, and s(t) = 0.1×xe(t,d) (following previous studies 49,81,82 ). Parameter sampling was performed using five Monte Carlo chains, each run for 50,000 iterations, the first 25,000 of which were considered burn-in and discarded, resulting in 125,000 parameter sets.…”

Programmatic modeling for biological systems

“…Lessons from e.g. hydrology and climate modeling have been highly influential toward addressing these issues [20*,42,43].…”