Single molecule electrophoresis of star polymers through nanopores: Simulations

Single-molecule resolution of macromolecules with nanopore devices

Nanopore-based electrical detection technology holds single-molecule resolution and combines the advantages of high sensitivity, high throughput, rapid analysis, and label-free detection. It is widely applied in the determination of organic and biological macromolecules, small molecules, and nanomaterials, as well as in nucleic acid and protein sequencing. There are a wide variety of organic polymers and biopolymers, and their chemical structures, and conformation in solution directly affect their ensemble properties. Currently, there is limited approach available for the analysis of single-molecule conformation and self-assembled topologies of polymers, dendrimers and biopolymers. Nanopore single-molecule platform offers unique advantages over other sensing technologies, particularly in molecular size differentiation of macromolecules and complex conformation analysis. In this review, the classification of nanopore devices, including solid-state nanopores (SSNs), biological nanopores, and hybrid nanopores is introduced. The recent developments and applications of nanopore devices are summarized, with a focus on the applications of nanopore platform in the resolution of the structures of synthetic polymer, including dendritic, star-shaped, block copolymers, as well as biopolymers, including polysaccharides, nucleic acids and proteins. The future prospects of nanopore sensing technique are ultimately discussed.

Langevin Dynamics Study on the Driven Translocation of Polymer Chains with a Hairpin Structure

The hairpin structure is a common and fundamental secondary structure in macromolecules. In this work, the process of the translocation of a model polymer chain with a hairpin structure is studied using Langevin dynamics simulations. The simulation results show that the dynamics of hairpin polymer translocation through a nanopore are influenced by the hairpin structure. Hairpin polymers can be classified into three categories, namely, linear-like, unsteady hairpin, and steady hairpin, according to the interaction with the stem structure. The translocation behavior of linear-like polymers is similar to that of a linear polymer chain. The time taken for the translocation of unsteady hairpin polymers is longer than that for a linear chain because it takes a long time to unfold the hairpin structure, and this time increases with stem interaction and decreases with the driving force. The translocation of steady hairpin polymers is distinct, especially under a weak driving force; the difficulty of unfolding the hairpin structure leads to a low translocation probability and a short translocation time. The translocation behavior of hairpin polymers can be explained by the theory of the free-energy landscape.

Escape dynamics of active ring polymers in a cylindrical nanochannel

We explore the escape dynamics of active ring polymers confined in a cylindrical nanochannel using Brownian dynamics. Our simulation results show that the escape time decreases with the increase of the Péclet number, which is not noticeable between the two stages of the escape process, based on whether the center of mass of the polymer is inside or outside the nanochannel. However, the monomer motion trajectory of the active polymer is very different from that of the passive polymer, similar to the snake-like motion with uniform velocity. The passive polymer, however, is in constant fugitive motion with increased velocity at the tail end of the escape. Our work is vital for understanding the escape dynamics of active ring polymers in the confined nanochannel, which provides new perspectives on their characterization and analysis.

Constraints on Knot Insertion, Not Internal Jamming, Control Polycatenane Translocation Dynamics through Crystalline Pores

The translocation of polymers through pores and channels is an archetypal process in biology and is widely studied and exploited for applications in bio- and nanotechnology. In recent times, the translocation of polymers of various different topologies has been studied both experimentally and by computer simulation. However, in some cases, a clear understanding of the precise mechanisms that drive their translocation dynamics can be challenging to derive. Experimental methods are able to provide statistical details of polymer translocation, but computer simulations are uniquely placed to uncover a finer level of mechanistic understanding. In this work, we use high-throughput molecular simulations to reveal the importance that knot insertion rates play in controlling translocation dynamics in the small pore limit, where unexpected nonpower law behavior emerges. This work both provides new predictive understanding of polycatenane translocation and shows the importance of carefully considering the role of the definition of translocation itself.

Translocation, Rejection and Trapping of Polyampholytes

Polyampholytes (PA) are a special class of polymers comprising both positive and negative monomers along their sequence. Most proteins have positive and negative residues and are PAs. Proteins have a well-defined sequence while synthetic PAs have a random charge sequence. We investigated the translocation behavior of random polyampholyte chains through a pore under the action of an electric field by means of Monte Carlo simulations. The simulations incorporated a realistic translocation potential profile along an extended asymmetric pore and translocation was studied for both directions of engagement. The study was conducted from the perspective of statistics for disordered systems. The translocation behavior (translocation vs. rejection) was recorded for all 220 sequences comprised of N = 20 charged monomers. The results were compared with those for 107 random sequences of N = 40 to better demonstrate asymptotic laws. At early times, rejection was mainly controlled by the charge sequence of the head part, but late translocation/rejection was governed by the escape from a trapped state over an antagonistic barrier built up along the sequence. The probability distribution of translocation times from all successful attempts revealed a power-law tail. At finite times, there was a population of trapped sequences that relaxed very slowly (logarithmically) with time. If a subensemble of sequences with prescribed net charge was considered the power-law decay was steeper for a more favorable net charge. Our findings were rationalized by theoretical arguments developed for long chains. We also provided operational criteria for the translocation behavior of a sequence, explaining the selection by the translocation process. From the perspective of protein translocation, our findings can help rationalize the behavior of intrinsically disordered proteins (IDPs), which can be modeled as polyampholytes. Most IDP sequences have a strong net charge favoring translocation. Even for sequences with those large net charges, the translocation times remained very dispersed and the translocation was highly sequence-selective.

Polymer Translocation Time

The force- and flow-induced translocation processes of linear and ring polymers are studied using a combination of multiparticle collision dynamics and molecular dynamics, focusing on the behavior of the polymer translocation time. We compare the force- and flow-induced translocations of linear and ring polymers. It is found that when the translocation time (τ*) is characterized by scaling exponents, δ, δ', and α, via the relations τ* ∼ f^δN^α and τ* ∼ J^δ'N^α, the scaling exponents are not constants. For long chains tested, α = 1.0 for both force- and flow-induced translocations. The difference between the force- and flow-induced translocations stems from different monomer crowding effects due to distinct flow patterns outside the channel. Furthermore, general relations for polymer translocation time are derived for these two translocation processes, which are in good agreement with the simulation results. Our results provide clear molecular pictures for the force- and flow-induced translocations, which shed light on the understanding of translocation dynamics and provide guidance for practical applications such as molecular sequencing and ultrafiltration.

Translocation of a looped polymer threading through a nanopore

Recent experiments reported that the complicated translocation dynamics of a looped DNA chain through a nanopore can be detected by ionic current blockade profiles. Inspired by the experimental results, we systematically study the translocation dynamics of a looped polymer, formed by three building blocks of a loop in the middle and two tails of the same length connected with the loop, by using Langevin dynamics simulations. Based on two entering modes (tail-leading and loop-leading) and three translocation orders (loop-tail-tail, tail-loop-tail, and tail-tail-loop), the translocation of the looped polymer is classified into six translocation pathways, corresponding to different current blockade profiles. The probabilities of the six translocation pathways are dependent on the loop length, polymer length, and pore radius. Moreover, the translocation times of the entire polymer and the loop are investigated. We find that the two translocation times show different dependencies on the translocation pathways and on the lengths of the loop and the entire polymer.

Polyelectrolyte Translocation through a Tortuous Nanopore

Although nanopores have shown tremendous promise for use in DNA sequencing, the rate of translocation through most pores studied previously is too rapid for the genetic information to be read accurately. In this study, dissipative particle dynamics simulations were employed to investigate the feasibility of using tortuous nanopores to control the rate of polyelectrolyte translocation. Unlike many previous studies, our simulation method incorporates the effects of hydrodynamic and electrostatic interactions and the spatial variation of electric field strength. The average translocation time, ⟨τ⟩, increases with the pore length and tortuosity but decreases as the pore width increases. For the longest pore investigated, the introduction of tortuosity results in ⟨τ⟩ increasing by as much as 187% as compared to a straight pore. The temporal variation of bond tension indicates that slower translocation in tortuous nanopores is caused by inhibition of tension propagation.

Flow-Induced Translocation of Star Polymers through a Nanopore

The flow-induced translocation of star polymers through a cylindrical nanopore has been studied using dissipative particle dynamics (DPD) simulations. The number of arms, f, was varied with the total number of monomers, N, kept constant. The effect of simulating the capture of the polymer into the pore upon the mean translocation time, <τ_t>, has been investigated by varying the chain's initial location. The results indicate that the incorporation of the capture process results in a reduction of <τ_t> by up to 15%. This is because the chain's initial location affects the extent of its stretching along the flow direction during translocation. <τ_t> exhibits nonmonotonic variation with f, in agreement with recently reported results for electric field-driven translocation of star polymers. Its value is larger and shows greater variation with f when the solvent quality is better. For the same value of f, the capture occurs faster in a good solvent. In addition, <τ_t> is greater for a semiflexible chain than its flexible counterpart as the time required for the branch point to enter the nanopore is longer in the former case.

Translocation of Star Polyelectrolytes through a Nanopore

The electric field driven translocation of charged star polymers through a cylindrical nanopore has been studied using dissipative particle dynamics simulations. The critical field strength required to induce translocation depends on both the number of arms and the number of beads per arm. It may therefore be possible to separate star polyelectrolytes of different arm lengths using electric field driven translocation through a nanopore. The average translocation time exhibits nonmonotonic variation with the number of arms for good solvent conditions. During translocation, a star polymer with many arms is stretched along the pore axis to a lesser extent as compared to its linear counterpart. Unlike a linear chain that shows tension propagation with large tensions for bonds about to enter the pore, a star has the tensest bonds closest to the branch point whose connectivity to multiple arms raises difficulty for its entry and passage.