Universality of the collapse transition of sticky polymers

Abstract

The universality of the swelling of the radius of gyration of a homopolymer relative to its value in the θ state, independent of polymer-solvent chemistry, in the crossover regime between θ and athermal solvent conditions, is well known. Here we study, by Brownian dynamics, a polymer model where a subset of monomers is labelled as "stickers". The mutual interaction of the stickers is more attractive than those of the other ("backbone") monomers, and has the additional important characteristic of "functionality" ϕ, i.e., the maximum number of stickers that can locally bind to a given sticker. A saturated bond formed in this manner remains bound until it breaks due to thermal fluctuations, a requirement which can be viewed as an additional Boolean degree of freedom that describes the bonding. This, in turn, makes the question of the order of the collapse transition a non-trivial one. Nevertheless, for the parameters that we have studied (in particular, ϕ = 1), we find a standard second-order θ collapse, using a renormalised solvent quality parameter that takes into account the increased average attraction due to the presence of stickers. We examine the swelling of the radius of gyration of such a sticky polymer relative to its value in the altered θ state, using a novel potential to model the various excluded volume interactions that occur between the monomers on the chain. We find that the swelling of such sticky polymers is identical to the universal swelling of homopolymers in the thermal crossover regime. Additionally, for our model, the Kuhn segment length under θ conditions is found to be the same for chains with and without stickers. arXiv:1903.07356v2 [cond-mat.soft]

Computing 3D Chromatin Configurations from Contact Probability Maps by Inverse Brownian Dynamics

“…The parameters a and b are determined by applying the two boundary conditions, namely, U SDK mn ¼ 0 at r à mn ¼ r à c and U SDK mn ¼ Àε mn at r à mn ¼ 2 1/6 s*. The appropriate choice of the cutoff radius r à c has been investigated extensively in a recent study (48), and it has been shown that a value of r à c ¼ 1.82s* leads to an accurate prediction of the static properties of a polymer chain in poor, theta, and good solvents. The same value is adopted in this study.…”

Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains

“…The repulsive part is unaffected by the choice of potential parameter ϵ ij . The attractive part is modelled by which, unlike the Lennard-Jones potential, smoothly reaches zero at the cut off radius r c = 1.82σ [44, 45]. We simulate the chromatin polymer using Brownian dynamics simulations, where the time evolution of the bead positions are governed by an Ito stochastic differential equation.…”

“…The repulsive part is unaffected by the choice of potential parameter ij . The attractive part is modelled by 1 2 ij cos (αr 2 ij + β) − 1 which, unlike the Lennard-Jones potential, smoothly reaches zero at the cut off radius r c = 1.82σ [44,45]. Parameters α and β control the r c (see supplementary information (SI)).…”

Wet and dry internal friction can be measured with the Jarzynski equality

“…constructed using a cosine function which smoothly approaches zero at the cut-off distance, r c . A detailed comparison of the SDK potential with the Lennard-Jones and the WCA potential has been performed recently [62].…”

“…Increasing the value of beyond zero results in a decrease in the solvent quality. A special feature of the SDK potential [62] is that modifying the value of allows one to tune the attractive interactions selectively, without affecting the repulsive branch of the pair-potential. This is in stark contrast to the more commonly used Lennard-Jones (LJ) potential for which changing the well-depth affects both the attractive and the repulsive branches.…”