Self-Avoiding Wormlike Chain Confined in a Cylindrical Tube: Scaling Behavior

Escape dynamics of active ring polymers in a cylindrical nanochannel

We explore the escape dynamics of active ring polymers confined in a cylindrical nanochannel using Brownian dynamics. Our simulation results show that the escape time decreases with the increase of the Péclet number, which is not noticeable between the two stages of the escape process, based on whether the center of mass of the polymer is inside or outside the nanochannel. However, the monomer motion trajectory of the active polymer is very different from that of the passive polymer, similar to the snake-like motion with uniform velocity. The passive polymer, however, is in constant fugitive motion with increased velocity at the tail end of the escape. Our work is vital for understanding the escape dynamics of active ring polymers in the confined nanochannel, which provides new perspectives on their characterization and analysis.

Entropic force of cone-tethered polymers interacting with a planar surface

Computer simulations are used to characterize the entropic force of one or more polymers tethered to the tip of a hard conical object that interact with a nearby hard flat surface. Pruned-enriched Rosenbluth method Monte Carlo simulations are used to calculate the variation of the conformational free energy F of a hard-sphere polymer with respect to the cone-tip-to-surface distance h from which the variation of the entropic force f≡|dF/dh| with h is determined. We consider the following cases: a single freely jointed tethered chain, a single semiflexible tethered chain, and several freely jointed chains of equal length each tethered to the cone tip. The simulation results are used to test the validity of a prediction by Maghrebi et al. [Maghrebi et al., Europhys. Lett. 96, 66002 (2011)0295-507510.1209/0295-5075/96/66002; Phys. Rev. E 86, 061801 (2012)1539-375510.1103/PhysRevE.86.061801] that f∝(γ_{∞}-γ_{0})h^{-1}, where γ_{0} and γ_{∞} are universal scaling exponents for the partition function of the tethered polymer for h=0 and h=∞, respectively. The measured functions f(h) are generally consistent with the predictions, with small quantitative discrepancies arising from the approximations employed in the theory. In the case of multiple tethered polymers, the entropic force per polymer is roughly constant, which is qualitatively inconsistent with the predictions.

Confinement anisotropy drives polar organization of two DNA molecules interacting in a nanoscale cavity

There is growing appreciation for the role phase transition based phenomena play in biological systems. In particular, self-avoiding polymer chains are predicted to undergo a unique confinement dependent demixing transition as the anisotropy of the confined space is increased. This phenomenon may be relevant for understanding how interactions between multiple dsDNA molecules can induce self-organized structure in prokaryotes. While recent in vivo experiments and Monte Carlo simulations have delivered essential insights into this phenomenon and its relation to bacteria, there are fundamental questions remaining concerning how segregated polymer states arise, the role of confinement anisotropy and the nature of the dynamics in the segregated states. To address these questions, we introduce an artificial nanofluidic model to quantify the interactions of multiple dsDNA molecules in cavities with controlled anisotropy. We find that two dsDNA molecules of equal size confined in an elliptical cavity will spontaneously demix and orient along the cavity poles as cavity eccentricity is increased; the two chains will then swap pole positions with a frequency that decreases with increasing cavity eccentricity. In addition, we explore a system consisting of a large dsDNA molecule and a plasmid molecule. We find that the plasmid is excluded from the larger molecule and will exhibit a preference for the ellipse poles, giving rise to a non-uniform spatial distribution in the cavity that may help explain the non-uniform plasmid distribution observed during in vivo imaging of high-copy number plasmids in bacteria.

Bistability of a helical filament confined on a cylinder

The natural configuration of an intrinsically curved and twisted filament is uniquely a helix so that it can be referred to as a helical filament. We find that confining a helical filament on a cylinder can create a bistable state. When c_{0}R=0.5, where c_{0} is the intrinsic curvature of filament and R is the radius of cylinder, the phase diagram for the stability of a helix contains three regimes. Regime I has a small intrinsic twisting rate (ITR) and exhibits a bistable state which consists of two isoenergic helices. In regime II, the filament has a moderate ITR and the bistable state consists of a metastable low-pitch helix and a stable nonhelix. In regime III, the helix is unstable, owing to a large ITR. A similar phenomenon occurs when c_{0}R∼0.5. Monte Carlo simulation confirms these conclusions and indicates further that there are bistable nonhelices in regime III. This bistable system offers a prospective green material since the wide range of parameters and distinctive configurations for bistable states favor its realization and application.

Crowding-induced polymer trapping in a channel

In this work, we report an intriguing phenomenon: crowding-induced polymer trapping in a channel. Using Langevin dynamics simulations and analytical calculations, we find that for a polymer confined in a channel, crowding particles can push a polymer into the channel corner through inducing an effective polymer-corner attraction due to the depletion effect. This phenomenon is referred to as polymer trapping. The occurrence of polymer trapping requires a minimum volume fraction of crowders, ϕ^{*}, which scales as ϕ^{*}∼(a_{c}/L_{p})^{1/3} for a_{c}≫a_{m} and ϕ^{*}∼(a_{c}/L_{p})^{1/3}(a_{c}/a_{m})^{1/2} for a_{c}≪a_{m}, where a_{c} is the crowder diameter, a_{m} is the monomer diameter, and L_{p} is the polymer persistence length. For DNA, ϕ^{*} is estimated to be around 0.25 for crowders with a_{c}=2nm. We find that ϕ^{*} also strongly depends on the shape of the channel cross section, and ϕ^{*} is much smaller for a triangle channel than a square channel. The polymer trapping leads to a nearly fully stretched polymer conformation along a channel corner, which may have practical applications, such as full stretching of DNA for the nanochannel-based genome mapping technology.

Single-molecule optical genome mapping in nanochannels: multidisciplinarity at the nanoscale

The human genome contains multiple layers of information that extend beyond the genetic sequence. In fact, identical genetics do not necessarily yield identical phenotypes as evident for the case of two different cell types in the human body. The great variation in structure and function displayed by cells with identical genetic background is attributed to additional genomic information content. This includes large-scale genetic aberrations, as well as diverse epigenetic patterns that are crucial for regulating specific cell functions. These genetic and epigenetic patterns operate in concert in order to maintain specific cellular functions in health and disease. Single-molecule optical genome mapping is a high-throughput genome analysis method that is based on imaging long chromosomal fragments stretched in nanochannel arrays. The access to long DNA molecules coupled with fluorescent tagging of various genomic information presents a unique opportunity to study genetic and epigenetic patterns in the genome at a single-molecule level over large genomic distances. Optical mapping entwines synergistically chemical, physical, and computational advancements, to uncover invaluable biological insights, inaccessible by sequencing technologies. Here we describe the method's basic principles of operation, and review the various available mechanisms to fluorescently tag genomic information. We present some of the recent biological and clinical impact enabled by optical mapping and present recent approaches for increasing the method's resolution and accuracy. Finally, we discuss how multiple layers of genomic information may be mapped simultaneously on the same DNA molecule, thus paving the way for characterizing multiple genomic observables on individual DNA molecules.

Bistability induced by a spontaneous twisting rate for a two-dimensional intrinsically curved filament

We find that a moderate intrinsic twisting rate (ITR) can induce a bistable state for a force-free two-dimensional intrinsically curved filament. There are two different configurations of equal energy in a bistable state so that the filament is clearly different from its three-dimensional counterpart. The smaller the ITR or the larger the intrinsic curvature (IC), the clearer the distinction between two isoenergetic configurations and the longer the filament. In bistable states, the relationship between length and ITR is approximately a hyperbola and relationship between IC and critical ITR is approximately linear. Thermal fluctuation can result in a shift between two isoenergetic configurations, but large bending and twisting rigidities can prevent the shift and maintain the filament in one of these two configurations. Moreover, a filament can have a metastable state and at a finite temperature such a filament has the similar property as that of a filament with bistable state.

Limitations of the equivalent neutral polymer assumption for theories describing nanochannel-confined DNA

The prevailing theories describing DNA confinement in a nanochannel are predicated on the assumption that wall-DNA electrostatic interactions are sufficiently short-ranged such that the problem can be mapped to an equivalent neutral polymer confined by hard walls with an appropriately reduced effective channel size. To determine when this hypothesis is valid, we leveraged a recently reported experimental data set for the fractional extension of DNA molecules in a 250-nm-wide poly(dimethyl siloxane) (PDMS) nanochannel with buffer ionic strengths between 0.075 and 48 mM. Evaluating these data in the context of the weakly correlated telegraph model of DNA confinement reveals that, at ionic strengths greater than 0.3 mM, the average fractional extension of the DNA molecules agree with theoretical predictions with a mean absolute error of 0.04. In contrast, experiments at ionic strengths below 0.3 mM produce average fractional extensions that are systematically smaller than the theoretical predictions with a larger mean absolute error of 0.15. The deviations between experiment and theory display a correlation coefficient of 0.82 with the decay length for the DNA-wall electrostatics, linking the deviations with a breakdown in approximating the DNA with an equivalent neutral polymer.

Extension distribution for DNA confined in a nanochannel near the Odijk regime

DNA confinement in a nanochannel typically is understood via mapping to the confinement of an equivalent neutral polymer by hard walls. This model has proven to be effective for confinement in relatively large channels where hairpin formation is frequent. An analysis of existing experimental data for Escherichia coli DNA extension in channels smaller than the persistence length, combined with an additional dataset for λ-DNA confined in a 34 nm wide channel, reveals a breakdown in this approach as the channel size approaches the Odijk regime of strong confinement. In particular, the predicted extension distribution obtained from the asymptotic solution to the weakly correlated telegraph model for a confined wormlike chain deviates significantly from the experimental distribution obtained for DNA confinement in the 34 nm channel, and the discrepancy cannot be resolved by treating the alignment fluctuations or the effective channel size as fitting parameters. We posit that the DNA-wall electrostatic interactions, which are sensible throughout a significant fraction of the channel cross section in the Odijk regime, are the source of the disagreement between theory and experiment. Dimensional analysis of the wormlike chain propagator in channel confinement reveals the importance of a dimensionless parameter, reflecting the magnitude of the DNA-wall electrostatic interactions relative to thermal energy, which has not been considered explicitly in the prevailing theories for DNA confinement in a nanochannel.

Orientationally ordered states of a wormlike chain in spherical confinement

One of the basic characteristics of a linear dsDNA molecule is its persistence length, typically of order 50 nm. The DNA chain inflicts a large energy penalty if it is bent sharply at that length scale. Viruses of bacteria, known as bacteriophages, typically have a dimension of a few tens of nanometers. Yet, it is known that a bacteriophage actively packages viral DNA inside the capsid and ejects it afterwards. Here, adopting a commonly used polymer model known as the wormlike chain, we answer an idealized question: Placing a linear DNA molecule inside a spherical cavity, what ordered states can we derive from known tools in statistical physics? Solving the model in a rigorous field-theory framework, we report a universal phase diagram for four orientationally ordered and disordered states, in terms of two relevant physical parameters.

Simulations corroborate telegraph model predictions for the extension distributions of nanochannel confined DNA

Hairpins in the conformation of DNA confined in nanochannels close to their persistence length cause the distribution of their fractional extensions to be heavily left skewed. A recent theory rationalizes these skewed distributions using a correlated telegraph process, which can be solved exactly in the asymptotic limit of small but frequent hairpin formation. Pruned-enriched Rosenbluth method simulations of the fractional extension distribution for a channel-confined wormlike chain confirm the predictions of the telegraph model. Remarkably, the asymptotic result of the telegraph model remains robust well outside the asymptotic limit. As a result, the approximations in the theory required to map it to the polymer model and solve it in the asymptotic limit are not the source of discrepancies between the predictions of the telegraph model and experimental distributions of the extensions of DNA during genome mapping. The agreement between theory and simulations motivates future work to determine the source of the remaining discrepancies between the predictions of the telegraph model and experimental distributions of the extensions of DNA in nanochannels used for genome mapping.

Nanochannel-Confined TAMRA-Polypyrrole Stained DNA Stretching by Varying the Ionic Strength from Micromolar to Millimolar Concentrations

Large DNA molecules have been utilized as a model system to investigate polymer physics. However, DNA visualization via intercalating dyes has generated equivocal results due to dye-induced structural deformation, particularly unwanted unwinding of the double helix. Thus, the contour length increases and the persistence length changes so unpredictably that there has been a controversy. In this paper, we used TAMRA-polypyrrole to stain single DNA molecules. Since this staining did not change the contour length of B-form DNA, we utilized TAMRA-polypyrrole stained DNA as a tool to measure the persistence length by changing the ionic strength. Then, we investigated DNA stretching in nanochannels by varying the ionic strength from 0.06 mM to 47 mM to evaluate several polymer physics theories proposed by Odijk, de Gennes and recent papers to deal with these regimes.

Measuring the wall depletion length of nanoconfined DNA

Efforts to study the polymer physics of DNA confined in nanochannels have been stymied by a lack of consensus regarding its wall depletion length. We have measured this quantity in 38 nm wide, square silicon dioxide nanochannels for five different ionic strengths between 15 mM and 75 mM. Experiments used the Bionano Genomics Irys platform for massively parallel data acquisition, attenuating the effect of the sequence-dependent persistence length and finite-length effects by using nick-labeled E. coli genomic DNA with contour length separations of at least 30 µm (88 325 base pairs) between nick pairs. Over 5 × 106 measurements of the fractional extension were obtained from 39 291 labeled DNA molecules. Analyzing the stretching via Odijk's theory for a strongly confined wormlike chain yielded a linear relationship between the depletion length and the Debye length. This simple linear fit to the experimental data exhibits the same qualitative trend as previously defined analytical models for the depletion length but now quantitatively captures the experimental data.