1924 Biophysical Journal Volume 107 October

Transcription

1 194 Biophysical Journal Volume 107 October Article Ejecting Phage DNA against Cellular Turgor Pressure Sanjin Marion 1, * and Antonio Siber 1 1 Institute of Physics, Zagreb, Croatia ABSTRACT We examine in vivo ejection of noncondensed DNA from tailed bacteriophages into bacteria. The ejection is dominantly governed by the physical conditions in the bacteria. The confinement of the DNA in the virus capsid only slightly helps the ejection, becoming completely irrelevant during its last stages. A simple calculation based on the premise of condensed DNA in the cell enables us to estimate the maximal bacterial turgor pressure against which the ejection can still be fully realized. The calculated pressure (~5 atm) shows that the ejection of DNA into Gram-negative bacteria could proceed spontaneously, i.e., without the need to invoke active mechanisms. INTRODUCTION The explanation of all the relevant (thermodynamical) forces guiding the ejection of dsdna from tailed bacteriophages into bacterial cells, despite 50 years of research, is still missing (1,). The ejection starts as a release of DNA from the fully packed capsid (protein coating of the bacteriophage). The DNA is packed to extreme densities and exerts a pressure of atm on the capsid (3,4). Models developed and tested in vitro (see Ponchon et al. (1) and Molineux and Panja () and references therein) predict that the ejecting force resulting from even such a dense packing is insufficient to completely transfer the DNA into the cell interior. Although cells have smaller turgor pressures than fully packed bacteriophages, the ejecting force (and pressure) in the capsid drops sharply as it empties (5). A recent single-molecule Hershey-Chase experiment (6) hints that in vivo ejection is controlled not by the amount of DNA left inside the capsid but by the amount ejected into the cell. This means that once the pressure built-up in the bacteriophage is spent on the DNA ejection, a cellular mechanism takes over, e.g., as seen in vivo for T5 (7) and T7 (8). There have been various proposed mechanisms for completion of the ejection: nanomotors or enzymes that ratchet in the stalled part of the DNA (8,9); a solvent flow through the semipermeable capsid and into the cell, simultaneously flushing the DNA trough the tail (10); osmotic pressure from proteins remaining in the capsid (11); and diffusion of DNA with assisted pulling from proteins in the cell (1). It appears that none of these models give a definite answer; experiments suggest a coexistence of many different mechanisms. We propose a scenario sufficient to explain ejection into some Gram-negative bacteria based only on thermodynamic considerations of DNA in a noncondensed state in the capsid and in a condensed state in the cell. Models of DNA ejection based on the continuum theory by Ubbink and Odijk (13,14) have been applied to explain the ejection when both the ejected and encapsidated DNA are condensed (15). In the model by Tzlil et al. (15), the capsid DNA, being confined, has a nonoptimal shape that generates the force ejecting it from the capsid until the ejection force is balanced by an external counterforce. If this in vitro model is applied to in vivo cases, it predicts that ejection stalls if the cellular turgor pressure is larger than ~0.5 atm. Typical cellular pressures are estimated to be atm for Gram-negative (16 18) and 15 5 atm for Gram-positive bacteria (17). In our approach, we maintain the hypothesis that DNA is condensed in the cell due to the presence of multivalent cations during ejection (the concept of a condensate in the cell during ejection is not incompatible with a high reaction rate (19), and accordingly, a high DNA transcription rate after complete ejection). Typical in vitro experiments cannot obtain a clean separation between two (three) compartments, one in the virus (immersed in the extracellular fluid), and the other in the cell (0), which is a key feature of our model. The results obtained in vitro may thus have a limited applicability in the in vivo conditions. Condensing agents are found in some capsids, but this is not a general situation and has been related to reduced infectivity (1). Many viral capsids are permeable to small ions, so we expect intracapsid solvent conditions analogous to those in their exterior. Because typical extracellular fluids have no condensing agents in sufficient concentrations, DNA ejection into bacteria is expected to proceed from a noncondensed state. MODELING THE DNA Submitted April 8, 014, and accepted for publication September 5, 014. *Correspondence: smarion@ifs.hr Editor: Jason Kahn. The phage DNA is one long strand of total length L 0 and persistence length L p z 50 nm () able to move between two compartments with different thermodynamic Ó 014 by the Biophysical Society /14/10/194/6 $.00

2 Ejecting Phage DNA against Cellular Turgor Pressure 195 conditions the virus and the cell. It is thus partitioned in two pieces, one of length L inside the virus head and tail, and the other of length L 0 L in the cell. Assuming no attractive DNA-DNA interaction in the capsid (noncondensed DNA), the total force on the DNA will tend to eject it from the capsid. There are two major contributions to the free energy F cell of the DNA in the cell. Because of the presence of multivalent ions and osmolytes, the ejected DNA will condense. The shape of the condensate is a result of a competition between DNA-DNA attractions mediated by the cations and unfavorable effects of DNA bending. On the other side, the major contribution to the free energy inside the capsid comes from bacteriophage tail confinement, bending from the capsid confinement, and excluded volume interactions between DNA segments. Ubbink-Odijk model for condensed DNA in the cell The state of the DNA inside the cell is described following the outline by Tzlil et al. (15). The volume V of the condensed DNA is proportional to the contour length L of the DNA V ¼ A 0 (L 0 L), with A 0 as the area per unit length of DNA. The DNA condensate is treated in the Ubbink-Odijk model (13,14) with free energy F tor ¼ gv þ ss þ E bend consisting of: a negative (attractive) part proportional to the DNA volume gv, where g is an effective condensing pressure; a positive contribution proportional to the surface area ss, where s reflects a lack of DNA neighbors at the condensate surface; and a bending energy of k B TL p /R per unit length of a circular loop of DNA wound at a radius R. The total free energy F tor is minimal when the condensate has the shape of a torus (15), but there are other possibilities like rods (3). A different shape would only change the surface and bending energy terms to some degree. The DNA in the condensate is assumed to be hexagonally packed with the area per unit length of DNA A 0 ¼ pd0 = p ffiffiffiffiffi 1 the packing constant for hexagonally packed cylinders (15). Here d 0 z.8 nm is the experimentally determined closest separation between DNA strands with added condensing agents (15,4). We note that Tzlil et al. (15) model the surface free energy contribution by assuming the loss of one-half of DNA nearest neighbors, whereas we assume the loss of one-third of the nearest neighbors as derived by Ubbink and Odijk (13). This constitutes a minor correction to the surface free energy term, but may become important if used to determine model parameters from toroid shapes. The approximation adopted by Tzlil et al. (15) assumes the toroid to have the shape of a thin torus. In this regime, the bending energy can be approximated as if the whole length of DNA is wound at a (mean) distance R to the center. We calculate the free energy of the condensate as for a torus with major radius R as F tor ðl; RÞ ¼ ga 0 ðl L 0 Þþ 3 þ k BT ðl L 0 ÞL p ; R A 0 g d pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffip ðl L 0 ÞA 0 and proceed with a variational approach. Minimizing Eq. 1 with respect to the major radius R yields optimal R, 3dLp 5ðL R ¼ p A 0 g ffiffiffiffiffiffiffi L0 Þ 1 5 ; () A 0 so that the minimal free energy of a toroidal condensate (15)is " # 1 A F tor ¼ e 0 ðl L 0 Þ 1 b 0 L 5 p ðl L d0e 4 0 Þ 5 ; (3) 0 where e 0 ¼ A 0 g is the condensation free energy per unit length in an ideally packed hexagonal lattice and b ¼ a constant. This e 0 was determined by Tzlil et al. (15) by fitting the toroid major axis to experimental data for l-dna in a solution of polylysine (3). Inasmuch as they chose a different surface term than that used by us and in the original model (14), the condensation free energy per unit length they obtain, g T, is smaller than ours, g UO. The shapes of toroids in the model of Ubbink-Odijk (14) are determined by a nondimensional parameter ffiffiffi R (1) a ¼ sv 1 3 d 0 L p ; (4) which depends on the choice of surface energy s. By using a different choice for the surface energy, the fits to toroid shapes would correspond to different values of the surface free energy s than in the original model. Because s is derived from g, we have that the Ubbink-Odijk model, adopted here also by us, has g UO ¼ 3/ g T ¼ 0.15 k B T/nm 3 (~6 atm). This difference arises solely from different treatment of the missing neighbors at the surface of the toroid (one-third missing in Ubbink and Odijk (14), g UO ~ 6 atm versus onehalf missing in Tzlil et al. (15), g T ~ 4 atm). An alternative approach to fitting is to obtain the condensing energy per unit length from osmotic force measurements (4), e.g., a 0 mm solution of Cobalt Hexamine corresponds to 0.04 k B T nm 3. The comparison of the two values obtained (0.15 k B T nm 3 vs k B T nm 3 ) shows that the thermodynamics of the condensed DNA importantly depends on the condensing agent. When most of the viral DNA is in the condensate, L << L 0, the surface and bending terms in Eq. 3 are negligible (15). This is because the largest contribution to the free energy of tori comes from the bulk term in F tor, as can be seen from 1 A 1:308 0 L 5 p d0 4e L 5 0 1: (5) 0

3 196 Marion and Siber Thus, the contributions of the surface and bending terms in the free energy are much smaller from the bulk contribution. To confirm the wider applicability of this statement, we performed a full minimization for toroids with different lengths of DNA, corresponding to different phage genomes. The minimization of free energy yielded the optimal shape together with the corresponding volume free energy F vol ¼ gv, surface free energy F surf ¼ ss, and the bending energy F bend. The optimal condensate shapes are shown in Fig. 1 with the corresponding energy contributions. We see that for all but the smallest phage, the correction to the turgor pressure that the phage can eject against, DP, is small. Excluded volume for two DNA rods in the capsid The dominant part of the excluded volume for two DNA molecules comes from screened electrostatic interactions. Due to its rigidity DNA behaves electrostatically as a long rod with a linear charge density. The interaction between two rods is obtained from the linearized Poisson-Boltzmann equation (Debye-Hückel). If the electrostatic screening length k and the DNA width are negligible with regard to the persistence length, then the interaction of two DNA rods that have their points of closest separation at r 1 and r, and a mutual angle of g between their directors n 1 and n, is given by (5,6) a b U 1 ðr 1 ; r Þ¼ pl ek e kjr 1 r j sin g ; (6) FIGURE 1 (a) Optimal DNA toroid shapes for several different lengths of DNA corresponding to well-known phages () (see text for details). Only the upper halves of the cross sections are shown. (b) The calculated correction to turgor pressure DP (see main text) for the DNA condensate in the cell as a function of genome length arising from the surface F surf and bending F bend free energies. where e is the solvent dielectric permittivity, and l is the effective linear charge density of the polyelectrolyte. The angle-dependent electrostatic excluded volume (second virial coefficient) is then calculated as (6,7) bðn 1 ; n Þ¼ 1 ZZ e U 1ðr 1 ;r Þ k B T 1 dr 1 dr ; (7) V which for the case of two charged rods of length L p at a closest separation of x ¼jr 1 r j evaluates to Z N bðn 1 ; n Þ¼L p 1 sin g e U 1 ðx;gþ k B T dx; (8) with D as the width of the DNA basepair. Based on the ambiguity in choosing the effective linear charge density attributable to correlations and screening, we opt to use an experimental fit that includes both hydration repulsion and screened electrostatic interactions between two parallel rods of DNA (8). We neglect the short-range hydration interactions and use only the electrostatic component of the fit D UðrÞ ¼b p e kr ffiffiffiffi: (9) kr We extend this potential to apply it to skewed rods following a scheme elaborated in Brenner and Parsegian (5) to obtain the angle-dependent potential rffiffiffi p Uðx; gþ ¼ b k e kx sin g : (10) From this, we calculate the angle-dependent excluded volume b(n 1, n ) with the value of b ¼ 85 k B T/nm for 100 mm added monovalent salt (8). The isotropic excluded volume v is obtained by integrating b(n 1, n ) over all possible angles for n 1 and n, assuming equal probability for finding them in any orientation. From v ¼ 1 ZZ bðn 1 ; n Þ 1 1 4p 4p du 1dU ; (11) we obtain the excluded volume of one DNA persistencelength segment (6) asv z L p DC 0, with D z.5 nm as the width of the DNA and C 0 z.87 as a constant depending on interaction specifics (100 mm monovalent salt concentration). A comparison with a viral capsid with R c ¼ 50 nm shows that v z V c /30; we see that 30 persistence-length segments of DNA are sufficient to exclude all the free volume, hinting at a large electrostatic intrastrand repulsion from confinement. The anisotropic excluded volume approximation is valid until the viral DNA transitions into a liquid crystal phase (6), which happens at roughly 30% packing (9) well beyond the encapsidated length for which the ejecting force is comparable to the opposing cellular force.

4 Ejecting Phage DNA against Cellular Turgor Pressure 197 RESULTS AND DISCUSSION We now study the balance of forces near the end of ejection so we can determine the maximal cellular turgor pressure that can be overcome. The crowded cellular interior exerts a turgor pressure P on the volume V of any foreign material to banish it from the cell (30) whereas the effect of condensation tries to draw more DNA into the condensate. The balance of free energy in the cell is thus F cell (L) ¼ (P g)a 0 (L 0 L), which is always negative if P < g. The DNA will tend to enter the cell despite the turgor pressure, due to the favorable condensation conditions. The value of g can be determined by fitting minimal toroid shapes to experimentally measured toroid parameters in the presence of condensing agents found in bacteria (15,3) (see previous section). Obtained estimates for g are between 4 and 6 atm, depending on the choice of model parameters. We take g z 4 atm for the turgor pressure that condensing agents in a cell could overcome on their own. Additional contributions, DP, come from unfavorable free energy contributions for DNA in the capsid F cap and corrections to surface and bending energy terms in the condensate (Fig. 1). The additional chemical potential for the DNA in the cell when the turgor pressure is increased by DP is DPA 0, and this should be matched by the chemical potential in the capsid to avoid the stalling of the ejection. When the two chemical potentials are equal, m ¼ vf cell vl ¼ DPA 0 ¼ vf cap vl ; (1) the ejection will stall at some length L*. From this we determine the maximal additional DP that can be overcome by the virus because the DNA is ejected from the cramped capsid. We now study effects due to the confinement in the capsid, which was the cause of the driving force in the early stage of ejection. Some bacteriophages have tails of considerable length t, so we examine whether they influence the ejection process. The entropic penalty for confining a semiflexible polymer in a tube of diameter w (31) is F tail k B T zt l ln L p l ; (13) where l ¼ w /3 L p 1/3 is the Odijk deflection length. For a tail with w z.75 nm (), the effective chemical potential from the tail is F tail /t z 0.7 k B T/nm. This is enough to oppose an additional 1.5 atm of turgor pressure in the cell and is independent of the length of the tail. However, the effect of the tail onsets only when the last DNA basepair exits the capsid and enters the tail (i.e., when L ¼ t). This suggests a barrier in the chemical potential that needs to be overcome for total ejection when 4 atm < P < 5.5 atm (as will be shown later). DNA is a charged polyelectrolyte with strong repulsive electrostatic interactions between any two points on its contour. Interactions between nearby parts of the contour act to give it its large persistence length comparable to the radius of the nearly spherical capsid L p ~ R c. When the DNA touches the capsid, any increase of length L will force the DNA to bend to conform to the shape of the capsid. The bending energy in such a situation can be approximated by that of a loop of DNA with radius R c, k B TL p /R c.a comparison with the previously neglected bending energy of the condensate (see Eq. 3) reveals that they are matched for R c z 50 nm according to the thin torus model (15). Smaller capsids could enhance the chemical potential; e.g., in the case of l-phage with R c z 30 nm (), the change of Dm z 0.03 k B T/nm is enough to compensate for z0.15 atm of turgor pressure. Any direct interactions between the DNA and capsid appear to be negligible; viral ejection experiments show no evidence of attractive forces (3), and because dsdna bacteriophage capsids have no considerable charge (33), only weak van der Waals interactions are possible. When the length of the DNA in the capsid is large enough, L >> L p, the DNA chain statistics resemble that of a random walk of n p ¼ L/L p persistence-length segments (34). This approximation is valid for steric interactions in bulk as long as L p % R c (35). In the case of strong spherical confinement and electrostatic interactions, the effective persistence length that governs correlations along the contour becomes as small as 10 nm (36), making the approach valid even in smaller capsids. The interaction energy between different parts of the DNA strand in confinement may be estimated on the basis of the excluded volume v between two segments (35). The corresponding Flory free energy of interaction in the capsid of volume V c is F v z k B Tn p v/v c. This contribution vanishes inasmuch as L / 0, so it cannot help the ejection in its latest stages. The excluded volume between two DNA segments can be approximated as that between two charged rods. This interaction is intrinsically anisotropic, but at low packing fractions (near the end of ejection) there is no order, and we can average this over all possible mutual angles between two cylinders (see Modeling the DNA). We obtain v ¼ L p DC 0, where D is the DNA diameter (D z.5 nm) and C 0 is a numeric constant. The excluded volume interactions will contribute to the total free energy as k B TC 0 DL /V c (but only in the regime when there are at least several persistence length segments inside the capsid). For R c z 50 nm in 100 mm monovalent salt, v/v c z 1/30, resulting in the effective chemical potential being an increasing function of length as vf v /vl z L/L p k B T/nm. If, say, 10 persistence lengths of DNA are in the capsid, the repulsive force is sufficient to oppose ~0.8 atm of turgor pressure. When the cellular turgor pressure P is larger than the effective condensing pressure g in the cell, the net driving pressure P g > 0 on the viral DNA will tend to repel it from the cell. When the net repulsive cellular pressure is smaller than the tail confinement penalty m tail ¼ F tail /t

5 198 Marion and Siber (corresponding to z1.5 atm), the DNA will be stuck in the virus; the tail of length t will be completely filled, and some length L t will reside inside the capsid. This stalling length is a result of all the repulsive interactions in the capsid canceling out with the net driving pressure in the cell. Note, however, that if the whole DNA from the capsid (of length L t) enters the tail, the additional asymmetry in the free energy of the two thermodynamic reservoirs onsets. This is due to the confinement penalty of the DNA in the tail. With a partially filled tail, the derivatives of the free energy per unit length (the chemical potential) in the virus (the DNA length increases in the virus) and in the cell (the DNA length increases in the cell) are not the same. The thermodynamical balance is thus broken, and the thermodynamical gradient necessary for the ejection is restored. Therefore, there exists a potential barrier that the DNA needs to overcome for its capsid-side end to enter the tail and be swiftly ejected. We now estimate whether the thermal fluctuations may overcome the free energy barrier. Because the capsid DNA is not condensed it is coupled to a solvent heat bath at temperature T. From the equipartition theorem, the encapsidated DNA will have ~1/ k B T thermal energy per degree of freedom. A semiflexible polymer of length L can be partitioned into a random walk of n p z L/L p steps with each step of length L p having two degrees of freedom (two angles) and the origin being at the tail entrance. The resulting DNA thermal energy is ~1 k B T/L p, or ~0.0 k B T/nm, which corresponds to a fluctuation in the maximal turgor pressure of ~0.1 atm. We can argue that the ejection can happen in a finite time if the barrier corresponds to up to, say, three standard deviations of ~0.3 atm. CONCLUSION We conclude that the mechanism of viral ejection from some dsdna bacteriophages into Gram-positive bacteria could be explained as a competition between the resisting turgor pressure and a free energy gain from condensation of the ejected part of the DNA. From the experimental data for the condensation energy of DNA, we estimate that a turgor pressure in excess of 4 atm can be overcome by unassisted ejection in line with recent moleculardynamics simulations (37). Our model does not exclude additional ejecting mechanisms such as the osmotic pressure from proteins remaining in the capsid (11) or pulling from proteins in the cell (1); such mechanisms can help to overcome even larger turgor pressures than obtained here. REFERENCES 1. Ponchon, L., S. Mangenot,., L. Letellier Encapsidation and transfer of phage DNA into host cells: from in vivo to single particles studies. Biochim. Biophys. Acta. 174: Molineux, I. J., and D. Panja Popping the cork: mechanisms of phage genome ejection. Nat. Rev. Microbiol. 11: Smith, D. E., S. J. Tans,., C. Bustamante The bacteriophage straight f9 portal motor can package DNA against a large internal force. Nature. 413: Siber, A., A. L. Bozic, and R. Podgornik. 01. Energies and pressures in viruses: contribution of nonspecific electrostatic interactions. Phys. Chem. Chem. Phys. 14: Purohit, P. K., M. M. Inamdar,., R. Phillips Forces during bacteriophage DNA packaging and ejection. Biophys. J. 88: van Valen, D., D. Wu,., R. Phillips. 01. A single-molecule Hershey-Chase experiment. Curr. Biol. : Mangenot, S., M. Hochrein,., L. Letellier Real-time imaging of DNA ejection from single phage particles. Curr. Biol. 15: Kemp, P., M. Gupta, and I. J. Molineux Bacteriophage T7 DNA ejection into cells is initiated by an enzyme-like mechanism. Mol. Microbiol. 53: Molineux, I. J No syringes please, ejection of phage T7 DNA from the virion is enzyme driven. Mol. Microbiol. 40: Panja, D., and I. J. Molineux Dynamics of bacteriophage genome ejection in vitro and in vivo. Phys. Biol. 7: Grayson, P., and I. J. Molineux Is phage DNA injected into cells biologists and physicists can agree. Curr. Opin. Microbiol. 10: Inamdar, M. M., W. M. Gelbart, and R. Phillips Dynamics of DNA ejection from bacteriophage. Biophys. J. 91: Ubbink, J., and T. Odijk Polymer- and salt-induced toroids of hexagonal DNA. Biophys. J. 68: Ubbink, J., and T. Odijk Deformation of toroidal DNA condensates under surface stress. Europhys. Lett. 33: Tzlil, S., J. T. Kindt,., A. Ben-Shaul Forces and pressures in DNA packaging and release from viral capsids. Biophys. J. 84: Arnoldi, M., M. Fritz,., A. Boulbitch Bacterial turgor pressure can be measured by atomic force microscopy. Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Topics. 6: São-José, C., M. de Frutos,., P. Tavares Pressure built by DNA packing inside virions: enough to drive DNA ejection in vitro, largely insufficient for delivery into the bacterial cytoplasm. J. Mol. Biol. 374: Deng, Y., M. Sun, and J. W. Shaevitz Direct measurement of cell wall stress stiffening and turgor pressure in live bacterial cells. Phys. Rev. Lett. 107: Zimmerman, S. B., and L. D. Murphy Macromolecular crowding and the mandatory condensation of DNA in bacteria. FEBS Lett. 390: Leforestier, A., A. Siber,., R. Podgornik Protein-DNA interactions determine the shapes of DNA toroids condensed in virus capsids. Biophys. J. 100: Köster, S., A. Evilevitch,., D. A. Weitz Influence of internal capsid pressure on viral infection by phage-l. Biophys. J. 97: R. Calendar, and S. T. Abedon, editors 006. The Bacteriophages, nd Ed. Oxford University Press, Oxford, United Kingdom. 3. Golan, R., L. I. Pietrasanta,., H. G. Hansma DNA toroids: stages in condensation. Biochemistry. 38: Rau, D. C., and V. A. Parsegian Direct measurement of the intermolecular forces between counterion-condensed DNA double helices. Evidence for long range attractive hydration forces. Biophys. J. 61: Brenner, S. L., and V. A. Parsegian A physical method for deriving the electrostatic interaction between rod-like polyions at all mutual angles. Biophys. J. 14: Odijk, T Theory of lyotropic polymer liquid crystals. Macromolecules. 19:

6 Ejecting Phage DNA against Cellular Turgor Pressure Onsager, L The effects of shape on the interaction of colloidal particles. Ann. N. Y. Acad. Sci. 51: Strey, H., V. Parsegian, and R. Podgornik Equation of state for DNA liquid crystals: fluctuation enhanced electrostatic double layer repulsion. Phys. Rev. Lett. 78: Leforestier, A Polymorphism of DNA conformation inside the bacteriophage capsid. J. Biol. Phys. 39: Castelnovo, M., R. K. Bowles,., W. M. Gelbart Osmotic force resisting chain insertion in a colloidal suspension. Eur. Phys. J. E Soft Matter. 10: Odijk, T On the statistics and dynamics of confined or entangled stiff polymers. Macromolecules. 16: Evilevitch, A., J. W. Gober,., W. M. Gelbart Measurements of DNA lengths remaining in a viral capsid after osmotically suppressed partial ejection. Biophys. J. 88: Losdorfer Bozic, A., A. Siber, and R. Podgornik. 01. How simple can a model of an empty viral capsid be? Charge distributions in viral capsids. J. Biol. Phys. 38: Khokhlov, A., and A. Semenov Liquid-crystalline ordering in the solution of long persistent chains. Phys. A. 108: Sakaue, T Semiflexible polymer confined in closed spaces. Macromolecules. 40: Cifra, P., Z. Benková, and T. Bleha Persistence lengths and structure factors of wormlike polymers under confinement. J. Phys. Chem. B. 11: Petrov, A. S., S. S. Douglas, and S. C. Harvey Effects of pulling forces, osmotic pressure, condensing agents and viscosity on the thermodynamics and kinetics of DNA ejection from bacteriophages to bacterial cells: a computational study. J. Phys. Condens. Matter. 5: