A systematic framework for the design of reduced-order models for signal transduction pathways from a control theoretic perspective

Systematic study of cellular signaling pathways facilitates improved understanding of processes including cell proliferation, metabolism and embryonic development. Key cell signaling pathway characteristics, such as transduction, amplification, feedback, and filtering display striking similarities to that of a control system. This leads us to believe that a control theoretic analysis of these pathways could enable a systems level understanding and help identify the role of individual modules in controlling the overall cellular behavior. Towards this end, this paper presents a framework with a step-by-step bottom-up methodology to guide the development of modular reduced-order signaling pathway components that collectively predict key observations and yet are simple. Critical steps of this iterative method include (1) modification of the pathway structure by addition and/or deletion of key nodes and/or arcs, (2) critical evaluation of multiple functional forms for fluxes and (3) estimation of the pathway model parameters. The parameter estimation minimizes the mismatch between the desired behavior and the predicted behavior using a hybrid procedure that involves a genetic algorithm to identify interesting regions in the parameter-space that are further explored using a local optimizer. The utility of this framework has been demonstrated by developing a reduced-order model for the mitogen-activated protein kinase (MAPK) pathway in mouse NIH-3T3 fibroblasts. The reduced-order model, consisting of five ordinary differential equations and 16 parameters, quantitatively predicts the bistable and proportional MAPK responses to the PDGF stimulus at different levels of MAP kinase phosphatase (MKP).

[1]  A. Nebreda,et al.  Regulation of the meiotic cell cycle in oocytes. , 2000, Current opinion in cell biology.

[2]  Venkat Venkatasubramanian,et al.  An Intelligent System for Reaction Kinetic Modeling and Catalyst Design , 2004 .

[3]  U. Bhalla,et al.  Emergent properties of networks of biological signaling pathways. , 1999, Science.

[4]  Prahlad T. Ram,et al.  MAP Kinase Phosphatase As a Locus of Flexibility in a Mitogen-Activated Protein Kinase Signaling Network , 2002, Science.

[5]  R. Fleischmann,et al.  Whole-genome random sequencing and assembly of Haemophilus influenzae Rd. , 1995, Science.

[6]  Chi-Ying F. Huang,et al.  Ultrasensitivity in the mitogen-activated protein kinase cascade. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[7]  D. Bar-Sagi,et al.  Ras effectors and their role in mitogenesis and oncogenesis , 1997, Journal of Molecular Medicine.

[8]  James E. Ferrell,et al.  c-Jun N-terminal Kinase Activation in Xenopus laevis Eggs and Embryos , 2001, The Journal of Biological Chemistry.

[9]  A. Thorburn,et al.  Activation of Ras by phorbol esters in cardiac myocytes , 1999, FEBS letters.

[10]  Ioannis P. Androulakis,et al.  Kinetic mechanism reduction based on an integer programming approach , 2000 .

[11]  P. Daoutidis,et al.  Nonlinear model reduction of chemical reaction systems , 2001 .

[12]  John R. Koza,et al.  Reverse Engineering of Metabolic Pathways from Observed Data Using Genetic Programming , 2000, Pacific Symposium on Biocomputing.

[13]  Ravi Iyengar,et al.  Robustness of the bistable behavior of a biological signaling feedback loop. , 2001, Chaos.

[14]  Venkat Venkatasubramanian,et al.  A hybrid genetic algorithm for efficient parameter estimation of large kinetic models , 2004, Comput. Chem. Eng..

[15]  David J. N. Limebeer,et al.  Linear Robust Control , 1994 .

[16]  George Stephanopoulos,et al.  Chemical Process Control: An Introduction to Theory and Practice , 1983 .

[17]  G. Stephanopoulos,et al.  Representation of process trends—Part I. A formal representation framework , 1990 .

[18]  K. Holmström,et al.  Model Reduction in Recursive Least Squares Identification , 1992 .

[19]  B. Séraphin,et al.  Positive feedback in eukaryotic gene networks: cell differentiation by graded to binary response conversion , 2001, The EMBO journal.

[20]  Paul I. Barton,et al.  Optimally-reduced kinetic models: reaction elimination in large-scale kinetic mechanisms , 2003 .

[21]  Raghunathan Rengaswamy,et al.  A Novel Interval-Halving Framework For Automated Identification of Process Trends , 2004 .

[22]  Yasuo Morishita,et al.  Mitogen-Activated Protein Kinase Cascade in Breast Cancer , 1999, Oncology.

[23]  J. Ferrell.  Self-perpetuating states in signal transduction: positive feedback, double-negative feedback and bistability. , 2002, Current opinion in cell biology.

[24]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[25]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[26]  J E Ferrell,et al.  Xenopus oocyte maturation: new lessons from a good egg. , 1999, BioEssays : news and reviews in molecular, cellular and developmental biology.

[27]  Venkat Venkatasubramanian,et al.  Automatic generation of qualitative descriptions of process trends for fault detection and diagnosis , 1991 .

[28]  D. Fell,et al.  Differential feedback regulation of the MAPK cascade underlies the quantitative differences in EGF and NGF signalling in PC12 cells , 2000, FEBS letters.

[29]  Vasilios Manousiouthakis,et al.  KINETIC MODEL REDUCTION USING GENETIC ALGORITHMS , 1998 .

[30]  H A KREBS,et al.  The tricarboxylic acid cycle. , 1948, Harvey lectures.

[31]  Venkat Venkatasubramanian,et al.  Catalyst design: knowledge extraction from high-throughput experimentation , 2003 .

[32]  B. Moore.  Principal component analysis in linear systems: Controllability, observability, and model reduction , 1981 .

[33]  E M Ross,et al.  GTPase-activating proteins for heterotrimeric G proteins: regulators of G protein signaling (RGS) and RGS-like proteins. , 2000, Annual review of biochemistry.

[34]  D. Lauffenburger,et al.  A Computational Study of Feedback Effects on Signal Dynamics in a Mitogen‐Activated Protein Kinase (MAPK) Pathway Model , 2001, Biotechnology progress.

[35]  J. Collins,et al.  Construction of a genetic toggle switch in Escherichia coli , 2000, Nature.

[36]  J. Ferrell.  Tripping the switch fantastic: how a protein kinase cascade can convert graded inputs into switch-like outputs. , 1996, Trends in biochemical sciences.

[37]  R. B. Newell,et al.  Control relevant model reduction: a reduced order model for ‘model IV’ fluid catalytic cracking units , 1994 .

[38]  M. Mavrovouniotis,et al.  Simplification of Mathematical Models of Chemical Reaction Systems. , 1998, Chemical reviews.

[39]  J. Doyle,et al.  Robust perfect adaptation in bacterial chemotaxis through integral feedback control. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[40]  L. Petzold,et al.  Model reduction for chemical kinetics: an optimization approach , 1999 .