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

Viral Genome Delivery Across Bacterial Cell Surfaces

In 1952, Hershey and Chase used bacteriophage T2 genome delivery inside Escherichia coli to demonstrate that DNA, not protein, is the genetic material. Over 70 years later, our understanding of bacteriophage structure has grown dramatically, mainly thanks to the cryogenic electron microscopy revolution. In stark contrast, phage genome delivery in prokaryotes remains poorly understood, mainly due to the inherent challenge of studying such a transient and complex process. Here, we review the current literature on viral genome delivery across bacterial cell surfaces. We focus on icosahedral bacterial viruses that we arbitrarily sort into three groups based on the presence and size of a tail apparatus. We inventory the building blocks implicated in genome delivery and critically analyze putative mechanisms of genome ejection. Bacteriophage genome delivery into bacteria is a topic of growing interest, given the renaissance of phage therapy in Western medicine as a therapeutic alternative to face the antibiotic resistance crisis.

Structure and the role of filling rate on model dsDNA packed in a phage capsid

The conformation of DNA inside bacteriophages is of paramount importance for understanding packaging and ejection mechanisms. Models describing the structure of the confined macromolecule have depicted highly ordered conformations, such as spooled or toroidal arrangements that focus on reproducing experimental results obtained by averaging over thousands of configurations. However, it has been seen that more disordered states, including DNA kinking and the presence of domains with different DNA orientation can also accurately reproduce many of the structural experiments. In this work we have compared the results obtained through different simulated filling rates. We find a rate dependence for the resulting constrained states showing different anisotropic configurations. We present a quantitative analysis of the density distribution and the DNA orientation across the capsid showing excellent agreement with structural experiments. Second, we have analyzed the correlations within the capsid, finding evidence of the presence of domains characterized by aligned segments of DNA characterized by the structure factor. Finally, we have measured the number and distribution of DNA defects such as the emergence of bubbles and kinks as function of the filling rate. We find the slower the rate the fewer kink defects that appear and they would be unlikely at experimental filling rates with our model parameters. DNA domains of various orientation get larger with slower rates.

Reduction of Salmonella in ground chicken using a bacteriophage

This study's goal was to ascertain the effectiveness of a commercially available Salmonella bacteriophage during ground chicken production focusing on: water source, different Salmonella serovars, and time. Salmonella-free boneless, skinless chicken meat was inoculated with 4.0 Log CFU/cm2 of either a cocktail of 3 Salmonella isolates derived from ground chicken (GC) or a cocktail of 3 Salmonella strains not isolated from ground chicken (non-GC). Bacteriophages were spread onto the chicken using sterile tap or filtered water for 30 min or 8 h. Salmonella was recovered using standard plating method. Greater Salmonella reduction was observed when the bacteriophage was diluted in sterile tap water than in sterile filtered water: 0.39 Log CFU/cm2 and 0.23 Log CFU/cm2 reduction after 30 min, respectively (P < 0.05). The non-GC isolates showed reductions of 0.71 Log CFU/cm2 and 0.90 Log CFU/cm2 after 30 min and 8 h, respectively (P < 0.05). The GC isolates were less sensitive to the bacteriophage: 0.39 Log CFU/cm2 and 0.67 Log CFU/cm2 reductions after 30 min and 8 h, respectively (P < 0.05). In conclusion, bacteriophage reduction was dependent on water used to dilute the bacteriophage, Salmonella's susceptibility to the bacteriophage, and treatment time.

Rigidity-induced scale invariance in polymer ejection from capsid

While the dynamics of a fully flexible polymer ejecting a capsid through a nanopore has been extensively studied, the ejection dynamics of semiflexible polymers has not been properly characterized. Here we report results from simulations of ejection dynamics of semiflexible polymers ejecting from spherical capsids. Ejections start from strongly confined polymer conformations of constant initial monomer density. We find that, unlike for fully flexible polymers, for semiflexible polymers the force measured at the pore does not show a direct relation to the instantaneous ejection velocity. The cumulative waiting time t(s), that is, the time at which a monomer s exits the capsid the last time, shows a clear change when increasing the polymer rigidity κ. The major part of an ejecting polymer is driven out of the capsid by internal pressure. At the final stage the polymer escapes the capsid by diffusion. For the driven part there is a crossover from essentially exponential growth of t with s of the fully flexible polymers to a scale-invariant form. In addition, a clear dependence of t on polymer length N_{0} was found. These findings combined give the dependence t(s)∝N_{0}^{0.55}s^{1.33} for the strongly rigid polymers. This crossover in dynamics where κ acts as a control parameter is reminiscent of a phase transition. This analogy is further enhanced by our finding a perfect data collapse of t for polymers of different N_{0} and any constant κ.

Phage-like packing structures with mean field sequence dependence

Packing of double-stranded DNA in phages must overcome both electrostatic repulsions and the problem of persistence length. We consider coarse-grained models with the ability to kink and with randomly generated disorder. We show that the introduction of kinking into configurations of the DNA polymer packaged within spherical confinement results in significant reductions of the overall energies and pressures. We use a kink model which has the ability to deform every 24 bp, close to the average length predicted from phage sequence. The introduction of such persistence length defects even with highly random packing models increases the local nematic ordering of the packed DNA polymer segments. Such local ordering allowed by kinking not only reduces the total bending energy of confined DNA due to nonlinear elasticity but also reduces the electrostatic component of the energy and pressure. We show that a broad ensemble of polymer configurations is consistent with the structural data. © 2016 Wiley Periodicals, Inc.

Computational virology: From the inside out

Viruses typically pack their genetic material within a protein capsid. Enveloped viruses also have an outer membrane made up of a lipid bilayer and membrane-spanning glycoproteins. X-ray diffraction and cryoelectron microscopy provide high resolution static views of viral structure. Molecular dynamics (MD) simulations may be used to provide dynamic insights into the structures of viruses and their components. There have been a number of simulations of viral capsids and (in some cases) of the inner core of RNA or DNA packaged within them. These simulations have generally focussed on the structural integrity and stability of the capsid and/or on the influence of the nucleic acid core on capsid stability. More recently there have been a number of simulation studies of enveloped viruses, including HIV-1, influenza A, and dengue virus. These have addressed the dynamic behaviour of the capsid, the matrix, and/or of the outer envelope. Analysis of the dynamics of the lipid bilayer components of the envelopes of influenza A and of dengue virus reveals a degree of biophysical robustness, which may contribute to the stability of virus particles in different environments. Significant computational challenges need to be addressed to aid simulation of complex viruses and their membranes, including the need to integrate structural data from a range of sources to enable us to move towards simulations of intact virions. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Ejecting phage DNA against cellular turgor pressure

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.

Calculation of a fluctuating entropic force by phase space sampling

A polymer chain pinned in space exerts a fluctuating force on the pin point in thermal equilibrium. The average of such fluctuating force is well understood from statistical mechanics as an entropic force, but little is known about the underlying force distribution. Here, we introduce two phase space sampling methods that can produce the equilibrium distribution of instantaneous forces exerted by a terminally pinned polymer. In these methods, both the positions and momenta of mass points representing a freely jointed chain are perturbed in accordance with the spatial constraints and the Boltzmann distribution of total energy. The constraint force for each conformation and momentum is calculated using Lagrangian dynamics. Using terminally pinned chains in space and on a surface, we show that the force distribution is highly asymmetric with both tensile and compressive forces. Most importantly, the mean of the distribution, which is equal to the entropic force, is not the most probable force even for long chains. Our work provides insights into the mechanistic origin of entropic forces, and an efficient computational tool for unbiased sampling of the phase space of a constrained system.

Length quantization of DNA partially expelled from heads of a bacteriophage T3 mutant

DNA packaging of phages phi29, T3 and T7 sometimes produces incompletely packaged DNA with quantized lengths, based on gel electrophoretic band formation. We discover here a packaging ATPase-free, in vitro model for packaged DNA length quantization. We use directed evolution to isolate a five-site T3 point mutant that hyper-produces tail-free capsids with mature DNA (heads). Three tail gene mutations, but no head gene mutations, are present. A variable-length DNA segment leaks from some mutant heads, based on DNase I-protection assay and electron microscopy. The protected DNA segment has quantized lengths, based on restriction endonuclease analysis: six sharp bands of DNA missing 3.7-12.3% of the last end packaged. Native gel electrophoresis confirms quantized DNA expulsion and, after removal of external DNA, provides evidence that capsid radius is the quantization-ruler. Capsid-based DNA length quantization possibly evolved via selection for stalling that provides time for feedback control during DNA packaging and injection.

RELATIONSHIP BETWEEN THE GENOME PACKING IN THE BACTERIOPHAGE CAPSID AND THE KINETICS OF DNA EJECTION

We present a general survey of experimental and theoretical observations of DNA structure and in vitro ejection kinetics for different bacteriophage species. In some species, like T5, the ejection may present pauses and arrests that have not been detected in others species like Lambda. We propose hypotheses to explain such differences and we discuss how the experimental conditions may be important for their detection. Our work highlights the role of DNA organization inside the bacteriophage capsid on the stochastic and out of equilibrium nature of the ejection process.