Model and computer simulations of the motion of DNA molecules during pulse field gel electrophoresis

A model is presented for the motion of individual molecules of DNA undergoing pulse field gel electrophoresis (PFGE). The molecule is represented by a chain of charged beads connected by entropic springs, and the gel is represented by a segmented tube surrounding the beads. This model differs from earlier reptation/tube models in that the tube is allowed to leak in certain places and the chain can double over and flow out of the side of the tube in kinks. It is found that these kinks often lead to the formation of U shapes, which are a major source of retardation in PFGE. The results of computer simulations using this model are compared with real DNA experimental results for the following cases: steady field motion as seen in fluorescence microscopy, mobility in steady fields, mobility in transverse field alternation gel electrophoresis (TFAGE), mobility in field inversion gel electrophoresis (FIGE), and linear dichroism (LD) of DNA in agarose gels during PFGE. Good agreement between the simulations and the experimental results is obtained.