Abstract:Spatio-temporal biochemical signaling in a large class of protein-protein interaction networks is well modeled by a reaction-diffusion system. The global existence of the solution to the reaction-diffusion system is determined by the reaction kinetics model and the protein network topology. We propose a novel reaction kinetics model that guarantees that the reaction-diffusion system with this model has a nonnegative invariant global classical solution for any network topology. We then present a computational method to identify the unknown parameters and initial values for a reaction-diffusion system with this reaction kinetics model. The identification approach solves an optimization problem that minimizes the cost function defined as the $L^2$-norm of the difference between the data and the solution of the reaction-diffusion system. We utilize an adjoint-based optimal control method to obtain the gradients of the cost function with respect to the parameters and initial values. The regularity of the global classical solutions of the reaction-diffusion system and its corresponding adjoint system avoids situations in which the gradients blow up, and therefore guarantees the success of the identification method for any network structure. Utilizing this gradient information, an efficient algorithm to solve the optimization problem is proposed and applied to estimate the mass diffusivities, rate constants and initial values of a reaction-diffusion system that models protein-protein interactions in a signaling network that regulates the actin cytoskeleton in a malignant breast cell.

A Computational Approach to Steady State Correspondence of Regular and Generalized Mass Action Systems

Network translation and steady state properties of chemical reaction systems

Translated Chemical Reaction Networks

A linear programming approach to dynamical equivalence, linear conjugacy, and the Deficiency One Theorem

A Decomposition-Based Approach for Deriving Positive Steady States of a Class of Chemical Reaction Networks with Non-Mass-Action Kinetics

Computational translation framework identifies biochemical reaction networks with special topologies and their long-term dynamics

Existence of Positive Steady States for Mass Conserving and Mass-Action Chemical Reaction Networks with a Single Terminal-Linkage Class

A Computational Approach to Multistationarity of Power-Law Kinetic Systems

Linear conjugacy of chemical kinetic systems

A generalization of Birch's theorem and vertex-balanced steady states for generalized mass-action systems

A Linear Programming Approach to Weak Reversibility and Linear Conjugacy of Chemical Reaction Networks

A computational approach to persistence, permanence, and endotacticity of biochemical reaction systems

Mass Action Dynamics of Coupled Reactions using Fluctuation Theory

Switching in mass action networks based on linear inequalities

Identification of Parameters and Initial Values for Reaction-Diffusion Systems in Protein Networks (Extended Version)

Computing weakly reversible linearly conjugate chemical reaction networks with minimal deficiency

Global stability of first order endotactic reaction systems

A network dynamics approach to chemical reaction networks

Stability Analysis of Biochemical Reaction Networks Linearly Conjugated to complex balanced Systems with Time Delays Added

Explicit exponential convergence to equilibrium for nonlinear reaction-diffusion systems with detailed balance condition