Description
Pluripotent stem cells (PSCs) exist in multiple stable states, each with specific cellular properties and molecular signatures. The process by which pluripotency is either maintained or destabilized to initiate specific developmental programs is poorly understood. We have developed a model to predict stabilized PSC gene regulatory network (GRN) states in response to combinations of input signals. While previous attempts to model PSC fate have been limited to static cell compositions, our approach enables simulations of dynamic heterogeneity by combining an Asynchronous Boolean Simulation (ABS) strategy with simulated single cell fate transitions using a Strongly Connected Components (SCCs). This computational framework was applied to a reverse-engineered and curated core GRN for mouse embryonic stem cells (mESCs) to simulate responses to LIF, Wnt/ß-catenin, FGF/ERK, BMP4, and Activin A/Nodal pathway activation. For these input signals, our simulations exhibit strong predictive power for gene expression patterns, cell population composition, and nodes controlling cell fate transitions. The model predictions extend into early PSC differentiation, demonstrating, for example, that a Cdx2-high/Oct4-low state can be efficiently generated from mESCs residing in a naïve and signal-receptive state sustained by combinations of signaling activators and inhibitors. Overall design: Examination of perturbed PSCs versus control PSCs and mesoderm progenitors Mouse pluripotent stem cells were grown on tissue culture plates for two days in serum-containing, feeder free medium supplemented with the following cytokines/small molecules: 2i = CHIR99021 (Reagents Direct 27-H76 – 3µM) & PD0325901 (Reagents Direct 39-C68 – 1µM) Jaki = JAK inhibitor (EMD Millipore 420097 – 2.0µM) BMP = BMP4 (R&D Systems 314-BP-010 – 10ng/ml) Alk5i = ALK5 inhibitor II (Cedarlane ALX-270-445 - 10µM)