RI-SINDy: Regulation-Informed Sparse Model Discovery in Biological Gene Networks

Abstract

Oscillatory and saturating dynamics in regulatory gene networks govern processes such as cell cycle timing, circadian rhythms, and inflammatory signaling, and sparse regression methods such as SINDy offer an interpretable method for recovering their governing equations from time-series data. However, these networks are typically governed by Hill-type nonlinearities that become nearly collinear with polynomial candidates over most bounds of an oscillating curve, and placing them in a shared library produces ill-conditioned regressions and unreliable coefficient estimates. Standard approaches such as ridge regularization, library normalization, and orthogonalization can marginally improve conditioning but cannot remove the structural multicollinearity of these systems, in turn limiting our ability to produce accurate parameter estimates. To address this issue, we propose the Regulation-Informed Sparse Identification of Nonlinear Dynamics (RI-SINDy) method, which removes known regulatory terms from the regression and recovers their coefficients directly from an asymptotic drift-balance condition. Because the time averaged derivative of a bounded trajectory approaches zero, the coefficients of regulatory nonlinearities can be computed independently by constraining them to offset the mean contribution of the remaining terms. We demonstrate that this method effectively reduces library collinearity and improves conditioning numbers by 3-45x, substantially improving coefficient recovery and regression accuracy. We test the RI-SINDy method on experimental Hes1 data and on synthetic Goodwin and NF-{kappa}B models. On Hes1, it achieves the lowest trajectory error among the compared methods (Traditional SINDy, SR3, and Nullcline Reconstruction). On the synthetic systems, it recovers the correct sparse structure with lower trajectory error than SR3 and Nullcline Reconstruction, and in a Goodwin noise sweep it recovers all active terms at noise levels up to 15%, where Traditional SINDy and E-SINDy fail to. In conclusion, the RI-SINDy method provides a robust framework to overcome the issue of multicollinearity for system identification in biological systems.