Semiflexible Polymers in Spherical Confinement: Bipolar Orientational Order Versus Tennis Ball States

Heterogeneous flexibility can contribute to chromatin segregation in the cell nucleus

The highly and slightly condensed forms of chromatin, heterochromatin and euchromatin, respectively, segregate in the cell nucleus. Heterochromatin is more abundant in the nucleus periphery. Here we study the mechanism of heterochromatin segregation by modeling interphase chromosomes as diblock ring copolymers confined in a rigid spherical shell using molecular dynamics simulations. In our model, heterochromatin and euchromatin are distinguished by their bending stiffnesses only, while an interaction potential between the spherical shell and chromatin is used to model lamin-associated proteins. Our simulations indicate that in the absence of attractive interactions between the nuclear shell and the chromatin, most heterochromatin segregates towards the nuclear interior due to the depletion of less flexible heterochromatin segments from the nuclear periphery. This inverted chromatin distribution,which is opposite to the conventional case with heterochromatin dominating at the periphery, is in accord with experimental observations in rod cells. This "inversion" is also found to be independent of the heterochromatin concentration and chromosome number. The chromatin distribution at the periphery found in vivo can be recovered by further increasing the bending stiffness of heterochromatin segments or by turning on attractive interactions between the nuclear shell and heterochromatin. Our results indicate that the bending stiffness of chromatin could be a contributor to chromosome organization along with differential effects of HP1α-driven phase segregation and of loop extruders and interactions with the nuclear envelope and topological constraints.

Effects of self-avoidance on the packing of stiff rods on ellipsoids

Using a statistical-mechanics approach, we study the effects of geometry and self-avoidance on the ordering of slender filaments inside nonisotropic containers, considering cortical microtubules in plant cells, and packing of genetic material inside viral capsids as concrete examples. Within a mean-field approximation, we show analytically how the shape of the container, together with self-avoidance, affects the ordering of the stiff rods. We find that the strength of the self-avoiding interaction plays a significant role in the preferred packing orientation, leading to a first-order transition for oblate cells, where the preferred orientation changes from azimuthal, along the equator, to a polar one, when self-avoidance is strong enough. While for prolate spheroids the ground state is always a polar-like order, strong self-avoidance results with a deep metastable state along the equator. We compute the critical surface describing the transition between azimuthal and polar ordering in the three-dimensional parameter space (persistence length, eccentricity, and self-avoidance) and show that the critical behavior of this system is in fact related to the butterfly catastrophe model. We calculate the pressure and shear stress applied by the filament on the surface, and the injection force needed to be applied on the filament in order to insert it into the volume. We compare these results to the pure mechanical study where self-avoidance is ignored, and discuss similarities and differences.

Kinetic trapping organizes actin filaments within liquid-like protein droplets

Several actin-binding proteins (ABPs) phase separate to form condensates capable of curating the actin network shapes. Here, we use computational modeling to understand the principles of actin network organization within VASP condensate droplets. Our simulations reveal that the different actin shapes, namely shells, rings, and mixture states are highly dependent on the kinetics of VASP-actin interactions, suggesting that they arise from kinetic trapping. Specifically, we show that reducing the residence time of VASP on actin filaments reduces degree of bundling, thereby promoting assembly of shells rather than rings. We validate the model predictions experimentally using a VASP-mutant with decreased bundling capability. Finally, we investigate the ring opening within deformed droplets and found that the sphere-to-ellipsoid transition is favored under a wide range of filament lengths while the ellipsoid-to-rod transition is only permitted when filaments have a specific range of lengths. Our findings highlight key mechanisms of actin organization within phase-separated ABPs.

Programmable Potentials Choreograph Defects in a Colloidal Crystal Shell

Crystallization on spherical surfaces is obliged by topology to induce lattice defects. But controlling the organization of such defects remains a great challenge due to the long-range constraints of the curved geometry. Here, we report on DNA-coated colloids whose programmable interaction potentials can be used to regulate the arrangement of defects and even achieve perfect icosahedral order on a sphere. Combined simulations and theoretical analysis show how the potential can be tuned by changing the temperature, thereby controlling the number of defects. An explicit expression for the effective potential is derived, allowing us to distinguish the effects of entropic repulsion and enthalpic attraction. Altogether, the present findings provide insights into the physics of crystallization on curved spaces and may be used for designing desired crystal geometries.

Mixing and demixing arising from compression of two semiflexible polymer chains in nanochannels

We use molecular dynamics simulation to probe the non-equilibrium physics of two nanochannel-confined semiflexible polymers in a homogeneous flow field. We find that for sufficiently stiff chains the internal organization of the two chains takes the form of interwoven folds and circular coils. This organization can lead to mixing or demixing depending on chain stiffness and flow speed. At low and intermediate flow, the two chains adopt a folded configuration, which favours mixing. At high flow, the two chains adopt a predominantly coiled configuration. The coiled configuration results in demixing when the chains are compressed from an initially demixed condition and mixing when the chains are compressed from an initially mixed condition. We find that the mixing/demixing behaviour is governed by the ratio of the number of folded segments of one chain relative to the other at low flow and by the degree of coiling in both chains at high flow. For decreasing stiffness, the chains start to aggregate locally instead of mixing smoothly at low and intermediate flow. In the limit of completely flexible chains, the two chains either completely segregate at low flow, or adopt a locally demixed configuration consisting of large aggregates of one chain relative to the other that undergo complex stochastic dynamics, diffusing, disintegrating, and reforming at intermediate flow. The transition from complete segregation to the aggregate-dominated configuration occurs when the linear intra-chain ordering breaks down.

Adherent Moving of Polymers in Spherical Confined Binary Semiflexible Ring Polymer Mixtures

Based on the coarse-grained model, we used molecular dynamics methods to calculate and simulate a semiflexible long ring–semiflexible short ring blended polymer system confined in a hard sphere. We systematically studied the distribution and motion characteristics of the long ring chain. The results show that when the short ring is short enough (Lshort < 20), the long ring (Llong = 50) is separated from the blend system and then distributed against the inner wall. As the length of the short ring increases (Lshort ≥ 20), the long ring can no longer be separated from the blending system. Moreover, we found that the long ring demonstrates a random direction of adherent walking behavior on the inner surface of the hard sphere. The velocity of the long ring decreases with the increase in the short ring length Lshort. Specifically for Lshort ≥ 20, the system does not undergo phase separation and the speed of the long ring decreases sharply along with the long ring distributed inside the confined bulk. This is related to the inner wall layer moving faster than the inside bulk of the restricted system. Our simulation results can help us to understand the distribution of macromolecules in biological systems in confined systems, including the restricted chromosome partitioning distribution and packing structure of circular DNA molecules.

Entropic control of nanoparticle self-assembly through confinement

Entropy can be the sole driving force for the construction and regulation of ordered structures of soft matter systems. Specifically, under confinement, the entropic penalty could induce enhanced entropic effects which potentially generate visually ordered structures. Therefore, spatial confinement or a crowding environment offers an important approach to control entropy effects in these systems. Here, we review how spatial confinement-mediated entropic effects accurately and even dynamically control the self-assembly of nanoscale objects into ordered structures, focusing on our efforts towards computer simulations and theoretical analysis. First, we introduce the basic principle of entropic ordering through confinement. We then introduce the applications of this concept to various systems containing nanoparticles, including polymer nanocomposites, biological macromolecular systems and macromolecular colloids. Finally, the future directions and challenges for tailoring nanoparticle organization through spatial confinement-mediated entropic effects are detailed. We expect that this review could stimulate further efforts in the fundamental research on the relationship between confinement and entropy and in the applications of this concept for designer nanomaterials.

RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates

RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.

Entropy-Driven Unconventional Crystallization of Spherical Colloidal Nanocrystals Confined in Wide Cylinders

The densest packings of identical spherical colloidal nanocrystals in a thin cylinder generally give rise to confinement-induced chiral ordering. Here, we demonstrate that entropy can invalidate Pauling's packing rules for the nanocrystals confined in wide cylinders and novel ordered phases, where chiral ordering is broken, emerge. The nucleation and growth of spherical colloidal nanocrystals in the wide cylinders exhibit unique mechanisms which are distinctly different from that of thin ones. Furthermore, theoretical models which capture the essential physics of the ordering transitions are developed to reproduce the achiral ordering and reveal that the ordered phases are thermodynamically stable and stabilized through confinement-mediated entropic effect. These findings demonstrate that entropy arising from thermal motion can invalidate Pauling's packing rules of spherical colloidal nanocrystals confined in cylinders, which provides new insights into confinement physics of colloidal particles and might inspire nonintuitive design rules for the fabrication of novel ordered phases through confinement.

Cylindrical confinement of solutions containing semiflexible macromolecules: surface-induced nematic order versus phase separation

Solutions of semiflexible polymers confined in cylindrical pores with repulsive walls are studied by Molecular Dynamics simulations for a wide range of polymer concentrations. Both the case where both lengths are of the same order and the case when the persistence length by far exceeds the contour length are considered, and the enhancement of nematic order along the cylinder axis is characterized. With increasing density the character of the surface effect changes from depletion to the formation of a layered structure. For binary 50 : 50 mixtures of the two types of polymers an interplay between surface enrichment of the stiffer component and the isotropic-nematic transition is found, and a phase separated structure with cylindrical symmetry occurs, with the isotropic phase located around the cylinder axis. For melt densities the mixed nematic phase forms at the wall a layer with a screw-like structure of a tilted smectic phase. The observed behavior is tentatively interpreted in terms of the competition of the chain orientational entropy with entropy of mixing and excluded volume due to the wall.

Enhanced Heterogeneous Diffusion of Nanoparticles in Semiflexible Networks

The transport of nanoparticles in semiflexible networks, which form diverse principal structural components throughout living systems, is important in biology and biomedical applications. By combining large-scale molecular simulations as well as theoretical analysis, we demonstrate here that nanoparticles in polymer networks with semiflexible strands possess enhanced heterogeneous diffusion characterized by more evident hopping dynamics. Particularly, the hopping energy barrier approximates to linear dependence on confinement parameters in the regime of moderate rigidity, in contrast to the quadratic dependence of both its soft and hard counterparts. This nonmonotonic feature can be attributed to the competition between the conformation entropy and the bending energy regulated by the chain rigidity, captured by developing an analytical model of a hopping energy barrier. Moreover, these theoretical results agree reasonably well with previous experiments. The findings bear significance in unraveling the fundamental physics of substance transport confined in network-topological environments and would provide an explanation for the dynamics diversity of nanoparticles within various networks, biological or synthetic.

Defects and defect engineering in Soft Matter

Soft matter covers a wide range of materials based on linear or branched polymers, gels and rubbers, amphiphilic (macro)molecules, colloids, and self-assembled structures. These materials have applications in various industries, all highly important for our daily life, and they control all biological functions; therefore, controlling and tailoring their properties is crucial. One way to approach this target is defect engineering, which aims to control defects in the material's structure, and/or to purposely add defects into it to trigger specific functions. While this approach has been a striking success story in crystalline inorganic hard matter, both for mechanical and electronic properties, and has also been applied to organic hard materials, defect engineering is rarely used in soft matter design. In this review, we present a survey on investigations on defects and/or defect engineering in nine classes of soft matter composed of liquid crystals, colloids, linear polymers with moderate degree of branching, hyperbranched polymers and dendrimers, conjugated polymers, polymeric networks, self-assembled amphiphiles and proteins, block copolymers and supramolecular polymers. This overview proposes a promising role of this approach for tuning the properties of soft matter.

Structural transformations in tetravalent nematic shells induced by a magnetic field

The role of applied fields on the structure of liquid crystals confined to shell geometries has been studied in past theoretical work, providing strategies to produce liquid crystal shells with controlled defect structure or valence. However, the predictions of such studies have not been experimentally explored yet. In this work, we study the structural transformations undergone by tetravalent nematic liquid crystal shells under a strong uniform magnetic field, using both experiments and simulations. We consider two different cases in terms of shell geometry and initial defect symmetry: (i) homogeneous shells with four s = +1/2 defects in a tetrahedral arrangement, and (ii) inhomogeneous shells with four s = +1/2 defects localized in their thinner parts. Consistently with previous theoretical results, we observe that the initial defect structure evolves into a bipolar one, in a process where the defects migrate towards the poles. Interestingly, we find that the defect trajectories and dynamics are controlled by curvature walls that connect the defects by pairs. Based on the angle between Bs, the local projection of the magnetic field on the shell surface, and n+½, a vector describing the defect orientations, we are able to predict the nature and shape of those inversion walls, and therefore, the trajectory and dynamics of the defects. This rule, based on symmetry arguments, is consistent with both experiments and simulations and applies for shells that are either homogeneous or inhomogeneous in thickness. By modifying the angle between Bs and n+½, we are able to induce, in controlled way, complex routes towards the final bipolar state. In the case of inhomogeneous shells, the specific symmetry of the shell allowed us to observe a hybrid splay-bend Helfrich wall for the first time.

Rodlike molecules in extreme confinement

A unique feature of colloid particles and biopolymers is the molecule's intrinsic rigidity characterized by a molecular-level length scale. Under extreme confinement conditions at cellular scales or in nanodevices, these molecules can display orientational ordering accompanied by severe density depletion. Conventional liquid-crystal theories, such as the Oseen-Frank or Landau-de Gennes theories, cannot capture the essential molecular-level properties: the boundary effects, which extend to a distance of the rigidity length scale, and the drastic variations of the inhomogeneous molecular density. Here we show, based on a simple interpretation of the Onsager model, that rodlike molecules in extreme annular confinement produce unusual liquid-crystal defect structures that are independent phases from the patterns usually seen in a weaker confinement environment.

Entropy-induced Separation of Binary Semiflexible Ring Polymer Mixtures in Spherical Confinement

Coarse-grained molecular dynamics simulations are used to investigate the conformations of binary semiflexible ring polymers (SRPs) of two different lengths confined in a hard sphere. Segregated structures of SRPs in binary mixtures are strongly dependent upon the number density of system (ρ), the bending energy of long SRPs (Kb, long), and the chain length ratio of long to short SRPs (α). With a low ρ or a weak Kb, long at a small ratio α, long SRPs are immersed randomly in the matrix of short SRPs. As ρ and bending energy of long SRPs (Kb, long) are increased up to a certain value for a large ratio α, a nearly complete segregation between long and short SRPs is observed, which can be further characterized by the ratio of tangential and radial components of long SRPs velocity. These explicit segregated structures of the two components in spherical confinement are induced by a delicate competition between the entropic excluded volume (depletion) effects and bending contributions.

Orientationally ordered states of a wormlike chain in spherical confinement

One of the basic characteristics of a linear dsDNA molecule is its persistence length, typically of order 50 nm. The DNA chain inflicts a large energy penalty if it is bent sharply at that length scale. Viruses of bacteria, known as bacteriophages, typically have a dimension of a few tens of nanometers. Yet, it is known that a bacteriophage actively packages viral DNA inside the capsid and ejects it afterwards. Here, adopting a commonly used polymer model known as the wormlike chain, we answer an idealized question: Placing a linear DNA molecule inside a spherical cavity, what ordered states can we derive from known tools in statistical physics? Solving the model in a rigorous field-theory framework, we report a universal phase diagram for four orientationally ordered and disordered states, in terms of two relevant physical parameters.

Effects of spherical confinement and backbone stiffness on flexible polymer jamming

We use molecular simulations to study jamming of a crumpled bead-spring model polymer in a finite container and compare to jamming of repulsive spheres. After proper constraint counting, the onset of rigidity is seen to occur isostatically as in the case of repulsive spheres. Despite this commonality, the presence of the curved container wall and polymer backbone bonds introduce new mechanical properties. Notably, these include additional bands in the vibrational density of states that reflect the material structure as well as oscillations in local contact number and density near the wall but with lower amplitude for polymers. Polymers have fewer boundary contacts, and this low-density surface layer strongly reduces the global bulk modulus. We further show that bulk-modulus dependence on backbone stiffness can be described by a model of stiffnesses in series and discuss potential experimental and biological applications.

Density fields for branching, stiff networks in rigid confining regions

We develop a formalism to describe the equilibrium distributions for segments of confined branched networks consisting of stiff filaments. This is applicable to certain situations of cytoskeleton in cells, such as for example actin filaments with branching due to the Arp2/3 complex. We develop a grand ensemble formalism that enables the computation of segment density and polarisation profiles within the confines of the cell. This is expressed in terms of the solution to nonlinear integral equations for auxiliary functions. We find three specific classes of behaviour depending on filament length, degree of branching and the ratio of persistence length to the dimensions of the geometry. Our method allows a numerical approach for semi-flexible filaments that are networked.

Emergent topological phenomena in active polymeric fluids

Polymeric fluids show a wealth of topological phenomena, from entanglement and reptation at microscales to orientational ordering and defect production at macroscales, which can be explained by statistical-mechanical theories. In the presence of activity, the latter must be augmented by forces that cause spontaneous chain motion and fluid flow. Here, using such augmented Langevin equations, we study active polymeric solutions and melts composed of chains of hydrodynamically interacting stresslets. In a spherical volume, contractile chains are unstable and self-knot into entangled melts at both low and high densities. Extensile chains in the same geometry form an unentangled reptating state at low densities and an entangled, coherently moving, non-reptating state at high densities. On a spherical surface, contractile chains show transitions, with increasing areal density, between isotropic, orientationally ordered and micro-phase separated states. Extensile chains in the same geometry show a transition between isotropic and nematic states. In both cases, defects in orientationally ordered states are produced athermally and without conserving topological charge. Our work reproduces the phenomenology of several recent experiments, highlights the importance of hydrodynamic interactions in active polymer fluids, and suggests non-equilibrium kinetic routes to topological structures that are otherwise difficult to obtain in equilibrium.

Self-assembly of semiflexible polymers confined to thin spherical shells

Confinement effects are critical for stiff macromolecules in biological cells, vesicles, and other systems in soft matter. For these molecules, the competition between the packing entropy and the enthalpic cost of bending is further shaped by strong confinement effects. Through coarse-grained molecular dynamics simulations, we explore the self-assembly of semiflexible polymers confined in thin spherical shells for various chain lengths, chain stiffnesses, and shell thicknesses. Here, we focus on the case where the contour and persistence length of the polymers are comparable to the radius of the confining cavity. The range of ordered structures is analyzed using several order parameters to elucidate the nature of orientational ordering in different parameter regimes. Previous simulations have revealed the emergence of bipolar and quadrupolar topological defects on the surface when the entire cavity was filled with a concentrated polymer solution [Phys. Rev. Lett., 2017, 118, 217803]. In contrast, spherical shell confinement restricts the appearance of a bipolar order. Instead, only the extent of the quadrupolar order changes with chain stiffness, as evidenced by the relative motion of topological defects. In the case of monolayers, we observe a nematic to smectic transition accompanied by a change in the nematic grain-size distribution as the contour length was decreased.

Quantized Self-Assembly of Discotic Rings in a Liquid Crystal Confined in Nanopores

Disklike molecules with aromatic cores spontaneously stack up in linear columns with high, one-dimensional charge carrier mobilities along the columnar axes, making them prominent model systems for functional, self-organized matter. We show by high-resolution optical birefringence and synchrotron-based x-ray diffraction that confining a thermotropic discotic liquid crystal in cylindrical nanopores induces a quantized formation of annular layers consisting of concentric circular bent columns, unknown in the bulk state. Starting from the walls this ring self-assembly propagates layer by layer towards the pore center in the supercooled domain of the bulk isotropic-columnar transition and thus allows one to switch on and off reversibly single, nanosized rings through small temperature variations. By establishing a Gibbs free energy phase diagram we trace the phase transition quantization to the discreteness of the layers' excess bend deformation energies in comparison to the thermal energy, even for this near room-temperature system. Monte Carlo simulations yielding spatially resolved nematic order parameters, density maps, and bond-orientational order parameters corroborate the universality and robustness of the confinement-induced columnar ring formation as well as its quantized nature.

Conformations and orientational ordering of semiflexible polymers in spherical confinement

Semiflexible polymers in lyotropic solution confined inside spherical nanoscopic "containers" with repulsive walls are studied by molecular dynamics simulations and density functional theory, as a first step to model confinement effects on stiff polymers inside of miniemulsions, vesicles, and cells. It is shown that the depletion effects caused by the monomer-wall repulsion depend distinctly on the radius R of the sphere. Further, nontrivial orientational effects occur when R, the persistence length ℓp, and the contour length L of the polymers are of similar magnitude. At intermediate densities, a "shell" of wall-attached chains is forming, such that the monomers belonging to those chains are in a layer at about the distance of one monomer from the container wall. At the same time, the density of the centers of mass of these chains is peaked somewhat further inside, but still near the wall. However, the arrangement of chains is such that the total monomer density is almost uniform in the sphere, apart from a small layering peak at the wall. It is shown that excluded volume effects among the monomers are crucial to account for this behavior, although they are negligible for comparable isolated single semiflexible chains of the same length.