We propose a hybrid molecular dynamics/multi-particle collision dynamics model to simulate a set of self-assembled semiflexible filaments and free monomers. Further, we introduce a Monte Carlo scheme to deal with single monomer addition (polymerization) or removal (depolymerization), satisfying the detailed balance condition within a proper statistical mechanical framework. This model of filaments, based on the wormlike chain, aims to represent equilibrium polymers with distinct reaction rates at both ends, such as self-assembled adenosine diphosphate-actin filaments in the absence of adenosine triphosphate (ATP) hydrolysis and other proteins. We report the distribution of filament lengths and the corresponding dynamical fluctuations on an equilibrium trajectory. Potential generalizations of this method to include irreversible steps like ATP-actin hydrolysis are discussed.
© 2012 American Institute of Physics