No Arabic abstract
We present results for a lattice model of bio-polymers where the type of $beta$-sheet formation can be controlled by different types of hydrogen bonds depending on the relative orientation of close segments of the polymer. Tuning these different interaction strengths leads to low-temperature structures with different types of orientational order. We perform simulations of this model and so present the phase diagram, ascertaining the nature of the phases and the order of the transitions between these phases.
We propose new polymer models for Monte Carlo simulation and apply them to a polymer chain confined in a relatively thin box which has both curved and flat sides, and show that either an ideal or an excluded-volume chain spends more time in the curved region than in the flat region. The ratio of the probability of finding a chain in the curved region and in flat region increases exponentially with increasing chain length. The results for ideal chains are quantitatively consistent with a previously published theory. We find that the same effect appears with excluded-volume chains and a similar scaling relation can be applied to them up to a certain length of the polymer.
It has become clear in recent years that the simple uniform wormlike chain model needs to be modified in order to account for more complex behavior which has been observed experimentally in some important biopolymers. For example, the large flexibility of short ds-DNA has been attributed to kink or hinge defects. In this paper, we calculate analytically, within the weak bending approximation, the force-extension relation of a wormlike chain with a permanent hinge defect along its contour. The defect is characterized by its bending energy (which can be zero, in the completely flexible case) and its position along the polymer contour. Besides the bending rigidity of the chain, these are the only parameters which describe our model. We show that a hinge defect causes a significant increase in the differential tensile compliance of a pre-stressed chain. In the small force limit, a hinge defect significantly increases the entropic elasticity. Our results apply to any pair of semiflexible segments connected by a hinge. As such, they may also be relevant to cytoskeletal filaments (F-actin, microtubules), where one may treat the cross-link connecting two filaments as a hinge defect.
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.