The Multistrand Simulator: Stochastic Simulation of the Kinetics of Multiple Interacting DNA Strands

Abstract

DNA has been used in vitro as a computational substrate due to programmable base-pairing interactions. This allows the construction of logic gates, self-assembled shapes, motors, walkers and other nanoscale devices comprised of DNA strands. These strands must be carefully designed in order to correctly perform their function, which requires effective models of the DNA system. The thermodynamics of multiple interacting DNA strands is a well-studied model which can make equilibrium predictions on these systems, but many devices operate far from equilibrium. The prediction of large non-equilibrium systems requires the simulation of a kinetics model, as these are both stochastic and computationally difficult to solve in an exact form. We developed the Multistrand kinetics simulator, which extends the previous work by expanding the kinetics model to handle multiple nucleic acid strands, as well as using optimized algorithms in order to speed up the simulations.