1
|
Upadhyay G, Kapri R, Chaudhuri A. Homopolymer and heteropolymer translocation through patterned pores under fluctuating forces. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2024; 47:23. [PMID: 38573533 DOI: 10.1140/epje/s10189-024-00417-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Accepted: 03/18/2024] [Indexed: 04/05/2024]
Abstract
We investigate the translocation of a semiflexible polymer through extended patterned pores using Langevin dynamics simulations, specifically focusing on the influence of a time-dependent driving force. Our findings reveal that, akin to its flexible counterpart, a rigid chain-like molecule translocates faster when subjected to an oscillating force than a constant force of equivalent average magnitude. The enhanced translocation is strongly correlated with the stiffness of the polymer and the stickiness of the pores. The arrangement of the pores plays a pivotal role in translocation dynamics, deeply influenced by the interplay between polymer stiffness and pore-polymer interactions. For heterogeneous polymers with periodically varying stiffness, the oscillating force introduces significant variations in the translocation time distributions based on segment sizes and orientations. On the basis of these insights, we propose a sequencing approach that harnesses distinct pore surface properties that are capable of accurately predicting sequences in heteropolymers with diverse bending rigidities.
Collapse
Affiliation(s)
- Gokul Upadhyay
- Department of Physical Sciences, Indian Institute of Science Education and Research Mohali, Sector 81, Knowledge City, S. A. S. Nagar, Manauli, 140306, India
| | - Rajeev Kapri
- Department of Physical Sciences, Indian Institute of Science Education and Research Mohali, Sector 81, Knowledge City, S. A. S. Nagar, Manauli, 140306, India
| | - Abhishek Chaudhuri
- Department of Physical Sciences, Indian Institute of Science Education and Research Mohali, Sector 81, Knowledge City, S. A. S. Nagar, Manauli, 140306, India.
| |
Collapse
|
2
|
Upadhyay G, Kapri R, Chaudhuri A. Gain reversal in the translocation dynamics of a semiflexible polymer through a flickering pore. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:185101. [PMID: 38262064 DOI: 10.1088/1361-648x/ad21a9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Accepted: 01/23/2024] [Indexed: 01/25/2024]
Abstract
We study the driven translocation of a semiflexible polymer through an attractive extended pore with a periodically oscillating width. Similar to its flexible counterpart, a stiff polymer translocates through an oscillating pore more quickly than a static pore whose width is equal to the oscillating pore's mean width. This efficiency quantified as a gain in the translocation time, highlights a considerable dependence of the translocation dynamics on the stiffness of the polymer and the attractive nature of the pore. The gain characteristics for various polymer stiffness exhibit a trend reversal when the stickiness of the pore is changed. The gain reduces with increasing stiffness for a lower attractive strength of the pore, whereas it increases with increasing stiffness for higher attractive strengths. Such a dependence leads to the possibility of a high degree of robust selectivity in the translocation process.
Collapse
Affiliation(s)
- Gokul Upadhyay
- Department of Physical Sciences, Indian Institute of Science Education and Research Mohali, Sector 81, Knowledge City, S. A. S. Nagar, Manauli 140306, India
| | - Rajeev Kapri
- Department of Physical Sciences, Indian Institute of Science Education and Research Mohali, Sector 81, Knowledge City, S. A. S. Nagar, Manauli 140306, India
| | - Abhishek Chaudhuri
- Department of Physical Sciences, Indian Institute of Science Education and Research Mohali, Sector 81, Knowledge City, S. A. S. Nagar, Manauli 140306, India
| |
Collapse
|
3
|
Singh SL, Chauhan K, Bharadwaj AS, Kishore V, Laux P, Luch A, Singh AV. Polymer Translocation and Nanopore Sequencing: A Review of Advances and Challenges. Int J Mol Sci 2023; 24:6153. [PMID: 37047125 PMCID: PMC10094227 DOI: 10.3390/ijms24076153] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Revised: 02/01/2023] [Accepted: 02/28/2023] [Indexed: 03/31/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Swarn Lata Singh
- Department of Physics, Mahila Mahavidyalaya (MMV), Banaras Hindu University, Varanasi 221005, UP, India
| | - Keerti Chauhan
- Department of Physics, Banaras Hindu University, Varanasi 221005, UP, India
| | - Atul S. Bharadwaj
- Department of Physics, CMP Degree College, University of Allahabad, Prayagraj 211002, UP, India
| | - Vimal Kishore
- Department of Physics, Banaras Hindu University, Varanasi 221005, UP, India
| | - Peter Laux
- Department of Chemical and Product Safety, German Federal Institute of Risk Assessment (BfR) Maxdohrnstrasse 8-10, 10589 Berlin, Germany
| | - Andreas Luch
- Department of Chemical and Product Safety, German Federal Institute of Risk Assessment (BfR) Maxdohrnstrasse 8-10, 10589 Berlin, Germany
| | - Ajay Vikram Singh
- Department of Chemical and Product Safety, German Federal Institute of Risk Assessment (BfR) Maxdohrnstrasse 8-10, 10589 Berlin, Germany
| |
Collapse
|
4
|
Lu LW, Wang ZH, Shi AC, Lu YY, An LJ. Polymer Translocation. CHINESE JOURNAL OF POLYMER SCIENCE 2023. [DOI: 10.1007/s10118-023-2975-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/29/2023]
|
5
|
Sharma A, Kapri R, Chaudhuri A. Driven translocation of a semiflexible polymer through a conical channel in the presence of attractive surface interactions. Sci Rep 2022; 12:19081. [PMID: 36351960 PMCID: PMC9646819 DOI: 10.1038/s41598-022-21845-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Accepted: 10/04/2022] [Indexed: 11/11/2022] Open
Abstract
We study the translocation of a semiflexible polymer through a conical channel with attractive surface interactions and a driving force which varies spatially inside the channel. Using the results of the translocation dynamics of a flexible polymer through an extended channel as control, we first show that the asymmetric shape of the channel gives rise to non-monotonic features in the total translocation time as a function of the apex angle of the channel. The waiting time distributions of individual monomer beads inside the channel show unique features strongly dependent on the driving force and the surface interactions. Polymer stiffness results in longer translocation times for all angles of the channel. Further, non-monotonic features in the translocation time as a function of the channel angle changes substantially as the polymer becomes stiffer, which is reflected in the changing features of the waiting time distributions. We construct a free energy description of the system incorporating entropic and energetic contributions in the low force regime to explain the simulation results.
Collapse
Affiliation(s)
- Andri Sharma
- grid.458435.b0000 0004 0406 1521Department of Physical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector 81, S. A. S. Nagar, Mohali, 140306 Punjab India
| | - Rajeev Kapri
- grid.458435.b0000 0004 0406 1521Department of Physical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector 81, S. A. S. Nagar, Mohali, 140306 Punjab India
| | - Abhishek Chaudhuri
- grid.458435.b0000 0004 0406 1521Department of Physical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector 81, S. A. S. Nagar, Mohali, 140306 Punjab India
| |
Collapse
|
6
|
Seth S, Bhattacharya A. How capture affects polymer translocation in a solitary nanopore. J Chem Phys 2022; 156:244902. [PMID: 35778106 PMCID: PMC9225749 DOI: 10.1063/5.0094221] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Accepted: 06/02/2022] [Indexed: 11/14/2022] Open
Abstract
DNA capture with high fidelity is an essential part of nanopore translocation. We report several important aspects of the capture process and subsequent translocation of a model DNA polymer through a solid-state nanopore in the presence of an extended electric field using the Brownian dynamics simulation that enables us to record statistics of the conformations at every stage of the translocation process. By releasing the equilibrated DNAs from different equipotentials, we observe that the capture time distribution depends on the initial starting point and follows a Poisson process. The field gradient elongates the DNA on its way toward the nanopore and favors a successful translocation even after multiple failed threading attempts. Even in the limit of an extremely narrow pore, a fully flexible chain has a finite probability of hairpin-loop capture, while this probability decreases for a stiffer chain and promotes single file translocation. Our in silico studies identify and differentiate characteristic distributions of the mean first passage time due to single file translocation from those due to translocation of different types of folds and provide direct evidence of the interpretation of the experimentally observed folds [M. Gershow and J. A. Golovchenko, Nat. Nanotechnol. 2, 775 (2007) and Mihovilovic et al., Phys. Rev. Lett. 110, 028102 (2013)] in a solitary nanopore. Finally, we show a new finding-that a charged tag attached at the 5' end of the DNA enhances both the multi-scan rate and the uni-directional translocation (5' → 3') probability that would benefit the genomic barcoding and sequencing experiments.
Collapse
Affiliation(s)
- Swarnadeep Seth
- Department of Physics, University of Central Florida, Orlando, Florida 32816-2385, USA
| | - Aniket Bhattacharya
- Department of Physics, University of Central Florida, Orlando, Florida 32816-2385, USA
| |
Collapse
|
7
|
Wang Z, Wang R, Lu Y, An L, Shi AC, Wang ZG. Mechanisms of Flow-Induced Polymer Translocation. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00288] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Zhenhua Wang
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, P. R. China
| | - Ruishu Wang
- Department of Mathematics, Jilin University, Changchun 130012, P. R. China
| | - Yuyuan Lu
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, P. R. China
| | - Lijia An
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, P. R. China
| | - An-Chang Shi
- Department of Physics and Astronomy, McMaster University, Hamilton, Ontario L8S 4M1, Canada
| | - Zhen-Gang Wang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
8
|
Bley M, Dzubiella J. Nonequilibrium free energy during polymer chain growth. J Chem Phys 2022; 156:084902. [DOI: 10.1063/5.0080786] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
During fast diffusion-influenced polymerization, nonequilibrium behavior of the polymer chains and the surrounding reactive monomers has been reported recently. Based on the laws of thermodynamics, the emerging nonequilibrium structures should be characterizable by some “extra free energy” (excess over the equilibrium Helmholtz free energy). Here, we study the nonequilibrium thermodynamics of chain-growth polymerization of ideal chains in a dispersion of free reactive monomers, using off-lattice, reactive Brownian dynamics computer simulations in conjunction with approximative statistical mechanics and relative entropy (Gibbs–Shannon and Kullback–Leibler) concepts. In the case of fast growing polymers, we indeed report increased nonequilibrium free energies Δ Fneq of several kB T compared to equilibrium and near-equilibrium, slowly growing chains. Interestingly, Δ Fneq is a non-monotonic function of the degree of polymerization and thus also of time. Our decomposition of the thermodynamic contributions shows that the initial dominant extra free energy is stored in the nonequilibrium inhomogeneous density profiles of the free monomer gas (showing density depletion and wakes) in the vicinity of the active center at the propagating polymer end. At later stages of the polymerization process, we report significant extra contributions stored in the nonequilibrium polymer conformations. Finally, our study implies a nontrivial relaxation kinetics and “restoring” of the extra free energy during the equilibration process after polymerization.
Collapse
Affiliation(s)
- Michael Bley
- Applied Theoretical Physics–Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Strasse 3, D-79104 Freiburg, Germany
| | - Joachim Dzubiella
- Applied Theoretical Physics–Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Strasse 3, D-79104 Freiburg, Germany
- Cluster of Excellence livMats@FIT–Freiburg Center for Interactive Materials and Bioinspired Technologies, Albert-Ludwigs-Universität Freiburg, Georges-Köhler-Allee 105, D-79110 Freiburg, Germany
| |
Collapse
|
9
|
Makhani EY, Zhang A, Haun JB. Quantifying and controlling bond multivalency for advanced nanoparticle targeting to cells. NANO CONVERGENCE 2021; 8:38. [PMID: 34846580 PMCID: PMC8633263 DOI: 10.1186/s40580-021-00288-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/11/2021] [Accepted: 11/07/2021] [Indexed: 06/13/2023]
Abstract
Nanoparticles have drawn intense interest as delivery agents for diagnosing and treating various cancers. Much of the early success was driven by passive targeting mechanisms such as the enhanced permeability and retention (EPR) effect, but this has failed to lead to the expected clinical successes. Active targeting involves binding interactions between the nanoparticle and cancer cells, which promotes tumor cell-specific accumulation and internalization. Furthermore, nanoparticles are large enough to facilitate multiple bond formation, which can improve adhesive properties substantially in comparison to the single bond case. While multivalent binding is universally believed to be an attribute of nanoparticles, it is a complex process that is still poorly understood and difficult to control. In this review, we will first discuss experimental studies that have elucidated roles for parameters such as nanoparticle size and shape, targeting ligand and target receptor densities, and monovalent binding kinetics on multivalent nanoparticle adhesion efficiency and cellular internalization. Although such experimental studies are very insightful, information is limited and confounded by numerous differences across experimental systems. Thus, we focus the second part of the review on theoretical aspects of binding, including kinetics, biomechanics, and transport physics. Finally, we discuss various computational and simulation studies of nanoparticle adhesion, including advanced treatments that compare directly to experimental results. Future work will ideally continue to combine experimental data and advanced computational studies to extend our knowledge of multivalent adhesion, as well as design the most powerful nanoparticle-based agents to treat cancer.
Collapse
Affiliation(s)
- Elliot Y Makhani
- Department of Materials Science and Engineering, University of California Irvine, Irvine, CA, 92697, USA
| | - Ailin Zhang
- Department of Biomedical Engineering, University of California Irvine, 3107 Natural Sciences II, Irvine, CA, 92697, USA
| | - Jered B Haun
- Department of Materials Science and Engineering, University of California Irvine, Irvine, CA, 92697, USA.
- Department of Biomedical Engineering, University of California Irvine, 3107 Natural Sciences II, Irvine, CA, 92697, USA.
- Department of Chemical and Biomolecular Engineering, University of California Irvine, Irvine, CA, 92697, USA.
- Chao Family Comprehensive Cancer Center, University of California Irvine, Irvine, CA, 92697, USA.
| |
Collapse
|
10
|
Bley M, Baul U, Dzubiella J. Controlling solvent quality by time: Self-avoiding sprints in nonequilibrium polymerization. Phys Rev E 2021; 104:034501. [PMID: 34654077 DOI: 10.1103/physreve.104.034501] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Accepted: 08/22/2021] [Indexed: 12/20/2022]
Abstract
A fundamental paradigm in polymer physics is that macromolecular conformations in equilibrium can be described by universal scaling laws, being key for structure, dynamics, and function of soft (biological) matter and in the materials sciences. Here we reveal that during diffusion-influenced, nonequilibrium chain-growth polymerization, scaling laws change qualitatively, in particular, the growing polymers exhibit a surprising self-avoiding walk behavior in poor and θ solvents. Our analysis, based on monomer-resolved, off-lattice reaction-diffusion computer simulations, demonstrates that this phenomenon is a result of (i) nonequilibrium monomer density depletion correlations around the active polymerization site, leading to a locally directed and self-avoiding growth, in conjunction with (ii) chain (Rouse) relaxation times larger than the competing polymerization reaction time. These intrinsic nonequilibrium mechanisms are facilitated by fast and persistent reaction-driven diffusion ("sprints") of the active site, with analogies to pseudochemotactic active Brownian particles. Our findings have implications for time-controlled structure formation in polymer processing, as in, e.g., reactive self-assembly, photocrosslinking, and three-dimensional printing.
Collapse
Affiliation(s)
- Michael Bley
- Applied Theoretical Physics-Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, D-79104 Freiburg, Germany
| | - Upayan Baul
- Applied Theoretical Physics-Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, D-79104 Freiburg, Germany
| | - Joachim Dzubiella
- Applied Theoretical Physics-Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, D-79104 Freiburg, Germany.,Cluster of Excellence livMats@FIT-Freiburg Center for Interactive Materials and Bioinspired Technologies, Albert-Ludwigs-Universität Freiburg, D-79110 Freiburg, Germany
| |
Collapse
|
11
|
Hoseinpoor SM, Nikoofard N, Ha BY. Characteristic time for the end monomers of a spherically confined polymer to find a nano-pore. J Chem Phys 2021; 154:114901. [PMID: 33752364 DOI: 10.1063/5.0040551] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Translocation of a polymer through a nano-pore is relevant in a variety of contexts such as passage of RNAs through a nuclear pore and transportation of proteins across a membrane. An essential step in polymer translocation is for the end monomers to search the pore. This process requires a characteristic time, referred to as the "attempt time" in this work. Here, we study the attempt time τ of a confined polymer inside a spherical surface by combining a scaling approach and Langevin dynamics simulations. For a moderately to strongly confined polymer, our results suggest that τ ∼ R3.67 for R > P and τ ∼ R2.67 for R < P, where R is the radius of the spherical surface and P is the persistence length of the polymer. All simulation data obtained for an intermediate range of the volume fraction of monomers ϕ(≲ 0.2) tend to collapse onto each other. This implies that τ does not explicitly depend on ϕ, in agreement with the theoretical predictions. These results will be useful for interpreting translocation as a two-step process: the initial attempt to find the pore and eventual pore crossing.
Collapse
Affiliation(s)
- S Mohammad Hoseinpoor
- Institute of Nanoscience and Nanotechnology, University of Kashan, Kashan 51167-87317, Iran
| | - Narges Nikoofard
- Institute of Nanoscience and Nanotechnology, University of Kashan, Kashan 51167-87317, Iran
| | - Bae-Yeun Ha
- Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
12
|
Sun LZ, Cao WP, Wang CH, Xu X. The translocation dynamics of the polymer through a conical pore: Non-stuck, weak-stuck, and strong-stuck modes. J Chem Phys 2021; 154:054903. [PMID: 33557527 DOI: 10.1063/5.0033689] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The external voltage-driven polymer translocation through a conical pore (with a large opening at the entry and a small tip at the exit) is studied by using the Langevin dynamics simulation in this paper. The entire translocation process is divided into an approaching stage and a threading stage. First, the approaching stage starts from the polymer entering the large opening and ends up at a terminal monomer reaching the pore tip. In this stage, the polymer will undergo the conformation adjustment to fit the narrowed cross-sectional area of the pore, leading to three approaching modes: the non-stuck mode with a terminal monomer arriving at the pore tip smoothly, the weak-stuck mode for the polymer stuck inside the pore for a short duration with minor conformational adjustments, and the strong-stuck mode with major conformational changes and a long duration. The approaching times (the duration of the approaching stage) of the three approaching modes show different behavior as a function of the pore apex angle. Second, the threading stage describes that the polymer threads through the pore tip with a linear fashion. In this stage, an increase in the apex angle causes the reduction of the threading time (the duration of the threading stage) due to the increase in the driving force with the apex angle at the tip. Moreover, we also find that with the increase in the apex angle or the polymer length, the polymer threading dynamics will change from the quasi-equilibrium state to the non-equilibrium state.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China
| | - Wei-Ping Cao
- Institute of Optoelectronic Technology, Lishui University, Lishui 323000, China
| | - Chang-Hui Wang
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China
| | - Xiaojun Xu
- Institute of Bioinformatics and Medical Engineering, Jiangsu University of Technology, Changzhou 213001, China
| |
Collapse
|
13
|
Qiao L, Ignacio M, Slater GW. An efficient kinetic Monte Carlo to study analyte capture by a nanopore: transients, boundary conditions and time-dependent fields. Phys Chem Chem Phys 2021; 23:1489-1499. [PMID: 33400742 DOI: 10.1039/d0cp03638b] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
To better understand the capture process by a nanopore, we introduce an efficient Kinetic Monte Carlo (KMC) algorithm that can simulate long times and large system sizes by mapping the dynamic of a point-like particle in a 3D spherically symmetric system onto the 1D biased random walk. Our algorithm recovers the steady-state analytical solution and allows us to study time-dependent processes such as transients. Simulation results show that the steady-state depletion zone near pore is barely larger than the pore radius and narrows at higher field intensities; as a result, the time to reach steady-state is much smaller than the time required to empty a zone of the size of the capture radius λe. When the sample reservoir has a finite size, a second depletion region propagates inward from the outer wall, and the capture rate starts decreasing when it reaches the capture radius λe. We also note that the flatness of the electric field near the pore, which is often neglected, induces a traffic jam that can increase the transient time by several orders of magnitude. Finally, we propose a new proof-of-concept scheme to separate two analytes of the same mobility but different diffusion coefficients using time-varying fields.
Collapse
Affiliation(s)
- Le Qiao
- Department of Physics, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada.
| | - Maxime Ignacio
- Department of Physics, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada.
| | - Gary W Slater
- Department of Physics, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada.
| |
Collapse
|
14
|
Haji Abdolvahab R, Niknam Hamidabad M. Pore shapes effects on polymer translocation. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2020; 43:76. [PMID: 33306147 DOI: 10.1140/epje/i2020-12001-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Accepted: 11/12/2020] [Indexed: 06/12/2023]
Abstract
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.
Collapse
|
15
|
Niknam Hamidabad M, Asgari S, Haji Abdolvahab R. Nanoparticle-assisted polymer translocation through a nanopore. POLYMER 2020. [DOI: 10.1016/j.polymer.2020.122847] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
16
|
Qiao L, Ignacio M, Slater GW. Voltage-driven translocation: Defining a capture radius. J Chem Phys 2019; 151:244902. [DOI: 10.1063/1.5134076] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Affiliation(s)
- Le Qiao
- Department of Physics, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Maxime Ignacio
- Department of Physics, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Gary W. Slater
- Department of Physics, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| |
Collapse
|
17
|
|
18
|
Abstract
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.
Collapse
|
19
|
Polson JM, Heckbert DR. Polymer translocation into cavities: Effects of confinement geometry, crowding, and bending rigidity on the free energy. Phys Rev E 2019; 100:012504. [PMID: 31499877 DOI: 10.1103/physreve.100.012504] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2019] [Indexed: 06/10/2023]
Abstract
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.
Collapse
Affiliation(s)
- James M Polson
- Department of Physics, University of Prince Edward Island, 550 University Avenue, Charlottetown, Prince Edward Island, Canada C1A 4P3
| | - David R Heckbert
- Department of Physics, University of Prince Edward Island, 550 University Avenue, Charlottetown, Prince Edward Island, Canada C1A 4P3
| |
Collapse
|
20
|
Kwon S, Lee S, Cho HW, Kim J, Kim JS, Sung BJ. The breakdown of the local thermal equilibrium approximation for a polymer chain during packaging. J Chem Phys 2019; 150:204901. [PMID: 31153198 DOI: 10.1063/1.5093946] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The conformational relaxation of a polymer chain often slows down in various biological and engineering processes. The polymer, then, may stay in nonequilibrium states throughout the process such that one may not invoke the local thermal equilibrium (LTE) approximation, which has been usually employed to describe the kinetics of various processes. In this work, motivated by recent single-molecule experiments on DNA packaging into a viral capsid, we investigate how the nonequilibrium conformations and the LTE approximation would affect the packaging of a polymer chain into small confinement. We employ a simple but generic coarse-grained model and Langevin dynamics simulations to investigate the packaging kinetics. The polymer segments (both inside and outside the confinement) stay away from equilibrium under strong external force. We devise a simulation scheme to invoke the LTE approximation during packaging and find that the relaxation of nonequilibrium conformations plays a critical role in regulating the packaging rate.
Collapse
Affiliation(s)
- Seulki Kwon
- Department of Chemistry, Sogang University, Seoul 04107, South Korea
| | - Seulgi Lee
- Department of Chemistry, Sogang University, Seoul 04107, South Korea
| | - Hyun Woo Cho
- Department of Chemistry, Sogang University, Seoul 04107, South Korea
| | - Jeongmin Kim
- Department of Chemistry, Sogang University, Seoul 04107, South Korea
| | - Jun Soo Kim
- Department of Chemistry and Nano Science, Ewha Womans University, Seoul 03760, South Korea
| | - Bong June Sung
- Department of Chemistry, Sogang University, Seoul 04107, South Korea
| |
Collapse
|
21
|
Sun LZ, Wang CH, Luo MB, Li H. Trapped and non-trapped polymer translocations through a spherical pore. J Chem Phys 2019; 150:024904. [PMID: 30646715 DOI: 10.1063/1.5063331] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The polymer translocation through a spherical pore is studied using the Langevin dynamics simulation. The translocation events are classified into two types: one is the trapped translocation in which the entire polymer is trapped in the pore and the other is the non-trapped translocation where the pore cannot hold the whole polymer. We find that the trapped translocation is favored at large spheres and small external voltages. However, the monomer-pore attraction would lead to the non-monotonic behavior of the trapped translocation possibility out of all translocation events. Moreover, both the trapped and non-trapped translocation times are dependent on the polymer length, pore size, external voltage, and the monomer-pore attraction. There exist two pathways for the polymer in the trapped translocation: an actively trapped pathway for the polymer trapped in the pore before the head monomer arrives at the pore exit, and a passively trapped pathway for the polymer trapped in the pore while the head monomer is struggling to move out of the pore. The studies of trapped pathways can provide a deep understanding of the polymer translocation behavior.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China
| | - Chang-Hui Wang
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China
| | - Meng-Bo Luo
- Department of Physics, Zhejiang University, Hangzhou 310027, China
| | - Haibin Li
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China
| |
Collapse
|
22
|
Kwon S, Sung BJ. Effects of solvent quality and non-equilibrium conformations on polymer translocation. J Chem Phys 2018; 149:244907. [PMID: 30599703 DOI: 10.1063/1.5048059] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The conformation and its relaxation of a single polymer depend on solvent quality in a polymer solution: a polymer collapses into a globule in a poor solvent, while the polymer swells in a good solvent. When one translocates a polymer through a narrow pore, a drastic conformational change occurs such that the kinetics of the translocation is expected to depend on the solvent quality. However, the effects of solvent quality on the translocation kinetics have been controversial. In this study, we employ a coarse-grained model for a polymer and perform Langevin dynamics simulations for the driven translocation of a polymer in various types of solvents. We estimate the free energy of polymer translocation using steered molecular dynamics simulations and Jarzynski's equality and find that the free energy barrier for the translocation increases as the solvent quality becomes poorer. The conformational entropy contributes most to the free energy barrier of the translocation in a good solvent, while a balance between entropy and energy matters in a poor solvent. Interestingly, contrary to what is expected from the free energy profile, the translocation kinetics is a non-monotonic function of the solvent quality. We find that for any type of solvent, the polymer conformation stays far away from the equilibrium conformation during translocation due to an external force and tension propagation. However, the degree of tension propagation differs depending on the solvent quality as well as the magnitude of the external force: the tension propagation is more significant in a good solvent than in a poor solvent. We illustrate that such differences in tension propagation and non-equilibrium conformations between good and poor solvents are responsible for the complicated non-monotonic effects of solvent quality on the translocation kinetics.
Collapse
Affiliation(s)
- Seulki Kwon
- Department of Chemistry, Sogang University, Seoul 04107, South Korea
| | - Bong June Sung
- Department of Chemistry, Sogang University, Seoul 04107, South Korea
| |
Collapse
|
23
|
Sun LZ, Li H, Xu X, Luo MB. Simulation study on the translocation of polyelectrolyte through conical nanopores. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:495101. [PMID: 30431017 DOI: 10.1088/1361-648x/aaeb19] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Experiments have suggested that the conical nanopore may be a reasonable sensor for the biopolymer analysis as it can provide high-resolution current signal. In this paper, we use Langevin dynamics simulation to study the translocation of charged polymer (polyelectrolyte) through three different conical nanopores, a single-conical nanopore with large entry and small exit (pore I), a single-conical nanopore with small entry and large exit (pore II), and a double-conical nanopore with the tip (narrowest place) at the middle (pore III). Simulation shows that the detailed translocation behaviors are of pore structure dependence. Pore I might be the most reasonable one for the polyelectrolyte analysis, especially at strong monomer-pore attraction, since it can efficiently reduce the polyelectrolyte speed at the tip. The simulation results are interpreted by the free energy profiles of the polyelectrolyte translocation through different pores and the time of individual monomer passing through the tips.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, People's Republic of China
| | | | | | | |
Collapse
|
24
|
Hsiao PY. Translocation of Charged Polymers through a Nanopore in Monovalent and Divalent Salt Solutions: A Scaling Study Exploring over the Entire Driving Force Regimes. Polymers (Basel) 2018; 10:E1229. [PMID: 30961154 PMCID: PMC6290626 DOI: 10.3390/polym10111229] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2018] [Revised: 10/31/2018] [Accepted: 11/02/2018] [Indexed: 11/16/2022] Open
Abstract
Langevin dynamics simulations are performed to study polyelectrolytes driven through a nanopore in monovalent and divalent salt solutions. The driving electric field E is applied inside the pore, and the strength is varied to cover the four characteristic force regimes depicted by a rederived scaling theory, namely the unbiased (UB) regime, the weakly-driven (WD) regime, the strongly-driven trumpet (SD(T)) regime and the strongly-driven isoflux (SD(I)) regime. By changing the chain length N, the mean translocation time is studied under the scaling form 〈 τ 〉 ∼ N α E - δ . The exponents α and δ are calculated in each force regime for the two studied salt cases. Both of them are found to vary with E and N and, hence, are not universal in the parameter's space. We further investigate the diffusion behavior of translocation. The subdiffusion exponent γ p is extracted. The three essential exponents ν s , q, z p are then obtained from the simulations. Together with γ p , the validness of the scaling theory is verified. Through a comparison with experiments, the location of a usual experimental condition on the scaling plot is pinpointed.
Collapse
Affiliation(s)
- Pai-Yi Hsiao
- Department of Engineering and System Science, National Tsing Hua University, Hsinchu 30013, Taiwan.
- Institute of Nuclear Engineering and Science, National Tsing Hua University, Hsinchu 30013, Taiwan.
| |
Collapse
|
25
|
Sun LZ, Luo MB, Cao WP, Li H. Theoretical study on the polymer translocation into an attractive sphere. J Chem Phys 2018; 149:024901. [PMID: 30007381 DOI: 10.1063/1.5025609] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
We report a non-sampling model, combining the blob method with the standard lattice-based approximation, to calculate the free energy for the polymer translocation into an attractive sphere (i.e., spherical confined trans side) through a small pore. The translocation time is then calculated by the Fokker-Planck equation based on the free energy profile. There is a competition between the confinement effect of the sphere and the polymer-sphere attraction. The translocation time is increased due to the confinement effect of the sphere, whereas it is reduced by the polymer-sphere attraction. The two effects offset each other at a special polymer-sphere attraction which is dependent on the sphere size, the polymer length, and the driving force. Moreover, the entire translocation process can be divided into an uncrowded stage where the polymer does not experience the confinement effect of the sphere and a crowded stage where the polymer is confined by the sphere. At the critical sphere radius, the durations of the two (uncrowded and crowded) stages are the same. The critical sphere radius R* has a scaling relation with the polymer length N as R* ∼ Nβ. The calculation results show that the current model can effectively treat the translocation of a three-dimensional self-avoiding polymer into the spherical confined trans side.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China
| | - Meng-Bo Luo
- Department of Physics, Zhejiang University, Hangzhou 310027, China
| | - Wei-Ping Cao
- Institute of Optoelectronic Technology, Lishui University, Lishui 323000, China
| | - Haibin Li
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China
| |
Collapse
|
26
|
Pastoriza-Gallego M, Thiébot B, Bacri L, Auvray L, Pelta J. Dynamics of a polyelectrolyte through aerolysin channel as a function of applied voltage and concentration ⋆. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2018; 41:58. [PMID: 29748865 DOI: 10.1140/epje/i2018-11661-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/23/2018] [Accepted: 04/12/2018] [Indexed: 06/08/2023]
Abstract
We describe the behaviour of a polyelectrolyte in confined geometry. The transport of a polyelectrolyte, dextran sulfate, through a recombinant protein channel, aerolysin, inserted into a planar lipid bilayer is studied as a function of applied voltage and polyelectrolyte concentration and chain length. The aerolysin pore has a weak geometry asymmetry, a high number of charged residues and the polyelectrolyte is strongly negatively charged. The resulting current blockades were characterized by short and long dwelling times. Their frequency varies exponentially as a function of applied voltage and linearly as a function of polyelectrolyte concentration. The long blockade duration decreases exponentially when the electrical force increases. The ratio of the population of short events to the one of long events decreases when the applied voltage increases and displays an exponential variation. The long residence time increases with the polyelectrolyte chain length. We measure a reduction of the effective charge of the polyelectrolyte at the pore entry and inside the channel. For a fixed applied voltage, + / - 100 mV, at both sides of the protein pore entrance, the events frequency is similar as a function of dextran sulfate concentration. The mean blockade durations are independent of polyelectrolyte concentration and are similar for both entrances of the pore and remain constant as a function of the electrical force.
Collapse
Affiliation(s)
| | - Bénédicte Thiébot
- LAMBE UMR 8587, Université Cergy-Pontoise, Université Paris-Seine, 95302, Cergy-Pontoise, France
| | - Laurent Bacri
- LAMBE UMR 8587, Université Evry, CNRS, CEA, Université Paris-Saclay, 91025, Evry, France
| | - Loïc Auvray
- Laboratoire Matière et Systèmes Complexes, Université Paris Diderot, F-75205, Paris, France
| | - Juan Pelta
- LAMBE UMR 8587, Université Evry, CNRS, CEA, Université Paris-Saclay, 91025, Evry, France.
| |
Collapse
|