Anomalous Dynamics of Forced Translocation

We consider the passage of long polymers of length N through a hole in a membrane. If the process is slow, it is in principle possible to focus on the dynamics of the number of monomers s on one side of the membrane, assuming that the two segments are in equilibrium. The dynamics of s(t) in such a limit would be diffusive, with a mean translocation time scaling as N^2 in the absence of a force, and proportional to N when a force is applied. We demonstrate that the assumption of equilibrium must break down for sufficiently long polymers (more easily when forced), and provide lower bounds for the translocation time by comparison to unimpeded motion of the polymer. These lower bounds exceed the time scales calculated on the basis of equilibrium, and point to anomalous (sub-diffusive) character of translocation dynamics. This is explicitly verified by numerical simulations of the unforced translocation of a self-avoiding polymer. Forced translocation times are shown to strongly depend on the method by which the force is applied. In particular, pulling the polymer by the end leads to much longer times than when a chemical potential difference is applied across the membrane. The bounds in these cases grow as N^2 and N^{1+\nu}, respectively, where \nu is the exponent that relates the scaling of the radius of gyration to N. Our simulations demonstrate that the actual translocation times scale in the same manner as the bounds, although influenced by strong finite size effects which persist even for the longest polymers that we considered (N=512).Comment: 13 pages, RevTeX4, 16 eps figure

Forced Imbibition - a Tool for Determining Laplace Pressure, Drag Force and Slip Length in Capillary Filling Experiments

When a very thin capillary is inserted into a liquid, the liquid is sucked into it: this imbibition process is controlled by a balance of capillary and drag forces, which are hard to quantify experimentally, in particularly considering flow on the nanoscale. By computer experiments using a generic coarse-grained model, it is shown that an analysis of imbibition forced by a controllable external pressure quantifies relevant physical parameter such as the Laplace pressure, Darcy's permeability, effective pore radius, effective viscosity, dynamic contact angle and slip length of the fluid flowing into the pore. In determining all these parameters independently, the consistency of our analysis of such forced imbibition processes is demonstrated.Comment: 4 pages, 5 figure

Criteria and geological setting for the generic geothermal underground research laboratory, GEOLAB

High flow rate injection and related hydromechanical interaction are the most important factors in reservoir development of Enhanced Geothermal Systems (EGS). GeoLaB, a new generic geothermal underground research laboratory (URL), is proposed for controlled high flow rate experiments (CHFE) to address limited comprehension of coupled processes connected to EGS reservoir flow conditions. As analogue for typical EGS development, CHFE require specific hydromechanical conditions including a connected fracture network in crystalline basement rock, sufficient hydraulic fracture transmissivities, a strike-slip to normal faulting tectonic regime, controllable hydraulic boundary conditions, and hydrothermal alteration fracture fillings that improve conditions for hydromechanical interaction. With the aim to identify most appropriate areas for future site selection, four criteria have been established based on the EGS reference site of Soultz. Two URLs in crystalline basement worldwide approximate the requirements of a new generic GeoLaB and may be used for accompanying experimentation. Besides favourable geological, hydraulic, and stress conditions, the vicinity to long-term EGS production favours the southern Black Forest as potential region for GeoLaB. Therefore, an exemplary site assessment has been carried out at “Wilhelminenstollen” in the southern Black Forest (Germany). New remote sensing, hydrochemical, and geophysical analyses as well as reactivation potential, and stress modelling were added to complement existing geological and hydrogeological information. At this site, reactivation potential analysis reveals two local maxima prone for shear reactivation as strike-slip faults. The highest lineament density is observed for the N110°E strike direction that is associated with both slip and dilation tendency maxima. Clay minerals occur in fractures and the matrix. Local, partly water-bearing fractures, when partly filled with ore minerals, were connected to veins in the tunnel using shallow geophysical methods. Hydrochemical data reveal infiltration of the tunnel water from at least 500 m above the tunnel. The results suggest a crystalline basement with a fracture network that is regionally connected and water-conducting. Hydraulic conductivity in the southern Black Forest granite is estimated to amount to about 4.5·10−8 m s−1 at 500 m depth. The hydraulic boundary conditions exclude unknown drainage. Analyses of the influence of topography on orientation and magnitude of the maximum stress indicate a minimum overburden of about 500 m for regional reactivation to be valid. In conclusion, the southern Black Forest and in particular “Wilhelminenstollen” offers favourable condition for CHFE. Final decision on the GeoLaB site is to be drawn from forthcoming exploration wells

Conformational distributions of isolated myosin motor domains encode their mechanochemical properties

Myosin motor domains perform an extraordinary diversity of biological functions despite sharing a common mechanochemical cycle. Motors are adapted to their function, in part, by tuning the thermodynamics and kinetics of steps in this cycle. However, it remains unclear how sequence encodes these differences, since biochemically distinct motors often have nearly indistinguishable crystal structures. We hypothesized that sequences produce distinct biochemical phenotypes by modulating the relative probabilities of an ensemble of conformations primed for different functional roles. To test this hypothesis, we modeled the distribution of conformations for 12 myosin motor domains by building Markov state models (MSMs) from an unprecedented two milliseconds of all-atom, explicit-solvent molecular dynamics simulations. Comparing motors reveals shifts in the balance between nucleotide-favorable and nucleotide-unfavorable P-loop conformations that predict experimentally measured duty ratios and ADP release rates better than sequence or individual structures. This result demonstrates the power of an ensemble perspective for interrogating sequence-function relationships

Self-energy limited ion transport in sub-nanometer channels

The current-voltage characteristics of the alpha-Hemolysin protein pore during the passage of single-stranded DNA under varying ionic strength, C, are studied experimentally. We observe strong blockage of the current, weak super-linear growth of the current as a function of voltage, and a minimum of the current as a function of C. These observations are interpreted as the result of the ion electrostatic self-energy barrier originating from the large difference in the dielectric constants of water and the lipid bilayer. The dependence of DNA capture rate on C also agrees with our model.Comment: more experimental material is added. 4 pages, 7 figure

Chaperone-assisted translocation of a polymer through a nanopore

Using Langevin dynamics simulations, we investigate the dynamics of chaperone-assisted translocation of a flexible polymer through a nanopore. We find that increasing the binding energy $\epsilon$ between the chaperone and the chain and the chaperone concentration $N_c$ can greatly improve the translocation probability. Particularly, with increasing the chaperone concentration a maximum translocation probability is observed for weak binding. For a fixed chaperone concentration, the histogram of translocation time $\tau$ has a transition from long-tailed distribution to Gaussian distribution with increasing $\epsilon$. $\tau$ rapidly decreases and then almost saturates with increasing binding energy for short chain, however, it has a minimum for longer chains at lower chaperone concentration. We also show that $\tau$ has a minimum as a function of the chaperone concentration. For different $\epsilon$, a nonuniversal dependence of $\tau$ on the chain length $N$ is also observed. These results can be interpreted by characteristic entropic effects for flexible polymers induced by either crowding effect from high chaperone concentration or the intersegmental binding for the high binding energy.Comment: 10 pages, to appear in J. Am. Chem. So

Does Young's equation hold on the nanoscale? A Monte Carlo test for the binary Lennard-Jones fluid

When a phase-separated binary ($A+B$) mixture is exposed to a wall, that preferentially attracts one of the components, interfaces between A-rich and B-rich domains in general meet the wall making a contact angle $\theta$. Young's equation describes this angle in terms of a balance between the $A-B$ interfacial tension $\gamma_{AB}$ and the surface tensions $\gamma_{wA}$, $\gamma_{wB}$ between, respectively, the $A$- and $B$-rich phases and the wall, $\gamma _{AB} \cos \theta =\gamma_{wA}-\gamma_{wB}$. By Monte Carlo simulations of bridges, formed by one of the components in a binary Lennard-Jones liquid, connecting the two walls of a nanoscopic slit pore, $\theta$ is estimated from the inclination of the interfaces, as a function of the wall-fluid interaction strength. The information on the surface tensions $\gamma_{wA}$, $\gamma_{wB}$ are obtained independently from a new thermodynamic integration method, while $\gamma_{AB}$ is found from the finite-size scaling analysis of the concentration distribution function. We show that Young's equation describes the contact angles of the actual nanoscale interfaces for this model rather accurately and location of the (first order) wetting transition is estimated.Comment: 6 pages, 6 figure

Anomalous Dynamics of Translocation

We study the dynamics of the passage of a polymer through a membrane pore (translocation), focusing on the scaling properties with the number of monomers $N$. The natural coordinate for translocation is the number of monomers on one side of the hole at a given time. Commonly used models which assume Brownian dynamics for this variable predict a mean (unforced) passage time $\tau$ that scales as $N^2$, even in the presence of an entropic barrier. However, the time it takes for a free polymer to diffuse a distance of the order of its radius by Rouse dynamics scales with an exponent larger than 2, and this should provide a lower bound to the translocation time. To resolve this discrepancy, we perform numerical simulations with Rouse dynamics for both phantom (in space dimensions $d=1$ and 2), and self-avoiding (in $d=2$) chains. The results indicate that for large $N$, translocation times scale in the same manner as diffusion times, but with a larger prefactor that depends on the size of the hole. Such scaling implies anomalous dynamics for the translocation process. In particular, the fluctuations in the monomer number at the hole are predicted to be non-diffusive at short times, while the average pulling velocity of the polymer in the presence of a chemical potential difference is predicted to depend on $N$.Comment: 9 pages, 9 figures. Submitted to Physical Review

Translocation of structured polynucleotides through nanopores

We investigate theoretically the translocation of structured RNA/DNA molecules through narrow pores which allow single but not double strands to pass. The unzipping of basepaired regions within the molecules presents significant kinetic barriers for the translocation process. We show that this circumstance may be exploited to determine the full basepairing pattern of polynucleotides, including RNA pseudoknots. The crucial requirement is that the translocation dynamics (i.e., the length of the translocated molecular segment) needs to be recorded as a function of time with a spatial resolution of a few nucleotides. This could be achieved, for instance, by applying a mechanical driving force for translocation and recording force-extension curves (FEC's) with a device such as an atomic force microscope or optical tweezers. Our analysis suggests that with this added spatial resolution, nanopores could be transformed into a powerful experimental tool to study the folding of nucleic acids.Comment: 9 pages, 5 figure

Driven polymer translocation through a nanopore: a manifestation of anomalous diffusion

We study the translocation dynamics of a polymer chain threaded through a nanopore by an external force. By means of diverse methods (scaling arguments, fractional calculus and Monte Carlo simulation) we show that the relevant dynamic variable, the translocated number of segments $s(t)$, displays an {\em anomalous} diffusive behavior even in the {\em presence} of an external force. The anomalous dynamics of the translocation process is governed by the same universal exponent $\alpha = 2/(2\nu +2 - \gamma_1)$, where $\nu$ is the Flory exponent and $\gamma_1$ - the surface exponent, which was established recently for the case of non-driven polymer chain threading through a nanopore. A closed analytic expression for the probability distribution function $W(s, t)$, which follows from the relevant {\em fractional} Fokker - Planck equation, is derived in terms of the polymer chain length $N$ and the applied drag force $f$. It is found that the average translocation time scales as $\tau \propto f^{-1}N^{\frac{2}{\alpha} -1}$. Also the corresponding time dependent statistical moments, $\propto t^{\alpha}$ and $\propto t^{2\alpha}$ reveal unambiguously the anomalous nature of the translocation dynamics and permit direct measurement of $\alpha$ in experiments. These findings are tested and found to be in perfect agreement with extensive Monte Carlo (MC) simulations.Comment: 6 pages, 4 figures, accepted to Europhys. Lett; some references were supplemented; typos were correcte