Dynamics of polymer packaging

Structure and dynamics of chemically active ring polymers: swelling to collapse

The ring structures are common in many synthetic or natural systems and experience both local and long-range forces by chemical sensing. This work is an effort to investigate the structural and dynamical properties of a chemically active ring in an explicit solvent bath utilizing hybrid molecular dynamics (MD) and multiparticle collision dynamics (MPCD) simulation techniques. We show that by tuning the chemical properties of the ring, it can be converted from a chemo-attractant to a chemo-repellent, thereby changing the steady state to be either collapsed or swelled as compared to its passive limit. We quantify these observations by comparing the scaling laws, local structures and the dynamics of active and passive rings. Furthermore, we show the impact of varying numbers of active sites by calculating the contact probability of the collapse state that highlights diverse structures. We also analyze the dynamics of the ring by finding the relaxation time and the mean square displacement of the centre of mass. A faster relaxation with enhanced diffusion is observed for the active rings.

Active particles in explicit solvent: Dynamics of clustering for alignment interaction

We study the dynamics of clustering in systems containing active particles that are immersed in an explicit solvent. For this, we have adopted a hybrid simulation method, consisting of molecular dynamics and multiparticle collision dynamics. In our model, the overlap-avoiding passive interaction of an active particle with another active particle or a solvent particle has been taken care of via variants of the Lennard-Jones potential. Dynamic interactions among the active particles have been incorporated via a Vicsek-like alignment rule in self-propulsion that facilitates clustering. We quantify the effects of activity and importance of hydrodynamics on the dynamics of clustering via variations of relevant system parameters. Results are obtained for low overall density of active particles, for which the state point is close to the vapor branch of the coexistence curve, and thus the morphology consists of disconnected clusters. In such a situation, the mechanism of growth switches among particle diffusion, diffusive coalescence, and ballistic aggregation, depending upon the presence or absence of active and hydrodynamic interactions providing different kinds of mobilities to the clusters. Corresponding growth laws have been quantified and discussed in the context of appropriate theoretical pictures. Our results suggest that multiparticle collision dynamics is an effective method for the investigation of hydrodynamic phenomena in phase-separating active matter systems.

Relaxation in a phase-separating two-dimensional active matter system with alignment interaction

Via computer simulations, we study kinetics of pattern formation in a two-dimensional active matter system. Self-propulsion in our model is incorporated via the Vicsek-like activity, i.e., particles have the tendency of aligning their velocities with the average directions of motion of their neighbors. In addition to this dynamic or active interaction, there exists passive inter-particle interaction in the model for which we have chosen the standard Lennard-Jones form. Following quenches of homogeneous configurations to a point deep inside the region of coexistence between high and low density phases, as the systems exhibit formation and evolution of particle-rich clusters, we investigate properties related to the morphology, growth, and aging. A focus of our study is on the understanding of the effects of structure on growth and aging. To quantify the latter, we use the two-time order-parameter autocorrelation function. This correlation, as well as the growth, is observed to follow power-law time dependence, qualitatively similar to the scaling behavior reported for passive systems. The values of the exponents have been estimated and discussed by comparing with the previously obtained numbers for other dimensions as well as with the new results for the passive limit of the considered model. We have also presented results on the effects of temperature on the activity mediated phase separation.

Quantitative Study of the Chiral Organization of the Phage Genome Induced by the Packaging Motor

Molecular motors that translocate DNA are ubiquitous in nature. During morphogenesis of double-stranded DNA bacteriophages, a molecular motor drives the viral genome inside a protein capsid. Several models have been proposed for the three-dimensional geometry of the packaged genome, but very little is known of the signature of the molecular packaging motor. For instance, biophysical experiments show that in some systems, DNA rotates during the packaging reaction, but most current biophysical models fail to incorporate this property. Furthermore, studies including rotation mechanisms have reached contradictory conclusions. In this study, we compare the geometrical signatures imposed by different possible mechanisms for the packaging motors: rotation, revolution, and rotation with revolution. We used a previously proposed kinetic Monte Carlo model of the motor, combined with Brownian dynamics simulations of DNA to simulate deterministic and stochastic motor models. We find that rotation is necessary for the accumulation of DNA writhe and for the chiral organization of the genome. We observe that although in the initial steps of the packaging reaction, the torsional strain of the genome is released by rotation of the molecule, in the later stages, it is released by the accumulation of writhe. We suggest that the molecular motor plays a key role in determining the final structure of the encapsidated genome in bacteriophages.

Packing of semiflexible polymers into viral capsid in crowded environments

We use coarse-grained Langevin dynamics simulations to study packing of semiflexible polymers into a spherical capsid, with and without a tail, inside a crowded cell. We use neutral and charged, but highly screened, polymers and compare packing rates of the two. Such packing conditions are relevant, for example, to λ DNA packing inside Escherichia coli bacterial cells, where the crowd particles are proteins, bacterial DNA, and salts. For a neutral polymer packing into a capsid with a tail, attractive interactions with the crowd particles make packing slightly harder at higher crowd densities, but repulsive interactions make it easier. Our results indicate that packing into a tailless capsid is less efficient at low crowd densities than into one with a long tail. However, this trend becomes opposite at higher densities. In addition, packing into a capsid with a long tail shows a highly variable waiting time before packing initiates, a feature absent for a tailless capsid. Electrical interactions at physiological conditions do not have much effect. Some bacterial cells, such as Pseudomonas chlororaphis, form a nucleuslike structure encapsulating the phage 201ϕ2-1 DNA. We also study here the packing dynamics with the nucleus present. We find packing is faster compared to the case of no-nucleus packing. We also observe knot formations but these knots untangle quickly while the polymer translocates. This knot formation is independent of polymer charge and presence of crowd particles.

Polymer translocation into cavities: Effects of confinement geometry, crowding, and bending rigidity on the free energy

Monte Carlo simulations are used to study the translocation of a polymer into a cavity. Modeling the polymer as a hard-sphere chain with a length up to N=601 monomers, we use a multiple-histogram method to measure the variation of the conformational free energy of the polymer with respect to the number of translocated monomers. The resulting free-energy functions are then used to obtain the confinement free energy for the translocated portion of the polymer. We characterize the confinement free energy for a flexible polymer in cavities with constant cross-sectional area A for various cavity shapes (cylindrical, rectangular, and triangular) as well as for tapered cavities with pyramidal and conical shape. The scaling of the free energy with cavity volume and translocated polymer subchain length is generally consistent with predictions from simple scaling arguments, with small deviations in the scaling exponents likely due to finite-size effects. The confinement free energy depends strongly on cavity shape anisometry and is a minimum for an isometric cavity shape with a length-to-width ratio of unity. Entropic depletion at the edges or vertices of the confining cavity are evident in the results for constant-A and pyramidal cavities. For translocation into infinitely long cones, the scaling of the free energy with taper angle is consistent with a theoretical prediction employing the blob model. We also examine the effects of polymer bending rigidity on the translocation free energy for cylindrical cavities. For isometric cavities, the observed scaling behavior is in partial agreement with theoretical predictions, with discrepancies arising from finite-size effects that prevent the emergence of well-defined scaling regimes. In addition, translocation into highly anisometric cylindrical cavities leads to a multistage folding process for stiff polymers. Finally, we examine the effects of crowding agents inside the cavity. We find that the confinement free energy increases with crowder density. At constant packing fraction the magnitude of this effect lessens with increasing crowder size for a crowder-to-monomer size ratio ≥1.

Effects of spherical confinement and backbone stiffness on flexible polymer jamming

We use molecular simulations to study jamming of a crumpled bead-spring model polymer in a finite container and compare to jamming of repulsive spheres. After proper constraint counting, the onset of rigidity is seen to occur isostatically as in the case of repulsive spheres. Despite this commonality, the presence of the curved container wall and polymer backbone bonds introduce new mechanical properties. Notably, these include additional bands in the vibrational density of states that reflect the material structure as well as oscillations in local contact number and density near the wall but with lower amplitude for polymers. Polymers have fewer boundary contacts, and this low-density surface layer strongly reduces the global bulk modulus. We further show that bulk-modulus dependence on backbone stiffness can be described by a model of stiffnesses in series and discuss potential experimental and biological applications.

Enhanced self-propulsion of a sphere-dimer in viscoelastic fluid

Micro-swimmers often have to encounter a medium that exhibits non-Newtonian behaviour. To understand the effect of complex environments on the propulsion dynamics of swimmers, here we have investigated a self-propelled sphere-dimer in a viscoelastic medium, using a coarse-grained hybrid mesoscopic simulation technique. We have shown that a viscoelastic fluid can result in the enhancement of swimming speed, as compared to the speed in a Newtonian fluid with the same viscosity. A non-linear response in the dimer velocity is seen for higher Péclet numbers in viscoelastic fluids. With help of various dynamical quantities, we have shown that the observed non-linear response of the directed velocity is associated with the micro-structural properties of the fluid. These include the alignment of the fluid elements and the density inhomogeneity around the moving dimer. The enhancement of self-propulsion velocity has been probed in detail, and the factors affecting the propulsion are identified.

A Molecular View of the Dynamics of dsDNA Packing Inside Viral Capsids in the Presence of Ions

Genome packing in viruses and prokaryotes relies on positively charged ions to reduce electrostatic repulsions, and induce attractions that can facilitate DNA condensation. Here we present molecular dynamics simulations spanning several microseconds of dsDNA packing inside nanometer-sized viral capsids. We use a detailed molecular model of DNA that accounts for molecular structure, basepairing, and explicit counterions. The size and shape of the capsids studied here are based on the 30-nanometer-diameter gene transfer agents of bacterium Rhodobacter capsulatus that transfer random 4.5-kbp (1.5 μm) DNA segments between bacterial cells. Multivalent cations such as spermidine and magnesium induce attraction between packaged DNA sites that can lead to DNA condensation. At high concentrations of spermidine, this condensation significantly increases the shear stresses on the packaged DNA while also reducing the pressure inside the capsid. These effects result in an increase in the packing velocity and the total amount of DNA that can be packaged inside the nanometer-sized capsids. In the simulation results presented here, high concentrations of spermidine³⁺ did not produce the premature stalling observed in experiments. However, a small increase in the heterogeneity of packing velocities was observed in the systems with magnesium and spermidine ions compared to the system with only salt. The results presented here indicate that the effect of multivalent cations and of spermidine, in particular, on the dynamics of DNA packing, increases with decreasing packing velocities.

DNA packaging in viral capsids with peptide arms

Strong chain rigidity and electrostatic self-repulsion of packed double-stranded DNA in viruses require a molecular motor to pull the DNA into the capsid. However, what is the role of electrostatic interactions between different charged components in the packaging process? Though various theories and computer simulation models were developed for the understanding of viral assembly and packaging dynamics of the genome, long-range electrostatic interactions and capsid structure have typically been neglected or oversimplified. By means of molecular dynamics simulations, we explore the effects of electrostatic interactions on the packaging dynamics of DNA based on a coarse-grained DNA and capsid model by explicitly including peptide arms (PAs), linked to the inner surface of the capsid, and counterions. Our results indicate that the electrostatic interactions between PAs, DNA, and counterions have a significant influence on the packaging dynamics. We also find that the packed DNA conformations are largely affected by the structure of the PA layer, but the packaging rate is insensitive to the layer structure.

The Semiflexible Polymer Translocation into Laterally Unbounded Region between Two Parallel Flat Membranes

Using the dynamic Monte Carlo method, we investigate dynamics of semiflexible polymer translocation through a nanopore into laterally unbounded region between two parallel flat membranes with separation R in presence of an electric field inside the pore. The average translocation time τ initially decreases rapidly with increase of R in the range of R < 10 and then almost keeps constant for R ≥ 10, and the decline range increases with increase of dimensionless bending stiffness κ. We mainly study the effect of chain length N, κ and electric field strength E on the translocation process for R = 5. The translocation dynamics is significantly altered in comparison to an unconfined environment. We find τ ~ N^α, where the exponent α increases with increase of E for small κ. α initially increases slowly with increase of E and then keeps constant for moderate κ. α decreases with increase of E for large κ. However, α decreases with increase of κ under various E. In addition, we find τ ~ κ^β. β decreases with increase of N under various E. These behaviors are interpreted in terms of the probability distribution of translocation time and the waiting time of an individual monomer segment passing through the pore during translocation.

Coarse-grained simulations of an active filament propelled by a self-generated solute gradient

A self-propelling semiflexible filament exhibits a variety of dynamical states depending on the flexibility and activity of the filament. Here we investigate the dynamics of such an active filament using a bead-spring model with the explicit hydrodynamic interactions. The activity in the filament is incorporated by inserting chemically active dimers at regular intervals along the chain. The chemical reactions at the catalytic bead of the dimer produces a self-generated concentration gradient and gives sufficient fuel to exhibit self-propulsion for the filament. Depending upon the rigidity and the configuration, the polymeric filament exhibits three distinct types of spontaneous motion, namely, rotational, snaking, and translational motion. The self-propulsion velocity of the filament for various rigidity and sizes has been calculated, and the factors affecting the propulsion are identified.

Nonequilibrium dynamics and ultraslow relaxation of confined DNA during viral packaging

Many viruses use molecular motors that generate large forces to package DNA to near-crystalline densities inside preformed viral proheads. Besides being a key step in viral assembly, this process is of interest as a model for understanding the physics of charged polymers under tight 3D confinement. A large number of theoretical studies have modeled DNA packaging, and the nature of the molecular dynamics and the forces resisting the tight confinement is a subject of wide debate. Here, we directly measure the packaging of single DNA molecules in bacteriophage phi29 with optical tweezers. Using a new technique in which we stall the motor and restart it after increasing waiting periods, we show that the DNA undergoes nonequilibrium conformational dynamics during packaging. We show that the relaxation time of the confined DNA is >10 min, which is longer than the time to package the viral genome and 60,000 times longer than that of the unconfined DNA in solution. Thus, the confined DNA molecule becomes kinetically constrained on the timescale of packaging, exhibiting glassy dynamics, which slows the motor, causes significant heterogeneity in packaging rates of individual viruses, and explains the frequent pausing observed in DNA translocation. These results support several recent hypotheses proposed based on polymer dynamics simulations and show that packaging cannot be fully understood by quasistatic thermodynamic models.

Confinement dynamics of a semiflexible chain inside nano-spheres

We study the conformations of a semiflexible chain, confined in nano-scaled spherical cavities, under two distinct processes of confinement. Radial contraction and packaging are employed as two confining procedures. The former method is performed by gradually decreasing the diameter of a spherical shell which envelopes a confined chain. The latter procedure is carried out by injecting the chain inside a spherical shell through a hole on the shell surface. The chain is modeled with a rigid body molecular dynamics simulation and its parameters are adjusted to DNA base-pair elasticity. Directional order parameter is employed to analyze and compare the confined chain and the conformations of the chain for two different sizes of the spheres are studied in both procedures. It is shown that for the confined chains in the sphere sizes of our study, they appear in spiral or tennis-ball structures, and the tennis-ball structure is more likely to be observed in more compact confinements. Our results also show that the dynamical procedure of confinement and the rate of the confinement are influential parameters of the structure of the chain inside spherical cavities.

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.

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

In this work, we report on simulations of double-stranded DNA (dsDNA) ejection from bacteriophage φ29 into a bacterial cell. The ejection was studied with a coarse-grained model, in which viral dsDNA was represented by beads on a torsion-less string. The bacteriophage's capsid and the bacterial cell were defined by sets of spherical constraints. To account for the effects of the viscous medium inside the bacterial cell, the simulations were carried out using a Langevin dynamics protocol. Our simplest simulations (involving constant viscosity and no external biasing forces) produced results compatible with the push-pull model of DNA ejection, with an ejection rate significantly higher in the first part of ejection than in the latter parts. Additionally, we performed more complicated simulations, in which we included additional factors such as external forces, osmotic pressure, condensing agents and ejection-dependent viscosity. The effects of these factors (independently and in combination) on the thermodynamics and kinetics of DNA ejection were studied. We found that, in general, the dependence of ejection forces and ejection rates on the amount of DNA ejected becomes more complex if the ejection is modeled with a broader, more realistic set of parameters and influences (such as variation in the solvent's viscosity and the application of an external force). However, certain combinations of factors and numerical parameters led to the opposition of some ejection-driving and ejection-inhibiting influences, ultimately causing an apparent simplification of the ejection profiles.

Effect of temperature and capsid tail on the packing and ejection of viral DNA

We use a simulation technique based on molecular dynamics and stochastic rotation model to present the effect of temperature and capsid tail on the packaging and ejection processes of semiflexible polymers. We consider two types of solvents, a good solvent, where the polymer is neutral and repulsion interactions among its various sections are favored, and one where the polymer is charged, giving rise to extra electrostatic reaction. For tailless capsids, we find that packing a neutral polymer is slightly slower at higher temperatures whereas its ejection is slightly slower at lower temperatures. We find the same trend for a charged polymer but the effect is noticeably larger. At a high enough temperature, we notice that packing a charged polymer can be stopped. On the other hand, at fixed temperature and regardless whether the polymer is charged, packing is much easier for a capsid with a tail whereas ejection is much slower. The effect of including the tail on the dynamics of a charged polymer, in particular, is rather significant: more packing fraction is facilitated at higher temperatures due to more ordered polymer configuration inside the capsid. In contrast, during ejection the tail traps the last remaining beads for quite some time before allowing full ejection. We interpret these results in terms of entropic and electrostatic forces.

A coarse-grain three-site-per-nucleotide model for DNA with explicit ions

The "three sites per nucleotide" (3SPN) model provides a coarse-grained representation of nucleic acids for simulation of molecular processes. Previously, this model has relied on an implicit representation of the surrounding ionic environment at the level of Debye-Hückel theory. In this work, we eliminate this limitation and present an explicit representation of ions, both monovalent and divalent. The coarse-grain ion-ion and ion-phosphate potential energy functions are inferred from all-atom simulations and parameterized to reproduce key features of the local structure and organization of ions in bulk water and in the presence of DNA. The resulting model, 3SPN.1-I, is capable of reproducing the local structure observed in detailed atomistic simulations, as well as the experimental melting temperature of DNA for a range of DNA oligonucleotide lengths, CG-content, Na(+) concentration, and Mg(2+) concentration.

Role of DNA-DNA interactions on the structure and thermodynamics of bacteriophages Lambda and P4

Electrostatic interactions play an important role in both packaging of DNA inside bacteriophages and its release into bacterial cells. While at physiological conditions DNA strands repel each other, the presence of polyvalent cations such as spermine and spermidine in solutions leads to the formation of DNA condensates. In this study, we discuss packaging of DNA into bacteriophages P4 and Lambda under repulsive and attractive conditions using a coarse-grained model of DNA and capsids. Packaging under repulsive conditions leads to the appearance of the coaxial spooling conformations; DNA occupies all available space inside the capsid. Under the attractive potential both packed systems reveal toroidal conformations, leaving the central part of the capsids empty. We also present a detailed thermodynamic analysis of packaging and show that the forces required to pack the genomes in the presence of polyamines are significantly lower than those observed under repulsive conditions. The analysis reveals that in both the repulsive and attractive regimes the entropic penalty of DNA confinement has a significant non-negligible contribution into the total energy of packaging. Additionally we report the results of simulations of DNA condensation inside partially packed Lambda. We found that at low densities DNA behaves as free unconfined polymer and condenses into the toroidal structures; at higher densities rearrangement of the genome into toroids becomes hindered, and condensation results in the formation of non-equilibrium structures. In all cases packaging in a specific conformation occurs as a result of interplay between bending stresses experienced by the confined polymer and interactions between the strands.

Polymer translocation into a fluidic channel through a nanopore

Using two-dimensional Langevin dynamics simulations, we investigate the dynamics of polymer translocation into a fluidic channel with diameter R through a nanopore under a driving force F . Due to the crowding effect induced by the partially translocated monomers, the translocation dynamics is significantly altered in comparison to an unconfined environment, namely, we observe a nonuniversal dependence of the translocation time τ on the chain length N . τ initially decreases rapidly and then saturates with increasing R , and a dependence of the scaling exponent α of τ with N on the channel width R is observed. The otherwise inverse linear scaling of τ with F breaks down and we observe a minimum of α as a function of F . These behaviors are interpreted in terms of the waiting time of an individual segment passing through the pore during translocation.

Biopolymer organization upon confinement

Biopolymers in vivo are typically subject to spatial restraints, either as a result of molecular crowding in the cellular medium or of direct spatial confinement. DNA in living organisms provides a prototypical example of a confined biopolymer. Confinement prompts a number of biophysics questions. For instance, how can the high level of packing be compatible with the necessity to access and process the genomic material? What mechanisms can be adopted in vivo to avoid the excessive geometrical and topological entanglement of dense phases of biopolymers? These and other fundamental questions have been addressed in recent years by both experimental and theoretical means. A review of the results, particularly of those obtained by numerical studies, is presented here. The review is mostly devoted to DNA packaging inside bacteriophages, which is the best studied example both experimentally and theoretically. Recent selected biophysical studies of the bacterial genome organization and of chromosome segregation in eukaryotes are also covered.

Macromolecular dynamics in crowded environments

The structural and dynamical properties of macromolecules in confining or crowded environments are different from those in simple bulk liquids. In this paper, both the conformational and diffusional dynamics of globular polymers are studied in solutions containing fixed spherical obstacles. These properties are studied as a function of obstacle volume fraction and size, as well as polymer chain length. The results are obtained using a hybrid scheme that combines multiparticle collision dynamics of the solvent with molecular dynamics that includes the interactions among the polymer monomers and between the polymer beads and obstacles and solvent molecules. The dynamics accounts for hydrodynamic interactions among the polymer beads and intermolecular forces with the solvent molecules. We consider polymers in poor solvents where the polymer chain adopts a compact globular structure in solution. Our results show that the collapse of the polymer chain to a compact globular state is strongly influenced by the obstacle array. A nonmonotonic variation in the radius of gyration with time is observed and the collapse time scale is much longer than that in simple solutions without obstacles. Hydrodynamic interactions are important at low obstacle volume fractions but are screened at high volume fractions. The diffusion of the globular polymer chain among the obstacles is subdiffusive in character on intermediate time scales where the dynamics explores the local structure of the heterogeneous environment. For large polymer chains in systems with high obstacle volume fractions, the chain adopts bloblike conformations that arise from trapping of portions of the chain in voids among the obstacles.

Viral assembly: a molecular modeling perspective

Icosahedral viruses are among the smallest and simplest of biological systems. The investigation of their structures represented the first step toward the establishment of molecular biophysics, over half a century ago. Many research groups are now pursuing investigations of viral assembly, a process that could offer new opportunities for the design of antiviral drugs and novel nanoparticles. A variety of experimental, theoretical and computational methods have been brought to bear on the study of virus structure and assembly. In this Perspective we review the contributions of theoretical and computational approaches to our understanding of the structure, energetics, thermodynamics and assembly of DNA bacteriophage and single-stranded icosahedral RNA viruses.

Characterization of DNA conformation inside bacterial viruses

In this study we develop a formalism to describe the organization of DNA inside bacteriophage capsids during genome packaging. We have previously shown that DNA inside bacteriophage phi29 (phi29) is organized into folded toroids [A. S. Petrov and S. C. Harvey, Structure 15, 21 (2007)], whereas epsilon15 (epsilon15) reveals the coaxial organization of the genetic material [A. S. Petrov, K. Lim-Hing, and S. C. Harvey, Structure 15, 807 (2007)]. We now show that each system undergoes two consecutive transitions. The first transition corresponds to the formation of global conformations and is analogous to a disorder-order conformational transition. The second transition is characterized by a significant loss of DNA mobility at the local level leading to glasslike dynamic behavior. Packing genetic material inside bacteriophages can be used as a general model to study the behavior of semiflexible chains inside confined spaces, and the proposed formalism developed here can be used to study other systems of linear polymer chains confined to closed spaces.

Knot-controlled ejection of a polymer from a virus capsid

We present a numerical study of the effect of knotting on the ejection of flexible and semiflexible polymers from a spherical, viruslike capsid. The polymer ejection rate is primarily controlled by the knot, which moves to the hole in the capsid and then acts as a ratchet. Polymers with more complex knots eject more slowly and, for large knots, the knot type, and not the flexibility of the polymer, determines the rate of ejection. We discuss the relation of our results to the ejection of DNA from viral capsids and conjecture that this process has the biological advantage of unknotting the DNA before it enters a cell.

Packings of a charged line on a sphere

We find equilibrium configurations of open and closed lines of charge on a sphere, and track them with respect to varying sphere radius. Closed lines transition from a circle to a spiral-like shape through two low-wave-number bifurcations-"baseball seam" and "twist"-which minimize Coulomb energy. The spiral shape is the unique stable equilibrium of the closed line. Other unstable equilibria arise through tip-splitting events. An open line transitions smoothly from an arc of a great circle to a spiral as the sphere radius decreases. Under repulsive potentials with faster-than-Coulomb power-law decay, the spiral is tighter in initial stages of sphere shrinkage, but at later stages of shrinkage the equilibria for all repulsive potentials converge on a spiral with uniform spacing between turns. Multiple stable equilibria of the open line are observed.

Simulations of knotting in confined circular DNA

The packing of DNA inside bacteriophages arguably yields the simplest example of genome organization in living organisms. As an assay of packing geometry, the DNA knot spectrum produced upon release of viral DNA from the P4 phage capsid has been analyzed, and compared to results of simulation of knots in confined volumes. We present new results from extensive stochastic sampling of confined self-avoiding and semiflexible circular chains with volume exclusion. The physical parameters of the chains (contour length, cross section, and bending rigidity) have been set to match those of P4 bacteriophage DNA. By using advanced sampling techniques, involving multiple Markov chain pressure-driven confinement combined with a thermodynamic reweighting technique, we establish the knot spectrum of the circular chains for increasing confinement up to the highest densities for which available algorithms can exactly classify the knots. Compactified configurations have an enclosing hull diameter approximately 2.5 times larger than the P4 caliper size. The results are discussed in relation to the recent experiments on DNA knotting inside the capsid of a P4 tailless mutant. Our investigation indicates that confinement favors chiral knots over achiral ones, as found in the experiments. However, no significant bias of torus over twist knots is found, contrary to the P4 results. The result poses a crucial question for future studies of DNA packaging in P4: is the discrepancy due to the insufficient confinement of the equilibrium simulation or does it indicate that out-of-equilibrium mechanisms (such as rotation by packaging motors) affect the genome organization, hence its knot spectrum in P4?

Shapes of semiflexible polymers in confined spaces

Biological macromolecules, living in the confines of a cell, often adopt conformations that are unlikely to occur in free space. In this paper, we investigate the effects of confinement on the shape of a semiflexible chain. Results of Monte Carlo simulations show the existence of a shape transition when the persistence length of the polymer becomes comparable to the dimensions of the box. An order parameter is introduced to quantify this behavior. A simple model is constructed to study the effect of the shape transition on the effective persistence length of the polymer.

Packaging double-helical DNA into viral capsids: structures, forces, and energetics

Small, icosahedral double-stranded DNA bacteriophage pack their genomes tightly into preformed protein capsids using an ATP-driven motor. Coarse-grain molecular-mechanics models provide a detailed picture of DNA packaging in bacteriophage, revealing how conformation depends on capsid size and shape, and the presence or absence of a protein core. The forces that oppose packaging have large contributions from both electrostatic repulsions and the entropic penalty of confining the DNA into the capsid, whereas elastic deformations make only a modest contribution. The elastic deformation energy is very sensitive to the final conformation, whereas the electrostatic and entropic penalties are not, so the packaged DNA favors conformations that minimize the bending energy.

Ejection dynamics of polymeric chains from viral capsids: effect of solvent quality

We present simulations investigating the effects of solvent quality on the dynamics of flexible (RNA-like) and semiflexible (DNA-like) polymers ejecting from spherical viral capsids. We find that the mean ejection time increases and the ejection time distributions are broadened as the solvent quality decreases. Our results thus suggest that DNA ejection may be very efficiently controlled by tuning the salt concentration in the environment, in agreement with recent experimental findings. We also observe random pauses in the ejection. These become extremely long for semiflexible polymers at lower solvent quality, and we interpret this as a signature of a low driving force for ejection. We find that, for most polymers, ejection is an all-or-nothing process at the solvent conditions we investigated: polymers normally completely eject once the process is initiated.

Portal motor velocity and internal force resisting viral DNA packaging in bacteriophage phi29

During the assembly of many viruses, a powerful molecular motor compacts the genome into a preassembled capsid. Here, we present measurements of viral DNA packaging in bacteriophage phi29 using an improved optical tweezers method that allows DNA translocation to be measured from initiation to completion. This method allowed us to study the previously uncharacterized early stages of packaging and facilitated more accurate measurement of the length of DNA packaged. We measured the motor velocity versus load at near-zero filling and developed a ramped DNA stretching technique that allowed us to measure the velocity versus capsid filling at near-zero load. These measurements reveal that the motor can generate significantly higher velocities and forces than detected previously. Toward the end of packaging, the internal force resisting DNA confinement rises steeply, consistent with the trend predicted by many theoretical models. However, the force rises to a higher magnitude, particularly during the early stages of packaging, than predicted by models that assume coaxial inverse spooling of the DNA. This finding suggests that the DNA is not arranged in that conformation during the early stages of packaging and indicates that internal force is available to drive complete genome ejection in vitro. The maximum force exceeds 100 pN, which is about one-half that predicted to rupture the capsid shell.

Ionic effects on viral DNA packaging and portal motor function in bacteriophage phi 29

In many viruses, DNA is confined at such high density that its bending rigidity and electrostatic self-repulsion present a strong energy barrier in viral assembly. Therefore, a powerful molecular motor is needed to package the DNA into the viral capsid. Here, we investigate the role of electrostatic repulsion on single DNA packaging dynamics in bacteriophage phi 29 via optical tweezers measurements. We show that ionic screening strongly affects the packing forces, confirming the importance of electrostatic repulsion. Separately, we find that ions affect the motor function. We separate these effects through constant force measurements and velocity versus load measurements at both low and high capsid filling. Regarding motor function, we find that eliminating free Mg(2+) blocks initiation of packaging. In contrast, Na(+) is not required, but it increases the motor velocity by up to 50% at low load. Regarding internal resistance, we find that the internal force was lowest when Mg(2+) was the dominant ion or with the addition of 1 mM Co(3+). Forces resisting DNA confinement were up to approximately 80% higher with Na(+) as the dominant counterion, and only approximately 90% of the genome length could be packaged in this condition. The observed trend of the packing forces is in accord with that predicted by DNA charge-screening theory. However, the forces are up to six times higher than predicted by models that assume coaxial spooling of the DNA and interaction potentials derived from DNA condensation experiments. The forces are also severalfold higher than ejection forces measured with bacteriophage lambda.

Stability diagram for dense suspensions of model colloidal Al2O3 particles in shear flow

In Al2O3 suspensions, depending on the experimental conditions, very different microstructures can be found, comprising fluidlike suspensions, a repulsive structure, and a clustered microstructure. For technical processing in ceramics, the knowledge of the microstructure is of importance, since it essentially determines the stability of a workpiece to be produced. To enlighten this topic, we investigate these suspensions under shear by means of simulations. We observe cluster formation on two different length scales: the distance of nearest neighbors and on the length scale of the system size. We find that the clustering behavior does not depend on the length scale of observation. If interparticle interactions are not attractive the particles form layers in the shear flow. The results are summarized in a stability diagram.

Biological consequences of tightly bent DNA: the other life of a macromolecular celebrity

The mechanical properties of DNA play a critical role in many biological functions. For example, DNA packing in viruses involves confining the viral genome in a volume (the viral capsid) with dimensions that are comparable to the DNA persistence length. Similarly, eukaryotic DNA is packed in DNA-protein complexes (nucleosomes), in which DNA is tightly bent around protein spools. DNA is also tightly bent by many proteins that regulate transcription, resulting in a variation in gene expression that is amenable to quantitative analysis. In these cases, DNA loops are formed with lengths that are comparable to or smaller than the DNA persistence length. The aim of this review is to describe the physical forces associated with tightly bent DNA in all of these settings and to explore the biological consequences of such bending, as increasingly accessible by single-molecule techniques.

Shear viscosity of claylike colloids in computer simulations and experiments

Dense suspensions of small strongly interacting particles are complex systems that are rarely understood on the microscopic level. We investigate properties of dense suspensions and sediments of small spherical Al2O3 particles in a shear cell by means of a combined molecular-dynamics and stochastic rotation dynamics simulation. We study structuring effects and the dependence of the suspension's viscosity on the shear rate and shear thinning for systems of varying salt concentration and pH value. To show the agreement of our results with experimental data, the relation between the bulk pH value and surface charge of spherical colloidal particles is modeled by Debye-Hückel theory in conjunction with a 2 pK charge regulation model.

Mesoscopic description of solvent effects on polymer dynamics

Solvent effects on polymer dynamics and structure are investigated using a mesoscopic solvent model that accounts for hydrodynamic interactions among the polymer beads. The simulation method combines molecular dynamics of the polymer chain, interacting with the solvent molecules through intermolecular forces, with mesoscopic multiparticle collision dynamics for the solvent molecules. Changes in the intermolecular forces between the polymer beads and mesoscopic solvent molecules are used to vary the solvent conditions from those for good to poor solvents. Polymer collapse and expansion dynamics following changes in solvent conditions are studied for homopolymer and block copolymer solutions. The frictional properties of polymers are also investigated.

Polymer packaging and ejection in viral capsids: shape matters

We use a mesoscale simulation approach to explore the impact of different capsid geometries on the packaging and ejection dynamics of polymers of different flexibility. We find that both packing and ejection times are faster for flexible polymers. For such polymers a sphere packs more quickly and ejects more slowly than an ellipsoid. For semiflexible polymers, however, the case relevant to DNA, a sphere both packs and ejects more easily. We interpret our results by considering both the thermodynamics and the relaxational dynamics of the polymers. The predictions could be tested with biomimetic experiments with synthetic polymers inside artificial vesicles. Our results suggest that phages may have evolved to be roughly spherical in shape to optimize the speed of genome ejection, which is the first stage in infection.

Polymer translocation: The effect of backflow

We investigate the effect of backflow on the translocation dynamics of short, flexible polymer chains threading through a small hole in a wall. We find that hydrodynamic interactions between polymer beads play an important role in determining the translocation time distribution: as a monomer moves through the hole it sets up a flow field which transfers momentum to neighboring monomers, thus helping them to move in the same direction. Translocation times are calculated by using the velocity-Verlet algorithm to solve the equations of motion of a polymer which moves in a fluid described by the stochastic rotation algorithm, a particle-based Navier-Stokes solver.