Stepwise unfolding of collapsed polymers

Self-Attractive Semiflexible Polymers under an External Force Field

The dynamical response of a tethered semiflexible polymer with self-attractive interactions and subjected to an external force field is numerically investigated by varying stiffness and self-interaction strength. The chain is confined in two spatial dimensions and placed in contact with a heat bath described by the Brownian multi-particle collision method. For strong self-attraction the equilibrium conformations range from compact structures to double-stranded chains, and to rods when increasing the stiffness. Under the external field at small rigidities, the initial close-packed chain is continuously unwound by the force before being completely elongated. For double-stranded conformations the transition from the folded state to the open one is sharp being steeper for larger stiffnesses. The discontinuity in the transition appears in the force-extension relation, as well as in the probability distribution function of the gyration radius. The relative deformation with respect to the equilibrium case along the direction normal to the force is found to decay as the inverse of the applied force.

Stretching semiflexible polymer chains: evidence for the importance of excluded volume effects from Monte Carlo simulation

Semiflexible macromolecules in dilute solution under very good solvent conditions are modeled by self-avoiding walks on the simple cubic lattice (d = 3 dimensions) and square lattice (d = 2 dimensions), varying chain stiffness by an energy penalty ε(b) for chain bending. In the absence of excluded volume interactions, the persistence length l(p) of the polymers would then simply be l(p) = l(b)(2d - 2)(-1)q(b) (-1) with q(b) = exp(-ε(b)/k(B)T), the bond length l(b) being the lattice spacing, and k(B)T is the thermal energy. Using Monte Carlo simulations applying the pruned-enriched Rosenbluth method (PERM), both q(b) and the chain length N are varied over a wide range (0.005 ≤ q(b) ≤ 1, N ≤ 50,000), and also a stretching force f is applied to one chain end (fixing the other end at the origin). In the absence of this force, in d = 2 a single crossover from rod-like behavior (for contour lengths less than l(p)) to swollen coils occurs, invalidating the Kratky-Porod model, while in d = 3 a double crossover occurs, from rods to Gaussian coils (as implied by the Kratky-Porod model) and then to coils that are swollen due to the excluded volume interaction. If the stretching force is applied, excluded volume interactions matter for the force versus extension relation irrespective of chain stiffness in d = 2, while theories based on the Kratky-Porod model are found to work in d = 3 for stiff chains in an intermediate regime of chain extensions. While for q(b) ≪ 1 in this model a persistence length can be estimated from the initial decay of bond-orientational correlations, it is argued that this is not possible for more complex wormlike chains (e.g., bottle-brush polymers). Consequences for the proper interpretation of experiments are briefly discussed.

Force-induced conformational transition in a system of interacting stiff polymers: application to unfolding

We consider a stiff polymer chain in poor solvent and apply a force at one end of the chain. We find that by varying the stiffness parameter, the polymer undergoes a transition from the globule state to the foldedlike state. The conformation of the folded state mimics the beta sheet as seen in the titin molecule. Using the exact enumeration technique, we study the extension-force and force-temperature diagrams of such a system. The force-temperature diagram shows the re-entrance behavior for a flexible chain. However, for a stiff chain this re-entrance behavior is absent and there is an enhancement in theta temperature with the rise of stiffness. We further propose that the internal information about the frozen structure of a polymer can be read from the distribution of end-to-end distance which shows a sawtoothlike behavior.

Nonlinear elasticity of single collapsed polyelectrolytes

Nonlinear elastic responses of short and stiff polyelectrolytes are investigated by dynamic simulations on a single-molecule level. When a polyelectrolyte condensate undergoes a mechanical unfolding, two types of force-extension curves--i.e., a force plateau and a stick-release pattern--are observed depending on the strength of the electrostatic interaction. We provide a physical interpretation of such force-extension behavior in terms of intramolecular structures of the condensates. We also describe charge distributions of counterions condensed onto a polyelectrolyte, which clarify formation of one-dimensional strongly correlated liquid at large Coulomb coupling regime. These findings may provide significant insights into the relationship between a molecular elasticity and a molecular mechanism of like-charge attractions observed in a wide range of charged biopolymer systems.

Force spectroscopy of single multidomain biopolymers: a master equation approach

Experiments using atomic force microscopy for unfolding single multidomain biopolymers cover a broad range of time scales from equilibrium to non-equilibrium. A master equation approach allows to identify and treat coherently three dynamical regimes for increasing linear ramp velocity: i) an equilibrium regime, ii) a transient regime where refolding events still occur, and iii) a saw-tooth regime without any refolding events. For each regime, analytical approximations are derived and compared to numerically investigated examples. We analyze in the framework of this model also a periodic experimental protocol instead of a linear ramp. In this case, a major simplification arises if the dynamics can be restricted to an effectively two-dimensional subspace. For transitions with an intermediate meta-stable state, like Immunoglobulin27, a refined model allows to extract previously unknown molecular parameters related to this meta-stable state.

Different pathways in mechanical unfolding/folding cycle of a single semiflexible polymer

Kinetics of conformational change of a semiflexible polymer under mechanical external field were investigated with Langevin dynamics simulations. It is found that a semiflexible polymer exhibits large hysteresis in mechanical folding/unfolding cycle even with a slow operation, whereas in a flexible polymer, the hysteresis almost disappears at a sufficiently slow operation. This suggests that the essential features of the structural transition of a semiflexible polymer should be interpreted at least on a two-dimensional phase space. The appearance of such large hysteresis is discussed in relation to different pathways in the loading and unloading processes. By using a minimal two-variable model, the hysteresis loop is described in terms of different pathways on the transition between two stable states.