From homogeneity to heterogeneity: Refining stochastic simulations of gene regulation

Cellular processes are intricately controlled through gene regulation, which is significantly influenced by intrinsic noise due to the small number of molecules involved. The Gillespie algorithm, a widely used stochastic simulation method, is pervasively employed to model these systems. However, this algorithm typically assumes that DNA is homogeneously distributed throughout the nucleus, which is not realistic. In this study, we evaluated whether stochastic simulations based on the assumption of spatial homogeneity can accurately capture the dynamics of gene regulation. Our findings indicate that when transcription factors diffuse slowly, these simulations fail to accurately capture gene expression, highlighting the necessity to account for spatial heterogeneity. However, incorporating spatial heterogeneity considerably increases computational time. To address this, we explored various stochastic quasi-steady-state approximations (QSSAs) that simplify the model and reduce simulation time. While both the stochastic total quasi-steady state approximation (stQSSA) and the stochastic low-state quasi-steady-state approximation (slQSSA) reduced simulation time, only the slQSSA provided an accurate model reduction. Our study underscores the importance of utilizing appropriate methods for efficient and accurate stochastic simulations of gene regulatory dynamics, especially when incorporating spatial heterogeneity.