No Arabic abstract
Stochastic simulations are used to characterize the knotting distributions of random ring polymers confined in spheres of various radii. The approach is based on the use of multiple Markov chains and reweighting techniques, combined with effective strategies for simplifying the geometrical complexity of ring conformations without altering their knot type. By these means we extend previous studies and characterize in detail how the probability to form a given prime or composite knot behaves in terms of the number of ring segments, $N$, and confining radius, $R$. For $ 50 le N le 450 $ we show that the probability of forming a composite knot rises significantly with the confinement, while the occurrence probability of prime knots are, in general, non-monotonic functions of 1/R. The dependence of other geometrical indicators, such as writhe and chirality, in terms of $R$ and $N$ is also characterized. It is found that the writhe distribution broadens as the confining sphere narrows.
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 30nm- 1mu m range and the confined DNA chains are 1 to 4mu m long. It is found that the knotting probability is maximum for slit widths in the 70-100nm 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 ~100nm 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.
We present the results of analytic calculations and numerical simulations of the behaviour of a new class of chain molecules which we call thick polymers. The concept of the thickness of such a polymer, viewed as a tube, is encapsulated by a special three body interaction and impacts on the behaviour both locally and non-locally. When thick polymers undergo compaction due to an attractive self-interaction, we find a new type of phase transition between a compact phase and a swollen phase at zero temperature on increasing the thickness. In the vicinity of this transition, short tubes form space filling helices and sheets as observed in protein native state structures. Upon increasing the chain length, or the number of chains, we numerically find a crossover from secondary structure motifs to a quite distinct class of structures akin to the semi-crystalline phase of polymers or amyloid fibers in polypeptides.
Using a recently developed bead-spring model for semiflexible polymers that takes into account their natural extensibility, we report an efficient algorithm to simulate the dynamics for polymers like double-stranded DNA (dsDNA) in the absence of hydrodynamic interactions. The dsDNA is modelled with one bead-spring element per basepair, and the polymer dynamics is described by the Langevin equation. The key to efficiency is that we describe the equations of motion for the polymer in terms of the amplitudes of the polymers fluctuation modes, as opposed to the use of the physical positions of the beads. We show that, within an accuracy tolerance level of $5%$ of several key observables, the model allows for single Langevin time steps of $approx1.6$, 8, 16 and 16 ps for a dsDNA model-chain consisting of 64, 128, 256 and 512 basepairs (i.e., chains of 0.55, 1.11, 2.24 and 4.48 persistence lengths) respectively. Correspondingly, in one hour, a standard desktop computer can simulate 0.23, 0.56, 0.56 and 0.26 ms of these dsDNA chains respectively. We compare our results to those obtained from other methods, in particular, the (inextensible discretised) WLC model. Importantly, we demonstrate that at the same level of discretisation, i.e., when each discretisation element is one basepair long, our algorithm gains about 5-6 orders of magnitude in the size of time steps over the inextensible WLC model. Further, we show that our model can be mapped one-on-one to a discretised version of the extensible WLC model; implying that the speed-up we achieve in our model must hold equally well for the latter. We also demonstrate the use of the method by simulating efficiently the tumbling behaviour of a dsDNA segment in a shear flow.
Biological molecules can form hydrogen bonds between nearby residues, leading to helical secondary structures. The associated reduction of configurational entropy leads to a temperature dependence of this effect: the helix-coil transition. Since the formation of helices implies a dramatic shortening of the polymer dimensions, an externally imposed end-to-end distance R affects the equilibrium helical fraction of the polymer and the resulting force- extension curves show anomalous plateau regimes. In this article, we investigate the behaviour of a cross-linked network of such helicogenic molecules, particularly, focusing on the coupling of the (average) helical content present in a network to the externally imposed strain. We show that both an elongation and compression can lead to an increase in helical domains under appropriate conditions.
We explore the effect of an attractive interaction between parallel-aligned polymers, which are perpendicularly grafted on a substrate. Such an attractive interaction could be due to, e.g., reversible cross-links. The competition between permanent grafting favoring a homogeneous state of the polymer brush and the attraction, which tends to induce in-plane collapse of the aligned polymers, gives rise to an instability of the homogeneous phase to a bundled state. In this latter state the in-plane translational symmetry is spontaneously broken and the density is modulated with a finite wavelength, which is set by the length scale of transverse fluctuations of the grafted polymers. We analyze the instability for two models of aligned polymers: directed polymers with a line tension and weakly bending chains with a bending stiffness.