| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |

* Mathematical and Statistical Computing Laboratory, Center for Information Technology, National Institutes of Health, Bethesda, Maryland 20892, and
Laboratory of Chemical Physics, National Institute of Diabetes, Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892
Correspondence: Address reprint requests to Irina V. Gopich, Laboratory of Chemical Physics, NIDDK, Bldg. 5/136, National Institutes of Health, Bethesda, MD 20892. Tel.: 301-496-1880; E-mail: irinag{at}speck.niddk.nih.gov.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
-hemolysin ion channels (Kasianowicz et al., 1996
The probability density found by Kasianowicz et al. (1996)
has three peaks. The authors interpreted the three-peak density as a superposition of two two-peak probability densities corresponding to different orientations of the DNA molecules in the channels. Position of the first (short time) peak very weakly depends on both the chain length and applied voltage. The authors attributed this peak to those molecules that did not translocate and escaped on the same side of the membrane where they had entered. The second peak located at larger times is due to the molecules that traverse the membrane. Its position depends on the orientation of the macromolecule. For both orientations the position linearly grows with the chain length and decreases, as the applied electric field, E, increases, as 1/E.
The Lubensky-Nelson theory (Lubensky and Nelson, 1999
) of driven polymer translocation through a narrow pore in a thin membrane exploits the idea that passage of a long polymer through the pore can be described in terms of one-dimensional diffusion of a point particle on interval (0, L), where L is the contour length of the polymer. Lubensky and Nelson calculated the probability density for the polymer lifetime in the channel by solving the corresponding diffusion equation with absorbing boundary conditions at the ends of the interval. These boundary conditions have an unpleasant feature: the particle cannot enter the interval through an absorbing end point because it is instantly trapped by the absorbing boundary. Lubensky and Nelson called this "a pathology of our model."
In this paper we generalize the Lubensky-Nelson theory by replacing absorbing boundary condition by radiation ones. Indeed, when a polymer inserts its head (or tail) into the channel, it can either escape on the same side of the membrane or make a couple of steps forward and, perhaps, eventually translocate the membrane. In terms of the one-dimensional diffusion model, such a behavior can be described using radiation (partially, but not perfectly absorbing) boundary conditions. As Lubensky and Nelson, we neglect the effect of the entropic barrier on the translocation, i.e., consider a rodlike polymer. For this model in the next section we find an exact solution for the probability density of the polymer lifetime in the channel. We also derive the translocation and return probabilities and average times that translocating and nontranslocating molecules spend in the channel as well as the probability densities for these times.
In "Results and Discussion", we analyze the dependence of the probability density of the polymer lifetime in the channel on the chain length and the driving force. In sufficiently strong fields, we found the same behavior of the probability density as was observed in the experiment by Kasianowicz et al. (1996)
. It is interesting that the probability density has two peaks only when the driving force is strong enough. For weaker forces the density has qualitatively the same behavior as in the absence of the force and monotonically decays with time. A critical driving force that separates monotonic and two-peak probability densities depends on the competition between the drift and diffusion as well as the trapping efficiency at the ends of interval. We analyze this question in "Results and Discussion, one peak or two peaks".
Finally, it is worth mentioning that the problem we are discussing in this paper is a fragment of a general problem of the polymer translocation through membrane pores. One can gain broader view on this subject from recent papers (Simon et al., 1992
; Peskin et al., 1993
; Baumgartner and Scolnick, 1994
, 1995
; Yoon and Deutsch, 1995
; Lee and Obukhov, 1996
; Schatz and Dobberstein, 1996
; Sung and Park, 1996
; Deutsch and Yoon, 1997
; Di Marzio and Mandell, 1997
; Carl, 1998
; Park and Sung, 1998a
,b
,c
; De Gennes, 1999a
,b
,c
; Di Marzio, 1999
; Han et al., 1999
; Muthukumar, 1999
, 2001
; Kumar and Sebastian, 2000
; Chern et al., 2001
; Chuang et al., 2002
) and references therein.
| THEORY |
|---|
|
|
|---|
![]() | (1) |
and kB and T are the Boltzmann constant and absolute temperature. Equation 1 describes diffusion of the particle in the linear potential U(x) = -Fx that gives rise to a constant drift velocity
![]() | (2) |
![]() | (3) |
(x) should be understood as
(x + 0), so that
The radiation boundary conditions imposed at the ends of the interval are
![]() | (4a) |
![]() | (4b) |
describes instantaneous escape into the bulk that corresponds to perfectly absorbing boundary conditions, the case studied by Lubensky and Nelson. k = 0 corresponds to reflecting boundary conditions that describe a situation when the macromolecule is anchored on both sides and therefore cannot escape from the membrane.
|
![]() | (5) |
![]() | (6) |
![]() | (7) |
describes macromolecules that do not translocate and escape the channel at time t on the same side of the membrane where they entered. The flux
gives the contribution into
(t) due to translocating macromolecules that escape the channel at time t.
These fluxes can be used to find probabilities that the macromolecule translocates,
, and does not translocate,
, through the membrane
![]() | (8) |
![]() | (9) |
![]() | (10) |
![]() | (11) |
monotonically decays with time whereas
is nonmonotonic; initially it grows with time, reaches a maximum, and then decays. Whether
(t) is monotonic or has two peaks depends on the competition between the two terms in the right-hand side of Eq. 11. This question is analyzed in detail in "Results and Discussion, one peak or two peaks".
Solutions
One can find exact solutions for the Laplace transforms of the fluxes and probability densities. The Laplace transform of a function F(t) denoted by
where s is the Laplace parameter, is defined as
![]() | (12) |
![]() | (13a) |
![]() | (13b) |
Then one can find the Laplace transform of the probability density
(t), which according to Eq. 7 is given by
![]() | (14) |
= 0 because the macromolecule anchored on both sides cannot escape from the membrane. In the opposite limiting case, when
,
and, hence,
, whereas
as it must be with absorbing boundary conditions.
The expressions in Eqs. 13 allows one to find the probabilities defined in Eq. 8. The translocation probability is given by
![]() | (15) |
![]() | (16) |
is finite. When
= 0, the fluxes in Eqs. 13 vanish and, hence,
as it must be when the macromolecule is anchored on both sides. When
,
, whereas
as it must be.
The Laplace transforms of the conditional densities defined in Eq. 9 are given by
![]() | (17) |
![]() | (18) |
![]() | (19) |
All the quantities considered above depend on the problem parameters. These dependencies are discussed in the next section. Before proceeding to this discussion, we note that the Laplace transform of the Lubensky-Nelson result for
is recovered from
in Eq. 17 in the limiting case of k
. Here, according to Eq. 15,
as it must be. However, the ratio in Eq. 17 remains finite because the flux in the numerator also vanishes as k
.
Note that the fluxes
and
can be derived using the eigenfunction expansion of the propagator. However, we found the Laplace transform formalism more convenient for our analysis.
Concluding this section, we indicate that the theory above should be modified in several respects to treat the case when two orientations of the macromolecule differ from one another. First, one has to use different rate constants in the radiation boundary conditions at x = 0 and x = L in Eqs. 4. Second, the diffusion constant in Eq. 1 may depend on the orientation of the macromolecule. Third, the probability density for the blockade time is a linear combination of the solutions found for different orientations. The coefficients in this linear combination are the entrance probabilities for macromolecules of different orientations. To find these probabilities, one has to analyze how macromolecules enter the channel, which seems to us an extremely difficult task.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
respectively. Using these scales, we introduce dimensionless time,
and will analyze the dependence of its density on the dimensionless driving force,
and the dimensionless length of the macromolecule,
The Laplace transform in Eq. 14 is numerically inverted using the Stehfest algorithm (Stehfest, 1970
and
Dependence on the field
To study how the density depends on
we inverted the dimensionless version of the Laplace transform in Eq. 14 numerically for
= 0, 1, 5, 10, 20, and 30 at constant
The densities are shown in Fig. 2. The density monotonically decreases with time when
whereas for all other values of
the density has two peaks. Position of the second peak,
moves to shorter times as
increases. Times
are given in Table 1. In this table we also give the product
One can see that the product is practically a constant when the field is strong enough. This means that in such fields
as was found in the experiment (Kasianowicz et al., 1996
).
|
|
where
is given by the dimensionless version of Eq. 19. To check this assumption, we give times
in Table 1. One can see that
and
are close in sufficiently strong fields that confirms the assumption. In weaker fields, the two times differ and, hence, the assumption is not valid.
Dependence on the length of the macromolecule
To study the
-dependence of the density, we inverted the dimensionless version of the Laplace transform in Eq. 14 numerically for
= 10, 20, and 30 at constant value of
The densities are shown in Fig. 3. As one might expect, the longer is the macromolecule, the larger is the time
Times
are given in Table 2, where we also give the ratio
and average translocation times,
One can see that the ratio is practically a constant. This means that
as was found experimentally (Kasianowicz et al., 1996
). One can also see good agreement between times
and
that supports the assumption about Gaussian shape of
for given values of
and
|
|
decreases, the second peak moves to larger times and its amplitude decreases. One might expect that the probability density
is monotonic only when
, whereas at any finite
it has two peaks. In fact, this is not true, and
remains a monotonic function even at finite
if it is sufficiently small.
To show this, we use Eq. 11 that gives
as a sum of the monotonically decaying term,
and the nonmonotonic term,
In Fig. 4 we show these two terms and their sum,
as well as first derivatives of these functions at
and
One can see that the derivative
is always negative. It has a local maximum at
, where
. When
increases, the monotonic term decreases whereas the growing part of the nonmonotonic term increases, as shown in Fig. 5. As a consequence, when the driving force increases,
approaches zero, then becomes equal to zero at a certain critical field,
and eventually becomes positive when
This means that at
, the growing part of
is sufficiently large to produce the second peak as shown in Fig. 6.
|
|
|
(t) is monotonic or has two peaks depends on these parameters also. For example, when k
0 and/or D
, the equilibration occurs much faster than escape and we have
![]() | (20) |
![]() | (21) |
= kßF coth(ßFL/2). One should not be confused by the absence of
in probabilities given in Eq. 20, which are the leading terms of the probabilities given in Eqs. 15 and 16 in the small
limit. Additional explanation of the dependence on
is given just below Eq. 16. | CONCLUDING REMARKS |
|---|
|
|
|---|
The translocation probability and Laplace transform of the residence time probability density were derived within the framework of our model. The Laplace transform was numerically inverted to study how the probability density depends on the length of the macromolecule and external driving force. It is interesting that the probability density changes its shape at a certain critical value of the driving force. The density has two peaks when the driving force is larger that the critical value and monotonically decreases with time when the force is smaller that this value. The dependencies predicted by the model in sufficiently strong fields are the same as those observed in experiments on driven translocation of single-stranded RNA and DNA molecules through single
-hemolysin ion channels by Kasianowicz et al. (1996)
.
The model analyzed in the paper is oversimplified because we assume that the polymer is a rigid rod. Single-stranded DNA and RNA molecules used in the experiment correspond more closely to a freely jointed polymer chain. For such polymers the entropy barrier has to be taken into consideration as first indicated by Sung and Park (1996)
. Neglecting the entropy potential, we considerably simplify the problem. This makes it possible to derive an exact solution for the Laplace transform of the probability density of the polymer lifetime in the channel, which is used in our analysis to study the shape of the probability density as a function of the driving force and the length of the polymer. It is impossible to find an analytical solution without this simplifying assumption, i.e., when the entropy barrier is included. The fact that our model predicts the behavior qualitatively similar to that observed in the experiment suggests that the contribution into the driving force due to the entropy potential is small compared to the contribution due to the external field when the field is strong enough. We suppose to analyze this question in the future.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
I. V. Gopich is on leave from the Institute of Chemical Kinetics and Combustion, Siberian Branch of the Russian Academy of Sciences, Novosibirsk 630090, Russia.
Submitted on June 13, 2002; accepted for publication October 2, 2002.
| REFERENCES |
|---|
|
|
|---|
Baumgartner, A., and J. Scolnick. 1994. Polymer electrophoresis across a model membrane. J. Phys. Chem. 98:1065510658.
Baumgartner, A., and J. Skolnick. 1995. Spontaneous translocation of a polymer across a curved membrane. Phys. Rev. Lett. 74:21422145.[Medline]
Bezrukov, S. M., I. Vodyanoy, and V. A. Parsegian. 1994. Counting polymers moving through a single ion channel. Nature. 370:279281.[Medline]
Bezrukov, S. M. 2000. Ion channels as molecular coulter counters to probe metabolite transport. J. Membr. Biol. 174:113.[Medline]
Carl, W. A. 1998. A polymer threading a membrane: model system for a molecular pump. J. Chem. Phys. 108:79217922.
Chern, S. S., A. E. Cardenas, and R. D. Coalson. 2001. Three-dimensional dynamic Monte Carlo simulations of driven polymer transport through a hole in a wall. J. Chem. Phys. 115:77727782.
Chuang, J., Y. Kantor, and M. Kardar. 2002. Anomalous dynamics of translocation. Phys. Rev. E. 65:118021118028.
De Gennes, P. G. 1999a. Problems of DNA entry into a cell. Physica A. 274:17.
De Gennes, P. G. 1999b. Flexible polymers in nanopores. Adv. Polym. Sci. 138:91105.
De Gennes, P. G. 1999c. Passive entry of a DNA molecule into a small pore. Proc. Natl. Acad. Sci. USA. 96:72627264.
Deutsch, J. M., and H. Yoon. 1997. Macromolecular separation through a porous surface. J. Chem. Phys. 106:93769381.
Di Marzio, E. A., and A. J. Mandell. 1997. Phase transition behavior of a linear macromolecule threading a membrane. J. Chem. Phys. 107:55105514.
Di Marzio, E. A. 1999. The ten classes of polymeric phase transitions: their use as models for self-assembly. Prog. Polym. Sci. 24:329377.
Han, J., S. W. Turner, and H. G. Craighead. 1999. Entropic trapping and escape of long DNA molecules at submicron size constriction. Phys. Rev. Lett. 83:16881691.
Henrickson, S. E., M. Misakian, B. Robertson, and J. J. Kasianowicz. 2000. Driven DNA transport into an asymmetric nanometer-scale pore. Phys. Rev. Lett. 85:30573060.[Medline]
Kasianowicz, J. J., E. Brandin, D. Branton, and D. W. Deamer. 1996. Characterization of individual polynucleotide molecules using a membrane channel. Proc. Natl. Acad. Sci. USA. 93:1377013773.
Kumar, K. K., and K. L. Sebastian. 2000. Adsorption-assisted translocation of a chain molecule through a pore. Phys. Rev. E. 62:75367539.
Lee, N., and S. Obukhov. 1996. Diffusion of a polymer chain through a thin membrane. J. Phys. II. 6:195204.
Lubensky, D. K., and D. R. Nelson. 1999. Driven polymer translocation through a narrow pore. Biophys. J. 77:18241838.
Meller, A., L. Nivon, and D. Branton. 2001. Voltage-driven DNA translocations through a nanopore. Phys. Rev. Lett. 86:34353433435348.
Muthukumar, M. 1999. Polymer translocation through a hole. J. Chem. Phys. 111:1037110374.
Muthukumar, M. 2001. Translocation of a confined polymer through a hole. Phys. Rev. Lett. 86:31883191.[Medline]
Park, P. J., and W. Sung. 1998a. A stochastic model of polymer translocation dynamics through biomembranes. Int. J. Bifurcat. Chaos. 8:927931.
Park, P. J., and W. Sung. 1998b. Polymer release out of a spherical vesicle through a pore. Phys. Rev. E. 57:730734.
Park, P. J., and W. Sung. 1998c. Polymer translocation induced by adsorption. J. Chem. Phys. 108:30133018.
Peskin, C. S., G. M. Odell, and G. F. Oster. 1993. Cellular motions and thermal fluctuations: the Brownian ratchet. Biophys. J. 65:316324.
Schatz, G., and B. Dobberstein. 1996. Common principles of protein translocation across membranes. Science. 271:15191526.[Abstract]
Simon, S. M., C. S. Peskin, and G. F. Oster. 1992. What drives the translocation of proteins? Proc. Natl. Acad. Sci. USA. 89:37703774.
Stehfest, H. 1970. Algorithm 368: Numerical inversion of Laplace transforms. Comm. ACM. 13:4749.
Sung, W., and P. J. Park. 1996. Polymer translocation through a pore in a membrane. Phys. Rev. Lett. 77:783786.[Medline]
Valko, P. P., and S. Vajda. 2002. Inversion of noise-free Laplace transforms: towards a standardized set of test problems. Inverse Problems in Engineering. 10:467483. (See also http://pumpjack.tamu.edu/
valko/).
Yoon, H. G., and J. M. Deutsch. 1995. Dynamics of a polymer in the presence of permeable membranes. J. Chem. Phys. 102:90909095.
This article has been cited by other articles:
![]() |
O. K. Dudko, J. Mathe, A. Szabo, A. Meller, and G. Hummer Extracting Kinetics from Single-Molecule Force Spectroscopy: Nanopore Unzipping of DNA Hairpins Biophys. J., June 15, 2007; 92(12): 4188 - 4195. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. R. Bauer and W. Nadler From the Cover: Molecular transport through channels and pores: Effects of in-channel interactions and blocking PNAS, August 1, 2006; 103(31): 11446 - 11451. [Abstract] [Full Text] [PDF] |
||||
![]() |
I.-C. Yeh and G. Hummer Nucleic acid transport through carbon nanotube membranes PNAS, August 17, 2004; 101(33): 12177 - 12182. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Flomenbom and J. Klafter Translocation of a Single-Stranded DNA Through a Conformationally Changing Nanopore Biophys. J., June 1, 2004; 86(6): 3576 - 3584. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Movileanu, S. Cheley, and H. Bayley Partitioning of Individual Flexible Polymers into a Nanoscopic Protein Pore Biophys. J., August 1, 2003; 85(2): 897 - 910. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |