Multiscale coupling of molecular dynamics and hydrodynamics: application to DNA translocation through a nanopore

We investigate the dynamics of DNA translocation through a nanopore driven by an external force using Langevin dynamics simulations in two dimensions (2D) to study how the translocation dynamics depend on the details of the DNA sequences. We consider a coarse-grained model of DNA built from two bases $A$ and $C$, having different base-pore interactions, {\textit e.g.}, a strong (weak) attractive force between the pore and the base $A$ ($C$) inside the pore. From a series of s...

The effect of the microscopic structure of a pore on polymer translocation is studied using Langevin dynamics simulation, and the consequence of introducing patterned stickiness inside the pore is investigated. It is found that the translocation process is extremely sensitive to the detailed structure of such patterns with faster than exponential dependence of translocation times on the stickiness of the pore. The stochastic nature of the translocation process leads to discer...

The advent of solid state nanodevices allows for interrogating the physico-chemical properties of a polyelectrolyte chain by electrophoretically driving it through a nanopore. Salient dynamical aspects of the translocation process have been recently characterized by theoretical and computational studies of model polymer chains free from self-entanglement. However, sufficiently long equilibrated chains are necessarily knotted. The impact of such topological "defects" on the tr...

We consider the influence of electrostatic forces on driven translocation dynamics of a flexible polyelectrolyte being pulled through a nanopore by an external force on the head monomer. To this end, we augment the iso-flux tension propagation (IFTP) theory with electrostatics for a negatively charged biopolymer pulled through a nanopore embedded in a similarly charged anionic membrane. We show that for the realistic case such as a single-stranded DNA, the translocation dynam...

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 form...

We investigate the dynamics of polymer translocation through nanopores under external driving by 3D Langevin Dynamics simulations, focusing on the scaling of the average translocation time $\tau$ versus the length of the polymer, $\tau\sim N^{\alpha}$. For slow translocation, i.e., under low driving force and/or high friction, we find $\alpha \approx 1+\nu \approx 1.588$ where $\nu$ denotes the Flory exponent. In contrast, $\alpha\approx 1.37$ is observed for fast translocati...

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 d...

We investigate the dynamics of polymer translocation through a nanopore using two-dimensional Langevin dynamics simulations. In the absence of external driving force, we consider a polymer which is initially placed in the middle of the pore and study the escape time $\tau_e$ required for the polymer to completely exit the pore on either side. The distribution of the escape times is wide and has a long tail. We find that $\tau_e$ scales with the chain length $N$ as $\tau_e \si...

We demonstrate the use of multiscale polymer modelling to quantitatively predict DNA linear dichroism (LD) in shear flow. LD is the difference in absorption of light polarised along two perpendicular axes, and has long been applied to study biopolymer structure and drug-biopolymer interactions. As LD is orientation-dependent, the sample must be aligned in order to measure a signal. Shear flow via a Couette cell can generate the required orientation, however it is challenging ...

We compare a newly developed hybrid simulation method which combines classical molecular dynamics (MD) and computational fluid dynamics (CFD) to a simulation consisting only of molecular dynamics. The hybrid code is composed of three regions: a classical MD region, a continuum domain where the dynamical equations are solved by standard CFD methods, and an overlap domain where transport information from the other two domains is exchanged. The exchange of information in the ove...