Highly Constrained Kinetic Models for Single-Cell Gene Expression Analysis

Advances in single-cell RNA sequencing (scRNA-seq) and high-resolution imaging techniques, such as single-molecule tracking (SMT) of RNA and transcription factors, allow researchers to quantitatively explore dynamics and variation but have never been integrated into a single coherent model. In this study, we propose a kinetic model that intakes multiple data types, including steady-state and time-resolved datasets, to simulate and fit stochastic models of gene transcription to experimental data. We find that 3-state models provide an essential improvement over the widely used 2-state model for most genes and have the property of kinetic proofreading, which we argue is advantageous in the cellular context. We further identify two dimensionless quantities derived from the rate equations which are broadly conserved across genes. Finally, we extend this model to scRNA-seq datasets to infer kinetic rates under defined perturbations and reveal biochemical insight into the mechanism of action of transcription factors.