Abstract. Human liver-type phosphofructokinase 1 (PFKL) has been shown to regulate glucose flux as a scaffolder arranging glycolytic and gluconeogenic enzymes into a multienzyme metabolic condensate, the glucosome. However, it has remained elusive of how phase separation of PFKL is governed and initiates glucosome formation in living cells, thus hampering to understand a mechanism of glucosome formation and its functional contribution to human cells. In this work, we developed a stochastic model in silico using the principle of Langevin dynamics to investigate how biological properties of PFKL contribute to the condensate formation. The significance of an intermolecular interaction between PFKLs, an effective concentration of PFKL at a region of interest, and its own self-assembled filaments in formation of PFKL condensates and control of their sizes were demonstrated by molecular dynamics simulation using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS). Such biological properties that define intracellular dynamics of PFKL appear to be essential for phase separation of PFKL, which may represent an initiation step for the formation of glucosome condensates. Collectively, our computational study provides mechanistic insights of glucosome formation, particularly an initial stage through the formation of PFKL condensates in living human cells.