Analytical and numerical studies of sequence dependence of passage times for translocation of heterobiopolymers through nanopores

Pore shapes effects on polymer translocation

We translocated polymers through pores of different shapes and interaction patterns in three dimensions by Langevin molecular dynamics. There were four simple cylindrical pores of the same length but with different diameters. The results showed that even though decreasing the pore diameter would always decrease the translocation velocity, it was strongly dependent on the shape of the increased pore diameter. Although increasing the pore diameter made the translocation faster in simple cylindrical pores, it was complicated in different pore shapes, e.g. increasing the diameter in the middle decreased the translocation velocity. Investigating polymer shapes through the translocation process and comparing the shapes by the cumulative waiting time for different pore structures reveals the non-equilibrium properties of translocation. Moreover, polymer shape parameters such as gyration radius, polymer center of mass, and average aspect ratio help us to distinguish different pore shapes and/or different polymers.

Translocation through a narrow pore under a pulling force

We employ a three-dimensional molecular dynamics to simulate a driven polymer translocation through a nanopore by applying an external force, for four pore diameters and two external forces. To see the polymer and pore interaction effects on translocation time, we studied nine interaction energies. Moreover, to better understand the simulation results, we investigate polymer center of mass, shape factor and the monomer spatial distribution through the translocation process. Our results reveal that increasing the polymer-pore interaction energy is accompanied by an increase in the translocation time and decrease in the process rate. Furthermore, for pores with greater diameter, the translocation becomes faster. The shape analysis of the polymer indicates that the polymer shape is highly sensitive to the interaction energy. In great interactions, the monomers come close to the pore from both sides. As a result, the translocation becomes fast at first and slows down at last. Overall, it can be concluded that the external force does not play a major role in the shape and distribution of translocated monomers. However, the interaction energy between monomer and nanopore has a major effect especially on the distribution of translocated monomers on the trans side.

Force-dependent diffusion coefficient of molecular Brownian ratchets

We study the mean velocity and diffusion constant in three related models of molecular Brownian ratchets. Brownian ratchets can be used to describe translocation of biopolymers like DNA through nanopores in cells in the presence of chaperones on the trans side of the pore. Chaperones can bind to the polymer and prevent it from sliding back through the pore. First, we study a simple model that describes the translocation in terms of an asymmetric random walk. It serves as an introductory example but already captures the main features of a Brownian ratchet. We then provide an analytical expression for the diffusion constant in the classical model of a translocation ratchet that was first proposed by Peskin et al. [C. S. Peskin, G. M. Odell, and G. F. Oster, Cellular motions and thermal fluctuations: The Brownian ratchet, Biophys. J. 65, 316 (1993)BIOJAU0006-349510.1016/S0006-3495(93)81035-X]. This model is based on the assumption that the binding and unbinding of the chaperones are much faster than the diffusion of the DNA strand. To remedy this shortcoming, we propose a modified model that is also applicable if the (un)binding rates are finite. We calculate the force-dependent mean velocity and diffusivity for this model and compare the results to the original one. Our analysis shows that for large pulling forces the predictions of both models can differ strongly even if the (un)binding rates are large in comparison to the diffusion timescale but still finite. Furthermore, implications of the thermodynamic uncertainty relation on the efficiency of Brownian ratchets are discussed.

Chaperone-driven polymer translocation through nanopore: Spatial distribution and binding energy

Chaperones are binding proteins working as a driving force in biopolymer translocation. They bind to the biopolymer near the pore and prevent its backsliding. Chaperones may have different spatial distributions. Recently, we showed the importance of their spatial distribution in translocation and its effects on the sequence dependency of the translocation time. Here we focus on homopolymers and exponential distribution. Because of the exponential distribution of chaperones, the energy dependency of the translocation time will change. Here we find a minimum in translocation time versus binding effective energy (EBE) curve. The same trend can be seen in the scaling exponent of time versus polymer length, [Formula: see text] ([Formula: see text]), when plotted against EBE. Interestingly in some special cases, e.g. chaperones of size [Formula: see text] and with an exponential distribution rate of [Formula: see text], the minimum even reaches to an amount of less than 1 ([Formula: see text]). We explain the possibility of this rare result. Moreover, based on a theoretical discussion we show that, by taking into account the velocity dependency of the translocation on polymer length, one can truly predict the value of this minimum.

Chaperone-assisted translocation of flexible polymers in three dimensions

Polymer translocation through a nanometer-scale pore assisted by chaperones binding to the polymer is a process encountered in vivo for proteins. Studying the relevant models by computer simulations is computationally demanding. Accordingly, previous studies are either for stiff polymers in three dimensions or flexible polymers in two dimensions. Here, we study chaperone-assisted translocation of flexible polymers in three dimensions using Langevin dynamics. We show that differences in binding mechanisms, more specifically, whether a chaperone can bind to a single site or multiple sites on the polymer, lead to substantial differences in translocation dynamics in three dimensions. We show that the single-binding mode leads to dynamics that is very much like that in the constant-force driven translocation and accordingly mainly determined by tension propagation on the cis side. We obtain β≈1.26 for the exponent for the scaling of the translocation time with polymer length. This fairly low value can be explained by the additional friction due to binding particles. The multiple-site binding leads to translocation the dynamics of which is mainly determined by the trans side. For this process we obtain β≈1.36. This value can be explained by our derivation of β=4/3 for constant-bias translocation, where translocated polymer segments form a globule on the trans side. Our results pave the way for understanding and utilizing chaperone-assisted translocation where variations in microscopic details lead to rich variations in the emerging dynamics.

Translocation of a semiflexible polymer through a nanopore in the presence of attractive binding particles

We study the translocation dynamics of a semiflexible polymer through a nanopore from the cis into the trans compartment containing attractive binding particles (BPs) using the Langevin dynamics simulation in two dimensions. The binding particles accelerate the threading process in two ways: (i) reducing the back-sliding of the translocated monomer, and (ii) providing the pulling force toward the translocation direction. We observe that for certain binding strength (ε_{c}) and concentration (ρ) of the BPs, the translocation is faster than the ideal ratcheting condition as elucidated by Simon, Peskin, and Oster [M. Simon, C. S. Peskin, and G. F. Oster, Proc. Natl. Acad. Sci. USA 89, 3770 (1992)PNASA60027-842410.1073/pnas.89.9.3770]. The asymmetry produced by the BPs at the trans-side leads to similarities of this process to that of a driven translocation with an applied force inside the pore manifested in various physical quantities. Furthermore, we provide an analytic expression for the force experienced by the translocating chain as well as for the scaled mean first passage time (MFPT), for which we observe that for various combinations of N, ε, and ρ the scaled MFPT (〈τ〉/N^{1.5}ρ^{0.8}) collapses onto the same master plot. Based on the analysis of our simulation data, we provide plausible arguments with regard to how the scaling theory of driven translocation can be generalized for such a directed diffusion process by replacing the externally applied force with an effective force.

Polymer translocation: the first two decades and the recent diversification

Probably no other field of statistical physics at the borderline of soft matter and biological physics has caused such a flurry of papers as polymer translocation since the 1994 landmark paper by Bezrukov, Vodyanoy, and Parsegian and the study of Kasianowicz in 1996. Experiments, simulations, and theoretical approaches are still contributing novel insights to date, while no universal consensus on the statistical understanding of polymer translocation has been reached. We here collect the published results, in particular, the famous-infamous debate on the scaling exponents governing the translocation process. We put these results into perspective and discuss where the field is going. In particular, we argue that the phenomenon of polymer translocation is non-universal and highly sensitive to the exact specifications of the models and experiments used towards its analysis.

Polymer translocation through a nanopore driven by binding particles: influence of chain rigidity

We investigate the influence of chain rigidity on the dynamics of polymer translocation in the presence of binding particles (BPs) through a nanopore using two-dimensional Langevin dynamics simulations. With increasing chain rigidity κ, the mean translocation time 〈τ〉 increases monotonically due to an increase in the radius of gyration and a decrease in the center of mass velocity. Particularly for weak binding, we further find that 〈τ〉 shows a power-law behavior with the persistence length lp. Analysis indicates a scaling relation between the average velocity of the center of mass of a chain 〈vc.m.〉 and lp. As the chain becomes stiffer, the distribution of the translocation time τ approximates the Gaussian distribution and gets broader with the peak position being shifted towards longer translocation time. The corresponding translocation coordinate smax of the maximum waiting time gets smaller with increasing chain rigidity. Finally, under an extremely low BP concentration, 〈τ〉 shows a minimum for small κ, while it decreases monotonically for large κ with increasing binding energy. Our results suggest a nontrivial effect of the intrinsic property of chains on the dynamics of polymer translocation driven by BPs.

Stochastic nature of series of waiting times

Although fluctuations in the waiting time series have been studied for a long time, some important issues such as its long-range memory and its stochastic features in the presence of nonstationarity have so far remained unstudied. Here we find that the "waiting times" series for a given increment level have long-range correlations with Hurst exponents belonging to the interval 1/2<H<1. We also study positive-negative level asymmetry of the waiting time distribution. We find that the logarithmic difference of waiting times series has a short-range correlation, and then we study its stochastic nature using the Markovian method and determine the corresponding Kramers-Moyal coefficients. As an example, we analyze the velocity fluctuations in high Reynolds number turbulence and determine the level dependence of Markov time scales, as well as the drift and diffusion coefficients. We show that the waiting time distributions exhibit power law tails, and we were able to model the distribution with a continuous time random walk.

A Brownian ratchet for protein translocation including dissociation of ratcheting sites

We study a model for the translocation of proteins across membranes through a nanopore using a ratcheting mechanism. When the protein enters the nanopore it diffuses in and out of the pore according to a Brownian motion. Moreover, it is bound by ratcheting molecules which hinder the diffusion of the protein out of the nanopore, i.e. the Brownian motion is reflected such that no ratcheting molecule exits the pore. New ratcheting molecules bind at rate γ. Extending our previous approach (Depperschmidt and Pfaffelhuber in Stoch Processes Appl 120:901-925, 2010) we allow the ratcheting molecules to dissociate (at rate δ) from the protein (Model I). We also provide an approximate model (Model II) which assumes a Poisson equilibrium of ratcheting molecules on one side of the current reflection boundary. Using analytical methods and simulations we show that the speeds of both models are approximately the same. Our analytical results on Model II give the speed of translocation by means of a solution of an ordinary differential equation. This speed gives an approximation for the time it takes to translocate a protein of given length.

Sequence dependence of the binding energy in chaperone-driven polymer translocation through a nanopore

We study the translocation of stiff polymers through a nanopore, driven by the chemical-potential gradient exerted by binding proteins (chaperones) on the trans side of the pore. Bound chaperones prevent backsliding through the pore and, therefore, partially rectify the polymer passage. We show that the sequence of chain monomers with different binding affinity for the chaperones significantly affects the translocation dynamics. In particular, we investigate the effect of the nearest-neighbor adjacency probability of the two monomer types. Depending on the magnitude of the involved binding energies, the translocation speed may either increase or decrease with the adjacency probability. We determine the mean first passage time and show that, by tuning the effective binding energy, the motion changes continuously from purely diffusive to ballistic translocation.

Molecular simulation of protein dynamics in nanopores. II. Diffusion

A novel combination of discontinuous molecular dynamics and the Langevin equation, together with an intermediate-resolution model of proteins, is used to carry out long (several microsecond) simulations in order to study transport of proteins in nanopores. We simulated single-domain proteins with the alpha-helical native structure. Both attractive and repulsive interaction potentials between the proteins and the pores' walls are considered. The diffusivity D of the proteins is computed not only under the bulk conditions but also as a function of their "length" (the number of the amino-acid groups), temperature T, pore size, and interaction potentials with the walls. Compared with the experimental data, the computed diffusivities under the bulk conditions are of the correct order of magnitude. The diffusivities both in the bulk and in the pores follow a power law in the length [script-l] of the proteins and are larger in pores with repulsive walls. D(+)/D(-), the ratio of the diffusivities in pores with attractive and repulsive walls, exhibits two local maxima in its dependence on the pore size h, which are attributed to the pore sizes and protein configurations that induce long-lasting simultaneous interactions with both walls of the pores. Far from the folding temperature T(f), D increases about linearly with T, but due to the thermal fluctuations and their effect on the proteins' structure near T(f), the dependence of D on T in this region is nonlinear. We propose a novel and general "phase diagram," consisting of four regions, that describes qualitatively the effect of h, T, and interaction potentials with the walls on the diffusivity D of a protein.