Michel G. Gauthier

A Monte Carlo algorithm to study polymer translocation through nanopores. I. Theory and numerical approach

The process during which a polymer translocates through a nanopore depends on many physical parameters and fundamental mechanisms. We propose a new one-dimensional lattice Monte Carlo algorithm that integrates various effects such as the entropic forces acting on the subchains that are outside the channel, the external forces that are pulling the polymer through the pore, and the frictional effects that involve the chain and its environment. Our novel approach allows us to study the polymer as a single Brownian particle diffusing while subjected to a position-dependent force that includes both the external driving forces and the internal entropic bias. Frictional effects outside and inside the pore are also considered. This Monte Carlo method is much more efficient than other simulation methods, and it can be used to obtain scaling laws for various polymer translocation regimes. In this first part, we derive the model and describe a subtle numerical approach that gives exact results for both the escape probability and the mean translocation time (and higher moments of its distribution). The scaling laws obtained from this model will be presented and discussed in the second part of this series.

The theory of DNA separation by capillary electrophoresis

The Human Genome has been sequenced in large part owing to the invention of capillary electrophoresis. Although this technology has matured enough to allow such amazing achievements, the physical mechanisms at play during separation have yet to be completely understood and optimized. Recently, new separation regimes and new physical mechanisms have been investigated. The use of free-flow electrophoresis and new modes of pulsed-field electrophoresis have been suggested, while we have observed a shift towards single nucleotide polymorphism analysis and microchip technologies. A strong theoretical basis remains essential for the efficient development of new methods.

A Monte Carlo algorithm to study polymer translocation through nanopores. II. Scaling laws

In the first paper of this series, we developed a new one-dimensional Monte Carlo approach for the study of flexible chains that are translocating through a small channel. We also presented a numerical scheme that can be used to obtain exact values for both the escape times and the escape probabilities given an initial pore-polymer configuration. We now present and discuss the fundamental scaling behaviors predicted by this Monte Carlo method. Our most important result is the fact that, in the presence of an external bias E, we observe a change in the scaling law for the translocation time tau as function of the polymer length N: In the general expression tau approximately N(beta)E, the exponent changes from beta=1 for moderately long chains to beta=1+nu or beta=2nu for very large values of N (for Rouse and Zimm dynamics, respectively). We also observe an increase in the effective diffusion coefficient due to the presence of entropic pulling on unbiased polymer chains.

Molecular Dynamics simulation of a polymer chain translocating through a nanoscopic pore - Hydrodynamic interactions versus pore radius

Abstract.The detection of linear polymers translocating through a nanoscopic pore is a promising idea for the development of new DNA analysis techniques. However, the physics of constrained macromolecules and the fluid that surrounds them at the nanoscopic scale is still not well understood. In fact, many theoretical models of polymer translocation neglect both excluded-volume and hydrodynamic effects. We use Molecular Dynamics simulations with explicit solvent to study the impact of hydrodynamic interactions on the translocation time of a polymer. The translocation time τ that we examine is the unbiased (no charge on the chain and no driving force) escape time of a polymer that is initially placed halfway through a pore perforated in a monolayer wall. In particular, we look at the effect of increasing the pore radius when only a small number of fluid particles can be located in the pore as the polymer undergoes translocation, and we compare our results to the theoretical predictions of Chuang et al. (Phys. Rev. E 65, 011802 (2001)). We observe that the scaling of the translocation time varies from τ ∼ N11/5 to τ ∼ N9/5 as the pore size increases (N is the number of monomers that goes up to 31 monomers). However, the scaling of the polymer relaxation time remains consistent with the 9/5 power law for all pore radii.

Nondriven polymer translocation through a nanopore : Computational evidence that the escape and relaxation processes are coupled

Most of the theoretical models describing the translocation of a polymer chain through a nanopore use the hypothesis that the polymer is always relaxed during the complete process. In other words, models generally assume that the characteristic relaxation time of the chain is small enough compared to the translocation time that nonequilibrium molecular conformations can be ignored. In this paper, we use molecular dynamics simulations to directly test this hypothesis by looking at the escape time of unbiased polymer chains starting with different initial conditions. We find that the translocation process is not quite in equilibrium for the systems studied, even though the translocation time tau is about 10 times larger than the relaxation time tau{r}. Our most striking result is the observation that the last half of the chain escapes in less than approximately 12% of the total escape time, which implies that there is a large acceleration of the chain at the end of its escape from the channel.

Regulation of DNA Replication within the Immunoglobulin Heavy-Chain Locus During B Cell Commitment

The temporal order of replication of mammalian chromosomes appears to be linked to their functional organization, but the process that establishes and modifies this order during cell differentiation remains largely unknown. Here, we studied how the replication of the Igh locus initiates, progresses, and terminates in bone marrow pro-B cells undergoing B cell commitment. We show that many aspects of DNA replication can be quantitatively explained by a mechanism involving the stochastic firing of origins (across the S phase and the Igh locus) and extensive variations in their firing rate (along the locus). The firing rate of origins shows a high degree of coordination across Igh domains that span tens to hundreds of kilobases, a phenomenon not observed in simple eukaryotes. Differences in domain sizes and firing rates determine the temporal order of replication. During B cell commitment, the expression of the B-cell-specific factor Pax5 sharply alters the temporal order of replication by modifying the rate of origin firing within various Igh domains (particularly those containing Pax5 binding sites). We propose that, within the Igh CH-3′RR domain, Pax5 is responsible for both establishing and maintaining high rates of origin firing, mostly by controlling events downstream of the assembly of pre-replication complexes.

Sequence effects on the forced translocation of heteropolymers through a small channel

By using a recently developed Monte Carlo algorithm and an exact numerical method, we calculate the translocation probability and the average translocation time for charged heterogeneous polymers driven through a nanopore by an external electric field. The heteropolymer chains are composed of two types of monomers (A and B) which differ only in terms of their electric charge. We present an exhaustive study of chains composed of eight monomers by calculating the average translocation time associated with the 256 possible arrangements for various ratios of the monomer charges (lambda(A)lambda(B)) and electric field intensities E. We find that each sequence leads to a unique value of the translocation probability and time. We also show that the distribution of translocation times is strongly dependent on the two forces felt by the monomers ( approximately lambda(A)E and approximately lambda(B)E). Finally, we present results that highlight the effect of having repetitive patterns by studying the translocation times of various block copolymer structures for a very long chain composed of N=2(18) monomers (all with the same number of A and B monomers).

Equation of state, stability, anisotropy and nonlinear elasticity of diamond-cubic (ZB) silicon by phonon imaging at high pressure

Experimental phonon imaging in diamond anvils cell is demonstrated to be an adequate tool to extract the complete set of elastic constants of single-crystalline silicon up to the ZB ->beta-Sn transition (10 GPa). Contrary to what was commonly admitted, we demonstrate that the development of the strain-energy density in terms of strains cannot be stopped, for silicon, after the terms containing the third-order elastic constants. Nonlinear elasticity, degree of anisotropy and pressure-induced mechanical stability of the cubic silicon structure are thus revisited and investigated in more detail.

Exactly solvable Ogston model of gel electrophoresis. IX. Generalizing the lattice model to treat high field intensities

Traditionally, the Ogston regime is studied solely in the limit of low field intensities. This explains why the theoretical discussion has focused until now on the relative roles of the fractional volume available to the analyte and the subtleties of the gel architecture. Over the past several years, we have developed a lattice model of gel electrophoresis that has allowed us to revisit the fundamental assumptions of the standard Ogston model. In particular, we demonstrated that the fractional free volume is not the relevant parameter for gel sieving. In this article, we continue the development of this model and we generalize our mathematical approach to treat nonvanishing electric field intensities. To do so, we must revisit the way biased random walks are normally modeled by stochastic processes. Straightforward generalizations based on standard Metropolis-like schemes fail at high field intensities. Moreover, our generalization requires the complete decoupling of the spatial directions parallel and pe...

Control of DNA replication by anomalous reaction-diffusion kinetics

We propose a simple model for the control of DNA replication in which the rate of initiation of replication origins is controlled by protein-DNA interactions. Analyzing recent data from Xenopus frog embryos, we find that the initiation rate is reaction limited until nearly the end of replication, when it becomes diffusion limited. Initiation of origins is suppressed when the diffusion-limited search time dominates. To fit the experimental data, we find that the interaction between DNA and the rate-limiting protein must be subdiffusive.