Abstract:We introduce a framework for analyzing ordinary differential equation (ODE) models of biological networks using statistical model checking (SMC). A key aspect of our work is the modeling of single-cell variability by assigning a probability distribution to intervals of initial concentration values and kinetic rate constants. We propagate this distribution through the system dynamics to obtain a distribution over the set of trajectories of the ODEs. This in turn opens the door for performing statistical analysis of the ODE system's behavior. To illustrate this we first encode quantitative data and qualitative trends as bounded linear time temporal logic (BLTL) formulas. Based on this we construct a parameter estimation method using an SMC-driven evaluation procedure applied to the stochastic version of the behavior of the ODE system. We then describe how this SMC framework can be generalized to hybrid automata by exploiting the given distribution over the initial states and the--much more sophisticated--system dynamics to associate a Markov chain with the hybrid automaton. We then establish a strong relationship between the behaviors of the hybrid automaton and its associated Markov chain. Consequently, we sample trajectories from the hybrid automaton in a way that mimics the sampling of the trajectories of the Markov chain. This enables us to verify approximately that the Markov chain meets a BLTL specification with high probability. We have applied these methods to ODE based models of Toll-like receptor signaling and the crosstalk between autophagy and apoptosis, as well as to systems exhibiting hybrid dynamics including the circadian clock pathway and cardiac cell physiology. We present an overview of these applications and summarize the main empirical results. These case studies demonstrate that the our methods can be applied in a variety of practical settings.

Dynamic sensitivity analysis of oscillating biochemical systems using modified collocation method

Sensitivity-based Adaptive Mesh Refinement Collocation Method for Dynamic Optimization of Chemical and Biochemical Processes

Analytical Sensitivity Computation Using Collocation Method with Non-Uniform Mesh Discretisation for Numerical Solutions of Optimal Control Problems

Dynamic Model of a Non-Linear Non-Ideal Multivariable-Coupling Hybrid System: Simulated Moving Bed Chromatography

A Novel Analytical Second‐order Sensitivity Calculation Approach Using the Finite Element Method for Chemical Engineering Problems

A New Sensitivity-Based Adaptive Control Vector Parameterization Approach for Dynamic Optimization of Bioprocesses

Application of Molecular Simulation in Biosensor Development

Sensitivity Analysis of Stoichiometric Networks: An Extension of Metabolic Control Analysis to Non-equilibrium Trajectories

Efficient Finite-Difference Estimation of Second-Order Parametric Sensitivities for Stochastic Discrete Biochemical Systems

A Bilevel NLP Sensitivity‐based Decomposition for Dynamic Optimization with Moving Finite Elements

A Computationally Efficient Method for Parameter Sensitivity Analysis of Microbially Explicit Biogeochemical Models Accounting for Long‐Term Behavior

Modeling and Analyzing Biological Oscillations in Molecular Networks

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

Global Versus Local Orthogonal Collocation in Simultaneous Approach

Sensitivity Embedded Homotopy-Based Backtracking Method for Chemical Process Simulation

Statistical Model Checking based Analysis of Biological Networks

Adjoint Sensitivity Analysis on Multi-Scale Bioprocess Stochastic Reaction Network

Transitory Sensitivity in Automatic Chemical Kinetic Mechanism Analysis

A Simple Comparison of Biochemical Systems Theory and Metabolic Control Analysis

Efficient computation of adjoint sensitivities at steady-state in ODE models of biochemical reaction networks