Abstract: | This paper describes inverse eigenvalue problems that arise in
studying qualitative dynamics in systems biology models.
An
algorithm based on lift-and-project iterations is proposed, where the lifting step entails solving a constrained matrix inverse
eigenvalue problem. In particular, prior to carrying out the iterative steps, $a$-$priori$ bounds on
the entries of the Jacobian matrix are computed by relying on the reaction network structure as well as the form of
the rate law expressions for the model under consideration. Numerical results on a number of models show that the proposed
algorithm can be used to computationally explore the possible dynamical scenarios while identifying the important
mechanisms via the use of sparsity-promoting regularization. |