DNA electrophoresis in confined, periodic geometries: a new lakes-straits model

We present a method to study the dynamics of long DNA molecules inside a cubic array of confining spheres, connected through narrow openings. Our method is based on the coarse-grained, lakes-straits model of Zimm and is therefore much faster than Brownian dynamics simulations. In contrast to Zimm's approach, our method uses a standard stochastic kinetic simulation to account for the mass transfer through the narrow straits and the formation of new lakes. The different rates, or propensities, of the reactions are obtained using first-passage time statistics and a Monte Carlo sampling to compute the total free energy of the chain. The total free energy takes into account the self-avoiding nature of the chain as well as confinement effects from the impenetrable spheres. The mobilities of various chains agree with biased reptation theory at low and high fields. At moderate fields, confinement effects lead to a new regime of reptation where the mobility is a linear function of molecular weight and the dispersion is minimal.