Translocation of a Polymer through a Crowded Channel under Electrical Force

The translocation of a polymer chain through a crowded cylindrical channel is studied using the Langevin dynamics simulations. The influences of the field strength F, the chain length N, and the crowding extent ρ on the translocation time are evaluated, respectively. Scaling relation τ ~ F−α is observed. With the crowding extent ρ increasing, the scaling exponent α becomes large. It is found that, for noncrowded channel, translocation probability drops when the field strength becomes large. However, for high-crowded channel, it is the opposite. Moreover, the translocation time and the average translocation time for all segments both have exponential growth with the crowding extent. The investigation of shape factor 〈δ〉 shows maximum value with increasing of the number of segments outside s. At last, the number of segments inside channel Nin in the process of translocation is calculated and a peak is observed. All the information from the study may benefit protein translocation.

Kasianowicz et al. [12] demonstrated that single-stranded DNA and RNA molecules can be driven through the waterfilled -hemolysin channel under electric field. And the passage of each molecule is signaled by a blockade in the channel current. Improving this technique, the nucleotide sequence of DNA or RNA can be read off. In addition, solidstate nanopores have been used instead of the -hemolysin channel in other experiments [13][14][15][16][17]. The scaling relationship of translocation time with chain length is ∝ 1.27 , which is different from that obtained by Kasianowicz et al., ∝ .
In actual fact, the effect of crowding on the translocation dynamics should not have been neglected, because in the cell cytoplasm crowding due to macromolecular species and structural obstacles can be as high as 50% by volume. Moreover, in biological processes, protein-protein interactions are fundamentally important [37]. And this crowding environment has considerable influence on reaction rates, protein folding rates, and equilibria in vivo.
In this paper, we investigate the dynamics of polymer translocation through a crowded channel. And the extent of crowding is considered. Our model and the simulation technique are described in Section 2. Simulation results and discussion are presented in Section 3 and Section 4 is a conclusion.

Method of Calculation
In the simulations, we use the Langevin equation to study the Brownian motion of particles where the equation of motion for each bead at position is described as where is the bead mass, is the friction coefficient, denotes the external force due to the applied voltage represented by =̂, and ( ) is the random force which satisfies the fluctuation-dissipation relation [38] ⟨ ( )⟩ = 0, The total interaction is as follows: The finitely extensible nonlinear elastic (FENE) [39] spring potential interaction between two successive beads is where is the spring constant, 0 is the maximum allowed separation between connected monomers, and here is the distance between consecutive monomers. And a repulsive Lennard-Jones (LJ) potential is applied between all bead pairs considering both excluded volume and Van der Waals interactions between beads: Here, is the diameter of a bead, is the depth of the potential, and is the distance between two beads. The model of translocation is illustrated in Figure 1. A voltage is applied across the pore. The wall is formed by columns of stationary particles. Between the bead-wall particle pairs and bead-pore particle pairs, there exists the same short range repulsive LJ interaction as described in  Figure 1: A schematic illustration of a polymer moving through a crowded channel. The channel length is . The channel radius is . The external electric force is applied in the -axis direction inside the channel. The black spherical obstacles represent crowding environments. And the diameter of the obstacle is . (5). Crowding is modeled by randomly distributed spherical obstacles. The diameter of the obstacles is . And the interactions between obstacles and polymers or other particles are all described as (5). The volume of an obstacle is = (4/3) ( /2) 2 , and the volume of the channel is = 2 ⋅ . We can describe the density of the obstacles in the channel = / . Here is the number of obstacles distributed inside the channel. In this work, the diameter and the LJ interaction strength fit the length units and the system energy. The time scale is LJ = ( 2 / ) 1/2 . The parameters are = 1, 0 = 2 , = 7 , = 0.7, and = 1.2 [30]. In the simulations, = 5, = 2, and = 1. Firstly, the first monomer of the polymer is placed in the entrance of the channel, and the remaining monomers are to obtain an equilibrium configuration undergoing thermal collisions described by the Langevin thermostat. Then, the first monomer is released, and under the external electrical force the polymer begins to cross through the channel. The translocation time is defined as the time interval between the first monomer in the channel and the last monomer out of the channel. 2000 independent runs are averaged in our simulations.

Translocation Time.
As a result of external electric force, the polymer can overcome the entropic barrier due to the loss of the number of the configurations during translocation. It is obvious that the translocation time decreases with increasing electric force. In Figure 2, we fit the translocation time on the external force and find the scaling relation ∼ − .
The exponent is 0.482 ± 0.002, 0.48 ± 0.01, and 0.477 ± 0.006 without crowding for = 128, 64, and 32, respectively. The same scaling behavior has been found using DPD simulation [40]. And the exponent is 0.48 ± 0.01, which is completely in accord with our results. However, no crowding environment is discussed in that work. From the figure, it is certainly found that as the crowding extent increases, the translocation time increases under the same channel and chain. That is because the more crowded the channel, the larger the entropic barrier existing during translocation. However, as the crowding extent increases, the scaling exponent also increases. The value of is 0.593 ± 0.003, 0.618 ± 0.009, and 0.597 ± 0.008 under = 0.4 for = 128, 64, and 32. This is due to the increased entropic resistance as increases. For = 0.6, the exponent apparently becomes larger. For = 128, 64, and 32, the value of is 0.80 ± 0.01, 0.78 ± 0.01, and 0.74 ± 0.02, respectively. It is shown that the dependence of the translocation time on the external force is stronger for longer chain under = 0.6, while the scaling exponent for different chain length is nearly unobvious under = 0 and = 0.4. That means the influence of chain length to the scaling relation only happens under large crowded channel.
In Figure 3, we show the translocation probability of successful events as a function of the external electric force for different crowding extent and chain length . First, when = 0, the translocation probability decreases with increasing translocation force . For short chain = 64 and 32, the change is slow. However, for = 128, the translocation probability is above 80% under = 10 and 15. Because, under weak electric force, the translocation velocity is relatively slow, the monomers already translocated can diffuse and let the rest of chain through the channel. However, the translocation probability quickly drops to 60% and 20% for = 20 and 25. And it is only 10% when = 30. This is due to the fact that, with increasing , the monomers inside channel will get larger acceleration speed and exit crowding happens easily. That causes the translocation probability drops. Above all, it shows that, for = 0, low electric field leads to more successful translocation. Considering the crowded extent of the channel, = 0.4 and 0.6 are both investigated. It shows that, for = 32 and 64, the translocation probability changes little with increasing. The value is around 30% and 60%, respectively. However, for longer chain, = 128, when = 10, the translocation probability is 17%, and it increases to 40% when changes to 15. The phenomenon is more obvious for = 0.6. We focus on the results of = 128. The translocation probability is only 3% under = 10. With increasing , translocation probability increases first and then approaches saturation when > 20. It is quite different from that of = 0. It can be well understood that, with increasing , entropic barrier for the chain translocation becomes greater; thus, the chain will be dragged back out of the channel if the electric field is weak, which led to the low success of translocation probability. On the contrary, with the electric force increasing, the chain's backward motion out of the channel will effectively be prohibited, and the translocation probability increases. Figure 4 shows the translocation time as a function of the crowding extent . We find that translocation time increases with increasing . It keeps appropriate exponential growth. As is shown, the translocation time of = 10 is higher than that of = 20. For the same length of polymer, with crowding extent increasing, the interval between = 10 and = 20 becomes larger.    decreases. Before = 0.7, the average translocation time of segments of = 128 is longer than that of = 64. When = 0.7, the average time reaches the same value. As = 0.8, the average translocation time of segments of = 64 is longer than that of = 128. The following can be interpreted: when the channel is much crowded, with the addition of weak electric field = 10, the longer chain is easier to adjust its conformation outside to let the rest of the monomers pass through the channel. Figure 6 shows the translocation time of bead ( ) for = 128 and crowding extent = 0, 0.4, and 0.6. For = 0.6, the translocation time of = 1 is the largest. The reason is that the crowding channel will prevent the chain entering into it, so a high entropy barrier exists for the first several monomers. This can interpret why the translocation probability of = 0.6 is smaller than = 0 in Figure 3. Behind quick decreasing, ( ) decreases slowly and approaches saturation. For the last several monomers, a short peak is shown. This is because the channel exit will be crowded with translocated monomers. For less crowding channel = 0.4, almost a similar trend is investigated. However, the saturation stage is longer. When the chain goes through the noncrowded channel, the translocation time is almost unchanged till = 120. It seems that the translocation is easier than crowded channel. Also a peak is shown at the end of the curve. The reason is the same as above.

Chain Size and Conformations.
The instantaneous shape of an individual configuration may be described by several ratios based on the principal components 2 1 ≤ 2 2 ≤ 2 3 of 2 = 2 1 + 2 2 + 2 3 , that is, the orthogonal components of the squared radius of gyration taken along the principal axes of inertia [41,42]. ⟨ ⟩ is obtained by combining the reduced components of 2 to a single quantity that varies between 0 (sphere) and 1 (rod) [43,44] ⟨ ⟩ = 1 − 3 ⟨ In Figure 7 we show average shape factor ⟨ ⟩ versus the number of segments out of channel . Different crowding extent and electric force are both considered. As shown in Figure 7, with the increasing of , ⟨ ⟩ shows a maximum with increasing . It represents that the shape of chain is most spread as ⟨ ⟩ reaches a maximum value. For = 0, the value of corresponding to the maximum of ⟨ ⟩ is larger than that of = 0.4 and 0.6. That means the incompact conformation in the translocation for = 0 appears later than crowded channel. Another finding is when ⟨ ⟩ goes through    the maximum, it drops more quickly for noncrowded channel than crowded channel.
The number of segments inside channel in versus the number of segments out of channel has been shown in Figure 8 for = 128, = 0, 0.4, and 0.6, and = 15 and 25. There is a peak that existed in each curve at = 120∼125.
When there is no crowded particle in the channel, it is shown that in decreases slightly and has a long plateau till = 120. It seems that the translocation through noncrowded channel is smooth. At last, in increases to a peak value and then it drops. The peak implies that the translocation velocity decreases at the last lap of translocation. It is due to the fact that translocated segments will cause the exit to be crowded. However, if the channel is crowded with many particles occupied, here = 0.4 and 0.6 are both considered. The value of in decreases relatively greater than = 0. It can be explained that the crowding particles inside of the channel will push the chain out of the channel. It leads to the translocation velocity increasing; thus in decreases. Another finding is that the peak of the curve for electric force = 15 is less clear than that for = 25. This is because larger electric force will accelerate the crowding nearby the exit. And some segment will enter into the channel to ensure a successful translocation.

Conclusions
Using the Langevin dynamics simulations, we investigate the translocation of a polymer chain through a crowded cylindrical channel. We first observe a scaling relation ∼ − between the translocation times on the external force . The scaling exponent increases as the crowding extent becomes larger. As an example, for = 128, the scaling exponent is 0.482 ± 0.002, 0.593 ± 0.003, and 0.80 ± 0.01 for = 0, 0.4, and 0.6, respectively. This shows that the more crowded the channel, the much larger the entropic barrier during translocation. The translocation probability under different field strength is also investigated. The translocation probability decreases when the field strength becomes large for noncrowded channel. However, for high-crowded channel, it is the opposite. At the same time, we find that the translocation time has an exponential growth with the crowding extent . Average translocation time for all segments as a function of the crowding extent has a similar law. Considering the translocation time of each segment under different crowding extent, the peak at the end of translocation shows the crowding of translocated segments near the exit. On the other hand, the shape factor ⟨ ⟩ is studied to give the shape change of chain during the translocation. With the number of segments outside increasing, a maximum value ⟨ ⟩ is shown. Behind the maximum value, ⟨ ⟩ decreases more quickly for noncrowded channel than that for crowded channel. At last, the number of segments inside channel in in the process of translocation is investigated. A peak existing at the end of translocation means the translocation velocity decreases at the last lap of translocation. These findings will shed light on the translocation dynamics of polymer.