Simulations of Viral DNA Packaging and Ejection: Geometrical Order and Topological Disorder

Abstract

The packing of DNA inside bacteriophages arguably yields the simplest example of genome organisation in living organisms. Cryo-em studies showed that DNA in bacteriophages epsilon-15 and phi-29 is neatly ordered in concentric shells close to the capsid wall, while an increasing level of disorder was measured when moving away from the capsid internal surface. On the other hand the detected spectrum of knots formed by DNA that is circularised inside the P4 viral capsid showed that DNA tends to be knotted with high probability, with a knot spectrum characterized by complex knots and biased towards torus knots and against achiral ones. Existing coarse-grain DNA models, while being capable of reproducing the salient physical aspects of free, unconstrained DNA, are not able to reproduce the experimentally observed features of packaged viral DNA. Here we show, using stochastic simulation techniques, that both the shell ordering and the knot spectrum can be reproduced quantitatively if one accounts for the preference of contacting DNA strands to juxtapose at a small twist angle, as in cholesteric liquid crystals. The DNA knots we observe are strongly delocalized and, intriguingly, this is shown not to interfere with genome ejection out of the phage.