Matrix stability and bifurcation analysis by a network-based approach

In this paper, we develop a network-based methodology to investigate the problems related to matrix stability and bifurcations in nonlinear dynamical systems. By matching a matrix with a network, i.e., interaction graph, we propose a new network-based matrix analysis method by proving a theorem about matrix determinant under which matrix stability can be considered in terms of feedback loops. Especially, the approach can tell us how a node, a path, or a feedback loop in the interaction graph affects matrix stability. In addition, the roles played by a node, a path, or a feedback loop in determining bifurcations in nonlinear dynamical systems can also be revealed. Therefore, the approach can help us to screen optimal node or node combinations. By perturbing them, unstable matrices can be stabilized more efficiently or bifurcations can be induced more easily to realize desired state transitions. To illustrate feasibility and efficiency of the approach, some simple matrices are used to show how single or combinatorial perturbations affect matrix stability and induce bifurcations. In addition, the main idea is also illustrated through a biological problem related to T cell development with three nodes: TCF-1, GATA3, and PU.1, which can be considered to be a three-variable nonlinear dynamical system. The approach is especially helpful in understanding crucial roles of single or molecule combinations in biomolecular networks. The approach presented here can be expected to analyze other biological networks related to cell fate transitions and systematic perturbation strategy selection.