Modeling the dynamics of hypoxia inducible factor-1α (HIF-1α) within single cells and 3D cell culture systems

•A mathematical model is developed to describe single-cell HIF-1α-hypoxic response.•Oxygen measurements are obtained in tumorspheres subject to a hypoxic shock.•A mathematical model for the oxygen and HIF-1α dynamics in spheres is developed.•High HIF-1α levels are predicted at surface of sphere subject to hypoxic shock.•Differential dynamics of transcriptional targets of HIF-1 are predicted.

HIF (hypoxia inducible factor) is an oxygen-regulated transcription factor that mediates the intracellular response to hypoxia in human cells. There is increasing evidence that cell signaling pathways encode temporal information, and thus cell fate may be determined by the dynamics of protein levels. We have developed a mathematical model to describe the transient dynamics of the HIF-1α protein measured in single cells subjected to hypoxic shock. The essential characteristics of these data are modeled with a system of differential equations describing the feedback inhibition between HIF-1α and prolyl hydroxylases (PHD) oxygen sensors. Heterogeneity in the single-cell data is accounted through parameter variation in the model. We previously identified the PHD2 isoform as the main PHD sensor responsible for controlling the HIF-1α transient response, and make here testable predictions regarding HIF-1α dynamics subject to repetitive hypoxic pulses. The model is further developed to describe the dynamics of HIF-1α in cells cultured as 3D spheroids, with oxygen dynamics parameterized using experimental measurements of oxygen within spheroids. We show that the dynamics of HIF-1α and transcriptional targets of HIF-1α display a non-monotone response to the oxygen dynamics. Specifically we demonstrate that the dynamic transient behavior of HIF-1α results in differential dynamics in transcriptional targets.