Arxiv:1907.00469V1 [Cond-Mat.Soft] 30 Jun 2019

Spontaneous Domain Formation in Spherically-Confined Elastic Filaments Tine Curk,1, 2, ∗ James Daniel Farrell,1, ∗ Jure Dobnikar,1, 3, 4, y and Rudolf Podgornik5, 1, z 1Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China. 2Faculty of Chemistry and Chemical Engineering, University of Maribor, 2000 Maribor, Slovenia 3Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China. 4Department of Chemistry, University of Cambridge, Cambridge, CB2 1EW United Kingdom. 5School of Physical Sciences and Kavli Institute for Theoretical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China. (Dated: July 2, 2019) Although the free energy of a genome packing into a virus is dominated by DNA-DNA interactions, ordering of the DNA inside the capsid is elasticity{driven, suggesting general solutions with DNA organized into spool-like domains. Using analytical calculations and computer simulations of a long elastic filament confined to a spherical container, we show that the ground state is not a single spool as assumed hitherto, but an ordering mosaic of multiple homogeneously-ordered domains. At low densities, we observe concentric spools, while at higher densities, other morphologies emerge, which resemble topological links. We discuss our results in the context of metallic wires, viral DNA, and flexible polymers. Introduction. Spatial organization of constrained elas- liquid-crystalline order of compacted DNA was found to tic filaments underlies a range of packing problems from play a crucial role promoting the formation of knots [19]. macroscopic wires in spherical cavities [1, 2] to genomes The slow dynamics of stiff chains in confinement in bacteriophages [3, 4]. A theoretical interpretation implies that the configurations observed in the non- of single molecule experiments on DNA packaging in equilibrium simulations crucially depend on the applied viruses [5] suggests single-domain spool-like structures protocol, e.g., on how the external force is applied in with axial symmetry as a prevailing morphology. How- packaging-motor driven DNA encapsidation in bacterio- ever, recent experiments report a range of different struc- phages [14]. It is therefore not surprising that, even in tures with inhomogeneously-ordered genomes, defects, cases where the models are very similar, a wide variety and phase transitions between the packings. Cryo- of behaviours is reported. Here, we focus on the na- electron microscopy of bacteriophages [6, 7] at vari- ture of the ground-state packing configurations. While ous physicochemical conditions indicates a range of pos- the present analysis thus answers a different question sible morphologies of DNA within the capsid: from from previous works, its results should nevertheless pro- isotropically disordered, cholesteric liquid crystalline or- vide a well-defined background scenario with which non- dering with a disordered core, to uniform, concentric equilibrium configurations can be compared and assessed. spools with local hexagonal order, reminiscent of the Theoretical Model. We consider a spherical enclosure nested spools first proposed in [8]. Measurements of C with radius R0 (volume V0), and a hard semi-flexible DNA dynamics during packaging [9], and intermittent cylinder with cross-section σ, length L, and a persistence ejection dynamics [10] suggest multi-domain structures. length lp = Kc=kBT , where Kc is the bending rigidity. Similarly, experiments and simulations of quasi-one- The chain cannot overlap with itself nor with the confin- dimensional elastoplastic filaments forced into spherical ing wall and it is torsionally relaxed; its elastic energy is confinement [1, 2], as well as the theoretical description of 1 R L 2 described by the integral Eel = 2 Kc 0 ds=R (s), with elastic filaments within the continuum polymer nematic R(s) the local radius of curvature. theory [11, 12], indicate a multidomain structures with a We explore the limit of a rigid, long, thin chain, where mixture of local disorder and nematic order. the elastic energy dominates over the configurational en- tropy. In this limit, the chain can be described by the Coarse-grained simulations of confined semiflexible arXiv:1907.00469v1 [cond-mat.soft] 30 Jun 2019 continuum theory of polymer nematics [11, 12], where polymers mimicking DNA in viral shells have been per- the polymer is described by a nematic-order vector field formed at various levels of detail. DNA interactions and or local polymer \current" t(r) [12] with ^t = t=t the unit entropic contributions need to be considered to explain tangent vector (director) to the chain. The total amount the energetics of confinement [13, 14], but the conforma- of material is V = R t dr, with t ≡ jtj the local poly- tion of the confined DNA is dictated by elasticity [15]. p C mer density. Due to the excluded volume of the chain, Packing into icosahedral capsids was shown to have a t is bounded by the densest (hexagonal) packing of rigid strong stochastic component [16], and structures without p m π 3 spool-like symmetry have been reported, including multi- disks in two dimensions: (t = 6 ≈ 0:907) [20]. The domain spool packing, observed when DNA is pushed ratio between the actual amount of polymer in thep capsid 3 3σ2 into the capsid through a portal [17]. The effect of elas- and its theoretical maximum, ν ≡ Vp=(tmV0) = 3 L, 8πR0 tic kinks on packing has been discussed in [18], while local is a key dimensionless parameter in the model. 2 The polymer current is only non-zero inside the capsid (a) 1 spool (b) 2 spools (c) 3 spools (d) Hopf link and has to be solenoidal (r · t = 0) except at the polymer end-points (the contribution of which vanishes in the long-chain limit). The lowest non-trivial order in the ex- pansion of the elastic free energy of the divergence-free z polymer current is given by [11, 12] (e) y Z 2 x 2Kc 1 ^ ^ E, 2 spools Eel = 2 t t × (r × t) dr ; (1) πσ V E, 3 spools where twist deformations (^t · (r ×^t)) are by assumption 1s n 0.9 excluded. The problem is to find the solenoidal field t(r) 0 0.2 0.4 0.6 0.8 11 E / within C which minimizes Eel. 0.8 Assuming axial symmetry, the solution is an spool [20]. w 0.6 w1, 2 s. However, relaxing the symmetry constraints might lead 0.8 w , 3 s. 1 0.4 to more complex structures such as those seen in experi- w2, 3 s. Elastic energy ments [7, 9], theory [11, 12], and simulations [15, 16, 18]. 0.2 Partitioning Indeed, as we show later, a single-domain spool is al- 0 most never the lowest elastic energy configuration. For 0.7 0 0.2 0.4 0.6 0.8 1 solutions containing more than one domain, t(r) is not Packing fraction n continuous in space, but it can still be solenoidal if at discontinuity interfaces S, the local field on both sides FIG. 1. Theoretical model. Visualization of the trial struc- is perpendicular to the normal, t · n(S) = 0. Multiple- tures: single spool (a) with polymer-occupied volume blue, domain structures should therefore be considered. Their double (b) and triple (c) spools (only the inner spools shown advantage is flexibility in arranging the material so as in green and orange, respectively), and the Hopf link (d). to avoid highly curved regions close to the capsid center, The structures shown are the solutions with optimal partitioning w at ν = 0:9. (e): The optimal ratio of elastic en- which cannot be done in single spool structures. 2s;3s 1s ergies Eel =Eel as a function of ν for double spools (green With this in mind, we assume that the field is parti- circles) and triple spools (red diamonds). Inset: The optimal Pnd tioned into nd domains, ν = 1 νi, with relative sizes partitioning as a function of ν for double spools (w1, green wi ≡ νi/ν. For nd ≤ 3, we construct candidate struc- stars) and triple spools (w1;2, red triangles). tures and search for the optimum partitioning fmig at several packing fractions 0 < ν < 1. Inspired by obser- q vations [1, 2, 15{17, 21], we consider four classes of trial 2=3 R0 1 − ν1 is composed of two small hemi-toroidal structures depicted in Fig. 1: single spools (a), double t caps of combined length proportional to ν2, connected spools (b) comprising two nested tori with perpendicular s t s by straight sections (ν2): ν2 = ν2 + ν2. The elastic en- axes, triply nested spools (c) with perpendicular axes, ergy is localized in the hemispheres and is determined by and|deviating from the nested-spool paradigm|Hopf νt alone (see section SI(B)): links (d), with a symmetric partitioning into two equiva- 2 lent, perpendicular, topologically-linked tori. We assume 2s 1s 1s t Eel (ν1; ν2;R0) = Eel (ν1;R0) + Eel (ν2;R1): (3) that the field density is tm within the defined domains 2s 1s and zero elsewhere and approximated (1) by a sum of in- Interestingly, at a given ν, the ratio Eel =Eel is com- tegrals over independent domains. We neglect the short pletely determined by the partitioning of the material segments connecting the domains, arguing that this con- w1. We calculate the optimal w1 for several values of ν tribution is negligible in the thin, long chain limit. by minimizing (3). The calculation for triple spools is 1s p 2 2 For a single spool, ^t = tm (−y; x; 0) = x + y , we analogous: ν = ν1 +ν2 +ν3, and the elastic energy a sum have [5]: over two outer spools (3) and an innermost spool with an 3s elaborate shape (see section SI(C)): Eel (ν1; ν2; ν3;R0) = 2 " # 3rdspool π KR 1=3 2s 1s p 0 1=3 1 1+ν Eel (ν1; ν2;R0) + Eel (ν1; ν2; ν3;R0).

Load more