Abstract
Particle-based simulators represent molecules of interest with point-like particles that diffuse and react in continuous space. These simulators are often used to investigate spatial or stochastic aspects of biochemical systems. This paper presents new particle-based simulation algorithms for modeling interactions between molecules and surfaces; they address irreversible and reversible molecular adsorption to, desorption from and transmission through membranes. Their central elements are: (i) relationships between adsorption, desorption and transmission coefficients on the one hand, and simulator interaction probabilities on the other, and (ii) probability densities for initial placements of desorbed molecules. These algorithms, which were implemented and tested in the Smoldyn simulator, are accurate, easy to implement and computationally efficient. They allow longer time steps and better address reversible processes than an algorithm that Erban and Chapman recently presented (Physical Biology 4:16-28, 2007). This paper also presents a method for simulating unbounded diffusion in a limited spatial domain using a partially absorbing boundary, as well as new solutions to the diffusion differential equation with reversible Robin boundary conditions.