Switching, Multiple Time-Scales and Geometric Blow-Up in a Low-Dimensional Gene Regulatory Network

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Jelbart, Samuel, Kristiansen, Kristian Uldall, Szmolyan, Peter
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866918156078940160
author Jelbart, Samuel
Kristiansen, Kristian Uldall
Szmolyan, Peter
author_facet Jelbart, Samuel
Kristiansen, Kristian Uldall
Szmolyan, Peter
contents ODE-based models for gene regulatory networks (GRNs) can often be formulated as smooth singular perturbation problems with multiple small parameters, some of which are related to time-scale separation, whereas others are related to 'switching' (proximity to a non-smooth singular limit). This motivates the study of reduced models obtained after (i) quasi-steady state reduction (QSSR), which utilises the time-scale separation, and (ii) piecewise-smooth approximations, which reduce the nonlinearity of the model by viewing highly nonlinear sigmoidal terms as singular perturbations of step functions. We investigate the interplay between the reduction methods (i)-(ii), in the context of a 4-dimensional GRN which has been used as a low-dimensional representative of an important class of (generally high-dimensional) GRN models in the literature. We begin by identifying a region in the small parameter plane for which this problem can be formulated as a smooth singularly perturbed system on a blown-up space, uniformly in the switching parameter. This allows us to apply Fenichel's coordinate-free theorems and obtain a rigorous reduction to a 2-dimensional system, that is a perturbation of the QSSR. Finally, we show that the reduced system features a Hopf bifurcation which does not appear in the QSSR system, due to the influence of higher order terms. Taken together, our findings suggest that the relative size of the small parameters is important for the validity of QSS reductions and the determination of qualitative dynamics in GRN models more generally. Although the focus is on the 4-dimensional GRN, our approach is applicable to higher dimensions.
format Preprint
id arxiv_https___arxiv_org_abs_2504_20467
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Switching, Multiple Time-Scales and Geometric Blow-Up in a Low-Dimensional Gene Regulatory Network
Jelbart, Samuel
Kristiansen, Kristian Uldall
Szmolyan, Peter
Dynamical Systems
34E13, 34E15, 34E10, 94C11, 37N25
ODE-based models for gene regulatory networks (GRNs) can often be formulated as smooth singular perturbation problems with multiple small parameters, some of which are related to time-scale separation, whereas others are related to 'switching' (proximity to a non-smooth singular limit). This motivates the study of reduced models obtained after (i) quasi-steady state reduction (QSSR), which utilises the time-scale separation, and (ii) piecewise-smooth approximations, which reduce the nonlinearity of the model by viewing highly nonlinear sigmoidal terms as singular perturbations of step functions. We investigate the interplay between the reduction methods (i)-(ii), in the context of a 4-dimensional GRN which has been used as a low-dimensional representative of an important class of (generally high-dimensional) GRN models in the literature. We begin by identifying a region in the small parameter plane for which this problem can be formulated as a smooth singularly perturbed system on a blown-up space, uniformly in the switching parameter. This allows us to apply Fenichel's coordinate-free theorems and obtain a rigorous reduction to a 2-dimensional system, that is a perturbation of the QSSR. Finally, we show that the reduced system features a Hopf bifurcation which does not appear in the QSSR system, due to the influence of higher order terms. Taken together, our findings suggest that the relative size of the small parameters is important for the validity of QSS reductions and the determination of qualitative dynamics in GRN models more generally. Although the focus is on the 4-dimensional GRN, our approach is applicable to higher dimensions.
title Switching, Multiple Time-Scales and Geometric Blow-Up in a Low-Dimensional Gene Regulatory Network
topic Dynamical Systems
34E13, 34E15, 34E10, 94C11, 37N25
url https://arxiv.org/abs/2504.20467