From toroidal to rod-like condensates of semiflexible polymers

Phase behaviour of semiflexible lattice polymers in poor-solvent solution: Mean-field theory and Monte Carlo simulations

We study a solution of interacting semiflexible polymers with curvature energy in poor-solvent conditions on the d-dimensional cubic lattice using mean-field theory and Monte Carlo computer simulations. Building upon past studies on a single chain, we construct a field-theory representation of the system and solve it within a mean-field approximation supported by Monte Carlo simulations in d = 3. A gas-liquid transition is found in the temperature-density plane that is then interpreted in terms of real systems. Interestingly, we find this transition to be independent of the bending rigidity. Past classical Flory-Huggins and Flory mean-field results are shown to be particular cases of this more general framework. Perspectives in terms of guiding experimental results towards optimal conditions are also proposed.

Energetic preference and topological constraint effects on the formation of DNA twisted toroidal bundles

DNA toroids are compact torus-shaped bundles formed by one or multiple DNA molecules being condensed from the solution due to various condensing agents. It has been shown that the DNA toroidal bundles are twisted. However, the global conformations of DNA inside these bundles are still not well understood. In this study, we investigate this issue by solving different models for the toroidal bundles and performing replica-exchange molecular dynamics (REMD) simulations for self-attractive stiff polymers of various chain lengths. We find that a moderate degree of twisting is energetically favorable for toroidal bundles, yielding optimal configurations of lower energies than for other bundles corresponding to spool-like and constant radius of curvature arrangements. The REMD simulations show that the ground states of the stiff polymers are twisted toroidal bundles with the average twist degrees close to those predicted by the theoretical model. Constant-temperature simulations show that twisted toroidal bundles can be formed through successive processes of nucleation, growth, quick tightening, and slow tightening of the toroid, with the two last processes facilitating the polymer threading through the toroid's hole. A relatively long chain of 512 beads has an increased dynamical difficulty to access the twisted bundle states due to the polymer's topological constraint. Interestingly, we also observed significantly twisted toroidal bundles with a sharp U-shaped region in the polymer conformation. It is suggested that this U-shaped region makes the formation of twisted bundles easier by effectively reducing the polymer length. This effect can be equivalent to having multiple chains in the toroid.

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.

Forming a Double-Helix Phase of Single Polymer Chains by the Cooperation between Local Structure and Nonlocal Attraction

Double-helix structures, such as DNA, are formed in nature to realize many unique functions. Inspired by this, researchers are pursuing strategies to design such structures from polymers. A key question is whether the double helix can be formed from the self-folding of a single polymer chain without specific interactions. Here, using Langevin dynamics simulation and theoretical analysis, we find that a stable double-helix phase can be achieved by the self-folding of single semiflexible polymers as a result of the cooperation between local structure and nonlocal attraction. The critical temperature of double-helix formation approximately follows T^{cri}∼ln(k_{θ}) and T^{cri}∼ln(k_{τ}), where k_{θ} and k_{τ} are the polymer bending and torsion stiffness, respectively. Furthermore, the double helix can exhibit major and minor grooves due to symmetric break for better packing. Our results provide a novel guide to the experimental design of the double helix.

Observations of three "re-entrant" twisted structures in double-stranded DNA dispersion particles

In this work we report on observations of new twisted (cholesteric-like) structures in liquid-crystalline dispersion particles with a hexagonal packing of double-stranded (ds) DNA molecules. Heating up to 80 °C of the DNA dispersion formed in a aqueous-salt solution with a high osmotic pressure (concentration) of poly(ethylene glycol) induces the formation of a new, optically active, spirally twisted structure of these molecules ("re-entrant" cholesteric structure (rest-A structure)). Cooling of this dispersion up to 22 °C is accompanied by the formation of an additional "re-entrant" cholesteric structure (rest-B). Modification of particles of the ds DNA dispersion (with rest-B structure) by replacing Na⁺ cations by multi-charged Gd³⁺ cations results in the third " re-entrant" structure (rest-C) despite a high density packing of ds nucleic acid molecules.

Supramolecular cyclization of semiflexible cylindrical micelles assembled from rod-coil graft copolymers

Uniform toroidal micelles can be constructed via the supramolecular cyclization of semiflexible cylindrical micelles, but revealing the conditions under which the cyclization occurs and the mechanism underlying the cyclization remains a challenge. In this study, we performed Brownian dynamics simulations of the supramolecular cyclization of semiflexible cylindrical micelles formed by rod-coil graft copolymers to obtain the cyclization conditions and understand the cyclization mechanism. It was found that the balance of the bending energy of the polymer backbones with the self-attraction energy between the pendant groups on the polymer backbones plays an important role in the cyclization process. A theoretical model based on this balance is developed to explain the cyclization mechanism, and the conditions required for realizing the supramolecular cyclization are obtained. The proposed mechanism is supported by our experimental findings regarding the supramolecular cyclization of polypeptide cylindrical micelles. The cyclization conditions and the revealed mechanism can guide further preparation of uniform toroidal micelles from semiflexible cylindrical micelles in an end-to-end closure manner.

Chain stiffness bridges conventional polymer and bio-molecular phases

Chain molecules play important roles in industry and in living cells. Our focus here is on distinct ways of modeling the stiffness inherent in a chain molecule. We consider three types of stiffnesses-one yielding an energy penalty for local bends (energetic stiffness) and the other two forbidding certain classes of chain conformations (entropic stiffness). Using detailed Wang-Landau microcanonical Monte Carlo simulations, we study the interplay between the nature of the stiffness and the ground state conformation of a self-attracting chain. We find a wide range of ground state conformations, including a coil, a globule, a toroid, rods, helices, and zig-zag strands resembling β-sheets, as well as knotted conformations allowing us to bridge conventional polymer phases and biomolecular phases. An analytical mapping is derived between the persistence lengths stemming from energetic and entropic stiffness. Our study shows unambiguously that different stiffnesses play different physical roles and have very distinct effects on the nature of the ground state of the conformation of a chain, even if they lead to identical persistence lengths.

Liquid-Crystalline Dispersions of Double-Stranded DNA

In this review, we compare the circular dichroism (CD) spectra of liquid-crystalline dispersion (LCD) particles formed in PEG-containing aqueous-salt solutions with the purpose of determining the packing of ds DNA molecules in these particles. Depending on the osmotic pressure of the solution, the phase exclusion of ds DNA molecules at room temperature results in the formation of LCD particles with the cholesteric or the hexagonal packing of molecules. The heating of dispersion particles with the hexagonal packing of the ds DNA molecules results in a new phase transition, accompanied by an appearance of a new optically active phase of ds DNA molecules. Our results are rationalized by way of a concept of orientationally ordered “quasinematic” layers formed by ds DNA molecules, with a parallel alignment in the hexagonal structure. These layers can adopt a twisted configuration with a temperature increase; and as a result of this process, a new, helicoidal structure of dispersion particle is formed (termed as the “re-entrant” cholesteric phase). To prove the cholesteric pattern of ds DNA molecules in this phase, the “liquid-like” state of the dispersion particles was transformed into its “rigid” counterpart.

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.

Self-assembly of bis-salphen compounds: from semiflexible chains to webs of nanorings

The recently-observed self-assembly of certain salphen-based compounds into neuron-like networks of microrings interconnected with nano-thin strings may suggest a new highly-potent tool for nanoscale patterning. However, the mechanism behind such phenomena needs to be clarified before they can be applied in materials design. Here we show that, in contrast with what was initially presumed, the emergence of a "rings-and-rods" pattern is unlikely to be explained by merging, collapse and piercing of vesicles as in previously reported cases of nanorings self-assembly via non-bonding interactions. We propose an alternative explanation: the compounds under study form a 1D coordination polymer, the fibres of which are elastic enough to fold into toroidal globules upon solvent evaporation, while being able to link separate chains into extended networks. This becomes possible because the structure of the compound's scaffold is found to adopt a very different conformation from that inferred in the original work. Based on ab initio and molecular dynamics calculations we propose a step-by-step description of self-assembly process of a supramolecular structure which explains all the observed phenomena in a simple and clear way. The individual roles of the compound' s scaffold structure, coordination centres, functional groups and solvent effects are also explained, opening a route to control the morphology of self-assembled networks and to synthesize new compounds exhibiting similar behaviour.

DNA like-charge attraction and overcharging by divalent counterions in the presence of divalent co-ions

Strongly correlated electrostatics of DNA systems has drawn the interest of many groups, especially the condensation and overcharging of DNA by multivalent counterions. By adding counterions of different valencies and shapes, one can enhance or reduce DNA overcharging. In this paper, we focus on the effect of multivalent co-ions, specifically divalent co-ions such as SO[Formula: see text]. A computational experiment of DNA condensation using Monte Carlo simulation in grand canonical ensemble is carried out where the DNA system is in equilibrium with a bulk solution containing a mixture of salt of different valency of co-ions. Compared to systems with purely monovalent co-ions, the influence of divalent co-ions shows up in multiple aspects. Divalent co-ions lead to an increase of monovalent salt in the DNA condensate. Because monovalent salts mostly participate in linear screening of electrostatic interactions in the system, more monovalent salt molecules enter the condensate leads to screening out of short-range DNA-DNA like charge attraction and weaker DNA condensation free energy. The overcharging of DNA by multivalent counterions is also reduced in the presence of divalent co-ions. Strong repulsions between DNA and divalent co-ions and among divalent co-ions themselves lead to a depletion of negative ions near the DNA surface as compared to the case without divalent co-ions. At large distances, the DNA-DNA repulsive interaction is stronger in the presence of divalent co-ions, suggesting that divalent co-ions' role is not only that of simple stronger linear screening.

Bundling in semiflexible polymers: A theoretical overview

Supramolecular assemblies of polymers are key modules to sustain the structure of cells and their function. The main elements of these assemblies are charged semiflexible polymers (polyelectrolytes) generally interacting via a long(er)-range repulsion and a short(er)-range attraction. The most common supramolecular structure formed by these polymers is the bundle. In the present paper, we critically review some recent theoretical and computational advances on the problem of bundle formation, and point a few promising directions for future work.

Structure and cohesive energy of dipolar helices

This paper deals with the investigation of cohesive energy in dipolar helices made up of hard spheres. Such tubular helical structures are ubiquitous objects in biological systems. We observe a complex dependence of cohesive energy on surface packing fraction and dipole moment distribution. As far as single helices are concerned, the lowest cohesive energy is achieved at the highest surface packing fraction. Besides, a striking non-monotonic behavior is reported for the cohesive energy as a function of the surface packing fraction. For multiple helices, we discover a new phase, exhibiting markedly higher cohesive energy. This phase is referred to as ZZ tube consisting of stacked crown rings (reminiscent of a pile of zig-zag rings), resulting in a local triangular arrangement with densely packed filaments parallel to the tube axis.

Theory and simulations of toroidal and rod-like structures in single-molecule DNA condensation

DNA condensation by multivalent cations plays a crucial role in genome packaging in viruses and sperm heads, and has been extensively studied using single-molecule experimental methods. In those experiments, the values of the critical condensation forces have been used to estimate the amplitude of the attractive DNA-DNA interactions. Here, to describe these experiments, we developed an analytical model and a rigid body Langevin dynamics assay to investigate the behavior of a polymer with self-interactions, in the presence of a traction force applied at its extremities. We model self-interactions using a pairwise attractive potential, thereby treating the counterions implicitly. The analytical model allows to accurately predict the equilibrium structures of toroidal and rod-like condensed structures, and the dependence of the critical condensation force on the DNA length. We find that the critical condensation force depends strongly on the length of the DNA, and finite-size effects are important for molecules of length up to 10(5)μm. Our Langevin dynamics simulations show that the force-extension behavior of the rod-like structures is very different from the toroidal ones, so that their presence in experiments should be easily detectable. In double-stranded DNA condensation experiments, the signature of the presence of rod-like structures was not unambiguously detected, suggesting that the polyamines used to condense DNA may protect it from bending sharply as needed in the rod-like structures.

Phase diagram of the ground states of DNA condensates

The phase diagram of the ground states of DNA in a bad solvent is studied for a semiflexible polymer model with a generalized local elastic bending potential characterized by a nonlinearity parameter x and effective self-attraction promoting compaction. x=1 corresponds to the wormlike chain model. Surprisingly, the phase diagram as well as the transition lines between the ground states are found to be a function of x. The model provides a simple explanation for the results of prior experimental and computational studies and makes predictions for the specific geometries of the ground states. The results underscore the impact of the form of the microscopic bending energy at macroscopic observable scales.

Effective interactions between fluid membranes

A self-consistent theory is proposed for the general problem of interacting undulating fluid membranes subject to the constraint that they do not interpenetrate. We implement the steric constraint via an exact functional integral representation and, through the use of a saddle-point approximation, transform it into a novel effective steric potential. The steric potential is found to consist of two contributions: one generated by zero-mode fluctuations of the membranes and the other by thermal bending fluctuations. For membranes of cross-sectional area S, we find that the bending fluctuation part scales with the intermembrane separation d as d-2 for d≪√S but crosses over to d-4 scaling for d≫√S, whereas the zero-mode part of the steric potential always scales as d-2. For membranes interacting exclusively via the steric potential, we obtain closed-form expressions for the effective interaction potential and for the rms undulation amplitude σ, which becomes small at low temperatures T and/or large bending stiffnesses κ. Moreover, σ scales as d for d≪√S but saturates at √kBTS/κ for d≫√S. In addition, using variational Gaussian theory, we apply our self-consistent treatment to study intermembrane interactions subject to different types of potentials: (i) the Moreira-Netz potential for a pair of strongly charged membranes with an intervening solution of multivalent counterions, (ii) an attractive square well, (iii) the Morse potential, and (iv) a combination of hydration and van der Waals interactions.