No Arabic abstract
We define a generalised model for three-stranded DNA consisting of two chains of one type and a third chain of a different type. The DNA strands are modelled by random walks on the three-dimensional cubic lattice with different interactions between two chains of the same type and two chains of different types. This model may be thought of as a classical analogue of the quantum three-body problem. In the quantum situation it is known that three identical quantum particles will form a triplet with an infinite tower of bound states at the point where any pair of particles would have zero binding energy. The phase diagram is mapped out, and the different phase transitions examined using finite-size scaling. We look particularly at the scaling of the DNA model at the equivalent Efimov point for chains up to 10000 steps in length. We find clear evidence of several bound states in the finite-size scaling. We compare these states with the expected Efimov behaviour.
We investigate the translocation of a single stranded DNA through a pore which fluctuates between two conformations, using coupled master equations. The probability density function of the first passage times (FPT) of the translocation process is calculated, displaying a triple, double or mono peaked behavior, depending on the interconversion rates between the conformations, the applied electric field, and the initial conditions. The cumulative probability function of the FPT, in a field-free environment, is shown to have two regimes, characterized by fast and slow timescales. An analytical expression for the mean first passage time of the translocation process is derived, and provides, in addition to the interconversion rates, an extensive characterization of the translocation process. Relationships to experimental observations are discussed.
We study voltage driven translocation of a single stranded (ss) DNA through a membrane channel. Our model, based on a master equation (ME) approach, investigates the probability density function (pdf) of the translocation times, and shows that it can be either double or mono-peaked, depending on the system parameters. We show that the most probable translocation time is proportional to the polymer length, and inversely proportional to the first or second power of the voltage, depending on the initial conditions. The model recovers experimental observations on hetro-polymers when using their properties inside the pore, such as stiffness and polymer-pore interaction.
We study the phase diagram of a system of $2times2times2$ hard cubes on a three dimensional cubic lattice. Using Monte Carlo simulations, we show that the system exhibits four different phases as the density of cubes is increased: disordered, layered, sublattice ordered, and columnar ordered. In the layered phase, the system spontaneously breaks up into parallel slabs of size $2times L times L$ where only a very small fraction cubes do not lie wholly within a slab. Within each slab, the cubes are disordered; translation symmetry is thus broken along exactly one principal axis. In the solid-like sublattice ordered phase, the hard cubes preferentially occupy one of eight sublattices of the cubic lattice, breaking translational symmetry along all three principal directions. In the columnar phase, the system spontaneously breaks up into weakly interacting parallel columns of size $2times 2times L$ where only a very small fraction cubes do not lie wholly within a column. Within each column, the system is disordered, and thus translational symmetry is broken only along two principal directions. Using finite size scaling, we show that the disordered-layered phase transition is continuous, while the layered-sublattice and sublattice-columnar transitions are discontinuous. We construct a Landau theory written in terms of the layering and columnar order parameters, which is able to describe the different phases that are observed in the simulations and the order of the transitions. Additionally, our results near the disordered-layered transition are consistent with the $O(3)$ universality class perturbed by cubic anisotropy as predicted by the Landau theory.
We present a new method for calculating internal forces in DNA structures using coarse-grained models and demonstrate its utility with the oxDNA model. The instantaneous forces on individual nucleotides are explored and related to model potentials, and using our framework, internal forces are calculated for two simple DNA systems and for a recently-published nanoscopic force clamp. Our results highlight some pitfalls associated with conventional methods for estimating internal forces, which are based on elastic polymer models, and emphasise the importance of carefully considering secondary structure and ionic conditions when modelling the elastic behaviour of single-stranded DNA. Beyond its relevance to the DNA nanotechnological community, we expect our approach to be broadly applicable to calculations of internal force in a variety of structures -- from DNA to protein -- and across other coarse-grained simulation models.
We study the dynamics of a double-stranded DNA (dsDNA) segment, as a semiflexible polymer, in a shear flow, the strength of which is customarily expressed in terms of the dimensionless Weissenberg number Wi. Polymer chains in shear flows are well-known to undergo tumbling motion. When the chain lengths are much smaller than the persistence length, one expects a (semiflexible) chain to tumble as a rigid rod. At low Wi, a polymer segment shorter than the persistence length does indeed tumble as a rigid rod. However, for higher Wi the chain does not tumble as a rigid rod, even if the polymer segment is shorter than the persistence length. In particular, from time to time the polymer segment may assume a buckled form, a phenomenon commonly known as Euler buckling. Using a bead-spring Hamiltonian model for extensible dsDNA fragments, we first analyze Euler buckling in terms of the oriented deterministic state (ODS), which is obtained as the steady-state solution of the dynamical equations by turning off the stochastic (thermal) forces at a fixed orientation of the chain. The ODS exhibits symmetry breaking at a critical Weissenberg number Wi$_{text c}$, analogous to a pitchfork bifurcation in dynamical systems. We then follow up the analysis with simulations and demonstrate symmetry breaking in computer experiments, characterized by a unimodal to bimodal transformation of the probability distribution of the second Rouse mode with increasing Wi. Our simulations reveal that shear can cause strong deformation for a chain that is shorter than its persistence length, similar to recent experimental observations.