Adjoint Systems for Models of Cell Signaling Pathways and their Application to Parameter Fitting

The paper concerns the problem of fitting mathematical models of cell signaling pathways. Such models frequently take the form of sets of nonlinear ordinary differential equations. While the model is continuous in time, the performance index used in the fitting procedure involves measurements taken at discrete time moments. Adjoint sensitivity analysis is a tool which can be used for finding the gradient of a performance index in the space of parameters of the model. In the paper, a structural formulation of adjoint sensitivity analysis called the generalized backpropagation through time (GBPTT) is used. The method is especially suited for hybrid, continuous-discrete time systems. As an example, we use the mathematical model NF-kappaB of the regulatory module, which plays a major role in the innate immune response in animals.

[1]  A. Ma,et al.  Failure to regulate TNF-induced NF-kappaB and cell death responses in A20-deficient mice. , 2000, Science.

[2]  Piet Hemker,et al.  Nonlinear parameter estimation in initial value problems , 1974 .

[3]  Marek Kimmel,et al.  Mathematical model of NF- κB regulatory module , 2004 .

[4]  K. Schittkowski Parameter estimation in systems of nonlinear equations , 1994 .

[5]  Katherine C. Chen,et al.  Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. , 2003, Current opinion in cell biology.

[6]  Farren J. Isaacs,et al.  Computational studies of gene regulatory networks: in numero molecular biology , 2001, Nature Reviews Genetics.

[7]  Krzysztof Fujarewicz,et al.  On fitting of mathematical models of cell signaling pathways using adjoint systems. , 2005, Mathematical biosciences and engineering : MBE.

[8]  A. Griewank,et al.  Automatic differentiation of algorithms : theory, implementation, and application , 1994 .

[9]  Ursula Klingmüller,et al.  Modeling the Nonlinear Dynamics of Cellular Signal Transduction , 2004, Int. J. Bifurc. Chaos.

[10]  Lennart Ljung,et al.  On global identifiability for arbitrary model parametrizations , 1994, Autom..

[11]  V. Belevitch Inter-Reciprocity Applied to Electrical Networks , 1956 .

[12]  A. Arkin,et al.  Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. , 1998, Genetics.

[13]  Paul J. Werbos,et al.  Backpropagation Through Time: What It Does and How to Do It , 1990, Proc. IEEE.

[14]  Hidde de Jong,et al.  Modeling and Simulation of Genetic Regulatory Systems: A Literature Review , 2002, J. Comput. Biol..

[15]  Eric A. Wan,et al.  Diagrammatic Derivation of Gradient Algorithms for Neural Networks , 1996, Neural Computation.

[16]  A.V. Oppenheim,et al.  Analysis of linear digital networks , 1975, Proceedings of the IEEE.

[17]  Stephen A. Dyer,et al.  Digital signal processing , 2018, 8th International Multitopic Conference, 2004. Proceedings of INMIC 2004..

[18]  H. Bock Numerical Treatment of Inverse Problems in Chemical Reaction Kinetics , 1981 .

[19]  D. A. Baxter,et al.  Modeling transcriptional control in gene networks—methods, recent results, and future directions , 2000, Bulletin of mathematical biology.

[20]  A. Hoffmann,et al.  The I (cid:1) B –NF-(cid:1) B Signaling Module: Temporal Control and Selective Gene Activation , 2022 .

[21]  S. Liberty,et al.  Linear Systems , 2010, Scientific Parallel Computing.

[22]  Krzysztof Fujarewicz,et al.  Generalized Backpropagation through Time for Continuous Time Neural Networks and Discrete Time Measurements , 2004, ICAISC.

[23]  U. R. Prasad,et al.  Back propagation through adjoints for the identification of nonlinear dynamic systems using recurrent neural models , 1994, IEEE Trans. Neural Networks.