The Hypothesis of an Undiscovered Master Regulator in the Pluripotency Gene Regulatory Network: Mathematical Modeling, Stochastic Simulations, and Nanotechnology-Enabled Discovery

Low reprogramming efficiency in induced pluripotent stem cells (iPSCs) generated with Yamanaka factors indicates potential missing elements in the core pluripotency gene regulatory network (GRN). This theoretical-computational study hypothesizes an undiscovered master regulator (X) that exerts strong positive feedback on Nanog, enhancing network stability and reprogramming kinetics. Building on established models, we develop a reproducible ordinary differential equation (ODE) framework with transient external drive to mimic factor overexpression. Deterministic simulations reveal elevated and more robust Nanog steady-states with X. Stochastic extensions using stochastic differential equations (SDEs) demonstrate noise reduction, addressing cellular heterogeneity. Global sensitivity analysis quantifies robustness, phase plane analysis shows expanded pluripotent basins of attraction, and bifurcation diagrams confirm enhanced bistability. Parameters are rigorously justified from literature-reported protein half-lives and nondimensionalized. A bioinformatic analysis of public RNA-seq datasets identifies candidates (e.g., SRSF3, Esrrb) with expression patterns matching model predictions. Full Python code ensures reproducibility. We propose concrete nanotechnology-enabled strategies (nanoparticle-based pull-down assays) for empirical discovery of X.