RI-SINDy: Regulation-Informed Sparse Model Discovery in Biological Gene Networks
- 1University of Missouri-Kansas City
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.
Related
- Novel two-stage deep learning-based approach applied to gene expression data pertaining to esophageal adenocarcinoma boosting biological knowledge discovery
- A Taxonomy-Informed Sparse DNA Foundation Model for Microbial Genomics
- Physiology-Informed Sparse Identification of Dynamic Crop Models
- Distributed AI agents enable comprehensive discovery of biological insights at scale
- Operationalizing open-ended biological discovery across single-cell representations
BibTeX
@misc{lucenti2026risindy,
title = {{RI-SINDy: Regulation-Informed Sparse Model Discovery in Biological Gene Networks}},
author = {N. J. Lucenti and M. Bani-Yaghoub},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.04.756519},
url = {https://www.biorxiv.org/content/10.64898/2026.10.04.756519v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - RI-SINDy: Regulation-Informed Sparse Model Discovery in Biological Gene Networks
AU - N. J. Lucenti
AU - M. Bani-Yaghoub
PY - 2026
DA - 2026/10/04
PB - bioRxiv
DO - 10.64898/2026.10.04.756519
UR - https://www.biorxiv.org/content/10.64898/2026.10.04.756519v1
AB - 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.
ER -