Kevin D. Dorfman

Microfluidic chemostat for measuring single cell dynamics in bacteria

We designed a microfluidic chemostat consisting of 600 sub-micron trapping/growth channels connected to two feeding channels. The microchemostat traps E. coli cells and forces them to grow in lines for over 50 generations. Excess cells, including the mother cells captured at the start of the process, are removed from both ends of the growth channels by the media flow. With the aid of time-lapse microscopy, we have monitored dynamic properties such as growth rate and GFP expression at the single-cell level for many generations while maintaining a population of bacteria of similar age. We also use the microchemostat to show how the population responds to dynamic changes in the environment. Since more than 100 individual bacterial cells are aligned and immobilized in a single field of view, the microchemostat is an ideal platform for high-throughput intracellular measurements. We demonstrate this capability by tracking with sub-diffraction resolution the movements of fluorescently tagged loci in more than one thousand cells on a single device. The device yields results comparable to conventional agar microscopy experiments with substantial increases in throughput and ease of analysis.

Electrophoretic separation of DNA in gels and nanostructures.

The development of nanostructure devices has opened the door to new DNA separation techniques and fundamental investigations. With advanced nanotechnologies, artificial gels (e.g. nanopillar arrays, nanofilters) can be manufactured with controlled and ordered geometries. This contrast with gels, where the pores are disordered and the range of available pore sizes is limited by the level of cross-linking and the mechanical properties of the gel. In this review, we recall the theories developed for free-solution and gel electrophoresis (extended Ogston model, biased reptation and entropic trapping) and from this perspective, suggestions for future concepts for fast DNA separation using nanostructures will be given.

Short-time movement of E. coli chromosomal loci depends on coordinate and subcellular localization

In bacteria, chromosomal architecture shows strong spatial and temporal organization, and regulates key cellular functions, such as transcription. Tracking the motion of chromosomal loci at short timescales provides information related to both the physical state of the nucleo-protein complex and its local environment, independent of large-scale motions related to genome segregation. Here we investigate the short-time (0.1-10 s) dynamics of fluorescently labelled chromosomal loci in Escherichia coli at different growth rates. At these timescales, we observe for the first time a dependence of the locis apparent diffusion on both their subcellular localization and chromosomal coordinate, and we provide evidence that the properties of the chromosome are similar in the tested growth conditions. Our results indicate that either non-equilibrium fluctuations due to enzyme activity or the organization of the genome as a polymer-protein complex vary as a function of the distance from the origin of replication.

Moving beyond Watson–Crick models of coarse grained DNA dynamics

DNA produces a wide range of structures in addition to the canonical B-form of double-stranded DNA. Some of these structures are stabilized by Hoogsteen bonds. We developed an experimentally parameterized, coarse-grained model that incorporates such bonds. The model reproduces many of the microscopic features of double-stranded DNA and captures the experimental melting curves for a number of short DNA hairpins, even when the open state forms complicated secondary structures. We demonstrate the utility of the model by simulating the folding of a thrombin aptamer, which contains G-quartets, and strand invasion during triplex formation. Our results highlight the importance of including Hoogsteen bonding in coarse-grained models of DNA.

Physical descriptions of the bacterial nucleoid at large scales, and their biological implications

Recent experimental and theoretical approaches have attempted to quantify the physical organization (compaction and geometry) of the bacterial chromosome with its complement of proteins (the nucleoid). The genomic DNA exists in a complex and dynamic protein-rich state, which is highly organized at various length scales. This has implications for modulating (when not directly enabling) the core biological processes of replication, transcription and segregation. We overview the progress in this area, driven in the last few years by new scientific ideas and new interdisciplinary experimental techniques, ranging from high space- and time-resolution microscopy to high-throughput genomics employing sequencing to map different aspects of the nucleoid-related interactome. The aim of this review is to present the wide spectrum of experimental and theoretical findings coherently, from a physics viewpoint. In particular, we highlight the role that statistical and soft condensed matter physics play in describing this system of fundamental biological importance, specifically reviewing classic and more modern tools from the theory of polymers. We also discuss some attempts toward unifying interpretations of the current results, pointing to possible directions for future investigation.

Interplay between chain stiffness and excluded volume of semiflexible polymers confined in nanochannels.

The properties of channel-confined semiflexible polymers are determined by a complicated interplay of chain stiffness and excluded volume effects. Using Pruned-Enriched Rosenbluth Method (PERM) simulations, we study the equilibrium properties of channel-confined polymers by systematically controlling chain stiffness and excluded volume. Our calculations of chain extension and confinement free energy for freely jointed chains with and without excluded volume show excellent agreement with theoretical predictions. For ideal wormlike chains, the extension is seen to crossover from Odijk behavior in strong confinement to zero-stretching, bulk-like behavior in weak confinement. In contrast, for self-avoiding wormlike chains, we always observe that the linear scaling of the extension with the contour length is valid in the long-chain limit irrespective of the regime of confinement, owing to the coexistence of stiffness and excluded volume effects. We further propose that the long-chain limit for the extension corresponds to chain lengths wherein the projection of the end-to-end distance along the axis of the channel is nearly equal to the mean span parallel to the axis. For DNA in nanochannels, this limit was identified using PERM simulations out to molecular weights of more than 1 megabase pairs; the molecular weight of λ-DNA is found to exhibit nearly asymptotic fractional extension for channels sizes used commonly in experiments.

Brownian dynamics simulations of single-stranded DNA hairpins

We present a Brownian dynamics model which we use to study the kinetics and thermodynamics of single-stranded DNA hairpins, gaining insights into the role of stem mismatches and the kinetics rates underlying the melting transition. The model is a base-backbone type in which the DNA bases and sugar-phosphate backbone are represented as single units (beads) in the context of the Brownian dynamics simulations. We employ a minimal number of bead-bead interactions, leading to a simple computational scheme. To demonstrate the veracity of our model for DNA hairpins, we show that the model correctly captures the effects of base stacking, hydrogen bonding, and temperature on both the thermodynamics and the kinetics of hairpin formation and melting. When cast in dimensionless form, the thermodynamic results obtained from the present model compare favorably with default predictions of the m-fold server, although the present model is not sufficiently robust to provide dimensional results. The kinetic data at low temperatures indicate frequent but short-lived opening events, consistent with the measured chain end-to-end probability distribution. The model is also used to study the effect of base mismatches in the stem of the hairpin. With the parameters used here, the model overpredicts the relative shift in the melting temperature due to mismatches. The melting transition can be primarily attributed to a rapid increase in the hairpin opening rate rather than an equivalent decrease in the closing rate, in agreement with single-molecule experimental data.

Modeling the relaxation time of DNA confined in a nanochannel

Using a mapping between a Rouse dumbbell model and fine-grained Monte Carlo simulations, we have computed the relaxation time of λ-DNA in a high ionic strength buffer confined in a nanochannel. The relaxation time thus obtained agrees quantitatively with experimental data [Reisner et al., Phys. Rev. Lett. 94, 196101 (2005)] using only a single O(1) fitting parameter to account for the uncertainty in model parameters. In addition to validating our mapping, this agreement supports our previous estimates of the friction coefficient of DNA confined in a nanochannel [Tree et al., Phys. Rev. Lett. 108, 228105 (2012)], which have been difficult to validate due to the lack of direct experimental data. Furthermore, the model calculation shows that as the channel size passes below approximately 100 nm (or roughly the Kuhn length of DNA) there is a dramatic drop in the relaxation time. Inasmuch as the chain friction rises with decreasing channel size, the reduction in the relaxation time can be solely attributed to the sharp decline in the fluctuations of the chain extension. Practically, the low variance in the observed DNA extension in such small channels has important implications for genome mapping.

Electrophoretic transport through channels of periodically varying cross section

We analyze the electrophoretic transport of a point-sized Brownian particle in a narrow channel of periodically varying cross section. The mean long-time transport coefficients (effective velocity and dispersivity) are calculated using generalized Taylor-Aris dispersion (macrotransport) theory for spatially periodic media. Analytical results for slowly varying geometries are obtained using a long-wavelength perturbation scheme for both O(1) and large Peclet numbers.

Distribution of distances between DNA barcode labels in nanochannels close to the persistence length

We obtained experimental extension data for barcoded E. coli genomic DNA molecules confined in nanochannels from 40 nm to 51 nm in width. The resulting data set consists of 1 627 779 measurements of the distance between fluorescent probes on 25 407 individual molecules. The probability density for the extension between labels is negatively skewed, and the magnitude of the skewness is relatively insensitive to the distance between labels. The two Odijk theories for DNA confinement bracket the mean extension and its variance, consistent with the scaling arguments underlying the theories. We also find that a harmonic approximation to the free energy, obtained directly from the probability density for the distance between barcode labels, leads to substantial quantitative error in the variance of the extension data. These results suggest that a theory for DNA confinement in such channels must account for the anharmonic nature of the free energy as a function of chain extension.