Equilibria and Stability of a Class of Positive Feedback Loops Mathematical Analysis and Its Application to Caspase-dependent Apoptosis

Positive feedback loops are common regulatory elements in metabolic and protein signalling pathways. The length of such feedback loops determines stability and sensitivity to network perturbations. Here we provide a mathematical analysis of arbitrary length positive feedback loops with protein production and degradation. These loops serve as an abstraction of typical regulation patterns in protein signalling pathways. We first perform a steady state analysis and, independently of the chain length, identify exactly two steady states that represent either biological activity or inactivity. We thereby provide two formulas for the steady state protein concentrations as a function of feedback length, strength of feedback, as well as protein production and degradation rates. Using a control theory approach, analysing the frequency response of the linearisation of the system and exploiting the Small Gain Theorem , we provide conditions for local stability for both steady states. Our results demonstrate that, under some parameter relationships, once a biological meaningful on steady state arises, it is stable, while the off steady state, where all proteins are inactive, becomes unstable. We apply our results to a three-tier feedback of caspase activation in apoptosis and demonstrate how an The authors would like to thank Ms. Niamh Connolly and the anonymous reviewers for thoroughly reviewing the manuscript and providing helpful ideas. intermediary protein in such a loop may be used as a signal amplifier within the cascade. Our results provide a rigorous mathematical analysis of positive feedback chains of arbitrary length, thereby relating pathway structure and stability.

[1]  U. Alon An introduction to systems biology : design principles of biological circuits , 2019 .

[2]  Míriam R. García,et al.  Positive feedback in the Akt/mTOR pathway and its implications for growth signal progression in skeletal muscle cells: an analytical study. , 2012, Journal of theoretical biology.

[3]  M. Hayden,et al.  Caspase-6 and neurodegeneration , 2011, Trends in Neurosciences.

[4]  Shinji Hara,et al.  Existence criteria of periodic oscillations in cyclic gene regulatory networks , 2011, Autom..

[5]  Jakub Wenus,et al.  Mathematical modelling of the mitochondrial apoptosis pathway. , 2011, Biochimica et biophysica acta.

[6]  E. Kuranaga Caspase signaling in animal development , 2011, Development, growth & differentiation.

[7]  Maike A. Laussmann,et al.  The Caspase-8 Dimerization/Dissociation Balance Is a Highly Potent Regulator of Caspase-8, -3, -6 Signaling* , 2010, The Journal of Biological Chemistry.

[8]  Markus Rehm,et al.  Diffusion is capable of translating anisotropic apoptosis initiation into a homogeneous execution of cell death , 2010, BMC Systems Biology.

[9]  D. Bernstein Matrix Mathematics: Theory, Facts, and Formulas , 2009 .

[10]  G Melino,et al.  Ordering of caspases in cells undergoing apoptosis by the intrinsic pathway , 2009, Cell Death and Differentiation.

[11]  R. E. Huber,et al.  Practical Considerations When Using Temperature to Obtain Rate Constants and Activation Thermodynamics of Enzymes with Two Catalytic Steps: Native and N460T-β-Galactosidase (E. coli) as Examples , 2009, The protein journal.

[12]  Marc Tessier-Lavigne,et al.  APP binds DR6 to trigger axon pruning and neuron death via distinct caspases , 2009, Nature.

[13]  D. Lauffenburger,et al.  Modeling a Snap-Action, Variable-Delay Switch Controlling Extrinsic Cell Death , 2008, PLoS biology.

[14]  Onn Brandman,et al.  Feedback Loops Shape Cellular Signals in Space and Time , 2008, Science.

[15]  M. Rehm,et al.  Intracellular signaling dynamics during apoptosis execution in the presence or absence of X-linked-inhibitor-of-apoptosis-protein. , 2008, Biochimica et biophysica acta.

[16]  D. Thieffry,et al.  Functional organisation of Escherichia coli transcriptional regulatory network , 2008, Journal of molecular biology.

[17]  Gheorghe Craciun,et al.  Homotopy methods for counting reaction network equilibria. , 2007, Mathematical biosciences.

[18]  Eduardo Sontag,et al.  A Passivity-Based Approach to Stability of Spatially Distributed Systems With a Cyclic Interconnection Structure , 2007, IEEE Transactions on Automatic Control.

[19]  Darren W. Williams,et al.  Local caspase activity directs engulfment of dendrites during pruning , 2006, Nature Neuroscience.

[20]  Jochen H M Prehn,et al.  Systems analysis of effector caspase activation and its control by X‐linked inhibitor of apoptosis protein , 2006, The EMBO journal.

[21]  Eduardo D. Sontag,et al.  Diagonal stability of a class of cyclic systems and its connection with the secant criterion , 2006, Autom..

[22]  Nils Blüthgen,et al.  Mathematical Modeling Identifies Inhibitors of Apoptosis as Mediators of Positive Feedback and Bistability , 2006, PLoS Comput. Biol..

[23]  R. Eils,et al.  Mathematical modeling reveals threshold mechanism in CD95-induced apoptosis , 2004, The Journal of cell biology.

[24]  F. Allgöwer,et al.  Bistability Analyses of a Caspase Activation Model for Receptor-induced Apoptosis* , 2004, Journal of Biological Chemistry.

[25]  D. Green,et al.  The Pathophysiology of Mitochondrial Cell Death , 2004, Science.

[26]  Chaohong Sun,et al.  Regulation of the Apaf-1/Caspase-9 Apoptosome by Caspase-3 and XIAP* , 2003, The Journal of Biological Chemistry.

[27]  J Downward,et al.  Caspase-6 is the direct activator of caspase-8 in the cytochrome c-induced apoptosis pathway: absolute requirement for removal of caspase-6 prodomain , 2002, Cell Death and Differentiation.

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

[29]  P. Paraskevopoulos Modern Control Engineering , 2001 .

[30]  Xiaodong Wang,et al.  Smac, a Mitochondrial Protein that Promotes Cytochrome c–Dependent Caspase Activation by Eliminating IAP Inhibition , 2000, Cell.

[31]  E. Slee,et al.  Cleavage of BID during cytotoxic drug and UV radiation-induced apoptosis occurs downstream of the point of Bcl-2 action and is catalysed by caspase-3: a potential feedback loop for amplification of apoptosis-associated mitochondrial cytochrome c release , 2000, Cell Death and Differentiation.

[32]  A. M. Turing,et al.  The chemical basis of morphogenesis , 1952, Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences.

[33]  Daniel Liberzon,et al.  Common Lyapunov functions for families of commuting nonlinear systems , 2005, Syst. Control. Lett..

[34]  Ian Postlethwaite,et al.  Multivariable Feedback Control: Analysis and Design , 1996 .