The study of hematopoietic stem cell (HSCs) maintenance and differentiation to supply the hematopoietic system presents unique challenges, given the complex regulation of the process and the difficulty in observing cellular interactions in the stem cell niche. Quantitative methods and tools have emerged as valuable mechanisms to address this issue
however, the stochasticity of HSCs presents significant challenges for mathematical modeling, especially when bridging the gap between theoretical models and experimental validation. In this work, we have built a flexible and user-friendly stochastic dynamical and spatial model for long-term HSCs (LT-HSCs) and short-term HSCs (ST-HSCs) that captures experimentally observed cellular variability and heterogeneity. Our model implements the behavior of LT-HSCs and ST-HSCs and predicts their homeostatic dynamics. Furthermore, our model can be modified to explore various biological scenarios, such as stress-induced perturbations mediated by apoptosis, and successfully implement these conditions. Finally, the model incorporates spatial dynamics, simulating cell behavior in a 2D environment by combining Brownian motion with spatially graded parameters.