Direct Lyapunov exponent analysis enables parametric study of transient signalling governing cell behaviour

Computational models aid in the quantitative understanding of cell signalling networks. One important goal is to ascertain how multiple network components work together to govern cellular responses, that is, to determine cell 'signal-response' relationships. Several methods exist to study steady-state signals in the context of differential equation-based models. However, many biological networks influence cell behaviour through time-varying signals operating during a transient activated state that ultimately returns to a basal steady-state. A computational approach adapted from dynamical systems analysis to discern how diverse transient signals relate to alternative cell fates is described. Direct finite-time Lyapunov exponents (DLEs) are employed to identify phase-space domains of high sensitivity to initial conditions. These domains delineate regions exhibiting qualitatively different transient activities that would be indistinguishable using steady-state analysis but which correspond to different outcomes. These methods are applied to a physicochemical model of molecular interactions among caspase-3, caspase-8 and X-linked inhibitor of apoptosis - proteins whose transient activation determines cell death against survival fates. DLE analysis enabled identification of a separatrix that quantitatively characterises network behaviour by defining initial conditions leading to apoptotic cell death. It is anticipated that DLE analysis will facilitate theoretical investigation of phenotypic outcomes in larger models of signalling networks.