Simulating bistable biochemical systems by means of reactive multiparticle collision dynamics

Based on a reactive multiple particle collision method, we construct a mesoscopic dynamics model to simulate chemical system. The validity of the reactive multiple particle collision method under various conditions in a double-feedback bi-stable chemical system is studied. Then, we extend it to simulate diffusion-limited reactions with fast reaction rate in cellular environment. Using the improved method, we observe bi-stable behavior with randomly distributed reactants and spatial domain separation of opposite phases. The particle-based mesoscopic method is computationally efficient, although hydrodynamic interactions and fluctuation are both properly accounted for. Stochastic effects shown to play dominant roles in biochemical dynamics are also considered. The improved method could be used to explore a variety of reactions with disparate scale of reaction rates.