Dynamics of polymer translocation through a nanopore induced by different sizes of crowding agents

Protocols for translocation processes of flexible polymers through a pore using LAMMPS

Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) is an open-source, powerful simulator with a customizable platform for extensive Langevin dynamics simulations. Here, we present a protocol for using LAMMPS to develop coarse-grained models of polymeric systems with macromolecular crowding, an integral part of any soft matter or biophysical system. We describe steps for installing software, using LAMMPS basic commands and code, and translocating polymers. This protocol has potential applications in developing mathematical models for DNA sequencing, controlled drug delivery, and cellular transport processes. For complete details on the use and execution of this protocol, please refer to Garg et al.1.

Macromolecular Crowding Facilitates ssDNA Capture within Biological Nanopores: Role of Size Variation and Solution Heterogeneity

Genetic sequencing is a vital process that requires the transport of charged nucleic acids through transmembrane nanopores. Single-molecule studies show that macromolecular bulk crowding facilitates the capture of these polymers, leading to a high throughput of nanopore sensors. Motivated by these observations, a minimal discrete-state stochastic framework was developed to describe the role of poly(ethylene glycol) (PEG) crowders in varying concentrations in the transport of ssDNA through α-hemolysin nanopores. This theory suggested that the cooperative partitioning of polycationic PEGs controls the capture of ssDNA due to underlying electrostatic interactions. Herein, we investigate the impact of the size variation of PEGs on the capture event. Even though larger crowders attract ssDNA strongly to enhance its capture, our results show that considerable cooperative partitioning of PEGs is also required to achieve high interevent frequency. The exact analytical results are supported by existing single-molecule studies. Since real cellular conditions are heterogeneous, its influence on the ssDNA capture rate is studied by introducing a binary mixture of crowders. Our results indicate that the "polymer-pushing-polymer" concept possibly affects the capture rate depending on the mixture composition. These new findings provide valuable insights into the microscopic mechanism of the capture process, which eventually allows for accurate genome sequencing in crowded solutions.

Forced and spontaneous translocation dynamics of a semiflexible active polymer in two dimensions

Polymer translocation is a fundamental topic in non-equilibrium physics and is crucially important to many biological processes in life. In the present work, we adopt two-dimensional Langevin dynamics simulations to study the forced and spontaneous translocation dynamics of an active filament. The influence of polymer stiffness on the underlying dynamics is explicitly analyzed. For the forced translocation, the results show a robust stiffness-induced inhibition, and the translocation time exhibits a dual-exponent scaling relationship with the bending modulus. Tension propagation (TP) is also examined, where we find prominent modifications in terms of both activity and stiffness. For spontaneous translocation into a pure solvent, the translocation time is almost independent of the polymer stiffness. However, when the polymer is translocated into a porous medium, an intriguing non-monotonic alteration of translocation time with increasing chain stiffness is demonstrated. The semiflexible chain is beneficial for translocation while the rigid chain is not conducive. Stiffness regulation on the diffusion dynamics of the polymer in porous media shows a consistent scenario. The interplay of activity, stiffness, and porous crowding provides a new mechanism for understanding the non-trivial translocation dynamics of an active filament in complex environments.

Translocation Dynamics of an Active Filament through a Long-Length Scale Channel

Active filament translocation through a confined space is crucial for diverse biological processes. By using Langevin dynamics simulations, we investigate the translocation dynamics of an axially self-propelled chain through a channel. First, results show a suggestive reciprocal scaling of translocation time versus active force. Second, in the case of a long channel, we demonstrate a very intriguing nonmonotonic change of translocation time with increasing channel width. The driving force shows a similar trend, providing a consistent picture to understand the unexpected channel width effect. In particular, in a moderately broad channel, the disordered chain conformation results in a loss of driving force and thus inhibits translocation dynamics. Chain adsorption might occur in a wide channel, which accounts for a facilitated translocation. Lastly, we connect the translocation process to tension propagation (TP). A modified TP picture is proposed to interpret the waiting time distribution. Our work highlights the new phenomenology owing to the crucial interplay of activity and spacial confinement, which drives the translocation dynamics, going beyond the traditional entropic barrier scenario.

Polymer Translocation and Nanopore Sequencing: A Review of Advances and Challenges

Various biological processes involve the translocation of macromolecules across nanopores; these pores are basically protein channels embedded in membranes. Understanding the mechanism of translocation is crucial to a range of technological applications, including DNA sequencing, single molecule detection, and controlled drug delivery. In this spirit, numerous efforts have been made to develop polymer translocation-based sequencing devices, these efforts include findings and insights from theoretical modeling, simulations, and experimental studies. As much as the past and ongoing studies have added to the knowledge, the practical realization of low-cost, high-throughput sequencing devices, however, has still not been realized. There are challenges, the foremost of which is controlling the speed of translocation at the single monomer level, which remain to be addressed in order to use polymer translocation-based methods for sensing applications. In this article, we review the recent studies aimed at developing control over the dynamics of polymer translocation through nanopores.

Monte Carlo Simulation of Static and Dynamic Properties of Linear Polymer in a Crowded Environment

In this paper, we investigate the static and dynamic properties of linear polymer in the presence of obstacles. A Monte Carlo (MC) simulation method in two dimensions with a bond fluctuation model (BFM) was used to achieve this goal. To overcome the entropic barrier, we put the middle monomer of the polymer in the middle of the pore, which is placed between ordered and disordered obstacles. We probed the static properties of the polymer by calculating the mean square of the radius of gyration and the mean square end-to-end distance of the polymer, and we found that the scaling exponents of both the mean square end-to-end distance $〈R^2〉$ and the mean square radius of gyration $〈R_g^2〉$ as a function of the polymer length $N$ vary with the area fraction of crowding agents, $\varphi$ . The dynamic properties have also been studied by exploring the translocation of the polymer. Our current research shows that the escape time $\tau$ increases as $\varphi$ increases. Moreover, in the power-law relation of escape time $\tau$ as a function of polymer length $N$ , the scaling exponent ( $\alpha$ ) changes with $\varphi$ . Furthermore, the study has shown that the translocation of the polymer favors the disordered barriers.

Polymer looping kinetics in active heterogeneous environments

A typical biological environment is usually featured by crowding and heterogeneity, leading to complex reaction kinetics of the immersed macromolecules. In the present work, we adopt Langevin dynamics simulations to systematically investigate polymer looping kinetics in active heterogeneous media crowded with a mixture of mobile active particles and immobile obstacles. For comparison, a parallel study is also performed in the passive heterogeneous media. We explicitly analyze the change of looping time and looping probability with the variation of obstacle ratio, volume fraction and crowder size. We reveal the novel phenomena of inhibition-facilitation transition of the looping rate induced by heterogeneity, crowdedness and activity. In addition, our results demonstrate a very non-trivial crowder size effect on the looping kinetics. The underlying mechanism is rationalized by the interplay of polymer diffusion, conformational change and looping free-energy barrier. The competing effect arising from active particles and obstacles on structural and dynamical properties of the polymer yields a consistent scenario for our observations. Lastly, the non-exponential kinetics of the looping process is also analyzed. We find that both activity and crowding can strengthen the heterogeneity degree of the looping kinetics.

Comparative study of polymer looping kinetics in passive and active environments

Intra-chain looping in complex environments is significant in advancing our understanding of biological processes in life. We adopt Langevin dynamics simulations to perform a comparative study of polymer looping kinetics in passive and active environments. From the analysis of looping quantities, including looping-unlooping times and looping probabilities, we unraveled the intriguing effects of active crowder size, activity and crowding. Firstly, we figured out the phase diagram involving a novel facilitation-inhibition transition in the parameter space of active crowder size and active force, and the two-fold roles of activity are clarified. In particular, we find that active particles of a size comparable to the polymer monomer are most favorable for facilitated looping, while those with a similar size to the polymer gyration radius impede the looping most seriously. Secondly, the underlying looping mechanisms in different active crowder size regimes are rationalized by the interplay among diffusion, polymer conformational change and the free-energy barrier. For small active crowders, activity significantly promotes end-to-end distance diffusion, which dominantly facilitates both looping and unlooping processes. In the case of moderate active crowders, the polymer chain suffers from prominent swelling, and thus inevitable inhibited looping will occur. For large active crowders, activity induces a counterintuitive non-cage effect on the looping kinetics, through yielding a higher effective temperature and larger unlooping free-energy barrier. This is in sharp contrast to the caging phenomena observed in passive media. Lastly, the volume-fraction dependence of the looping quantities in an active bath demonstrates dramatic discrepancies from that in a passive bath, which highlights the contrasting effects of activity and crowding.

Effects of active crowder size and activity-crowding coupling on polymer translocation

Polymer translocation in complex environments is crucially important to many biological processes in life. In the present work, we adopted two-dimensional Langevin dynamics simulations to study the forced and unbiased polymer translocation dynamics in active and crowded media. The translocation time and probability are analyzed in terms of active force F_a, volume fraction φ and also the crowder size. The non-trivial active crowder size effect and activity-crowding coupling effect as well as the novel mechanism of unbiased translocation between two active environments with different active particle sizes are clarified. Firstly, for forced translocation, we reveal an intriguing non-monotonic dependence of the translocation time on the crowder size in the case of large activity. In particular, crowders of intermediate size similar to the polymer segment are proven to be the most favorable for translocation. Moreover, a facilitation-inhibition crossover of the translocation time with increasing volume fraction is observed, indicating a crucial activity-crowding coupling effect. Secondly, for unbiased translocation driven by different active crowder sizes, the translocation probability demonstrates a novel turnover phenomenon, implying the appearance of an opposite directional preference as the active force exceeds a critical value. The translocation time in both directions decreases monotonically with the active force. The asymmetric activity effect together with the entropic driving scenario provides a reasonable picture for the peculiar behavior observed in unbiased translocation.

Non-monotonic dependence of polymer chain dynamics on active crowder size

Configuration dynamics of flexible polymer chains is of ubiquitous importance in many biological processes. Here, we investigate a polymer chain immersed in a bath of size-changed active particles in two dimensional space using Langevin dynamics simulations. Particular attention is paid to how the radius of gyration R_g of the polymer chain depends on the size σ_c of active crowders. We find that R_g shows nontrivial non-monotonic dependence on σ_c: The chain first swells upon increasing σ_c, reaching a fully expanded state with maximum R_g, and then, R_g decreases until the chain collapses to a compact coil state if the crowder is large enough. Interestingly, the chain may oscillate between a collapse state and a stretched state at moderate crowder size. Analysis shows that it is the competition between two effects of active particles, one stretching the chain from inside due to persistence motion and the other compressing the chain from outside, that leads to the non-monotonic dependence. Besides, the diffusion of the polymer chain also shows nontrivial non-monotonic dependence on σ_c. Our results demonstrate the important interplay between particle activity and size associated with polymer configurations in active crowding environments.

Polymer translocation into cavities: Effects of confinement geometry, crowding, and bending rigidity on the free energy

Monte Carlo simulations are used to study the translocation of a polymer into a cavity. Modeling the polymer as a hard-sphere chain with a length up to N=601 monomers, we use a multiple-histogram method to measure the variation of the conformational free energy of the polymer with respect to the number of translocated monomers. The resulting free-energy functions are then used to obtain the confinement free energy for the translocated portion of the polymer. We characterize the confinement free energy for a flexible polymer in cavities with constant cross-sectional area A for various cavity shapes (cylindrical, rectangular, and triangular) as well as for tapered cavities with pyramidal and conical shape. The scaling of the free energy with cavity volume and translocated polymer subchain length is generally consistent with predictions from simple scaling arguments, with small deviations in the scaling exponents likely due to finite-size effects. The confinement free energy depends strongly on cavity shape anisometry and is a minimum for an isometric cavity shape with a length-to-width ratio of unity. Entropic depletion at the edges or vertices of the confining cavity are evident in the results for constant-A and pyramidal cavities. For translocation into infinitely long cones, the scaling of the free energy with taper angle is consistent with a theoretical prediction employing the blob model. We also examine the effects of polymer bending rigidity on the translocation free energy for cylindrical cavities. For isometric cavities, the observed scaling behavior is in partial agreement with theoretical predictions, with discrepancies arising from finite-size effects that prevent the emergence of well-defined scaling regimes. In addition, translocation into highly anisometric cylindrical cavities leads to a multistage folding process for stiff polymers. Finally, we examine the effects of crowding agents inside the cavity. We find that the confinement free energy increases with crowder density. At constant packing fraction the magnitude of this effect lessens with increasing crowder size for a crowder-to-monomer size ratio ≥1.

Interactions of a Polypeptide with a Protein Nanopore Under Crowding Conditions

Molecular crowding, a ubiquitous feature of the cellular environment, has significant implications in the kinetics and equilibrium of biopolymer interactions. In this study, a single charged polypeptide is exposed to competing forces that drive it into a transmembrane protein pore versus forces that pull it outside. Using single-molecule electrophysiology, we provide compelling experimental evidence that the kinetic details of the polypeptide-pore interactions are substantially affected by high concentrations of less-penetrating polyethylene glycols (PEGs). At a polymer concentration above a critical value, the presence of these neutral macromolecular crowders increases the rate constant of association but decreases the rate constant of dissociation, resulting in a stronger polypeptide-pore interaction. Moreover, a larger-molecular weight PEG exhibits a lower rate constant of association but a higher rate constant of dissociation than those values corresponding to a smaller-molecular weight PEG. These outcomes are in accord with a lower diffusion constant of the polypeptide and higher depletion-attraction forces between the polypeptide and transmembrane protein pore under crowding and confinement conditions.

DNA size in confined environments

For short DNA molecules in crowded environments, we evaluate macroscopic parameters such as the average end-to-end distance and the twist conformation by tuning the strength of the site specific confinement driven by the crowders.

Polymer translocation through nanopore into active bath

Polymer translocation through nanopores into a crowded environment is of ubiquitous importance in many biological processes. Here we investigate polymer translocation through a nanopore into an active bath of self-propelled particles in two-dimensional space using Langevin dynamics simulations. Interestingly, we find that the mean translocation time τ can show a bell-shape dependence on the particle activity F_a at a fixed volume fraction ϕ, indicating that the translocation process may become slower for small activity compared to the case of the passive media, and only when the particle activity becomes large enough can the translocation process be accelerated. In addition, we also find that τ can show a minimum as a function of ϕ if the particle activity is large enough, implying that an intermediate volume fraction of active particles is most favorable for the polymer translocation. Detailed analysis reveals that such nontrivial behaviors result from the two-fold effect of active bath: one that active particles tend to accumulate near the pore, providing an extra pressure hindering the translocation, and the other that they also aggregate along the polymer chain, generating an effective pulling force accelerating the translocation. Such results demonstrate that active bath plays rather subtle roles on the polymer translocation process.

Translocation of polymers into crowded media with dynamic attractive nanoparticles

The translocation of polymers through a small pore into crowded media with dynamic attractive nanoparticles is simulated. Results show that the nanoparticles at the trans side can affect the translocation by influencing the free-energy landscape and the diffusion of polymers. Thus the translocation time τ is dependent on the polymer-nanoparticle attraction strength ɛ and the mobility of nanoparticles V. We observe a power-law relation of τ with V, but the exponent is dependent on ɛ and nanoparticle concentration. In addition, we find that the effect of attractive dynamic nanoparticles on the dynamics of polymers is dependent on the time scale. At a short time scale, subnormal diffusion is observed at strong attraction and the diffusion is slowed down by the dynamic nanoparticles. However, the diffusion of polymers is normal at a long time scale and the diffusion constant increases with the increase in V.

Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size

The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.

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.