Biopolymer organization upon confinement

Self-Organized States from Solutions of Active Ring Polymers in Bulk and under Confinement

In the present work, we study, by means of numerical simulations, the structural and dynamical behavior of a suspension of active ring polymers in bulk and under lateral confinement. At high activity, when changing the distance between the confining planes and the polymers' density, we identify the emergence of a self-organized dynamical state, characterized by the coexistence of slowly diffusing clusters of rotating disks and faster rings moving in between them. We further assess that self-organization is robust in a range of polymer sizes, and we identify a critical value of the activity, necessary to trigger cluster formation. This system has distinctive features resembling at the same time polymers, liquid crystals, and active systems, where the interplay between activity, topology, and confinement leads to a spontaneous segregation in an initially one-component solution.

Mixing and demixing arising from compression of two semiflexible polymer chains in nanochannels

We use molecular dynamics simulation to probe the non-equilibrium physics of two nanochannel-confined semiflexible polymers in a homogeneous flow field. We find that for sufficiently stiff chains the internal organization of the two chains takes the form of interwoven folds and circular coils. This organization can lead to mixing or demixing depending on chain stiffness and flow speed. At low and intermediate flow, the two chains adopt a folded configuration, which favours mixing. At high flow, the two chains adopt a predominantly coiled configuration. The coiled configuration results in demixing when the chains are compressed from an initially demixed condition and mixing when the chains are compressed from an initially mixed condition. We find that the mixing/demixing behaviour is governed by the ratio of the number of folded segments of one chain relative to the other at low flow and by the degree of coiling in both chains at high flow. For decreasing stiffness, the chains start to aggregate locally instead of mixing smoothly at low and intermediate flow. In the limit of completely flexible chains, the two chains either completely segregate at low flow, or adopt a locally demixed configuration consisting of large aggregates of one chain relative to the other that undergo complex stochastic dynamics, diffusing, disintegrating, and reforming at intermediate flow. The transition from complete segregation to the aggregate-dominated configuration occurs when the linear intra-chain ordering breaks down.

Adherent Moving of Polymers in Spherical Confined Binary Semiflexible Ring Polymer Mixtures

Based on the coarse-grained model, we used molecular dynamics methods to calculate and simulate a semiflexible long ring–semiflexible short ring blended polymer system confined in a hard sphere. We systematically studied the distribution and motion characteristics of the long ring chain. The results show that when the short ring is short enough (Lshort < 20), the long ring (Llong = 50) is separated from the blend system and then distributed against the inner wall. As the length of the short ring increases (Lshort ≥ 20), the long ring can no longer be separated from the blending system. Moreover, we found that the long ring demonstrates a random direction of adherent walking behavior on the inner surface of the hard sphere. The velocity of the long ring decreases with the increase in the short ring length Lshort. Specifically for Lshort ≥ 20, the system does not undergo phase separation and the speed of the long ring decreases sharply along with the long ring distributed inside the confined bulk. This is related to the inner wall layer moving faster than the inside bulk of the restricted system. Our simulation results can help us to understand the distribution of macromolecules in biological systems in confined systems, including the restricted chromosome partitioning distribution and packing structure of circular DNA molecules.

Knot Factories with Helical Geometry Enhance Knotting and Induce Handedness to Knots

We performed molecular dynamics simulations of DNA polymer chains confined in helical nano-channels under compression in order to explore the potential of knot-factories with helical geometry to produce knots with a preferred handedness. In our simulations, we explore mutual effect of the confinement strength and compressive forces in a range covering weak, intermediate and strong confinement together with weak and strong compressive forces. The results find that while the common metrics of polymer chain in cylindrical and helical channels are very similar, the DNA in helical channels exhibits greatly different topology in terms of chain knottedness, writhe and handedness of knots. The results show that knots with a preferred chirality in terms of average writhe can be produced by using channels with a chosen handedness.

Geometric Predictors of Knotted and Linked Arcs

Inspired by how certain proteins “sense” knots and entanglements in DNA molecules, here, we ask if local geometric features that may be used as a readout of the underlying topology of generic polymers exist. We perform molecular simulations of knotted and linked semiflexible polymers and study four geometric measures to predict topological entanglements: local curvature, local density, local 1D writhe, and nonlocal 3D writhe. We discover that local curvature is a poor predictor of entanglements. In contrast, segments with maximum local density or writhe correlate as much as 90% of the time with the shortest knotted and linked arcs. We find that this accuracy is preserved across different knot types and also under significant spherical confinement, which is known to delocalize essential crossings in knotted polymers. We further discover that nonlocal 3D writhe is the best geometric readout of the knot location. Finally, we discuss how these geometric features may be used to computationally analyze entanglements in generic polymer melts and gels.

Organized states arising from compression of single semiflexible polymer chains in nanochannels

We use molecular dynamics simulation to probe the nonequilibrium physics of single nanochannel-confined semiflexible polymers in a homogeneous flow field. The flow field compresses the polymer against the end of the nanochannel, simulating an experiment of a nanochannel confined chain compressed against a slit barrier. The flow-based compression gives rise to a packing of the chain against the channel end that possesses a striking organization, consisting of interweaving of folds and circular coils. For stiff chains at low flow, we find that the organization is dominated by repeated hairpin folds. For stiff chains at higher flow, we observe that circular coils arise along with the folds, with folding and coiling domains becoming interwoven at the highest flow speeds. Chain organization is retained even when the chain persistence length is on order of the channel width. We show that the global polymer organization, consisting of a number of defined folds and coiled loops, arises from the minimization of the total chain free energy.

A statistical approach to knot confinement via persistent homology

In this paper, we study how randomly generated knots occupy a volume of space using topological methods. To this end, we consider the evolution of the first homology of an immersed metric neighbourhood of a knot's embedding for growing radii. Specifically, we extract features from the persistent homology (PH) of the Vietoris-Rips complexes built from point clouds associated with knots. Statistical analysis of our data shows the existence of increasing correlations between geometric quantities associated with the embedding and PH-based features, as a function of the knots' lengths. We further study the variation of these correlations for different knot types. Finally, this framework also allows us to define a simple notion of deviation from ideal configurations of knots.

Knot formation of dsDNA pushed inside a nanochannel

Recent experiments demonstrated that knots in single molecule dsDNA can be formed by compression in a nanochannel. In this manuscript, we further elucidate the underlying molecular mechanisms by carrying out a compression experiment in silico, where an equilibrated coarse-grained double-stranded DNA confined in a square channel is pushed by a piston. The probability of forming knots is a non-monotonic function of the persistence length and can be enhanced significantly by increasing the piston speed. Under compression knots are abundant and delocalized due to a backfolding mechanism from which chain-spanning loops emerge, while knots are less frequent and only weakly localized in equilibrium. Our in silico study thus provides insights into the formation, origin and control of DNA knots in nanopores.

Equilibrium Conformation of Catenated DNA Networks in Slitlike Confinement

A kinetoplast is a planar network of catenated DNA rings with topology that resembles that of chain mail armor. In this work, we use single-molecule experiments to probe the conformation of kinetoplasts confined to slits. We find that the in-plane size of kinetoplasts increases with degree of confinement, akin to the slitlike confinement of linear DNA. The change in kinetoplast size with channel height is consistent with the scaling prediction from a Flory-type approach for a 2D polymer. With an increase in extent of confinement, the kinetoplasts appear to unfold and take on more uniform circular shapes, in contrast to the broad range of conformations observed for kinetoplasts in bulk.

Mesoscale Simulation of Bacterial Chromosome and Cytoplasmic Nanoparticles in Confinement

In this study we investigated, using a simple polymer model of bacterial chromosome, the subdiffusive behaviors of both cytoplasmic particles and various loci in different cell wall confinements. Non-Gaussian subdiffusion of cytoplasmic particles as well as loci were obtained in our Langevin dynamic simulations, which agrees with fluorescence microscope observations. The effects of cytoplasmic particle size, locus position, confinement geometry, and density on motions of particles and loci were examined systematically. It is demonstrated that the cytoplasmic subdiffusion can largely be attributed to the mechanical properties of bacterial chromosomes rather than the viscoelasticity of cytoplasm. Due to the randomly positioned bacterial chromosome segments, the surrounding environment for both particle and loci is heterogeneous. Therefore, the exponent characterizing the subdiffusion of cytoplasmic particle/loci as well as Laplace displacement distributions of particle/loci can be reproduced by this simple model. Nevertheless, this bacterial chromosome model cannot explain the different responses of cytoplasmic particles and loci to external compression exerted on the bacterial cell wall, which suggests that the nonequilibrium activity, e.g., metabolic reactions, play an important role in cytoplasmic subdiffusion.

Packing of stiff rods on ellipsoids: Geometry

We suggest a geometrical mechanism for the ordering of slender filaments inside nonisotropic containers, using cortical microtubules in plant cells and the packing of viral genetic material inside capsids as concrete examples. We show analytically how the shape of the cell affects the ordering of phantom elastic rods that are not self-avoiding (i.e., self-crossing is allowed). We find that for oblate cells, the preferred orientation is along the equator, while for prolate spheroids with an aspect ratio close to 1, the orientation is along the principal (long axis). Surprisingly, at a high enough aspect ratio, a configurational phase transition occurs and the rods no longer point along the principal axis, but at an angle to it, due to high curvature at the poles. We discuss some of the possible effects of self-avoidance using energy considerations. These results are relevant to other packing problems as well, such as the spooling of filament in the industry or spider silk inside water droplets.

Effects of Topological Constraints on Penetration Structures of Semi-Flexible Ring Polymers

The effects of topological constraints on penetration structures of semi-flexible ring polymers in a melt are investigated using molecular dynamics simulations, considering simultaneously the effects of the chain stiffness. Three topology types of rings are considered: 01-knot (the unknotted), 31-knot and 61-knot ring polymers, respectively. With the improved algorithm to detect and quantify the inter-ring penetration (or inter-ring threading), the degree of ring threading does not increase monotonously with the chain stiffness, existing a peak value at the intermediate stiffness. It indicates that rings interpenetrate most at intermediate stiffness where there is a balance between coil expansion (favoring penetrations) and stiffness (inhibiting penetrations). Meanwhile, the inter-ring penetration would be suppressed with the knot complexity of the rings. The analysis of effective potential between the rings provides a better understanding for this non-monotonous behavior in inter-ring penetration.

Polymer in wedge-shaped confinement: Effect on the θ temperature

The equilibrium properties of a finite-length linear polymer chain confined in an infinite wedge composed of two perfectly reflecting hard walls meeting at a variable apex angle (α) are presented. One end of the polymer is anchored a distance y from the apex on the conical axis of symmetry, while the other end is free. We report here, the nonmonotonic behavior of θ temperature as a function of y for a finite-length chain. Data collapse for different chain lengths indicates that such behavior will exist for all finite lengths. We delineate the origin of such nonmonotonic behavior, which may have potential applications in understanding the cellular process occurring in nanoconfined geometries.

Packing structure of semiflexible rings

Unraveling the packing structure of dense assemblies of semiflexible rings is not only fundamental for the dynamical description of polymer rings, but also key to understand biopackaging, such as observed in circular DNA of viruses or genome folding. Here we use X-ray tomography to study the geometrical and topological features of disordered packings of rubber bands in a cylindrical container. Assemblies of short bands assume a liquid-like disordered structure, with short-range orientational order, and reveal only minor influence of the container. In the case of longer bands, the confinement causes folded configurations and the bands interpenetrate and entangle. Most of the systems are found to display a threading network which percolates the system. Surprisingly, for long bands whose diameter is more than twice the diameter of the container, we found that all bands interpenetrate each other, in a complex fully entangled structure.

Transcriptional supercoiling boosts topoisomerase II-mediated knotting of intracellular DNA

Recent studies have revealed that the DNA cross-inversion mechanism of topoisomerase II (topo II) not only removes DNA supercoils and DNA replication intertwines, but also produces small amounts of DNA knots within the clusters of nucleosomes that conform to eukaryotic chromatin. Here, we examine how transcriptional supercoiling of intracellular DNA affects the occurrence of these knots. We show that although (−) supercoiling does not change the basal DNA knotting probability, (+) supercoiling of DNA generated in front of the transcribing complexes increases DNA knot formation over 25-fold. The increase of topo II-mediated DNA knotting occurs both upon accumulation of (+) supercoiling in topoisomerase-deficient cells and during normal transcriptional supercoiling of DNA in TOP1 TOP2 cells. We also show that the high knotting probability (P^kn ≥ 0.5) of (+) supercoiled DNA reflects a 5-fold volume compaction of the nucleosomal fibers in vivo. Our findings indicate that topo II-mediated DNA knotting could be inherent to transcriptional supercoiling of DNA and other chromatin condensation processes and establish, therefore, a new crucial role of topoisomerase II in resetting the knotting–unknotting homeostasis of DNA during chromatin dynamics.

Self-assembly of core-corona particles confined in a circular box

Using Monte Carlo simulations, we study the assembly of colloidal particles interacting via isotropic core-corona potentials in two dimensions and confined in a circular box. We explore the structural variety at low temperatures as function of the number of particles (N) and the size of the confining box and find a rich variety of patterns that are not observed in unconfined flat space. For a small number of particles [Formula: see text], we identify the zero-temperature minimal energy configurations at a given box size. When the number of particles is large ([Formula: see text]), we distinguish different regimes that appear in route towards close packing configurations as the box size decreases. These regimes are characterized by the increase in the number of branching points and their coordination number. Interestingly, we obtain anisotropic open structures with unexpected variety of rotational symmetries that can be controlled by changing the model parameters, and some of the structures have chirality, in spite of the isotropy of the interactions and of the confining box. For arbitrary temperatures, we employ Monte Carlo integration to obtain the average energy and the configurational entropy of the system, which are then used to construct a phase diagram as function of temperature and box radius. Our findings show that confined core-corona particles can be a suitable system to engineer particles with highly complex internal structure that may serve as building blocks in hierarchical assembly.

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.

Shape and Size Control of Artificial Cells for Bottom-Up Biology

Bottom-up biology is an expanding research field that aims to understand the mechanisms underlying biological processes via in vitro assembly of their essential components in synthetic cells. As encapsulation and controlled manipulation of these elements is a crucial step in the recreation of such cell-like objects, microfluidics is increasingly used for the production of minimal artificial containers such as single-emulsion droplets, double-emulsion droplets, and liposomes. Despite the importance of cell morphology on cellular dynamics, current synthetic-cell studies mainly use spherical containers, and methods to actively shape manipulate these have been lacking. In this paper, we describe a microfluidic platform to deform the shape of artificial cells into a variety of shapes (rods and discs) with adjustable cell-like dimensions below 5 μm, thereby mimicking realistic cell morphologies. To illustrate the potential of our method, we reconstitute three biologically relevant protein systems (FtsZ, microtubules, collagen) inside rod-shaped containers and study the arrangement of the protein networks inside these synthetic containers with physiologically relevant morphologies resembling those found in living cells.

Topologically Linked Chains in Confinement

We examine how channel confinement affects the equilibrium properties of topologically linked ring polymers and, by contrast, of equivalent unlinked rings, too. By performing extensive simulations of semiflexible rings of different chain length, N, and channel diameter, D, we discover three notable properties purely due to linking. First, upon entering the weak confinement regime, the length of the physically linked portion, lLKThe, becomes independent of chain length. Next, even when confinement is strong enough to pull apart and segregate unlinked rings, lLK stays much larger than in the highly stretched limit. Finally, at fixed N, lLK varies approximately as D^0.5, and we provide a simple scaling argument for this power-law behavior. These properties, which may hold for different link topologies, can be tested by current experimental setups on DNA rings confined in microchannels. Moreover, they could be relevant for the efficient in vivo unlinking of newly replicated bacterial chromosomes.

Entropic elasticity and dynamics of the bacterial chromosome: A simulation study

We study the compression and extension dynamics of a DNA-like polymer interacting with non-DNA binding and DNA-binding proteins, by means of computer simulations. The geometry we consider is inspired by recent experiments probing the compressional elasticity of the bacterial nucleoid (DNA plus associated proteins), where DNA is confined into a cylindrical container and subjected to the action of a "piston"-a spherical bead to which an external force is applied. We quantify the effect of steric interactions (excluded volume) on the force-extension curves as the polymer is compressed. We find that non-DNA-binding proteins, even at low densities, exert an osmotic force which can be a lot larger than the entropic force exerted by the compressed DNA. The trends we observe are qualitatively robust with respect to changes in protein sizes and are similar for neutral and charged proteins (and DNA). We also quantify the dynamics of DNA expansion following removal of the "piston": while the expansion is well fitted by power laws, the apparent exponent depends on protein concentration and protein-DNA interaction in a significant way. We further highlight an interesting kinetic process which we observe during the expansion of DNA interacting with DNA-binding proteins when the interaction strength is intermediate: the proteins bind while the DNA is packaged by the compression force, but they "pop-off" one-by-one as the force is removed, leading to a slow unzipping kinetics. Finally, we quantify the importance of supercoiling, which is an important feature of bacterial DNA in vivo.

Relevance of packing to colloidal self-assembly

Since the 1920s, packing arguments have been used to rationalize crystal structures in systems ranging from atomic mixtures to colloidal crystals. Packing arguments have recently been applied to complex nanoparticle structures, where they often, but not always, work. We examine when, if ever, packing is a causal mechanism in hard particle approximations of colloidal crystals. We investigate three crystal structures composed of their ideal packing shapes. We show that, contrary to expectations, the ordering mechanism cannot be packing, even when the thermodynamically self-assembled structure is the same as that of the densest packing. We also show that the best particle shapes for hard particle colloidal crystals at any finite pressure are imperfect versions of the ideal packing shape.

Conformations and orientational ordering of semiflexible polymers in spherical confinement

Semiflexible polymers in lyotropic solution confined inside spherical nanoscopic "containers" with repulsive walls are studied by molecular dynamics simulations and density functional theory, as a first step to model confinement effects on stiff polymers inside of miniemulsions, vesicles, and cells. It is shown that the depletion effects caused by the monomer-wall repulsion depend distinctly on the radius R of the sphere. Further, nontrivial orientational effects occur when R, the persistence length ℓp, and the contour length L of the polymers are of similar magnitude. At intermediate densities, a "shell" of wall-attached chains is forming, such that the monomers belonging to those chains are in a layer at about the distance of one monomer from the container wall. At the same time, the density of the centers of mass of these chains is peaked somewhat further inside, but still near the wall. However, the arrangement of chains is such that the total monomer density is almost uniform in the sphere, apart from a small layering peak at the wall. It is shown that excluded volume effects among the monomers are crucial to account for this behavior, although they are negligible for comparable isolated single semiflexible chains of the same length.

Simulations of knotting of DNA during genome mapping

Genome mapping involves the confinement of long DNA molecules, in excess of 150 kilobase pairs, in nanochannels near the circa 50 nm persistence length of DNA. The fidelity of the map relies on the assumption that the DNA is linearized by channel confinement, which assumes the absence of knots. We have computed the probability of forming different knot types and the size of these knots for long chains (approximately 164 kilobase pairs) via pruned-enriched Rosenbluth method simulations of a discrete wormlike chain model of DNA in channel sizes ranging from 35 nm to 60 nm. Compared to prior simulations of short DNA in similar confinement, these long molecules exhibit both complex knots, with up to seven crossings, and multiple knots per chain. The knotting probability is a very strong function of channel size, ranging from 0.3% to 60%, and rationalized in the context of Odijk's theory for confined semiflexible chains. Overall, the knotting probability and knot size obtained from these equilibrium measurements are not consistent with experimental measurements of the properties of anomalously bright regions along the DNA backbone during genome mapping experiments. This result suggests that these events in experiments are either knots formed during the processing of the DNA prior to injection into the nanochannel or regions of locally high DNA concentration without a topological constraint. If so, knots during genome mapping are not an intrinsic problem for genome mapping technology.

Translocation of a granular chain in a horizontally vibrated saw-tooth channel

We study the translocation mechanism of a granular chain in a horizontally vibrated saw-tooth channel using MD simulations and macro-scale experiments and show that the translocation speed is independent of the chain length as long as the chain length is larger than the spatial period of the saw-tooth. With the help of simulation, we explore the effect of geometry of the container and frequency and amplitude of vibration as well as chain flexibility on the chain drift speed. We observe that the most efficient transport is achieved when one of the channel walls is shifted with respect to the other wall by an amount equal to half the spatial period of the saw-tooth. We define a persistence length for the chain and show that the translocation speed depends on the ratio of persistence length over the spatial period of the saw-tooth. The optimum translocation occurs when this ratio is about 0.4. We also determine the optimum saw-tooth angle for the translocation of the chain as well as the optimum distance between the two walls. Some properties of this system are similar to those of polymer systems.

Synthesis of semicrystalline nanocapsular structures obtained by Thermally Induced Phase Separation in nanoconfinement

Phase separation of a polymer solution exhibits a peculiar behavior when induced in a nanoconfinement. The energetic constraints introduce additional interactions between the polymer segments that reduce the number of available configurations. In our work, this effect is exploited in a one-step strategy called nanoconfined-Thermally Induced Phase Separation (nc-TIPS) to promote the crystallization of polymer chains into nanocapsular structures of controlled size and shell thickness. This is accomplished by performing a quench step of a low-concentrated PLLA-dioxane-water solution included in emulsions of mean droplet size <500 nm acting as nanodomains. The control of nanoconfinement conditions enables not only the production of nanocapsules with a minimum mean particle diameter of 70 nm but also the tunability of shell thickness and its crystallinity degree. The specific properties of the developed nanocapsular architectures have important implications on release mechanism and loading capability of hydrophilic and lipophilic payload compounds.

Packaging stiff polymers in small containers: A molecular dynamics study

The question of how stiff polymers are able to pack into small containers is particularly relevant to the study of DNA packaging in viruses. A reduced version of the problem based on coarse-grained representations of the main components of the system-the DNA polymer and the spherical viral capsid-has been studied by molecular dynamics simulation. The results, involving longer polymers than in earlier work, show that as polymers become more rigid there is an increasing tendency to self-organize as spools that wrap from the inside out, rather than the inverse direction seen previously. In the final state, a substantial part of the polymer is packed into one or more coaxial spools, concentrically layered with different orientations, a form of packaging achievable without twisting the polymer.

Clusters of polyhedra in spherical confinement

Dense particle packing in a confining volume remains a rich, largely unexplored problem, despite applications in blood clotting, plasmonics, industrial packaging and transport, colloidal molecule design, and information storage. Here, we report densest found clusters of the Platonic solids in spherical confinement, for up to [Formula: see text] constituent polyhedral particles. We examine the interplay between anisotropic particle shape and isotropic 3D confinement. Densest clusters exhibit a wide variety of symmetry point groups and form in up to three layers at higher N. For many N values, icosahedra and dodecahedra form clusters that resemble sphere clusters. These common structures are layers of optimal spherical codes in most cases, a surprising fact given the significant faceting of the icosahedron and dodecahedron. We also investigate cluster density as a function of N for each particle shape. We find that, in contrast to what happens in bulk, polyhedra often pack less densely than spheres. We also find especially dense clusters at so-called magic numbers of constituent particles. Our results showcase the structural diversity and experimental utility of families of solutions to the packing in confinement problem.

Pore Translocation of Knotted Polymer Chains: How Friction Depends on Knot Complexity

Knots can affect the capability of polymers to translocate through narrow pores in complex and counterintuitive ways that are still relatively unexplored. We report here on a systematic theoretical and computational investigation of the driven translocation of flexible chains accommodating a large repertoire of knots trapped at the pore entrance. These include composite knots, which are the most common form of spontaneous entanglement in long polymers. Two unexpected results emerge from this study. First, the high force translocation compliance does not decrease systematically with knot complexity. Second, the response of composite knots is so dependent on the order of their factor knots, that their hindrance can even be lower than some of their prime components. We show that the resulting rich and seemingly disparate phenomenology can be captured in a seamless framework based on the mechanism by which the tractive force is propagated along and past the knots. The quantitative scheme can be viably used for predictive purposes and, hence, ought to be useful in applicative contexts, too.

Scaling behavior of topologically constrained polymer rings in a melt

Large scale molecular dynamics simulations on graphic processing units (GPUs) are employed to study the scaling behavior of ring polymers with various topological constraints in melts. Typical sizes of rings containing 3(1), 5(1) knots and catenanes made up of two unknotted rings scale like N(1/3) in the limit of large ring sizes N. This is consistent with the crumpled globule model and similar findings for unknotted rings. For small ring lengths knots occupy a significant fraction of the ring. The scaling of typical ring sizes for small N thus depends on the particular knot type and the exponent is generally larger than 0.4.

Knotting dynamics of DNA chains of different length confined in nanochannels

Langevin dynamics simulations are used to characterize the typical mechanisms governing the spontaneous tying, untying and the dynamical evolution of knots in coarse-grained models of DNA chains confined in nanochannels. In particular we focus on how these mechanisms depend on the chain contour length, Lc, at a fixed channel width D = 56 nm corresponding to the onset of the Odijk scaling regime where chain backfoldings and hence knots are disfavoured but not suppressed altogether. We find that the lifetime of knots grows significantly with Lc, while that of unknots varies to a lesser extent. The underlying kinetic mechanisms are clarified by analysing the evolution of the knot position along the chain. At the considered confinement, in fact, knots are typically tied by local backfoldings of the chain termini where they are eventually untied after a stochastic motion along the chain. Consequently, the lifetime of unknots is mostly controlled by backfoldings events at the chain ends, which is largely independent of Lc. The lifetime of knots, instead, increases significantly with Lc because knots can, on average, travel farther along the chain before being untied. The observed interplay of knots and unknots lifetimes underpins the growth of the equilibrium knotting probability of longer and longer chains at fixed channel confinement.

An Anisotropic Effective Model for the Simulation of Semiflexible Ring Polymers

We derive and introduce anisotropic effective pair potentials to coarse-grain solutions of semiflexible ring polymers of various lengths. The system has been recently investigated by means of full monomer-resolved computer simulations, revealing a host of unusual features and structure formation, which, however, cannot be captured by a rotationally averaged effective pair potential between the rings' centers of mass [Bernabei M.; Soft Matter2013, 9, 1287]. Our new coarse-graining strategy is to picture each ring as a soft, penetrable disk. We demonstrate that for the short- and intermediate-length rings the new model is quite capable of capturing the physics in a quantitative fashion, whereas for the largest rings, which resemble flexible ones, it fails at high densities. Our work opens the way for the physical justification of general, anisotropic penetrable interaction potentials.

Cluster Glasses of Semiflexible Ring Polymers

We present computer simulations of concentrated solutions of unknotted nonconcatenated semiflexible ring polymers. Unlike in their flexible counterparts, shrinking involves a strong energetic penalty, favoring interpenetration and clustering of the rings. We investigate the slow dynamics of the centers-of-mass of the rings in the amorphous cluster phase, consisting of disordered columns of oblate rings penetrated by bundles of prolate ones. Scattering functions reveal a striking decoupling of self- and collective motions. Correlations between centers-of-mass exhibit slow relaxation, as expected for an incipient glass transition, indicating the dynamic arrest of the cluster positions. However, self-correlations decay at much shorter time scales. This feature is a manifestation of the fast, continuous exchange and diffusion of the individual rings over the matrix of clusters. Our results reveal a novel scenario of glass formation in a simple monodisperse system, characterized by self-collective decoupling, soft caging, and mild dynamic heterogeneity.

How molecular knots can pass through each other

We propose a mechanism in which two molecular knots pass through each other and swap positions along a polymer strand. Associated free energy barriers in our simulations only amount to a few kBT, which may enable the interchange of knots on a single DNA strand.

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.

Colocalization of coregulated genes: a steered molecular dynamics study of human chromosome 19

The connection between chromatin nuclear organization and gene activity is vividly illustrated by the observation that transcriptional coregulation of certain genes appears to be directly influenced by their spatial proximity. This fact poses the more general question of whether it is at all feasible that the numerous genes that are coregulated on a given chromosome, especially those at large genomic distances, might become proximate inside the nucleus. This problem is studied here using steered molecular dynamics simulations in order to enforce the colocalization of thousands of knowledge-based gene sequences on a model for the gene-rich human chromosome 19. Remarkably, it is found that most (≈ 88%) gene pairs can be brought simultaneously into contact. This is made possible by the low degree of intra-chromosome entanglement and the large number of cliques in the gene coregulatory network. A clique is a set of genes coregulated all together as a group. The constrained conformations for the model chromosome 19 are further shown to be organized in spatial macrodomains that are similar to those inferred from recent HiC measurements. The findings indicate that gene coregulation and colocalization are largely compatible and that this relationship can be exploited to draft the overall spatial organization of the chromosome in vivo. The more general validity and implications of these findings could be investigated by applying to other eukaryotic chromosomes the general and transferable computational strategy introduced here.

Islands of conformational stability for filopodia

Filopodia are long, thin protrusions formed when bundles of fibers grow outwardly from a cell surface while remaining closed in a membrane tube. We study the subtle issue of the mechanical stability of such filopodia and how this depends on the deformation of the membrane that arises when the fiber bundle adopts a helical configuration. We calculate the ground state conformation of such filopodia, taking into account the steric interaction between the membrane and the enclosed semiflexible fiber bundle. For typical filopodia we find that a minimum number of fibers is required for filopodium stability. Our calculation elucidates how experimentally observed filopodia can obviate the classical Euler buckling condition and remain stable up to several tens of . We briefly discuss how experimental observation of the results obtained in this work for the helical-like deformations of enclosing membrane tubes in filopodia could possibly be observed in the acrosomal reactions of the sea cucumber Thyone, and the horseshoe crab Limulus. Any realistic future theories for filopodium stability are likely to rely on an accurate treatment of such steric effects, as analysed in this work.

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.

Knotting of linear DNA in nano-slits and nano-channels: a numerical study

The amount and type of self-entanglement of DNA filaments is significantly affected by spatial confinement, which is ubiquitous in biological systems. Motivated by recent advancements in single DNA molecule experiments based on nanofluidic devices and by the introduction of algorithms capable of detecting knots in open chains, we investigate numerically the entanglement of linear, open DNA chains confined inside nano-slits. The results regard the abundance, type, and length of occurring knots and are compared with recent findings for DNA inside nano-channels. In both cases, the width of the confining region, D, spans the 30 nm-1 μm range and the confined DNA chains are 1-4 μm long. It is found that the knotting probability is maximum for slit widths in the 70-100 nm range. However, over the considered DNA contour lengths, the maximum incidence of knots remains below 20%, while for channel confinement it tops 50%. Further differences of the entanglement are seen for the average contour length of the knotted region, which drops significantly below D ~100 nm for channel-confinement, while it stays approximately constant for slit-like confinement. These properties ought to reverberate in different kinetic properties of linear DNA depending on confinement and could be detectable experimentally or exploitable in nano-technological applications.