Translocation of a polymer chain driven by a dichotomous noise

Translocation of a polymer chain driven by a dichotomous noise

Alessandro Fiasconaro, Juan José Mazo, Fernando Falo Departamento de Física de la Materia Condensada, Universidad de Zaragoza. 50009 Zaragoza, Spain Instituto de Ciencia de Materiales de Aragón, CSIC-Universidad de Zaragoza. 50009 Zaragoza, Spain Instituto de Biocomputación y Física de Sistemas Complejos, Universidad de Zaragoza. 50018 Zaragoza, Spain
August 2, 2019

We consider the translocation of a one-dimensional polymer through a pore channel helped by a motor driven by a dichotomous noise with time exponential correlation. We are interested in the study of the translocation time, mean velocity and stall force of the system as a function of the mean driving frequency. We find a monotonous translocation time, in contrast with the mean velocity which shows a pronounced maximum at a given frequency. Interestingly, the stall force shows a nonmonotonic behavior with the presence of a minimum. The influence of the spring elastic constant to the mean translocation times and velocities is also presented.

1 Introduction

Translocation features of polymers through natural and artificial pores is a current active research topic in biophysics and nanotechnology [1, 2, 3]. Motivated by many broad interest experimental results, different models have been introduced to describe and study in a simple way this and related problems. For instance, single barrier potentials [4], as well as flashing ratchet models [5], have been studied to describe the polymer translocation and polymer transport dynamics. The passage of small molecules through passive cell channels can be also modeled by stochastic and rachetlike forces [6]. In some cases the transport phenomena involves not translocation through pores, but also molecular motors, whose complex action has been recently addressed at high attention [7, 8]. In addition, nanotechnological applications try to emulate the complex biological process related to the translocation dynamics [9, 10].

Recently, we have studied different models for the 1d translocation of a spring-bead polymer helped by a motor using a sinusoidal force [11]. The introduction of a time dependent driving force imposes a new time scale on the system, and provides new and richer phenomenology: for sinusoidal driving, the translocation time shows an oscillatory behavior as a function of the frequency.

In order to introduce stochasticity in the motor action and motivated by the relevant role played by dichotomous noise in biological problems, in this manuscript we consider the case of a polymer driven by a two-state force: constant force which pushes the polymer chain in one direction during the activity of the motor, and zero force which leaves the polymer to diffuse freely otherwise. This pure dichotomous mechanism constitutes a first approach in describing a machine working dichotomously between two on-off states [12, 13].

The motor modeled in [12, 13] acts during a fixed time, while the waiting times are exponentially distributed with a mean time depending on the ATP concentration. In the present work a simpler dichotomous mechanism which can well point out, by contrast, the specific behavior of the ATP based machines is studied.

On the other hand, pure dichotomous driving makes sense in the nanotechnological context as well as in the biophysical one. In the first case the passage of a polymer can be induced through a graphene pore or solid state channeling [14, 15] by applying a dichotomous force between the two sides of the layer. In the second case, the model can describe the translocation of a linear molecule through a cell membrane gate having a chemical potential difference between their two sides. The driving is in this case induced by the typical open/close mechanism which follow the purely dichotomous switching largely used in literature [6, 16, 17].

The purpose of our work is to model phenomenologically the possible physical systems described above. We want to stress here the qualitative specific results connected to the purely dichotomous driving.

Thus, differently from the sinusoidal case, no special behavior is observed in the mean translocation time of the polymer for the case here studied. However, for this problem, another observable parameter can be studied. In fact, single molecule experiments are able to detect and use the instantaneous velocity in order to quantify the translocation process in forced systems [7]. Remarkably, we find a non trivial behavior of the polymer translocation velocity as a function of the mean frequency of the driving with the presence of a maximum, even if the translocation time shows only a monotonic behavior. This difference reveals the importance of dealing with several measures to explore the complex behavior of the polymer translocation.

The dependence of the stall force of the machine is also calculated. We find, again, a strong nonmonotonic behavior of with the frequency, similar to the one found in [11].

The paper is organized as follows: first we present the model for polymer and the properties of the stochastic driving force. The main properties of the translocation process are then calculated: translocation time, mean velocity and stall force. Finally, we analyze the dependence of the above properties with the chain stiffness.

2 The model

Figure 1: Scheme of a linear chain driven by a dichotomous force restricted to a small space region (width ). is the mean time during which the force maintains a same value .

The polymer is modeled as a unidimensional chain of dimensionless monomers connected by harmonic springs [20].

1d models are suitable in order to describe the dynamics of polymers constrained to move in confined channels [15]. Also, in many experimental situations [7] the polymer is stretched, thus removing the dimensionality dependence of the measured quantities. Moreover, in this work, we want to fix our attention to the motor activity in the translocation more than the delay given by other effects, such as entropic contributions.

The total potential energy is


where is the elastic constant, the position of the -th particle, and the equilibrium distance between adjacent monomers.

The translocation is helped by the presence of a motor which is activated dichotomously. The machine has a spatial working width and the position represents the right edge of its action (see Fig. 1). Thus the monomers such that experience a force made by the motor. We define to represent the dichotomous force, which fluctuates between two values (no force) and . Thus


Here gives the mean residence time in each state. With respect to the spatial dependence


The dynamics of the monomer of the chain is then described by the following overdamped Langevin equations:


where the viscosity parameter for each monomer is included in the normalized time units. stands for Gaussian uncorrelated thermal fluctuation and follows the usual statistical properties and .

3 Results

We performed a set of numerical experiment with a stochastic Runge-Kutta algorithm, using a time step of . The polymer is compound by monomers and starts with all the spring at the rest length (), and the last monomer of the chain lies at (), just in the final action range of the dichotomous force. The noise intensity is held fixed at the value , , and . The choice of the number of monomers , or equivalently the length , is arbitrary and this small number has been used for computational convenience. We note that in a previous work [11], also with 1d chain, it was found that scales with . Similarly we find that scales with .

In this first part, the elastic constant is held equal to 1, a meaningful choice that corresponds to a not too rigid approximation for the polymer. We will study the main observables of the system as a function of the mean frequency transition .

3.1 Translocation times

Figure 2: Polymer translocation time as a function of the frequency of the fluctuating force. The solid line shows the theoretical prediction of Eq. 6. In the inset is shown the standard deviation of the exit times distribution.

The translocation time is computed as a mean first passage time of the center of mass of the polymer: the average over the realizations of the time spent by the center of mass of the chain to reach the position . In Fig. 2 we see the value of plotted as function of the mean frequency of the driving, and, in the inset, the standard deviation whose values are of the same order of magnitude than the mean time values, as expected. We find that is a monotonic function of . This result is different from that for a periodic force (sinusoidal or square wave ones) where it is observed a minimum in the translocation for and an oscillating behavior for higher frequencies [11].

In contrast with the behavior of the translocation time, as we will see, the velocity is not a monotonic function of and a maximum is found in this function for (see Fig. 3). Both effects (minimum translocation time [11] or maximum mean velocity) reveal some interesting similarities with the resonant activation phenomenon [18, 19].

We can make a simple analytical prediction for in the low frequency region which however is found to be valid in a broad frequency range (see solid line in Fig. 2). Let be the value of the exit time when a constant force is applied during all the dynamics. In the limit we have to distinguish between two cases depending on the initial value of the force, or . In the first case the translocation time is corrected in a first approximation by a long waiting time if the system switches to the off state before , which occurs with a probability . This correction gives a contribution of to the total time. In the second case there is an additional time in the off state for escaping. Thus the total translocation time is


Since this equation is derived in the low frequency limit where we have


The intermediate frequency region is characterized by the presence of the constant force alternated by the absence of the force (diffusive dynamics) with an average time ratio between them different for different values of the mean frequency. Surprisingly, Eq. (6) also describes in a good way that frequency region. The third region is instead characterized by a high frequency switching rate between the two force states. There the translocation time is much smaller than , the polymer experiences a mean force , and . A careful observation of our numerical results show that in this high frequency regime


not observable in Fig.2.

Figure 3: Mean velocity and number of monomers inside the motor while in its active state (inset) as a function of the mean frequency of the fluctuating force.
Figure 4: First passage time (left) and velocity (right) probability distributions at four frequency values.

3.2 Mean velocity

Fig. 3 shows the mean velocity of the polymer as a function of . The inset of the figure shows the average number of monomers inside the motor during the active states, . We can see that this number is not constant for different values of the mean frequency , at least for the value of the elastic constant used in those calculations.

The main result in the velocity curves is the presence of a well pronounced maximum, which put in evidence the qualitative difference between the calculation of the mean first passage time and the mean velocity. In fact, the velocity is computed as , where is the escape time in the -th realization.

As visible in Fig. 4, the exit times distribution changes in shape by changing the mean driving frequency . For high values of , the time distribution is very narrow around its mean value . The corresponding probability distribution function for the velocity is also a narrow function. Decreasing the value of , the distributions are more asymmetric and the width increases. The maximum of the time distribution moves toward lower values of time, but the asymmetry changes and higher and higher values of translocation times are involved. That’s why in that region the mean first passage time increases, although the time of the ’s maximum decreases. The distribution of the velocity change as well; but because of the increased width of the time distribution, the mean value of the velocity in that region does not follow the relation and increases with respect to the high frequency limit, in opposite direction as the one expected from the time behavior. The reason of this effect is that in the average, the smaller times have a higher weight in the inverse than bigger ones. Thus the mean velocity increases up to a maximum. Decreasing in the low frequency region, the average of the times continues rising up, because the distribution involves higher and higher times. The velocity, however, now decreases since the very high times escapes do not contribute importantly to the mean velocity.

In a first approach the translocation velocity in the high frequency limit is given by


a fraction of monomers given by experience a force . Then, the corresponding translocation time is (remind that is the escape time if the motor is always working). On the contrary, in the low frequency limit one half of the realizations give a very long escape time (and velocity goes to zero) and another half give . Thus for low frequencies we obtain the same value of the velocity that for high frequencies.

However, we can see in Fig. 3 that the low frequency limit of the mean velocity does not satisfy the relationship just derived, being lower than the high frequency value . This happens because the force exerted on the polymer is affected by the number of monomers inside the motor which, as shown in the inset of the figure, also depends on . We will see below a confirmation of the given relation by using a strong elastic constant between the monomers, which guarantees a constant number of monomers inside the motor nevertheless the dynamical conditions are (see Fig. 7).

From Eq. (4) it is easy to derive the following equation for the mean velocity


where, for each experiment , is the average number of particles inside the motor when the motor is , is the total motor working time, and is the translocation time of each realization.

At high frequency, . However, decreasing the frequency for most of the cases, , during the translocation the motor spends more time activated that deactivated since most translocations happen during the activation stage of the motor. Thus both, the translocation time and the mean velocity increase111This is not the case at low values of , where the strong change in with the frequency dominates the overall behavior and suppress the velocity maximum as shown in the inset of Fig. 7 for .. This behavior changes when . Then remains constant in Eq. (9), increases when decreases and the velocity also decreases towards the expected value moderate by the mean number of monomers in the motor in the low frequency limit. This explains the presence of the maximum in the velocity.

A rough estimation of is given by the fixed value , corresponding to the distribution of monomer inside the motor in the case that they maintain the same relative distance, equal to the rest, over all the dynamics. This condition will be completely satisfied for high values of the elastic constant (rigid chain limit), when becomes independent on . As we will see below, both the high and low frequency limits for the mean velocity take in that case the same value (see the inset of Fig. 7)

which is slightly higher than the limit value shown in the inset of Fig. 2 because .

3.3 Stall Force

The stall force is the force that we need to apply against the motor in order to stop the polymer translocation. It is a measure of the strength of the motor and, in this model, it depends on the frequency of the driving.

Figure 5: A force pull is applied at the first monomer to measure the motor stall force.

A set of simulations have been performed by applying a pull force (see Fig. 5) on the left extremum of the chain, in opposite direction to the motor driving force.  The initial condition for the chain has been fixed with the polymer center of mass in the center of the motor. Then, the velocity of the center of mass is measured waiting for the exit on the left or on the right of the motor region. That way, the force for which the mean velocity is zero gives .

Fig. 6 shows the stall force as a function of the frequency. As shown in the lower inset, for a given frequency the mean velocity decreases linearly with . The upper inset, shows that for pull forces of the order of the stall force the velocity presents a minimum, contrary to the behavior at (Fig. 3). Then the stall force, which presents a similar trend, shows a clear minimum in the same frequency region.

Figure 6: Stall force as a function of the frequency of the driving. The upper inset shows the mean velocity as a function of for the three pull force values . The bottom  inset shows the linear behavior of the mean velocity as a function of the pull force for the three frequency values . The other parameters are the same of Fig. 2.

As in the oscillating case [11], the scale variation of the stall force is small (around ), and an experimental verification of its behavior with the mean frequency the minimum could be not immediately simple to perform.

3.4 Elastic constant dependence

Finally, we investigate the dependence of translocation time and velocity on the elastic constant of the polymer. A magnitude that strongly depends on is the mean number of monomer inside the motor during the pushing cycle, . This number modules the velocity as it is show in Eq (9). Results are plotted in Fig. 7, where the translocation time and velocity (in the inset) are presented for different value of .

Figure 7: Translocation time and mean velocity (in the inset) of the polymer chain for different values of the elastic constant .

We can see that notable differences (especially visible in the mean velocity plot) are evident by changing the value of . For the case the velocity looses the maximum, which is always present for the higher values of . As expected, a clear saturating behavior of the whole curve is evident by increasing when the chain behaves like a rigid bar. As announced before in the text, in this limit the mean velocity in the cases of both high and low switching frequency gives the expected value given above. This limit is already fulfilled for .

4 Conclusions

The interest in the introduction of simple models is that they can capture the more relevant features of different processes. In that way, they can result to be very useful for a coarse-grain description of different systems.

The model described here studies the translocation process of a polymer driven by a simple motor which exerts a dichotomous force. We analyze the dependence of the translocation time with the mean frequency of the driving field, and find an analytical expression for the low frequency regime. In spite of the monotonic behavior of the translocation time, the velocity presents a clear maximum at a resonant value of the mean frequency. We argue that this maximum comes from the optimization of the ”on states” duration of the driving forces with the corresponding translocation time. The detection of this maximum, (also seen in the periodic case) could be tackled with the recent single molecule experimental techniques.

The stall force able to stop the polymer translocation against the motor has been also evaluated, finding in our calculations results very close to the oscillating driving, previously studied. The stall force show a very clear minimum at a resonant mean frequency of the driving.

The model can have application in artificial nanotechnological devices driven by dichotomously fluctuating fields, as well as biological pore membrane with intrinsic noise.

This work has been supported by the project FIS2008-01240 of the Spain MICINN.


  • [1] J. J. Kasianowicz, E. Brandin, D. Branton, and D. W. Deamer, Proc. Natl. Acad. Sci. USA 93, 13770 (1996)
  • [2] M.Zwolak and M. Di Ventra, Rev. Mod. Phys 80, 141 (2008)
  • [3] G. F. Schneider et al., Nano Lett. 10 3163 (2010).
  • [4] N. Pizzolato, A. Fiasconaro, B. Spagnolo, Int. J. Bifurc. Chaos 18, 2871 (2008); N. Pizzolato, A. Fiasconaro, and B. Spagnolo, J. Stat. Mech: Theory and Exp. P01011 (2009).
  • [5] M.T. Downton, M.J.Zuckermann, E.M.Craig, M.Plischke, H.Linke, Phys. Rev. E 73, 011909 (2006); E.M.Craig, M.J.Zuckermann, H.Linke, Phys. Rev. E 73, 051106 (2006)
  • [6] I. Kosztin and K. Schulten. Phys. Rev. Lett. 93, 238102 (2004)
  • [7] D. E. Smith, S. J. Tans, S. B. Smith, S. Grimes, D. L. Anderson, and C. Bustamante, Nature 413, 748 (2001)
  • [8] J. R. Moffitt, Y. R. Chemla, K.Aathavan, S. Grimes, P. J. Jardine, D. L. Anderson, and C. Bustamante, Nature 457, 446 (2009).
  • [9] M. Mickler, E. Schleiff, and T. Hugel Chem. Phys. Chem. 9, 1503 (2008).
  • [10] E. B. Starikov, D. Henning, H. Yamada, R. Gutierrez, B. Norden, and G. Cuniberti, Biophys. Rev. and Lett. 4(3), 209 (2009)
  • [11] A. Fiasconaro, F. Falo, and J.J. Mazo, Phys. Rev. E 82, 031803 (2010).
  • [12] A. Gomez-Marin and J. M. Sancho, Eur. Phys. Lett. 86, 40002 (2009)
  • [13] A. Fiasconaro, F. Falo, and J.J. Mazo. New Journal of Physics 14, 023004 (2012).
  • [14] J. Han, S. W. Turner, and H. G. Craighead, Phys. Rev. Lett. 83, 1688 (1999).
  • [15] B. Luan, H. Peng, S. Polonsky, S. Rossnagel, G. Stolovitzky, and G. Martyna, Phys. Rev. Lett. 104, 238103 (2010)
  • [16] M. M. Millonas,and D. R. Chialvo, Phys. Rev. Lett. 76, 550 (1996).
  • [17] A. Kargol and K. Kabza, Physical Biology 5, 026003 (2008).
  • [18] N. Pizzolato, A. Fiasconaro, D. Persano Adorno, B. Spagnolo, Physical Biology 7, 034001 (2010).
  • [19] C.R. Doering, and J. C. Gadoua, Phys. Rev. Lett. 69, 2318 (1992); M. Bier, and R.D. Astumian, Phys. Rev. Lett. 71, 1649 (1993).
  • [20] P. E. J. Rouse, J. Chem. Phys. 21, 1272 (1953).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description