Langevin dynamics simulation of DNA ejection from a phage

We have performed Langevin dynamics simulations of a coarse-grained model of ejection of dsDNA from Φ29 phage. Our simulation results show significant variations in the local ejection speed, consistent with experimental observations reported in the literature for both in vivo and in vitro systems. In efforts to understand the origin of such variations in the local speed of ejection, we have investigated the correlations between the local ejection kinetics and the packaged structures created at various motor forces and chain flexibility. At lower motor forces, the packaged DNA length is shorter with better organization. On the other hand, at higher motor forces typical of realistic situations, the DNA organization inside the capsid suffers from significant orientational disorder, but yet with long orientational correlation times. This in turn leads to lack of registry between the direction of the DNA segments just to be ejected and the direction of exit. As a result, a significant amount of momentum transfer is required locally for successful exit. Consequently, the DNA ejection temporarily slows down exhibiting pauses. This slowing down occurs at random times during the ejection process, completely determined by the particular starting conformation created by prescribed motor forces. In order to augment our inference, we have additionally investigated the ejection of chains with deliberately changed persistence length. For less inflexible chains, the demand on the occurrence of large momentum transfer for successful ejection is weaker, resulting in more uniform ejection kinetics. While being consistent with experimental observations, our results show the nonergodic nature of the ejection kinetics and call for better theoretical models to portray the kinetics of genome ejection from phages.