No Arabic abstract
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.
Motivated by the structure of networks of cross-linked cytoskeletal biopolymers, we study the orientationally ordered phases in two-dimensional networks of randomly cross-linked semiflexible polymers. We consider permanent cross-links which prescribe a finite angle and treat them as quenched disorder in a semi-microscopic replica field theory. Starting from a fluid of un-cross-linked polymers and small polymer clusters (sol) and increasing the cross-link density, a continuous gelation transition occurs. In the resulting gel, the semiflexible chains either display long range orientational order or are frozen in random directions depending on the value of the crossing angle, the crosslink concentration and the stiffness of the polymers. A crossing angle $thetasim 2pi/M$ leads to long range $M$-fold orientational order, e.g., hexatic or tetratic for $theta=60^{circ}$ or $90^{circ}$, respectively. The transition is discontinuous and the critical cross-link density depends on the bending stiffness of the polymers and the cross-link geometry: the higher the stiffness and the lower $M$, the lower the critical number of cross-links. In between the sol and the long range ordered state, we always observe a gel which is a statistically isotropic amorphous solid (SIAS) with random positional and random orientational localization of the participating polymers.
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.
We calculate the scattering intensity of dilute and semi-dilute solutions of star polymers. The star conformation is described by a model introduced by Daoud and Cotton. In this model, a single star is regarded as a spherical region of a semi-dilute polymer solution with a local, position dependent screening length. For high enough concentrations, the outer sections of the arms overlap and build a semi-dilute solution (a sea of blobs) where the inner parts of the actual stars are embedded. The scattering function is evaluated following a method introduced by Auvray and de Gennes. In the dilute regime there are three regions in the scattering function: the Guinier region (low wave vectors, q R << 1) from where the radius of the star can be extracted; the intermediate region (1 << q R << f^(2/5)) that carries the signature of the form factor of a star with f arms: I(q) ~ q^(-10/3); and a high wavevector zone (q R >> f^(2/5)) where the local swollen structure of the polymers gives rise to the usual q^(-5/3) decay. In the semi-dilute regime the different stars interact strongly, and the scattered intensity acquires two new features: a liquid peak that develops at a reciprocal position corresponding to the star-star distances; and a new large wavevector contribution of the form q^(-5/3) originating from the sea of blobs.
We report Monte Carlo simulations of the self-assembly of supramolecular polymers based on a model of patchy particles. We find a first-order phase transition, characterized by hysteresis and nucleation, toward a solid bundle of polymers, of length much greater than the average gas phase length. We argue that the bundling transition is the supramolecular equivalent of the sublimation transition, that results from a weak chain-chain interaction. We provide a qualitative equation of state that gives physical insight beyond the specific values of the parameters used in our simulations.