Next Article in Journal / Special Issue
Plasmonic Biosensors Based on Deformed Graphene
Previous Article in Journal
Characterisation of Amyloid Aggregation and Inhibition by Diffusion-Based Single-Molecule Fluorescence Techniques
Previous Article in Special Issue
RBL-2H3 Mast Cell Receptor Dynamics in the Immunological Synapse
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

College of Physics and Optoelectronic Engineering, Harbin Engineering University, 150001 Harbin, China
*
Author to whom correspondence should be addressed.
Biophysica 2022, 2(4), 525-537; https://doi.org/10.3390/biophysica2040044
Submission received: 24 October 2022 / Revised: 15 November 2022 / Accepted: 22 November 2022 / Published: 24 November 2022
(This article belongs to the Collection Feature Papers in Biophysics)

Abstract

:
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.

1. Introduction

Cyclic polymer-restricted systems can be found in biological systems, such as the circular chromosome [1,2], circular DNA of viruses [3,4,5,6] and cyclic peptides [7]. The physical effects of topological constraints and restrictions play an important role in the structure and function of individual genetic material [8,9,10,11,12], such as the shape of chromosomes in elongated bacterial cells [13], the self-organization of DNA (or DNA-actin filament mixing systems) in cell confinement [14], and the drug delivery of spherical vesicles [15,16].
Polymer simulation is an important research method for studying the biosystem, such as the partitioning of chromatin in the nucleus [17], observing the circular DNA of viruses [18,19] and genome folding. Halverson et al. reviewed the research progress and reliability of modeling chromatin folding in eukaryotic cells using topological constraints on the topological constraints of polymers by comparing the latest theories related to simple chains (untangled and unlinked) or chain systems with the latest experiments related to genome folding. Due to the no-ends feature and topological constraints of ring chains, there is territorial segregation in a simulation system composed of 200 unconnected uncommitted uncombusted rings with a chain length of 1600 monomers [18]. This is consistent with experimental results for a partitioning phenomenon in the chromosomal region of fibroblast nuclei [20].
Phase separation is an important phenomenon in biological systems, such as membraneless organelles formed by liquid–liquid phase separation [21], the formation of the bacterial nucleoid due to the segregative phase separation [22] and compartmentalization of conventional nuclei [17]. Abnormal phase-separation mechanisms may even cause diseases, including neurodegenerative diseases, tumors, aging and infectious diseases, etc. Joyeux et al. performed a numerical simulation of the DNA chain inside the confining sphere using a coarse-grained method. They showed that the method has high sensitivity to the dissymmetry of DNA–DNA, DNA–crowder, crowder–crowder repulsive interactions, especially for dumbbells and octahedra crowders. Additionally, the demixing of the DNA coil and non-binding globular macromolecules present in the cytoplasm leads to nucleoid formation [22]. Falk et al. combined Hi-C analysis of inverted rod nuclei with microscopy and polymer simulations. They found that attractions between heterochromatic regions are crucial for establishing both compartmentalization and the concentric shells of pericentromeric heterochromatin, facultative heterochromatin and euchromatin in the inverted nucleus [17].
The semiflexible ring polymer (SRP) is an important model for simulating the study of circular biomolecules, such as circular chromosomes, eukaryotic mitochondrial DNA, circular plasmids and circular structure monosaccharides of E. coli. By introducing the stiffness of the chain, the ring chains stretch into a disk-like ring shape [23,24], similar to molecular structural characteristics. SRPs have been extensively studied theoretically and experimentally [13,22,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40]. The folding of SRPs in confinement and their interspersing with other SRPs are consistent with the observed biological systems phenomenon. Ostermeir et al. employed a simulation method to investigate the conformations of a single SRP confined in a sphere. By analyzing shape parameters calculated from average simulation data and theoretical analysis, they concluded that the construction of polymer rings causes buckling due to elasticity and entropy [25]. They also discovered and quantified a SRP’s conformational transition to an eight-shaped ring chain confined in a hard sphere [26]. Jeong et al. studied the morphology of SRPs when they were adsorbed on a topographically or chemically structured substrate surface using computational simulation. They found four equilibrium shapes: a round toroidal and a confined elongated shape, as well as two shapes containing bulges. Additionally, there was a cascade of transitions between elongated shapes [27]. Fritsche et al. found that different geometrical constraints can shape the spatial organization of SRP in confined conditions. Elongated, rod-like geometries reduce the number of chain overcrossings and induce ordering. Additionally, there is no preferred orientational axis in the case of spherical confinement. During the simulation process, SRP migrates from the center of the accessible space to the surrounding confined surface, forming a spool-like structure similar to DNA condensation in the capsid [28]. Experimentally, Gómez et al. used X-ray tomography to study the geometrical and topological features of disordered packages of rubber bands in a cylindrical container. For short SRPs, there was a liquid-like disordered structure with a short-range orientational order. However, for longer SRPs, the confinement caused folded configurations, and the bands were interpenetrated and entangled [36].
Biological systems are often complex and diverse, and this paper examines the kinetic behavior of semiflexible circular polymers confined in spheres of blending systems. This study helps us to understand the microscopic mechanisms of folding and movement of biomolecules in confined systems.
In this paper, we constructed a simulation model with a long SRP–short SRP blending system limited to the hard ball and studied the distribution of the long SRP in the binary mixture. We found an instance of phase separation between the long and short SRP, and the long ring is adhered to the wall.
This article consists of the following parts: In Section 2, the model and simulation details are provided. In Section 3, our results on the long ring polymer in spherically confined binary semiflexible ring polymer mixtures are given, and in Section 4, the conclusion is presented. Vesicle-restricted mixed macromolecules can deliver active macromolecules, such as DNA, peptides and some synthetic macromolecules. Our spherical shells can be regarded as vesicles. Furthermore, the simulation results in this paper can help us understand the physical mechanisms of interactions in vesicle-restricted mixed macromolecules.

2. Model and Method

During the simulation process, all ring chains, including long and short chains, are untangled and unconnected structures. We use the classical bead spring model to model the ring polymer chain, which is widely used in polymer simulations, and the simulation results of this model are consistent with the experimental results. Each ring chain is composed of L spherical monomers with a diameter of σ and a mass of m, and the length of the long ring in this article is fixed as Llong = 50, and the length of the short ring ranges from Lshort = 5 to Lshort = 40. The corresponding potential energy of each chain is given by
U = U bond + U angle + U LJ
where adjacent monomers on the chain are connected to each other by finite stretchable nonlinear elasticity (FENE) potentials [41]:
U bond ( r ) = K R 0 2 2 ln [ 1 ( r R 0 ) 2 ] , r < R 0
where r is the distance between two adjacent monomers on a ring chain. K = 30kBT/σ2 is the spring coefficient, and the finite ductility correlation parameter R0 = 1.5σ is used to avoid chain crossing.
To describe the rigid property of the semiflexible chains, we introduce the bond angle potential energy between adjacent bonds, as follows [42,43]:
U angle = K b e n d i n g [ 1 cos ( θ θ 0 ) ]
where θ is the angle between two adjacent bonds and Kbending is the bending energy. Additionally, θ0 is the equilibrium value of the angle, for different length of rings, θ0 = [π × (L − 2)]/L, such as for Llong, θ0 = [π × (50 − 2)]/50 = 0.96π. When Kbending is larger, the chain is more difficult to bend, that is, the chain is rigid. In this paper, the bending energy of the long ring is Kb-long = 100, and the bending energy of the short ring is Kb-short = 50.
To prevent overlap between all monomers, all bonded and non-bonded monomer interactions use the cutoff Lennard–Jones (LJ) potential:
U L J ( r ) = 4 ε [ ( σ r ) 12 ( σ r ) 6 + 1 4 ] , r 2 1 / 6 σ 0 , r > 2 1 / 6 σ
where r is the distance between two monomers, the truncation radius is rc = 21/6σ—that is, long ring–long ring, short ring–short ring, long ring–short ring, long ring–hard ball and short ring–hard ball interactions are all pure repulsion, and the interaction strength ε = 1.0kBT.
During the simulation, the binary semiflexible ring hybrid system is confined to a hard ball with a radius of R = 20σ. The spherical surface consists of LJ particles with the same diameter size σ and mass m as the polymer monomer. The velocity Verlet algorithm is used to integrate Newtonian equations of motion with a time step of Δt = 0.006τ0, where τ0 = (mσ2/kBT)1/2 is the inherent MD unit of time. The units of reduction kBT = 1, σ = 1, m = 1, τ0 = (mσ2/kBT)1/2 = 1 are units of energy, length, mass and time, respectively. The total particle number density of the binary mixed system is defined by ρ = (Llong ∗ Nlong + Lshort ∗ Nshort)/V, where Nlong and Nshort are the number of long rings and the number of short rings, Llong and Lshort are the long ring chain lengths and short ring lengths, respectively, and V is the volume of the hard ball. According to the results of a previous study, the necessary condition for phase separation in the long-short ring blending system is the high particle number density, and to prevent semiflexible ring deformation, the particle number density in this paper is fixed to ρ = 0.5σ−3. In this paper, Llong = 50, Nlong = 1 is fixed, Nshort can be calculated according to the number of different particles Lshort, and the specific values are given in Table 1 of Section 3.1. Unless otherwise stated, the relevant result is Llong = 50, Lshort = 1, Lshort = 10.
First, we place all ring chains completely randomly inside the hard ball. All rings remain untangled and unconnected at all times in the simulation with other semiflexible rings. Then, the system is balanced in the NVT ensemble, the total running time is 108 steps, where the first step running time Δt1 = 107 is used to ensure balance. After Δt1, data are recorded every 104 steps as a sample. Using a Nosé–Hoover thermostat, the reduction temperature is T* = 1.0 in ε/kB. All The simulations were performed in the open source molecular dynamics software package LAMMPS [44] (Sandia National Lab, Albuquerque, NM, USA).

3. Results and Discussion

3.1. Distribution of the Binary Mixture

In our system, the radius of the hard ball outside is fixed at R = 20. Since there is a strongly repulsive force for the distance between the monomer of the long SRP and a sphere surface less than d < 1, the distance between the monomer of the long SRP and the sphere center ranges of d is d = 0–19. The length of the long SRP is fixed at Llong = 50, which is commonly used to calculate long SRPs in simulation studies [23,24]. Additionally, the length of the short ring in the bulk involved in this article is in the range of Lshort = 5–40, and the corresponding number of short rings can be calculated, and the specific values are given in the following Table 1. Calculation details of the simulation parameters in Table 1 are given in Supplementary Information. A screenshot of the simulation system is shown in Figure 1. The snapshots are after long simulation steps and the system is stable. In general, the system is stable after 107 steps, and 107–109 steps are used for statistics, if we take 104 steps as a data point, there are about 99,000 points per-sample. Additionally, there are Nsample = 50 samples. All calculated results ensemble averaging over many conformations and all SRPs involved.
The rings are initially set as unknotted and unlinked. All chains cannot cross and are preserved during the whole simulation process. To examine the topological state of the system, we calculate the mean square radius of gyration Rg2 and the prolateness parameter p, defined, respectively, as
R g 2 = 1 N i = 1 N ( r i r c ) 2
Additionally, we analyze the shape of SRPs by calculating the prolateness parameters, which are defined from the radius of the gyration tensor:
T α β = 1 N i = 1 N ( r i α r α c m ) ( r i β r β c m )     ( α , β = x , y , z )
where N is the chain length, and α and β represent the Cartesian components. The vectors ri and ricm are the positions of the monomers and the center-of-mass of the SRPs, respectively. The prolateness parameter p is defined as
p = ( 2 λ 1 λ 2 λ 3 ) ( 2 λ 2 λ 1 λ 3 ) ( 2 λ 3 λ 1 λ 2 ) 2 ( λ 1 2 + λ 2 2 + λ 3 2 λ 1 λ 2 λ 1 λ 3 λ 2 λ 3 ) 3 2
where λ1, λ2, and λ3 are the three eigenvalues of the tensor with λ1 ≤ λ2 ≤ λ3. For perfectly oblate objects ( λ 1 < λ 2 = λ 3 ), the prolateness is p = −1. In addition, by taking the average of the stable structures of SRP, we can calculate the average prolateness of the ring chains; the prolateness parameter for a 50-bond ring polymer chain with Kb-long = 100 is <p> = −0.8, and the bracket represents the ensemble averaging over all conformations. For our system, using 107 steps to 108 steps, every 10,000 steps to sample points, and using 50 samples to average, we used a total of 450,000 sample points. Therefore, SRPs can be regarded as disk-like structures and the topological properties of SRPs are completely different from those of flexible ones. In our simulation, < R g > = 7.9 for the long SRP, which is close to the radius of the 50-bond circle (r = Llong/2/π = 7.958). Here, b is the bond length, and our simulation result is b = 0.97σ = 0.97, which is widely used in polymer simulation [23,24,38,43]. Thus, the SRPs with Kb,ring = 100 are completely extended and have stiff structures.
As shown in the snapshot in Figure 2, the system reaches equilibrium; the long SRP gradually detaches from the bulk. To avoid the effect of position placement of the long SRP at the initial moment, the long SRP is placed near the center of the ball, as shown in Figure 2. We recorded a typical simulation process in which the rings are distributed at different timesteps, screenshots are shown in Figure 2. As the simulation time increases, the long SRP gradually detaches from the bulk, eventually adheres to the wall distribution and reaches equilibrium.
After reaching equilibrium, we calculated the distribution of long ring beads within the spherical shell, demonstrating that the ring chain is distributed on the inner surface.
To verify the existence of separation structures in a binary hybrid system of short SRPs and a long SRP, we calculated the probability of the occurrence of the long SRP being located in different shells as follows:
p ( d ) = M l a y e r ( d ) L l o n g
where Mlayer is the monomer number located in the d-th shell, as shown in Figure 3a. As shown in Figure 3b, for Llong > 20, the p(d) curve has no peak and even shows a slight decrease. The distribution of the long SRP is almost uniform throughout the sphere, with a small peak near the wall (d ≈ 19). As Llong decreases to Llong = 5 and Llong = 10, there is a significant peak near the wall; for this condition, the long SRP adheres to the inner surface of the hard sphere, and there is a structure of long-short ring separation in the binary hybrid system with the long-short SRP binary mixture. For Llong > 20 binary hybrid systems, the long SRP and the short SRPs are fully mixed and evenly distributed throughout the confined space, while for the system where Llong is reduced to 5 and 10, the long and short rings are almost completely separated. Experimentally, there exists phase separation for large and small blended particles confined within the sphere.
We also calculated the average distance between the monomer of the long SRP and the sphere center <dis>, defined as:
< d i s > = i = 1 N s a m p l e j = 1 L l o n g d i s i j N s a m p l e L l o n g
Figure 3c shows the relationship between the average distance of the long-ring monomer from the center of the sphere <dis> with Lshort under different Lshort. <dis> decreases as Lshort increases. When Llong = 5, 10, <dis> the value reaches above 18, as shown by the dotted line. The results show that most monomers of the long SRP are distributed in the layer closest to the inner surface of the hard sphere, while the short semiflexible ring constitutes the anisotropic matrix inside the sphere.
In fact, in the case of spherical restriction, there is entropy attraction between the long SRP and the spherical shell. This strong entropy attraction causes the long SRP to move towards the spherical shell, eventually forming a separation structure of the long and short rings in the binary semiflexible ring hybrid system.
Yodhet et al. [45,46] studied the separation behavior of binary colloids of different sizes limited to rigid vesicles, and the results showed that this separation structure was formed at a sufficiently high colloid volume fraction. Our system is similar to a binary hardball hybrid system. In a binary hard-ball hybrid system, two larger balls close to each other do not change their interaction energy but increase the space reachable by other particles, causing the system to increase entropy, decrease free energy, and satisfy (3/2) α′∙φs∙kBT [46,47]. Here, α′ is the ratio of the radius of the large ball to the radius of the ball (RL/RS), and φs is the volume fraction of the ball. Although our binary semiflexible ring system is more complex than the binary hard ball system, there is a similar emptying interaction between the long SRP and the wall, and this entropy attraction also depends on the chain length ratio of the long SRP to the short semiflexible ring, i.e., α = Llong/Lshort. The results showed that the shorter the Lshort (the larger the α), the stronger the attraction between the long SRP and the wall. If α ≤ 2.5 (i.e., Lshort ≥ 20), no significant separation structure can be observed in a binary semiflexible ring hybrid system for all particle densities and long SRP bending energies.

3.2. Adherent Moving of the Long SRP

We found adherent moving of the long SRP in spherically confined binary semiflexible ring polymer mixtures. Figure 4 shows all the regions that the monomers of the long SRP can cover at different times. The time interval is 104 steps. Additionally, the total number of trajectories is 4000. As shown in Figure 4a, for the Llong = 30 simulation system, beads of the long SRP have a high range of motion throughout the space inside the sphere. As can be seen in Figure 4b, for Llong = 10, the range of activity of beads of the long SRP is limited to the inner shell area of the adherent, and the ball formed by the trajectory points is hollow. Therefore, the long SRP-shaped polymer walks against the wall in the blending system.
We further calculated the time-dependent distance S(t) and the average velocity v of long SRPs after equilibrium as follows:
S ( t ) = < 0 t s ( t ) d t >
The distance S(t) we calculated is the total length of the trajectory of the center of mass of the long SRP. For example, when the center of mass of the long SRP moves from t = 0 to position t, St(t) is the length of the trajectory that the long SRP passes through during this motion. The distance St(t) is the sum of St(t) = d1 + d2 + d3 + … + dt. To investigate the speed of the chain motion in the restricted system, we computed the scalar of the distance and calculated the average velocity v based on the slope of S(t), as shown in Figure 5 and Figure 6. Here, we discuss the distance S(t) instead of the displacement r ( t ) and the bracket <…> represents ensemble averaging over many conformations. Additionally, velocity V is the ratio of the distance ∆S and the time ∆t, the slope of S(t). As Lshort < 20, nearly complete segregation occurs between long and short SRPs, where long SRPs are attached to the sphere surface and only move along the sphere surface. When ρ = 0.5σ−3, the average velocity of long SRPs v decreases with the increase in Lshort, as shown in Figure 6. As Lshort < 20, long SRPs diffuse in a wall-attached shell region with phase separation. With the increase in Lshort, i.e., Lshort = 20 and 30, long SRPs with low content are immersed in the matrix of short SRPs. Additionally, the inside bulk of the restricted system moves faster than the inner wall layer, and the motion velocity of long SRPs is dominated by the short SRPs matrix and significantly slows down as Lshort increases.
For polymer-restricted systems, there is a decrease in system density at the inner layer near the confined surface. There are density oscillations (“layering”) near the inner wall of the sphere, well known for dense fluids near hard walls [48,49,50,51,52,53,54,55,56]. As shown in Figure 7, we calculate reduced polymer number density ρmonomerbulk, where ρmonomer is the number density of different shell and ρbulk is the number density of the inside bulk. This shows that the oscillation of monomer densities fades away for d < 10 by forming layers. The ordered monomers in the nearest inner layer (d = 19) of the sphere help the monomers to be ordered in the second layer (d = 18), and in turn, the third (d = 17) and fourth layers (d = 16) are also formed in a layer-by-layer manner leading to the density oscillations. Additionally, ordering becomes weak as d decreases. Such a layering structure is observed both in polymer systems and in the various confined liquid systems, such as ionic liquids in graphene double layers [57].
Due to the number of possible conformations of polymer chains being limited by an impenetrable object, polymer chains tend to move away from the wall. The monomer density in the closest adherent layer d ≥ 19 is relatively low (region colored in red). Additionally, as the distance from the wall increases, monomers significantly aggregate in the interface layer (d ≈ 18.8), and then the monomer density gradually decreases. Inside part of the sphere with d < 10, the density is stable, ρmonomer and ρbulk are almost equal, and ρmonomerbulk≈1. Chains are affected by the wall’s limitation and the system’s uniform density. That is, the lower density of the polymer chain nearest the wall inevitably increases the chain density of the central area of the bulk, and the increase in the density of the central area of the bulk will in turn produce a force that pushes the chains towards the direction of the wall. Therefore, a lower density layer of is formed nearest the wall surface. The movement speed in the close layer is faster than the internal bulk. When the long ring separates from the blending system, it locates in this adherent layer, where the density is lower and the chain movement speed is faster. Moreover, there is a peak in density immediately left to the minimum, as shown in Figure 7. Additionally, because of the existence of highly dense areas, once the long SRP is distributed in the inner surface layer d≈19, the high-density layer blocks the long SRP from entering the interior. It is difficult for the long SRP to enter the interior d < 19 and then move throughout the inner surface. When Lshort ≥ 20, the motion is slower as the ring length increases, and the long ring cannot be separated into the adherent layer and stay in the bulk. Because of the increased length of short chains, the ring movement slows down naturally. Moreover, the long ring cannot be separated from the inside bulk. The long ring distributes both in the interface layer and in the bulk. Combining the above two factors, the average velocity of the long ring v drops significantly for Lshort ≥ 20.
All SRPs and the hard ball have repulsive volume, i.e., where other individuals cannot access it. As shown in Figure 8, the shades of red and blue represent the excluded regions of the long SRP and the wall, respectively. When the number of small rings increases, corresponding to the ρ = 0.5σ−3 in this article, the large ring is required to adhere to the wall; that is, the exclusion volume of the large ring coincides with the exclusion volume of the spherical monomer, giving the small ring more room for activity so that the entropy of the small ring increases, and the system is oriented to the separation of the two phases of the large ring and the small ring.
When the long SRP moves to the wall, the excluded regions of the long SRP and the wall (shades of red and blue) overlap. The short SRPs’ entropy therefore increases. The reduction in free energy produces an “entropic force” that pushes the long spheres to the inside surface. When the long SRP is moved to a flat wall, moreover, the overlap volume and the free-energy loss are approximately doubled [58]. This phenomenon of entropy leading to the separation of mixed phases is often found in other systems. In binary hard-sphere mixtures, these effects are known to drive the crystallization of large spheres in the bulk [13,59,60] and on flat surfaces.
Furthermore, the shape of the wall can lead to entropic forces in a specific direction along the wall. For example, the larger spheres are locally repelled from an edge cut into the wall [61] and attracted to a corner (i.e., where the “wall” meets the “floor”) [62]. If the wall has a constantly changing radius of curvature, these forces are predicted to act everywhere along it [61]. As shown in Figure 8, when the long SRP is near the wall, the overlap volume depends on the wall’s curvature radius. The large sphere therefore moves in the direction of increasing curvature to minimize the small SRPs’ excluded volume. For our system, the outer surface is spherical, i.e., the curvature of the inner surface is equal everywhere, so the long ring moves in a random direction along the inner surface. The long SRP moves along the wall to maximize the size of the overlap region, as indicated by the arrow.

4. Conclusions

We simulated a long semiflexible ring–short semiflexible ring blending system confined in a sphere, and by changing the length of the short ring, the phase separation results were obtained when the short ring length Lshort < 20. The phase separation, related to entropy, coincides with the exclusion volume of the wall when the large ring is attached to the wall, providing more activity space for the small ring, and the entropy of the small ring increases. Additionally, after the large ring is attached to the wall, it moves randomly along the wall in a random direction. By calculating the distance and speed of the long ring, after comparison, we find that as the length of the short ring chain increases, the movement of the long ring gradually slows down, and when there is no more phase separation, the movement of the long ring slows down rapidly. As phase separation disappears, long rings are distributed within the system rather than against the wall. In the restricted sphere, due to the relatively low particle density of the adherent layer, the range of motion space is large, and the internal motion is slower than the motion of the adherent layer. As a result, the movement of the long ring slows down rapidly. Our research helps to understand the separation of macromolecular blending systems in confined biological systems, such as chromosomal partitioning and the self-organization of DNA–actin filament mixing systems with cell confinement.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biophysica2040044/s1.

Author Contributions

Conceptualization, X.Z.; methodology, X.Z.; resources, X.Z.; data curation, X.Z.; W.W.; supervision, X.Z.; funding acquisition, X.Z.; writing, original draft preparation, X.Z.; writing, review and editing, X.Z.; W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Fundamental Research Funds for the Central Universities (grant no. 3072021CFT2508).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author, [Xiaolin Zhou], upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, F.B.; Japaridze, A.; Zheng, X.; Wiktor, J.; Kerssemakers, J.W.J.; Dekker, C. Direct imaging of the circular chromosome in a live bacterium. Nat. Commun. 2019, 10, 2194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Shao, Y.Y.; Lu, N.; Cai, C.; Zhou, F.; Wang, S.S.; Zhao, Z.H.; Zhao, G.P.; Zhou, J.Q.; Xue, X.L.; Qin, Z.J. A single circular chromosome yeast. Cell Res. 2019, 29, 87–89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jacob, F.; Wollman, E.L. Genetic and physical determinations of chromosomal segments in Escherichia coli. Sym. Soc. Exp. Biol. 1958, 12, 75–92. [Google Scholar]
  4. Freifelder, D.; Kleinschmidt, A.K.; Sinsheimer, R.L. Electron microscopy of single-stranded DNA: Circularity of DNA of bacteriophage ϕX174. Science 1964, 146, 254–255. [Google Scholar] [CrossRef] [PubMed]
  5. Wasserman, S.A.; Cozzarelli, N.R. Biochemical topology: Applications to DNA recombination and replication. Science 1986, 232, 951–960. [Google Scholar] [CrossRef] [PubMed]
  6. Kundukad, B.; Yan, J.; Doyle, P.S. Efect of YOYO-1 on the mechanical properties of DNA. Soft Matter 2014, 10, 9721–9728. [Google Scholar] [CrossRef]
  7. Ireland, D.C.; Clark, R.J.; Daly, N.L.; Craik, D.J. Isolation, sequencing, and structure–activity relationships of cyclotides. J. Nat. Prod. 2010, 73, 1610–1622. [Google Scholar] [CrossRef] [PubMed]
  8. Parisi, D.; Truzzolillo, D.; Deepak, V.D. Transition from confined to bulk dynamics in symmetric star–linear polymer mixtures. Macromolecules 2019, 52, 5872–5883. [Google Scholar] [CrossRef] [Green Version]
  9. Li, S.J.; Ding, M.M.; Shi, T.F. Spatial distribution of entanglements and dynamics in polymer films confined by smooth walls. Polymer 2019, 172, 365–371. [Google Scholar] [CrossRef]
  10. Song, J.; Kahraman, R.; Collinson, D.W. Temperature effects on the nanoindentation characterization of stiffness gradients in confined polymers. Soft Matter 2019, 15, 359–370. [Google Scholar] [CrossRef]
  11. Raffaele, P.; Alessio, D.; Mose, C. Influence of wall heterogeneity on nanoscopically confined polymers. Phys. Chem. Chem. Phys. 2019, 21, 772–779. [Google Scholar]
  12. Zuo, B.; Zhou, H.; Davis, M.J.B. Effect of local chain conformation in adsorbed nanolayers on confined polymer molecular mobility. Phys. Rev. Lett. 2019, 122, 217801. [Google Scholar] [CrossRef] [PubMed]
  13. Jun, S.; Mulder, B. Entropy-driven spatial organization of highly confined polymers: Lessons for the bacterial chromosome. Proc. Natl. Acad. Sci. USA 2006, 103, 12388–12393. [Google Scholar] [CrossRef] [Green Version]
  14. e Silva, M.S.; Alvarado, J.; Nguyen, J.; Georgoulia, N.; Mulder, B.; Koenderink, G. Self-organized patterns of actin filaments in cell-sized confinement. Soft Matter 2011, 7, 10631–10641. [Google Scholar] [CrossRef]
  15. Qi, W.; Yan, X.; Fei, J.; Wang, A.; Cui, Y.; Li, J. Triggered release of insulin from glucose-sensitive enzyme multilayer shells. Biomaterials 2009, 30, 2799–2806. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Y.; Wang, L.; Li, B.; Cheng, Y.; Zhou, D.; Chen, X.; Jing, X.; Huang, Y. Compact vesicles self-assembled from binary graft copolymers with high hydrophilic fraction for potential drug/protein delivery. ACS Macro Lett. 2017, 6, 1186–1190. [Google Scholar] [CrossRef] [PubMed]
  17. Falk, M.; Feodorova, Y.; Naumova, N.; Imakaev, M.; Lajoie, B.; Leonhardt, H.; Joffe, B.; Dekker, J.; Fudenberg, G.; Solovei, I.; et al. Heterochromatin drives compartmentalization of inverted and conventional nuclei. Nature 2019, 510, 395. [Google Scholar] [CrossRef] [PubMed]
  18. Halverson, J.D.; Smrek, J.; Kremer, K.; Grosberg, A.Y. From a melt of rings to chromosome territories: The role of topological constraints in genome folding. Rep. Prog. Phys. 2014, 77, 022601. [Google Scholar] [CrossRef] [Green Version]
  19. Marenduzzo, D.; Micheletti, C.; Orlandini, E. Biopolymer organization upon confinement. J. Phys. Condens. Matter 2010, 22, 283102. [Google Scholar] [CrossRef] [PubMed]
  20. Bolzer, A.; Kreth, G.; Solovei, I.; Koehler, D.; Saracoglu, K.; Fauth, C.; Eils, S.R.; Cremer, C.; Speicher, M.R.; Cremer, T. Three-dimensional maps of all chromosomes in human male fibroblast nuclei and prometaphase rosettes. PLoS Biol. 2005, 3, 157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Brangwynne, C.P.; Eckmann, C.R.; Courson, D.S.; Rybarska, A.; Hoege, C.; Gharakhani, J.; Julicher, F.; Hyman, A.A. Germline P Granules Are Liquid Droplets That Localize by Controlled Dissolution/Condensation. Science 2009, 324, 1729–1732. [Google Scholar] [CrossRef]
  22. Joyeux, M. A segregative phase separation scenario of the formation of the bacterial nucleoid. Soft Matter 2018, 14, 7368–7381. [Google Scholar] [CrossRef] [PubMed]
  23. Bernabei, M.; Bacova, P.; Moreno, A.J. Fluids of semiflexible ring polymers: Effective potentials and clustering. Soft Matter 2013, 9, 1287–1300. [Google Scholar] [CrossRef] [Green Version]
  24. Poier, P.; Bacova, P.; Moreno, A.J.; Likos, C.N.; Blaak, R. Anisotropic effective interactions and stack formation in mixtures of semiflexible ring polymers. Soft Matter 2016, 12, 4805–4820. [Google Scholar] [CrossRef] [PubMed]
  25. Ostermeir, K.; Alim, K.; Frey, E. Buckling of stiff polymer rings in weak spherical confinement. Phys. Rev. E 2010, 81, 061802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Ostermeir, K.; Alim, K.; Frey, E. Confinement induces conformational transition of semiflexible polymer rings to figure eight form. Soft Matter 2010, 6, 3467–3471. [Google Scholar] [CrossRef] [Green Version]
  27. Jeong, S.; Cho, S.; Kim, J.M.; Baig, C. Interfacial Molecular structure and dynamics of confined ring polymer melts under shear Flow. Macromolecules 2018, 51, 4670–4677. [Google Scholar] [CrossRef]
  28. Fritsche, M.; Heermann, D.W. Confinement driven spatial organization of semiflexible ring polymers: Implications for biopolymer packaging. Soft Matter 2011, 7, 6906–6913. [Google Scholar] [CrossRef]
  29. Tubiana, L.; Orlandini, E.; Micheletti, C. Multiscale Entanglement in ring polymers under spherical confinement. Phys. Rev. Lett. 2011, 107, 188302. [Google Scholar] [CrossRef] [PubMed]
  30. Micheletti, C.; Marenduzzo, D.; Orlandini, E.; Sumners, D.W. Knotting of random ring polymers in confined spaces. J. Chem. Phys. 2006, 124, 064903. [Google Scholar] [CrossRef] [Green Version]
  31. Ding, M.M.; Duan, X.Z.; Lu, Y.Y.; Shi, T.F. Flow-induced ring polymer translocation through nanopores. Macromolecules 2015, 48, 6002–6007. [Google Scholar] [CrossRef]
  32. Gutjahr, P.; Lipowsky, R.; Kierfeld, J. Semiflflexible polymer rings on topographically and chemically structured Surfaces. Soft Matter 2010, 6, 5461–5475. [Google Scholar] [CrossRef]
  33. Polson, J.M.; Kerry, D.R.M. Segregation of polymers under cylindrical confinement: Effects of polymer topology and crowding. Soft Matter 2018, 14, 6360–6373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Minton, A.P. The influence of macromolecular crowding and macromolecular confinement on biochemical reactions in physiological media. J. Biol. Chem. 2001, 276, 10577–10580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Gaitho, F.M.; Tsige, M.; Mola, G.T.; Pellicane, G. Surface segregation of cyclic chains in binary melts of thin polymer films: The Influence of Constituent Concentration. Polymers 2018, 10, 324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Gómez, L.R.; García, N.A.; Pöschel, T. Packing structure of semiflexible rings. Proc. Natl. Acad. Sci. USA 2020, 117, 3382–3387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Zhou, Y.C.; Hsiao, K.W.; Regan, K.E.; Kong, D.J.; McKenna, G.B.; Robertson-Anderson, R.M.; Schroeder, C.M. Effect of molecular architecture on ring polymer dynamics in semidilute linear polymer solutions. Nat. Commun. 2019, 10, 1753. [Google Scholar] [CrossRef] [Green Version]
  38. Zhou, X.L.; Li, K.; Guo, F.C.; Zhang, L.X. Ordered aggregation structures of semiflexible ring polymers in ring-linear blends. Polymer 2019, 175, 129–136. [Google Scholar] [CrossRef]
  39. Witz, G.; Rechendorff, K.; Adamcik, J.; Dietler, G. Conformation of ring polymers in 2D constrained environments. Phys. Rev. Lett. 2011, 106, 248301. [Google Scholar] [CrossRef] [Green Version]
  40. Zhou, X.L.; Liu, L.; Chen, J.M.; Zhang, L.X. Unusual conformations of semiflexible ring polymers confined in two parallel surfaces. Polymer 2018, 157, 180–189. [Google Scholar] [CrossRef]
  41. Grest, G.S.; Kremer, K. Molecular dynamics simulation for polymers in the presence of a heat bath. Phys. Rev. A 1986, 33, 3628–3631. [Google Scholar] [CrossRef] [PubMed]
  42. Cifra, P. Channel confinement of flexible and semiflexible macromolecules. J. Chem. Phys. 2009, 131, 224903. [Google Scholar] [CrossRef] [PubMed]
  43. Narros, A.; Moreno, A.J.; Likos, C.N. Influence of topology on effective potentials: Coarse-graining ring polymers. Soft Matter 2010, 6, 2435–2441. [Google Scholar] [CrossRef]
  44. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  45. Dinsmore, A.D.; Yodh, A.G. Phase diagrams of nearly-hard-sphere binary colloids. Phys. Rev. E 1995, 52, 4045–4057. [Google Scholar] [CrossRef]
  46. Dinsmore, A.D.; Wong, D.T.; Nelson, P.; Yodh, A.G. Hard spheres in vesicles: Curvature-induced forces and particle-induced curvature. Phys. Rev. Lett. 1998, 80, 409–412. [Google Scholar] [CrossRef] [Green Version]
  47. Asakura, S.; Oosawa, F. Surface tension of high-polymer solutions. J. Chem. Phys. 1954, 22, 1255. [Google Scholar] [CrossRef]
  48. Pais, A.A.C.C.; Miguel, M.G.; Linse, P.; Lindman, B. Polyelectrolytes confined to spherical cavities. J. Chem. Phys. 2002, 117, 1385. [Google Scholar] [CrossRef]
  49. Milchev, A.; Egorov, S.A.; Nikoubashman, A.; Binder, K. Conformations and orientational ordering of semiflexible polymers in spherical confinement. J. Chem. Phys. 2017, 146, 194907. [Google Scholar] [CrossRef] [Green Version]
  50. Sarabadani, J.; Milchev, A.; Vilgis, T.A. Structure and dynamics of polymer melt confined between two solid surfaces: A molecular dynamics study. J. Chem. Phys. 2014, 141, 044907. [Google Scholar] [CrossRef]
  51. Milchev, A.; Egorov, S.A.; Nikoubashman, A.; Binder, K. Adsorption and structure formation of semiflflexible polymers on spherical surfaces. Polymer 2018, 145, 463–472. [Google Scholar] [CrossRef]
  52. Zhou, Y.X.; Milner, S.T. Short-time dynamics reveals Tg suppression in simulated polystyrene thin films. Macromolecules 2017, 50, 5599–5610. [Google Scholar] [CrossRef]
  53. Milchev, A.; Egorov, S.A.; Vega, D.A.; Binder, K.; Nikoubashman, A. Densely Packed semiflflexible macromolecules in a rigid spherical capsule. Macromolecules 2018, 51, 2002–2016. [Google Scholar] [CrossRef]
  54. Nikoubashman, A.; Vega, D.A.; Binder, K.; Milchev, A. Semiflexible polymers in spherical confinement: Bipolar orientational order versus tennis ball states. Phys. Rev. Lett. 2017, 118, 217803. [Google Scholar] [CrossRef]
  55. Vacatello, M. Monte Carlo Simulations of the interface between polymer melts and solids. effects of chain stiffness. Macromol. Theory Simul. 2001, 10, 187–195. [Google Scholar] [CrossRef]
  56. Lee, E.; Jung, Y.J. Segregated structures of ring polymer melts near the surface: A molecular dynamics simulation study. Soft Matter 2015, 11, 6018–6028. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. DeYoung, A.D.; Park, S.W.; Dhumal, N.R.; Shim, Y.; Jung, Y.; Kim, H.J. Graphene oxide supercapacitors: A computer simulation study. J. Phys. Chem. C 2014, 118, 18472–18480. [Google Scholar] [CrossRef]
  58. Morrison, G.; Thirumalai, D. Semiflexible chains in confined spaces. Phys. Rev. E 2009, 79, 011924. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Jones, R.L.; Kumar, S.K.; Ho, D.L.; Briber, R.M.; Russell, T.P. Chain conformation in ultrathin polymer films. Nature 1999, 400, 146. [Google Scholar] [CrossRef]
  60. Riggleman, R.A.; Yoshimoto, K.; Douglas, J.F.; de Pablo, J.J. Influence of confinement on the fragility of antiplasticized and pure polymer films. Phys. Rev. Lett. 2006, 97, 045502. [Google Scholar] [CrossRef]
  61. Micheletti, C.; Marenduzzo, D.; Orlandini, E. Polymers with spatial or topological constraints: Theoretical and computational results. Phys. Rep. 2011, 504, 1–73. [Google Scholar] [CrossRef] [Green Version]
  62. Kato, A.; Shindo, E.; Sakaue, T.; Tsuji, A.; Yoshikawa, K. Conformational transition of giant DNA in a confined space surrounded by a phospholipid membrane. Biophys. J. 2009, 97, 1678. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Screenshot of simulation system. The long SRP is highlighted in red, and short SRPs are shown in different color for clarity. Here, Lshort = 10.
Figure 1. Screenshot of simulation system. The long SRP is highlighted in red, and short SRPs are shown in different color for clarity. Here, Lshort = 10.
Biophysica 02 00044 g001
Figure 2. Typical snapshots of the long SRP confined in a spherical confinement at different time points. Here, Lshort = 10.
Figure 2. Typical snapshots of the long SRP confined in a spherical confinement at different time points. Here, Lshort = 10.
Biophysica 02 00044 g002
Figure 3. (a) Schematic of shell partitions and monomers belonging to the layer. (b) Probability distribution function p(d) of long SRP with different lengths of short rings Lshort. (c) Average distance <dis> between the center of a sphere and monomers of the long SRP as a function of Lshort.
Figure 3. (a) Schematic of shell partitions and monomers belonging to the layer. (b) Probability distribution function p(d) of long SRP with different lengths of short rings Lshort. (c) Average distance <dis> between the center of a sphere and monomers of the long SRP as a function of Lshort.
Biophysica 02 00044 g003
Figure 4. The trajectories of the long SRP in the cross section of confined sphere with different Lshort for Lshort = 30 (a), 10 (b) in the interval of 104 steps. Additionally, the total number of trajectories is 4000.
Figure 4. The trajectories of the long SRP in the cross section of confined sphere with different Lshort for Lshort = 30 (a), 10 (b) in the interval of 104 steps. Additionally, the total number of trajectories is 4000.
Biophysica 02 00044 g004
Figure 5. Distance S(t) of long SRPs as a function of time for different chain length Lshort.
Figure 5. Distance S(t) of long SRPs as a function of time for different chain length Lshort.
Biophysica 02 00044 g005
Figure 6. Speed of long SRPs v as a function of Lshort. The red rings in the illustration indicate possible locations at different times.
Figure 6. Speed of long SRPs v as a function of Lshort. The red rings in the illustration indicate possible locations at different times.
Biophysica 02 00044 g006
Figure 7. Reduced polymer monomer number density distribution of polymer blends in spherical confinement.
Figure 7. Reduced polymer monomer number density distribution of polymer blends in spherical confinement.
Biophysica 02 00044 g007
Figure 8. Description of depletion forces in a binary hard-sphere mixture. When the long SRP moves to the wall, the excluded regions of long SRP and the wall (shades of red and blue) overlap. The short SRPs’ entropy therefore increases. The long SRP moves along the wall to maximize the size of the overlap region, as indicated by the arrow.
Figure 8. Description of depletion forces in a binary hard-sphere mixture. When the long SRP moves to the wall, the excluded regions of long SRP and the wall (shades of red and blue) overlap. The short SRPs’ entropy therefore increases. The long SRP moves along the wall to maximize the size of the overlap region, as indicated by the arrow.
Biophysica 02 00044 g008
Table 1. The detailed simulation parameters covered in this paper.
Table 1. The detailed simulation parameters covered in this paper.
Radius of SphereLlongNlongLshrotNshortLlong ∗ Nlong + Lshort ∗ Nshort
205015286314,365
101432
20716
30477
40358
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, X.; Wei, W. Adherent Moving of Polymers in Spherical Confined Binary Semiflexible Ring Polymer Mixtures. Biophysica 2022, 2, 525-537. https://doi.org/10.3390/biophysica2040044

AMA Style

Zhou X, Wei W. Adherent Moving of Polymers in Spherical Confined Binary Semiflexible Ring Polymer Mixtures. Biophysica. 2022; 2(4):525-537. https://doi.org/10.3390/biophysica2040044

Chicago/Turabian Style

Zhou, Xiaolin, and Wan Wei. 2022. "Adherent Moving of Polymers in Spherical Confined Binary Semiflexible Ring Polymer Mixtures" Biophysica 2, no. 4: 525-537. https://doi.org/10.3390/biophysica2040044

APA Style

Zhou, X., & Wei, W. (2022). Adherent Moving of Polymers in Spherical Confined Binary Semiflexible Ring Polymer Mixtures. Biophysica, 2(4), 525-537. https://doi.org/10.3390/biophysica2040044

Article Metrics

Back to TopTop