A computational framework for a Lyapunov-enabled analysis of biochemical reaction networks

Persistence and stability of generalized ribosome flow models with time-varying transition rates

In this paper some important qualitative dynamical properties of generalized ribosome flow models are studied. Ribosome flow models known from the literature are generalized by allowing an arbitrary directed network structure between compartments, and by assuming general time-varying rate functions corresponding to the transitions. Persistence of the dynamics is shown using the chemical reaction network (CRN) representation of the system where the state variables correspond to ribosome density and the amount of free space in the compartments. The L1 contractivity of solutions is also proved in the case of periodic reaction rates having the same period. Further we prove the stability of different compartmental structures including strongly connected ones with entropy-like logarithmic Lyapunov functions through embedding the model into a weakly reversible CRN with time-varying reaction rates in a reduced state space. Moreover, it is shown that different Lyapunov functions may be assigned to the same model depending on the non-unique factorization of the reaction rates. The results are illustrated through several examples with biological meaning including the classical ribosome flow model on a ring.

A mathematical model exhibiting the effect of DNA methylation on the stability boundary in cell-fate networks

Cell-fate networks are traditionally studied within the framework of gene regulatory networks. This paradigm considers only interactions of genes through expressed transcription factors and does not incorporate chromatin modification processes. This paper introduces a mathematical model that seamlessly combines gene regulatory networks and DNA methylation (DNAm), with the goal of quantitatively characterizing the contribution of epigenetic regulation to gene silencing. The 'Basin of Attraction percentage' is introduced as a metric to quantify gene silencing abilities. As a case study, a computational and theoretical analysis is carried out for a model of the pluripotent stem cell circuit as well as a simplified self-activating gene model. The results confirm that the methodology quantitatively captures the key role that DNAm plays in enhancing the stability of the silenced gene state.