Deterministic and stochastic aspects of VEGF-A production and the cooperative behavior of tumoral cell colony

A model is proposed to study the process of hypoxia-induced angiogenesis in cancer cells. The model accounts for the role played by the vascular endothelial growth factor (VEGF)-A in regulating the oxygen intake. VEGF-A is dynamically controlled by the HIF-1alpha concentration. If not degraded, HIF-1alpha can bind to the subunit termed HIF-1beta and so experience translocation to the nucleus, to exert its proper transcriptional activity. The delicate balance between these opposing tendencies translates into the emergence of distinct macroscopic behaviors in terms of the associated molecular concentrations that we here trace back to normoxia, hypoxia and death regimes. These aspects are firstly analyzed with reference to the ideal mean-field scenario. Stochastic fluctuations are also briefly discussed and shown to seed a cooperative interaction among cellular units, competing for the same oxygen reservoir.

💡 Research Summary

The paper presents a comprehensive mathematical framework for the hypoxia‑driven angiogenic response of cancer cells, focusing on the interplay between HIF‑1α, HIF‑1β, and VEGF‑A. The authors first construct a set of nine elementary reactions (Eqs. 1‑9) that capture oxygen diffusion, hydroxylase activation, HIF‑α degradation, HIF‑α/HIF‑β dimerisation, VEGF‑A synthesis, VEGF export, VEGF‑mediated oxygen recruitment, and degradation of both VEGF and HIF species. By introducing “empty” compartments (E₀, Eᵢ) they enforce a strict conservation of total molecular count, a device that facilitates stochastic treatment.

Using the van Kampen system‑size expansion, the stochastic master equation is reduced to deterministic mean‑field ordinary differential equations (Eqs. 11‑18) for the concentrations of O₂, inactive/active hydroxylases (W, Wₐ), HIF‑α, HIF‑β, intracellular and extracellular VEGF (Vi, V₀), and empty sites. A global conservation law (Eq. 19) reduces the dimensionality. Fixed‑point analysis of this deterministic system yields three biologically relevant steady states: (i) a normoxic equilibrium with low HIF‑α and VEGF, (ii) a hypoxic equilibrium where HIF‑α accumulates, drives VEGF production, and VEGF₀ recruits additional oxygen, and (iii) a death state where oxygen is exhausted and all species collapse to empty sites. Linear stability analysis shows how the basins of attraction depend on kinetic parameters a‑f (reaction rates) and the synthesis rates gα, gβ.

To capture finite‑size effects, the authors implement Gillespie stochastic simulations. These reveal that demographic noise can qualitatively alter system behavior, especially in small cell clusters. While the deterministic model predicts a single cell either remains in a normoxic or hypoxic steady state, stochastic runs show that occasional VEGF₀ release by one cell can create a shared oxygen pool that benefits neighboring cells. This “cooperative interaction” lowers HIF‑α levels across the cluster, suppresses further VEGF synthesis, and increases overall survival probability. The cooperative effect disappears in the continuum limit (N → ∞), highlighting the essential role of intrinsic noise in generating multistability and collective adaptation.

The authors discuss the biological implications: tumor microenvironments are intrinsically heterogeneous, and stochastic fluctuations in molecular numbers can foster cooperative angiogenic responses that enable tumor cells to survive under severe hypoxia. Consequently, anti‑VEGF therapies that target only the average VEGF production may be insufficient; therapeutic strategies must consider the stochastic, cooperative dynamics of cell populations.

In summary, the study integrates deterministic mean‑field analysis with stochastic simulations to elucidate how VEGF‑A production is regulated under hypoxia and how intrinsic noise can induce cooperative behavior among tumor cells competing for limited oxygen. The findings provide a quantitative basis for re‑evaluating anti‑angiogenic treatment approaches, emphasizing the need to account for population‑level stochasticity and inter‑cellular cooperation in cancer therapy design.