Folding and unfolding kinetics of a single semiflexible polymer

Expansion of Single Chains Released from a Spherical Cavity

A two-stage model is developed to explain the phenomena of chain expansion, released from a confining cavity. In the first stage, the chain is assumed to expand as a sphere, while in the second stage it expands like a coil. The kinetic equations for the variation of chain size are derived in the two stages by balancing the rate of the free energy change with the rate of the energy dissipation. Langevin dynamics simulations are then performed to examine the theory. We find that the expansion process is dominated by the second stage and the evolution of chain size follows, mainly, the predicted curve for coil expansion, which depends on the chain length and is not sensitive to the confining volume fraction. It permits to define the expansion time for the process. Further study reveals that the chain does undergo a spherical expansion in the first stage with the characteristic time much shorter than the one for the second stage. As a consequence, the first-stage variation of chain size can be regarded as an add-on to the principal curve of expansion designated by the second stage. The scaling behaviors and the associated scaling exponents are analyzed in details. The simulation results well support the theory.

Competing single-chain folding and multi-chain aggregation pathways control solution-phase aggregate morphology of organic semiconducting polymers

Understanding the solution-phase behaviour of organic semiconducting polymers is important for systematically improving the performance of devices based on solution-processed thin films of these molecules. Conventional polymer theory predicts that polymer conformations become more compact as solvent quality decreases, but recent experiments have shown the high-performance organic-semiconducting polymer P(NDI2OD-T2) to form extended rod-like aggregates much larger than a single chain in poor solvents, with the formation of these extended aggregates correlated with enhanced electron mobility in films deposited from these solutions. We explain the unexpected formation of extended aggregates using a novel coarse-grained simulation model of P(NDI2OD-T2) that we have developed to study the effect of solvent quality on its solution-phase behaviour. In poor solvents, we find that aggregation through only a few monomers gives effectively inseparable chains, leading to the formation of extended structures of partially overlapping chains via non-equilibrium assembly. This behaviour requires that multi-chain aggregation occurs faster than chain folding, which we show is the case for the chain lengths and concentrations shown experimentally to form rod-like aggregates. This kinetically controlled process introduces a dependence of aggregate structure on concentration, chain length, and chain flexibility, which we show is able to reconcile experimental findings and is generalisable to the solution-phase assembly of other semiflexible polymers.

Organization of associating or crosslinked actin filaments in confinement

A key factor of actin cytoskeleton organization in cells is the interplay between the dynamical properties of actin filaments and cell geometry, which restricts, confines and directs their orientation. Crosslinking interactions among actin filaments, together with geometrical cues and regulatory proteins can give rise to contractile rings in dividing cells and actin rings in neurons. Motivated by recent in vitro experiments, in this work we performed computer simulations to study basic aspects of the interplay between confinement and attractive interactions between actin filaments. We used a spring-bead model and Brownian dynamics to simulate semiflexible actin filaments that polymerize in a confining sphere with a rate proportional to the monomer concentration. We model crosslinking, or attraction through the depletion interaction, implicitly as an attractive short-range potential between filament beads. In confining geometries smaller than the persistence length of actin filaments, we show rings can form by curving of filaments of length comparable to, or longer than the confinement diameter. Rings form for optimal ranges of attractive interactions that exist in between open bundles, irregular loops, aggregated, and unbundled morphologies. The probability of ring formation is promoted by attraction to the confining sphere boundary and decreases for large radii and initial monomer concentrations, in agreement with prior experimental data. The model reproduces ring formation along the flat plane of oblate ellipsoids.

The folding pathways and thermodynamics of semiflexible polymers

Inspired by the protein folding and DNA packing, we have systematically studied the thermodynamic and kinetic behaviors of single semiflexible homopolymers by Langevin dynamics simulations. In line with experiments, a rich variety of folding products, such as rod-like bundles, hairpins, toroids, and a mixture of them, are observed in the complete diagram of states. Moreover, knotted structures with a significant population are found in a certain range of bending stiffness in thermal equilibrium. As the solvent quality becomes poorer, the population of the intermediate occurring in the folding process increases, which leads to a severe chevron rollover for the folding arm. However, the population of the intermediates in the unfolding process is very low, insufficient to induce unfolding arm rollover. The total types of folding pathways from the coil state to the toroidal state for a semiflexible polymer chain remain unchanged by varying the solvent quality or temperature, whereas the kinetic partitioning into different folding events can be tuned significantly. In the process of knotting, three types of mechanisms, namely, plugging, slipknotting, and sliding, are discovered. Along the folding evolution, a semiflexible homopolymer chain can knot at any stage of folding upon leaving the extended coil state, and the probability to find a knot increases with chain compactness. In addition, we find rich types of knotted topologies during the folding of a semiflexible homopolymer chain. This study should be helpful in gaining insight into the general principles of biopolymer folding.

Nanoconfinement greatly speeds up the nucleation and the annealing in single-DNA collapse

Manipulating and measuring single-molecule dynamics and reactions in nanofluidics is a rapidly growing field with broad applications in developing new biotechnologies, understanding nanoconfinement effects in vivo, and exploring new phenomena in confinement. In this work, we investigate the kinetics of DNA collapse in nanoslits using single T4-DNA (165.6 kbp) and λ-DNA (48.5 kbp), with particular focus on the measurement of the nucleation and annealing times. Fixing the ethanol concentration at 35% and varying the slit height from 2000 to 31 nm, the nucleation time dramatically decreases from more than 1 hour to a few minutes or less. The increased collapsed rate results from the larger free energy experienced by coiled DNA in confinement relative to compacted DNA. Our results also shed light on other conformational transitions in confinement, such as protein folding.

Collapse kinetics of vibrated granular chains

The kinetics of the collapse of the coil state into condensed states is studied with vibrated granular chain composed of N metal beads partially immersed in water. The radius of gyration of the chain, R(g) is measured. For short chains (N < 140), disk-like condensed state is formed and R(g) decreases with time such that the function ΔR(g)(2) (≡ R(g)(2) - R(g)(2)(∞)) = A e(-t/τ), where the relaxation time τ follows a power-law dependence on the chain length N with an exponent γ = 1.9 ± 0.2. For the chains with length N ≥ 300, rod-like clusters are observed during the initial stage of collapse and R(g)(2) = R(g)(2)(0) - Bt(β), with β = 0.6 ± 0.1. In the coarsening stage, the exponential dependence of ΔR(g)(2) on time still holds, however, the relaxation time τ fluctuates and has no simple dependence on N. Furthermore, the time dependence of the averaged radius of gyration of the individual clusters, R(g,cl) can be described by the theory of Lifshitz and Slyozov. A peak in the structure function of long chains is observed in the initial stage of the collapse transition. The collapse transition in the bead chains is a first order phase transition. However, features of the spinodal decomposition are also observed.

Compression and self-entanglement of single DNA molecules under uniform electric field

We experimentally study the effects of a uniform electric field on the conformation of single DNA molecules. We demonstrate that a moderate electric field (~200 V/cm) strongly compresses isolated DNA polymer coils into isotropic globules. Insight into the nature of these compressed states is gained by following the expansion of the molecules back to equilibrium after halting the electric field. We observe two distinct types of expansion modes: a continuous molecular expansion analogous to a compressed spring expanding, and a much slower expansion characterized by two long-lived metastable states. Fluorescence microscopy and stretching experiments reveal that the metastable states are the result of intramolecular self-entanglements induced by the electric field. These results have broad importance in DNA separations and single molecule genomics, polymer rheology, and DNA-based nanofabrication.

Effect of confinement on the collapsing mechanism of a flexible polymer chain

In this paper, Brownian dynamics simulation (BDS) studies are executed to demonstrate the distinctive influences of the extent of confinement on the collapsing mechanism and kinetics of a flexible hydrophobic polymer chain in a poor solvent. The collapsing behavior is quantified by the time of collapse, which below a critical dimension of the confinement (h(c)), encounters a drastic reduction with a further strengthening in the degree of confinement. For dimensions greater than this critical one, the collapse occurs through the well-known hydrodynamic interaction (HI) controlled multiple-globule-mediated mechanisms. However, for channel dimensions less than this critical one, the collapse mechanism is drastically altered. Under such circumstances, the collapse gets predominantly controlled by the confinement effects (with negligible contribution of the HIs) and occurs via the formation of a single central globule. This central globule rapidly engulfs the noncondensed polymer segments, and in the process largely hastens up the collapsing event. Under such circumstances, the collapse time is found to decrease linearly with decrements in the channel height. On the contrary, for channel heights greater than h(c), the multiple-globule-mediated collapse is characterized by a collapse time that shows an exponential dependence on the channel height, rapidly attaining a state in which the confinement effect becomes inconsequential and HIs dictate the entire collapsing behavior. We further propose detailed arguments based on physical reasoning as well as free energy estimations to conclusively support the qualitative and quantitative nature of influences of the confinement on the polymer collapse.

Refolding dynamics of stretched biopolymers upon force quench

Single-molecule force spectroscopy methods can be used to generate folding trajectories of biopolymers from arbitrary regions of the folding landscape. We illustrate the complexity of the folding kinetics and generic aspects of the collapse of RNA and proteins upon force quench by using simulations of an RNA hairpin and theory based on the de Gennes model for homopolymer collapse. The folding time, tau(F), depends asymmetrically on deltaf(S) = f (S) - f (m) and deltaf (Q) = f (m) - f (Q) where f (S) (f (Q)) is the stretch (quench) force and f (m) is the transition midforce of the RNA hairpin. In accord with experiments, the relaxation kinetics of the molecular extension, R(t), occurs in three stages: A rapid initial decrease in the extension is followed by a plateau and finally, an abrupt reduction in R(t) occurs as the native state is approached. The duration of the plateau increases as lambda = tau (Q)/tau (F) decreases (where tau (Q) is the time in which the force is reduced from f (S) to f (Q)). Variations in the mechanisms of force-quench relaxation as lambda is altered are reflected in the experimentally measurable time-dependent entropy, which is computed directly from the folding trajectories. An analytical solution of the de Gennes model under tension reproduces the multistage stage kinetics in R(t). The prediction that the initial stages of collapse should also be a generic feature of polymers is validated by simulation of the kinetics of toroid (globule) formation in semiflexible (flexible) homopolymers in poor solvents upon quenching the force from a fully stretched state. Our findings give a unified explanation for multiple disparate experimental observations of protein folding.

Dynamics of polymer decompression: expansion, unfolding, and ejection

The dynamics of polymer decompression from a compact state to swollen conformation can be formally described as nonlinear diffusion. We discuss two basic examples: (i) the expansion, or unfolding from a compact state, and (ii) the ejection of a compressed polymer through a pore. The problem can be solved exactly for case (i), but not for case (ii). Even in such situations, a scheme called uniform approximation is shown to be useful to get a physical insight involved. Its application to case (ii) is able to account for conflicting numerical data in a consistent way.