High-Fidelity Capture, Threading, and Infinite-Depth Sequencing of Single DNA Molecules with a Double-Nanopore System

Machine learning empowered next generation DNA sequencing: perspective and prospectus

The pursuit of ultra-rapid, cost-effective, and accurate DNA sequencing is a highly sought after aspect of personalized medicine development. With recent advancements, mainstream machine learning (ML) algorithms hold immense promise for high throughput DNA sequencing at the single nucleotide level. While ML has revolutionized multiple domains of nanoscience and nanotechnology, its implementation in DNA sequencing is still in its preliminary stages. ML-aided DNA sequencing is especially appealing, as ML has the potential to decipher complex patterns and extract knowledge from complex datasets. Herein, we present a holistic framework of ML-aided next-generation DNA sequencing with domain knowledge to set directions toward the development of artificially intelligent DNA sequencers. This perspective focuses on the current state-of-the-art ML-aided DNA sequencing, exploring the opportunities as well as the future challenges in this field. In addition, we provide our personal viewpoints on the critical issues that require attention in the context of ML-aided DNA sequencing.

Nanopore DNA sequencing technologies and their applications towards single-molecule proteomics

Sequencing of nucleic acids with nanopores has emerged as a powerful tool offering rapid readout, high accuracy, low cost and portability. This label-free method for sequencing at the single-molecule level is an achievement on its own. However, nanopores also show promise for the technologically even more challenging sequencing of polypeptides, something that could considerably benefit biological discovery, clinical diagnostics and homeland security, as current techniques lack portability and speed. Here we survey the biochemical innovations underpinning commercial and academic nanopore DNA/RNA sequencing techniques, and explore how these advances can fuel developments in future protein sequencing with nanopores.

DNA Barcodes Using a Dual Nanopore Device

We report a novel method based on the current blockade (CB) characteristics obtained from a dual nanopore device that can determine DNA barcodes with near-perfect accuracy using a Brownian dynamics simulation strategy. The method supersedes our previously reported velocity correction algorithm (S. Seth and A. Bhattacharya, RSC Advances, 11:20781-20787, 2021), taking advantage of the better measurement of the time-of-flight (TOF) protocol offered by the dual nanopore setup. We demonstrate the efficacy of the method by comparing our simulation data from a coarse-grained model of a polymer chain consisting of 2048 excluded volume beads of diameter 𝜎 = 24 bp using with those obtained from experimental CB data from a 48,500 bp λ-phage DNA, providing a 48500 2400 ≅ 24 base pair resolution in simulation. The simulation time scale is compared to the experimental time scale by matching the simulated time-of-flight (TOF) velocity distributions with those obtained experimentally (Rand et al., ACS Nano, 16:5258-5273, 2022). We then use the evolving coordinates of the dsDNA and the molecular features to reconstruct the current blockade characteristics on the fly using a volumetric model based on the effective van der Waal radii of the species inside and in the immediate vicinity of the pore. Our BD simulation mimics the control-zoom-in-logic to understand the origin of the TOF distributions due to the relaxation of the out-of-equilibrium conformations followed by a reversal of the electric fields. The simulation algorithm is quite general and can be applied to differentiate DNA barcodes from different species.

DNA Carrier-Assisted Molecular Ping-Pong in an Asymmetric Nanopore

Nanopore analysis relies on ensemble averaging of translocation signals obtained from numerous molecules, requiring a relatively high sample concentration and a long turnaround time from the sample to results. The recapture and subsequent re-reading of the same molecule is a promising alternative that enriches the signal information from a single molecule. Here, we describe how an asymmetric nanopore improves molecular ping-pong by promoting the recapture of the molecule in the trans reservoir. We also demonstrate that the molecular recapture could be improved by linking the target molecule to a long DNA carrier to reduce the diffusion, thereby achieving over 100 recapture events. Using this ping-pong methodology, we demonstrate its use in accurately resolving nanostructure motifs along a DNA scaffold through repeated detection. Our method offers novel insights into the control of DNA polymer dynamics within nanopore confinement and opens avenues for the development of a high-fidelity DNA detection platform.

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

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.

Nonequilibrium Thermodynamics of DNA Nanopore Unzipping

Using theory and simulations, we carried out a first systematic characterization of DNA unzipping via nanopore translocation. Starting from partially unzipped states, we found three dynamical regimes depending on the applied force f: (i) heterogeneous DNA retraction and rezipping (f<17  pN), (ii) normal (17  pN<f<60  pN), and (iii) anomalous (f>60  pN) drift-diffusive behavior. We show that the normal drift-diffusion regime can be effectively modeled as a one-dimensional stochastic process in a tilted periodic potential. We use the theory of stochastic processes to recover the potential from nonequilibrium unzipping trajectories and show that it corresponds to the free-energy landscape for single-base-pair unzipping. Applying this general approach to other single-molecule systems with periodic potentials ought to yield detailed free-energy landscapes from out-of-equilibrium trajectories.

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.

How capture affects polymer translocation in a solitary nanopore

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.

Atomistic Simulation of Molecules Interacting with Biological Nanopores: From Current Understanding to Future Directions

Biological nanopores have been at the focus of numerous studies due to their role in many biological processes as well as their (prospective) technological applications. Among many other topics, recent studies on nanopores have addressed two key areas: antibiotic permeation through bacterial channels and sensing of analytes. Although the two areas are quite far apart in terms of their objectives, in both cases atomistic simulations attempt to understand the solute dynamics and the solute-protein interactions within the channel lumen. While decades of studies on various channels have culminated in an improved understanding of the key molecular factors and led to practical applications in some cases, successful utilization is limited. In this Perspective we summarize recent progress in understanding key issues in molecular simulations of antibiotic translocation and in the development of nanopore sensors. Moreover, we comment on possible advancements in computational algorithms that can potentially resolve some of the issues.

Single-molecule biophysics experiments in silico: Toward a physical model of a replisome

The interpretation of single-molecule experiments is frequently aided by computational modeling of biomolecular dynamics. The growth of computing power and ongoing validation of computational models suggest that it soon may be possible to replace some experiments outright with computational mimics. Here, we offer a blueprint for performing single-molecule studies in silico using a DNA-binding protein as a test bed. We demonstrate how atomistic simulations, typically limited to sub-millisecond durations and zeptoliter volumes, can guide development of a coarse-grained model for use in simulations that mimic single-molecule experiments. We apply the model to recapitulate, in silico, force-extension characterization of protein binding to single-stranded DNA and protein and DNA replacement assays, providing a detailed portrait of the underlying mechanics. Finally, we use the model to simulate the trombone loop of a replication fork, a large complex of proteins and DNA.

•

Coarse-grained model derived from all-atom simulation recapitulates experiments

•

Model reproduces the elastic response to force and exchange dynamics

•

Model reveals structure of intermediate states usually inaccessible to experiment

•

Model applied to viral replisome with trombone loop containing tens of SSB proteins

Single Molecule DNA Analysis Based on Atomic-Controllable Nanopores in Covalent Organic Frameworks

We explored the application of two-dimensional covalent organic frameworks (2D COFs) in single molecule DNA analysis. Two ultrathin COF nanosheets were exfoliated with pore sizes of 1.1 nm (COF-1.1) and 1.3 nm (COF-1.3) and covered closely on a quartz nanopipette with an orifice of 20 ± 5 nm. COF nanopores exhibited high size selectivity for fluorescent dyes and DNA molecules. The transport of long (calf thymus DNA) and short (DNA-80) DNA molecules through the COF nanopores was studied. Because of the strong interaction between DNA bases and the organic backbones of COFs, the DNA-80 was transported through the COF-1.1 nanopore at a speed of 270 μs/base, which is the slowest speed ever observed compared with 2D inorganic nanomaterials. This study shows that the COF nanosheet can work individually as a nanopore monomer with controllable pore size like its biological counterparts.

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.

Nanoconfinement Effect for Signal Amplification in Electrochemical Analysis and Sensing

Owing to the urgent need for electrochemical analysis and sensing of trace target molecules in various fields such as medical diagnosis, agriculture and food safety, and environmental monitoring, signal amplification is key to promoting analysis and sensing performance. The nanoconfinement effect, derived from nanoconfined spaces and interfaces with sizes approaching those of target molecules, has witnessed rapid development for ultra-sensitive analyzing and sensing. In this review, the two main types of nanoconfinement systems - confined nanochannels and planes - are assessed and recent progress is highlighted. The merits of each nanoconfinement system, the nanoconfinement effect mechanisms, and applications for electrochemical analysis and sensing are summarized and discussed. This review aims to help deepen the understanding of nanoconfinement devices and their effects in order to develop new analysis and sensing applications for researchers in various fields.

A general theory of polymer ejection tested in a quasi two-dimensional space

A general ejection theory of polymer is developed in a two- and three-dimensional space. A polymer is confined initially in a cavity and ejects spontaneously to the outer space through a nanopore channel without the help of any external stimulus. A reflective wall boundary is set at the pore entrance to prevent the falling of the head monomer of chain into the cavity. Three stages are distinguished in a process: (1) an entering stage, in which the head monomer enters the pore to search for a way to traverse the pore channel, (2) a main ejection stage, in which the chain body is transported from the cavity to the outer space, (3) a leaving stage, in which the tail monomer passes through and leaves the pore channel. Depending on the number of the monomers remaining in the cavity, the main ejection stage can be divided into the confined and the non-confined stages. The non-confined stage can be further split into the thermal escape and the entropic pulling stages. The Onsager’s variational principle is applied to derive the kinetics equation of ejection. The escape time is calculated from the corresponding Kramers’ escape problem. Extensive molecular dynamics simulations are then performed in a quasi two-dimensional space to verify the theory. The variation of the ejection speed is carefully examined. The decreasing behavior of the number of monomers in the cavity is studied in details. The scaling properties of the spending time at each processing stage are investigated systematically by varying the chain length, the cavity diameter, and the initial volume fraction of chain. The results of simulation support firmly the predictions of the theory, cross-checked in the studies of various topics. In combining with the previous investigations in the three-dimensional space, the generalized theory is very robust, able to explain the two seemly different phenomena, polymer ejection and polymer translocation, together under the same theoretical framework in the two space dimensions.

Protein Transport through Nanopores Illuminated by Long-Time-Scale Simulations

The transport of molecules through nanoscale confined space is relevant in biology, biosensing, and industrial filtration. Microscopically modeling transport through nanopores is required for a fundamental understanding and guiding engineering, but the short duration and low replica number of existing simulation approaches limit statistically relevant insight. Here we explore protein transport in nanopores with a high-throughput computational method that realistically simulates hundreds of up to seconds-long protein trajectories by combining Brownian dynamics and continuum simulation and integrating both driving forces of electroosmosis and electrophoresis. Ionic current traces are computed to enable experimental comparison. By examining three biological and synthetic nanopores, our study answers questions about the kinetics and mechanism of protein transport and additionally reveals insight that is inaccessible from experiments yet relevant for pore design. The discovery of extremely frequent unhindered passage can guide the improvement of biosensor pores to enhance desired biomolecular recognition by pore-tethered receptors. Similarly, experimentally invisible nontarget adsorption to pore walls highlights how to improve recently developed DNA nanopores. Our work can be expanded to pressure-driven flow to model industrial nanofiltration processes.

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.