Dynamic Gene Regulatory Network Inference with Interpretable, Biophysically-Motivated Neural ODEs

Abstract

Gene Regulatory Networks (GRNs) are complex dynamical systems that modulate gene expression and drive transitions between phenotypic cell states. Determining these networks is crucial in understanding how gene dysregulation can lead to phenotypic variation and perturbation responses. We present a novel biophysically-motivated neural ordinary differential equation (ODE) model framework with a biologically interpretable deep learning architecture that leverages dynamic single-cell data: in-CAHOOTTS (gene regulatory network Inference with Context Aware Hybrid neural-ODEs on Transcriptional Time-series Systems). Our approach combines accurate prediction with mechanistic interpretability, decomposing gene expression dynamics into fundamental biophysical processes (mRNA transcription and degradation) while inferring regulatory network structure. We validate in-CAHOOTTS by learning the Saccharomyces cerevisiae response to rapamycin treatment, and learning the dynamic cell cycle progression. The trained model accurately predicts gene expression trajectories and identifies biologically relevant rapamycin response regulators. The framework achieves stable long-term predictions with time frames extending 30 times beyond training data, maintaining realistic oscillatory dynamics for over 40 hours in cell cycle modeling, demonstrating that neural ODEs can capture true biological attractors rather than merely fitting data. We argue that interpretable neural ODEs can successfully model complex biological dynamics while revealing mechanistic insights essential for understanding living systems. By bridging extensible machine learning frameworks and mechanistic biology, in-CAHOOTTS represents a significant advance toward deep learning systems that both predict biological outcomes and explain the mechanisms driving them.