Tug of war in a double-nanopore system

Two Detection Modes of Nanoslit Sensing Based on Planar Heterostructure of Graphene/Hexagonal Boron Nitride

Solid-state nanopore sequencing is now confronted with problems of stochastic pore clogging and too fast speed during the DNA permeation through a nanopore, although this technique is revolutionary with long readability and high efficiency. These two problems are related to controlling molecular transportation during sequencing. To control the DNA motion and identify the four bases, we propose nanoslit sensing based on the planar heterostructure of two-dimensional graphene and hexagonal boron nitride. Molecular dynamics simulations are performed on investigating the motion of DNA molecules on the heterostructure with a nanoslit sensor. Results show that the DNA molecules are confined within the hexagonal boron nitride (HBN) domain of the heterostructure. And the confinement effects of the heterostructure can be optimized by tailoring the stripe length. Besides, there are two ways of DNA permeation through nanoslits: the DNA can cross or translocate the nanoslit under applied voltages along the y and z directions. The two detection modes are named cross-slit and trans-slit, respectively. In both modes, the ionic current drops can be observed when the nanoslit is occupied by the DNA. And the ionic currents and dwell times can be simultaneously detected to identify the four different DNA bases. This study can shed light on the sensing mechanism based on the nanoslit sensor of a planar heterostructure and provide theoretical guidance on designing devices controlling molecular transportation during nanopore sequencing.

Discriminating protein tags on a dsDNA construct using a Dual Nanopore Device

We report Brownian dynamics simulation results with the specific goal to identify key parameters controlling the experimentally measurable characteristics of protein tags on a dsDNA construct translocating through a double nanopore setup. First, we validate the simulation scheme in silico by reproducing and explaining the physical origin of the asymmetric experimental dwell time distributions of the oligonucleotide flap markers on a 48 kbp long dsDNA at the left and the right pore. We study the effect of the electric field inside and beyond the pores, critical to discriminate the protein tags based on their effective charges and masses revealed through a generic power-law dependence of the average dwell time at each pore. The simulation protocols monitor piecewise dynamics at a sub-nanometer length scale and explain the disparate velocity using the concepts of nonequilibrium tension propagation theory. We further justify the model and the chosen simulation parameters by calculating the Péclet number which is in close agreement with the experiment. We demonstrate that our carefully chosen simulation strategies can serve as a powerful tool to discriminate different types of neutral and charged tags of different origins on a dsDNA construct in terms of their physical characteristics and can provide insights to increase both the efficiency and accuracy of an experimental dual-nanopore setup.

Ion Transport in pH-Regulated Double-Barreled Nanopores

Single-molecule detection and characterization with nanopores is a powerful technique that does not require labeling. Multinanopore systems, especially double nanopores, have attracted wide attention and have been applied in many fields. However, theoretical studies of electrokinetic ion transport in nanopores mainly focus on single nanopores. In this paper, for the first time, a theoretical study of pH-regulated double-barreled nanopores is conducted using three-dimensional Poisson-Nernst-Planck equations and Navier-Stokes equations. Four ionic species and the surface chemistry on the walls of the nanopores are included. The results demonstrate that the properties of the bulk salt solution significantly affect nanopore conductivity and ion transport phenomena in nanopores. There are two ion-enriched zones and two ion-depleted zones in double-barreled nanopores. Due to the symmetry of the double-barreled nanopore structure and surface charge density, there is no ionic rectification effect in double-barreled nanopores. The ion selectivity is similar to that of conventional single pH-regulated nanopores.

Multi-resolution simulation of DNA transport through large synthetic nanostructures

Modeling and simulation has become an invaluable partner in development of nanopore sensing systems. The key advantage of the nanopore sensing method - the ability to rapidly detect individual biomolecules as a transient reduction of the ionic current flowing through the nanopore - is also its key deficiency, as the current signal itself rarely provides direct information about the chemical structure of the biomolecule. Complementing experimental calibration of the nanopore sensor readout, coarse-grained and all-atom molecular dynamics simulations have been used extensively to characterize the nanopore translocation process and to connect the microscopic events taking place inside the nanopore to the experimentally measured ionic current blockades. Traditional coarse-grained simulations, however, lack the precision needed to predict ionic current blockades with atomic resolution whereas traditional all-atom simulations are limited by the length and time scales amenable to the method. Here, we describe a multi-resolution framework for modeling electric field-driven passage of DNA molecules and nanostructures through to-scale models of synthetic nanopore systems. We illustrate the method by simulating translocation of double-stranded DNA through a solid-state nanopore and a micron-scale slit, capture and translocation of single-stranded DNA in a double nanopore system, and modeling ionic current readout from a DNA origami nanostructure passage through a nanocapillary. We expect our multi-resolution simulation framework to aid development of the nanopore field by providing accurate, to-scale modeling capability to research laboratories that do not have access to leadership supercomputer facilities.

DNA barcode by flossing through a cylindrical nanopore

We report an accurate method to determine DNA barcodes from the dwell time measurement of protein tags (barcodes) along the DNA backbone using Brownian dynamics simulation of a model DNA and use a recursive theoretical scheme which improves the measurements to almost 100% accuracy. The heavier protein tags along the DNA backbone introduce a large speed variation in the chain that can be understood using the idea of non-equilibrium tension propagation theory. However, from an initial rough characterization of velocities into “fast” (nucleotides) and “slow” (protein tags) domains, we introduce a physically motivated interpolation scheme that enables us to determine the barcode velocities rather accurately. Our theoretical analysis of the motion of the DNA through a cylindrical nanopore opens up the possibility of its experimental realization and carries over to multi-nanopore devices used for barcoding.

We report a method for DNA barcoding from the dwell time measurement of protein tags (barcodes) along the DNA backbone using Brownian dynamics simulation of a model DNA and use a recursive scheme to improve the measurements to almost 100% accuracy.

DNA barcodes using a double nanopore system

The potential of a double nanopore system to determine DNA barcodes has been demonstrated experimentally. By carrying out Brownian dynamics simulation on a coarse-grained model DNA with protein tag (barcodes) at known locations along the chain backbone, we demonstrate that due to large variation of velocities of the chain segments between the tags, it is inevitable to under/overestimate the genetic lengths from the experimental current blockade and time of flight data. We demonstrate that it is the tension propagation along the chain's backbone that governs the motion of the entire chain and is the key element to explain the non uniformity and disparate velocities of the tags and DNA monomers under translocation that introduce errors in measurement of the length segments between protein tags. Using simulation data we further demonstrate that it is important to consider the dynamics of the entire chain and suggest methods to accurately decipher barcodes. We introduce and validate an interpolation scheme using simulation data for a broad distribution of tag separations and suggest how to implement the scheme experimentally.

Polymer escape through a three dimensional double-nanopore system

We study the escape dynamics of a double-stranded DNA (dsDNA) through an idealized double nanopore geometry subject to two equal and opposite forces (tug-of-war) using Brownian dynamics (BD) simulation. In addition to the geometrical restrictions imposed on the cocaptured dsDNA segment in between the pores, the presence of tug-of-war forces at each pore results in a variation of the local chain stiffness for the segment of the chain in between the pores, which increases the overall stiffness of the chain. We use the BD simulation results to understand how the intrinsic chain stiffness and the tug-of-war forces affect the escape dynamics by monitoring the local chain persistence length ℓ_p, the residence time of the individual monomers W(m) in the nanopores, and the chain length dependence of the escape time ⟨τ⟩ and its distribution. Finally, we generalize the scaling theory for the unbiased single nanopore translocation for a fully flexible chain for the escape of a semi-flexible chain through a double nanopore in the presence of tug-of-war forces. We establish that the stiffness dependent part of the escape time is approximately independent of the translocation mechanism so that ⟨τ⟩∼ℓ_p ^2/D+2, and therefore, the generalized escape time for a semi-flexible chain can be written as ⟨τ⟩=AN^αℓ_p ^2/D+2. We use the BD simulation results to compare the predictions of the scaling theory. Our numerical studies supplemented by scaling analysis provide fundamental insights to design new experiments where a dsDNA moves slowly through a series of graphene nanopores.

An efficient kinetic Monte Carlo to study analyte capture by a nanopore: transients, boundary conditions and time-dependent fields

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.