Source author record

Aniket Bhattacharya

Aniket Bhattacharya appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

Discriminating protein tags on dsDNA constructs using a dual Nanopore device

We report a novel simulation strategy that enables us to identify key parameters controlling the experimentally measurable characteristics of structural protein tags on dsDNA construct translocating through a double nanopore setup. First, we validate the scheme in silico by reproducing and explaining the physical origin of the experimental dwell time distributions of the Streptavidin markers on a 48 kbp long dsDNA. These studies reveal the important differences in the characteristics of the protein tags compared to the dynamics of dsDNA segments, immediately providing clues on how to improve the measurement protocols to decipher the unknown genomic lengths accurately. Of particular importance is the in silico studies on the effect of electric field inside and beyond the pores which we find is critical to discriminate 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 enable to monitor piecewise dynamics of the individual monomers at a sub-nanometer length scale and provide an explanation of the disparate velocity variation from one tag to the other using the nonequilibrium tension propagation theory, - a key element to decipher genomic lengths accurately. We further justify the model and the chosen simulation parameters by calculating the Peclet number which is in close agreement with the experiment. Analysis of our simulation results from the CG model has the capability to refine the accuracy of the experimentally obtained genomic lengths and 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.

preprint2022arXiv

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 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 towards 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 evidences of the interpretation of the experimentally observed folds [M. Gershow et al., Nat. Nanotech. 2, 775 (2007) and M. Mihovilovic et al. Phys. Rev. Letts. 110, 028102 (2013)] in a solitary nanopore. Finally, we show a new finding, - that a charged tag attached at the $5^{\prime}$ end of the DNA enhances both the multi-scan rate as well as the uni-directional translocation ($5^{\prime} \rightarrow 3^{\prime}$) probability that would benefit the genomic barcoding and sequencing experiments.

preprint2021arXiv

DNA Barcodes using 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.

preprint2021arXiv

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.

preprint2020arXiv

Polymer escape through a three dimensional Double-Nanopore System

We study escape dynamics of a double-stranded DNA (dsDNA) through an idealized double nanopore (DNP) 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 BD simulation results to understand how the intrinsic chain stiffness and the TOW forces affect the escape dynamics by monitoring the local chain persistence length $\ell_p$, the residence time of the individual monomers $W(m)$ in the nanopores, and the chain length dependence of the escape time $\langle τ\rangle$ 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 DNP in presence of TOW forces. We establish that the stiffness dependent part of the escape time is approximately independent of the translocation mechanism so that $\langle τ\rangle \sim \ell_p^{2/D+2}$, and therefore the generalized escape time for a semi-flexible chain can be written as $\langle τ\rangle = AN^α\ell_p^{2/D+2}$. We use 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.

preprint2020arXiv

Tug-of-War in a Double-Nanopore System

We simulate a tug-of-war (TOW) scenario for a model double-stranded DNA threading through a double nanopore (DNP) system. The DNA, simultaneously captured at both pores is subject to two equal and opposite forces $-\vec{f}_L= \vec{f}_R$ (TOW), where $\vec{f}_L$ and $\vec{f}_R$ are the forces applied to the left and the right pore respectively. Even though the net force on the DNA polymer $Δ\vec{f}_{LR}=\vec{f}_L+ \vec{f}_R=0$, the mean first passage time (MFPT) $\langle τ\rangle$ depends on the magnitude of the TOW forces $ \left | f_L \right | = \left |f_R \right | = f_{LR}$. We qualitatively explain this dependence of $\langle τ\rangle$ on $f_{LR}$ from the known results for the single-pore translocation of a triblock copolymer. We demonstrate that the time of flight (TOF) of a monomer with index $m$ ($\langle τ_{LR}(m) \rangle$) from one pore to the other exhibits quasi-periodic structure commensurate with the distance between the pores $d_{LR}$. Finally, we study the case $Δ\vec{f}_{LR}=\vec{f}_L+ \vec{f}_R \ne 0$, and qualitatively reproduce the experimental result of the dependence of the MFPT on $Δ\vec{f}_{LR}$. For a moderate bias, the MFPT for the DNP system for a chain length $N$ follows the same scaling ansatz as that of for the single nanopore, $\langle τ\rangle = \left( AN^{1+ν} + η_{pore}N \right) \left(Δf_{LR}\right)^{-1}$, where $η_{pore}$ is the pore friction, which enables us to estimate $\langle τ\rangle $ for a long chain. Our Brownian dynamics simulation studies provide fundamental insights and valuable information about the details of the translocation speed obtained from $\langle τ_{LR}(m) \rangle$, and accuracy of the translation of the data obtained in the time-domain to units of genomic distances.

preprint2015arXiv

Deconvoluting chain heterogeneities from driven translocation through a nano-pore

We study translocation dynamics of a driven compressible semi-flexible chain consisting of alternate blocks of stiff ($S$) and flexible ($F$) segments of size $m$ and $n$ respectively for different chain length $N$ in two dimension (2D). The free parameters in the model are the bending rigidity $κ_b$ which controls the three body interaction term, the elastic constant $k_F$ in the FENE (bond) potential between successive monomers, as well as the segmental lengths $m$ and $n$ and the repeat unit $p$ ($N=m_pn_p$) and the solvent viscosity $γ$. We demonstrate that due to the change in entropic barrier and the inhomogeneous viscous drag on the chain backbone a variety of scenarios are possible amply manifested in the waiting time distribution of the translocating chain. These information can be deconvoluted to extract the mechanical properties of the chain at various length scales and thus can be used to nanopore based methods to probe bio-molecules, such as DNA, RNA and proteins.

preprint2014arXiv

Conformations, Transverse Fluctuations and Crossover Dynamics of a Semi-Flexible Chain in Two Dimensions

We present a unified scaling description for the dynamics of monomers of a semiflexible chain under good solvent condition in the free draining limit. We consider both the cases where the contour length $L$ is comparable to the persistence length $\ell_p$ and the case $L\gg \ell_p$. Our theory captures the early time monomer dynamics of a stiff chain characterized by $t^{3/4}$ dependence for the mean square displacement(MSD) of the monomers, but predicts a first crossover to the Rouse regime of $t^{2ν/{1+2ν}}$ for $τ_1 \sim \ell_p^3$, and a second crossover to the purely diffusive dynamics for the entire chain at $τ_2 \sim L^{5/2}$. We confirm the predictions of this scaling description by studying monomer dynamics of dilute solution of semi-flexible chains under good solvent conditions obtained from our Brownian dynamics (BD) simulation studies for a large choice of chain lengths with number of monomers per chain N = 16 - 2048 and persistence length $\ell_p = 1 - 500$ Lennard-Jones (LJ) units. These BD simulation results further confirm the absence of Gaussian regime for a 2d swollen chain from the slope of the plot of $\langle R_N^2 \rangle/2L \ell_p \sim L/\ell_p$ which around $L/\ell_p \sim 1$ changes suddenly from $\left(L/\ell_p \right) \rightarrow \left(L/\ell_p \right)^{0.5} $, also manifested in the power law decay for the bond autocorrelation function disproving the validity of the WLC in 2d. We further observe that the normalized transverse fluctuations of the semiflexible chains for different stiffness $\sqrt{\langle l_{\bot}^2\rangle}/L$ as a function of renormalized contour length $L/\ell_p$ collapse on the same master plot and exhibits power law scaling $\sqrt{\langle l_{\bot}^2\rangle}/L \sim (L/\ell_p)^η$ at extreme limits, where $η= 0.5$ for extremely stiff chains ($L/\ell_p \gg 1$), and $η= -0.25$ for fully flexible chains.

preprint2014arXiv

DNA confined in a two-dimensional strip geometry

Semiflexible polymers characterized by the contour length $L$ and persistent length $\ell_p$ confined in a spatial region $D$ have been described as a series of ``{\em spherical blobs}'' and ``{\em deflecting lines}'' by de Gennes and Odjik for $\ell_p < D$ and $\ell_p \gg D$ respectively. Recently new intermediate regimes ({\em extended de Gennes} and {\em Gauss-de Gennes}) have been investigated by Tree {\em et al.} [Phys. Rev. Lett. {\bf 110}, 208103 (2013)]. In this letter we derive scaling relations to characterize these transitions in terms of universal scaled fluctuations in $d$-dimension as a function of $L,\ell_p$, and $D$, and show that the Gauss-de Gennes regime is absent and extended de Gennes regime is vanishingly small for polymers confined in a 2D strip. We validate our claim by extensive Brownian dynamics (BD) simulation which also reveals that the prefactor $A$ used to describe the chain extension in the Odjik limit is independent of physical dimension $d$ and is the same as previously found by Yang {\em et al.}[Y. Yang, T. W. Burkhardt, G. Gompper, Phys. Rev. E {\bf 76}, 011804 (2007)]. Our studies are relevant for optical maps of DNA stretched inside a nano-strip.

preprint2013arXiv

Universal monomer dynamics of a two dimensional semi-flexible chain

We present a unified scaling theory for the dynamics of monomers for dilute solutions of semiflexible polymers under good solvent conditions in the free draining limit. Our theory encompasses the well-known regimes of mean square displacements (MSDs) of stiff chains growing like t^{3/4} with time due to bending motions, and the Rouse-like regime t^{2 ν/ (1+ 2ν)} where νis the Flory exponent describing the radius R of a swollen flexible coil. We identify how the prefactors of these laws scale with the persistence length l_p, and show that a crossover from stiff to flexible behavior occurs at a MSD of order l^2_p (at a time proportional to l^3_p). A second crossover (to diffusive motion) occurs when the MSD is of order R^2. Large scale Molecular Dynamics simulations of a bead-spring model with a bond bending potential (allowing to vary l_p from 1 to 200 Lennard-Jones units) provide compelling evidence for the theory, in D=2 dimensions where ν=3/4. Our results should be valuable for understanding the dynamics of DNA (and other semiflexible biopolymers) adsorbed on substrates.

preprint2012arXiv

Influence of non-universal effects on dynamical scaling in driven polymer translocation

We study the dynamics of driven polymer translocation using both molecular dynamics (MD) simulations and a theoretical model based on the non-equilibrium tension propagation on the {\it cis} side subchain. We present theoretical and numerical evidence that the non-universal behavior observed in experiments and simulations are due to finite chain length effects that persist well beyond the relevant experimental and simulation regimes. In particular, we consider the influence of the pore-polymer interactions and show that they give a major contribution to the non-universal effects. In addition, we present comparisons between the theory and MD simulations for several quantities, showing extremely good agreement in the relevant parameter regimes. Finally, we discuss the potential limitations of the present theories.

preprint2012arXiv

Scaling theory of driven polymer translocation

We present a theoretical argument to derive a scaling law between the mean translocation time $τ$ and the chain length $N$ for driven polymer translocation. This scaling law explicitly takes into account the pore-polymer interactions, which appear as a correction term to asymptotic scaling and are responsible for the dominant finite size effects in the process. By eliminating the correction-to-scaling term we introduce a rescaled translocation time and show, by employing both the Brownian Dynamics Tension Propagation theory [Ikonen {\it et al.}, Phys. Rev. E {\bf 85}, 051803 (2012)] and molecular dynamics simulations that the rescaled exponent reaches the asymptotic limit in a range of chain lengths that is easily accessible to simulations and experiments. The rescaling procedure can also be used to quantitatively estimate the magnitude of the pore-polymer interaction from simulations or experimental data. Finally, we also consider the case of driven translocation with hydrodynamic interactions (HIs). We show that by augmenting the BDTP theory with HIs one reaches a good agreement between the theory and previous simulation results found in the literature. Our results suggest that the scaling relation between $τ$ and $N$ is retained even in this case.

preprint2012arXiv

Unifying model of driven polymer translocation

We present a Brownian dynamics model of driven polymer translocation, in which non-equilibrium memory effects arising from tension propagation (TP) along the cis side subchain are incorporated as a time-dependent friction. To solve the effective friction, we develop a finite chain length TP formalism, expanding on the work of Sakaue [Sakaue, PRE 76, 021803 (2007)]. The model, solved numerically, yields results in excellent agreement with molecular dynamics simulations in a wide range of parameters. Our results show that non-equilibrium TP along the cis side subchain dominates the dynamics of driven translocation. In addition, the model explains the different scaling of translocation time w.r.t chain length observed both in experiments and simulations as a combined effect of finite chain length and pore-polymer interactions.

preprint2010arXiv

Out of Equilibrium Characteristics of a Forced Translocating Chain through a Nanopore

Polymer translocation through a nano-pore in a thin membrane is studied using a coarse-grained bead-spring model and Langevin dynamics simulation with a particular emphasis to explore out of equilibrium characteristics of the translocating chain. We analyze the out of equilibrium chain conformations both at the $cis$ and the $trans$ side separately either as a function of the time during the translocation process or as as function of the monomer index $m$ inside the pore. A detailed picture of translocation emerges by monitoring the center of mass of the translocating chain, longitudinal and transverse components of the gyration radii and the end to end vector. We observe that polymer configurations at the $cis$ side are distinctly different from those at the $trans$ side. During the translocation, and immediately afterwards, the chain is clearly out of equilibrium, as different parts of the chain are characterized by a series of effective Flory exponents. We further notice that immediately after the translocation the last set of beads that have just translocated take a relatively compact structure compared to the first set of beads that translocated earlier, and the chain immediately after translocation is described by an effective Flory exponent $0.45 \pm 0.01$. The analysis of these results is further strengthened by looking at the conformations of chain segments of equal length as they cross from the $cis$ to the $trans$ side, We discuss implications of these results to the theoretical estimates and numerical simulation studies of the translocation exponent reported by various groups.

preprint2010arXiv

Polymer Translocation Induced by a Bad Solvent

We employ 3D Langevin Dynamics simulations to study the dynamics of polymer chains translocating through a nanopore in presence of asymmetric solvent conditions. Initially a large fraction ($>$ 50%) of the chain is placed at the \textit{cis} side in a good solvent while the $trans$ segments are placed in a bad solvent that causes the chain to collapse and promotes translocation from the $cis$ to the $trans$ side. In particular, we study the ratcheting effect of a globule formed at the \textit{trans} side created by the translocated segment, and how this ratchet drives the system towards faster translocation. Unlike in the case of unbiased or externally forced translocation where the mean first passage time $\langle τ\rangle $ is often characterized by algebraic scaling as a function of the chain length $N$ with a single scaling exponent $α$, and the histogram of the mean first passage time $P(τ/\langleτ\rangle)$ exhibits scaling, we find that scaling is not well obeyed. For relatively long chains we find $\langle τ\rangle \sim N^α$ where $α\approx 1$ for $\varepsilon/k_{B}T > 1$. In this limit, we also find that translocation proceeds with a nearly constant velocity of the individual beads(monomers), which is attributed to the coiling of the globule. We provide an approximate theory assuming rotat ional motion restricted on a 2D disc to demonstrate that there is a crossover from diffusive behavior of the center of mass for short chains to a single file translocation for long chains, where the average translocation time scales linearly with the chain length $N$.

preprint2008arXiv

Scaling exponents of Forced Polymer Translocation through a nano-pore

We investigate several scaling properties of a translocating homopolymer through a thin pore driven by an external field present inside the pore only using Langevin Dynamics (LD) simulation in three dimension (3D). Specifically motivated by several recent theoretical and numerical studies that are apparently at odds with each other, we determine the chain length dependence of the scaling exponents of the average translocation time, the average velocity of the center of mass, $<v_{CM}>$, the effective radius of gyration during the translocation process, and the scaling exponent of the translocation coordinate ($s$-coordinate) as a function of the translocation time. We further discuss the possibility that in the case of driven translocation the finite pore size and its geometry could be responsible that the veclocity scaling exponent is less than unity and discuss the dependence of the scaling exponents on the pore geometry for the range of $N$ studied here.