David K. Lubensky

Driven polymer translocation through a narrow pore

Motivated by experiments in which a polynucleotide is driven through a proteinaceous pore by an electric field, we study the diffusive motion of a polymer threaded through a narrow channel with which it may have strong interactions. We show that there is a range of polymer lengths in which the system is approximately translationally invariant, and we develop a coarse-grained description of this regime. From this description, general features of the distribution of times for the polymer to pass through the pore may be deduced. We also introduce a more microscopic model. This model provides a physically reasonable scenario in which, as in experiments, the polymers speed depends sensitively on its chemical composition, and even on its orientation in the channel. Finally, we point out that the experimental distribution of times for the polymer to pass through the pore is much broader than expected from simple estimates, and speculate on why this might be.

Unzipping kinetics of double-stranded DNA in a nanopore.

We studied the unzipping of single molecules of double-stranded DNA by pulling one of their two strands through a narrow protein pore. Polymerase chain reaction analysis yielded the first direct proof of DNA unzipping in such a system. The time to unzip each molecule was inferred from the ionic current signature of DNA traversal. The distribution of times to unzip under various experimental conditions fit a simple kinetic model. Using this model, we estimated the enthalpy barriers to unzipping and the effective charge of a nucleotide in the pore, which was considerably smaller than previously assumed.

DNA unzipped under a constant force exhibits multiple metastable intermediates.

Single molecule studies, at constant force, of the separation of double-stranded DNA into two separated single strands may provide information relevant to the dynamics of DNA replication. At constant applied force, theory predicts that the unzipped length as a function of time is characterized by jumps during which the strands separate rapidly, followed by long pauses where the number of separated base pairs remains constant. Here, we report previously uncharacterized observations of this striking behavior carried out on a number of identical single molecules simultaneously. When several single λ phage molecules are subject to the same applied force, the pause positions are reproducible in each. This reproducibility shows that the positions and durations of the pauses in unzipping provide a sequence-dependent molecular fingerprint. For small forces, the DNA remains in a partially unzipped state for at least several hours. For larger forces, the separation is still characterized by jumps and pauses, but the double-stranded DNA will completely unzip in less than 30 min.

An allosteric model of circadian KaiC phosphorylation

In a recent series of ground-breaking experiments, Nakajima et al. [Nakajima M, Imai K, Ito H, Nishiwaki T, Murayama Y, Iwasaki H, Oyama T, Kondo T (2005) Science 308:414–415] showed that the three cyanobacterial clock proteins KaiA, KaiB, and KaiC are sufficient in vitro to generate circadian phosphorylation of KaiC. Here, we present a mathematical model of the Kai system. At its heart is the assumption that KaiC can exist in two conformational states, one favoring phosphorylation and the other dephosphorylation. Each individual KaiC hexamer then has a propensity to be phosphorylated in a cyclic manner. To generate macroscopic oscillations, however, the phosphorylation cycles of the different hexamers must be synchronized. We propose a novel synchronization mechanism based on differential affinity: KaiA stimulates KaiC phosphorylation, but the limited supply of KaiA dimers binds preferentially to those KaiC hexamers that are falling behind in the oscillation. KaiB sequesters KaiA and stabilizes the dephosphorylating KaiC state. We show that our model can reproduce a wide range of published data, including the observed insensitivity of the oscillation period to variations in temperature, and that it makes nontrivial predictions about the effects of varying the concentrations of the Kai proteins.

Single molecule statistics and the polynucleotide unzipping transition.

We present an extensive theoretical investigation of the mechanical unzipping of double-stranded DNA under the influence of an applied force. In the limit of long polymers, there is a thermodynamic unzipping transition at a critical force value of order 10 pN, with different critical behavior for homopolymers and for random heteropolymers. We extend results on the disorder-averaged behavior of DNAs with random sequences [D. K. Lubensky and D. R. Nelson, Phys. Rev. Lett. 85, 1572 (2000)] to the more experimentally accessible problem of unzipping a single DNA molecule. As the applied force approaches the critical value, the double-stranded DNA unravels in a series of discrete, sequence-dependent steps that allow it to reach successively deeper energy minima. Plots of extension versus force thus take the striking form of a series of plateaus separated by sharp jumps. Similar qualitative features should reappear in micromanipulation experiments on proteins and on folded RNA molecules. Despite their unusual form, the extension versus force curves for single molecules still reveal remnants of the disorder-averaged critical behavior. Above the transition, the dynamics of the unzipping fork is related to that of a particle diffusing in a random force field; anomalous, disorder-dominated behavior is expected until the applied force exceeds the critical value for unzipping by roughly 5 pN.

Dynamics of Molecular Motors and Polymer Translocation with Sequence Heterogeneity

The effect of sequence heterogeneity on polynucleotide translocation across a pore and on simple models of molecular motors such as helicases, DNA polymerase/exonuclease, and RNA polymerase is studied in detail. Pore translocation of RNA or DNA is biased due to the different chemical environments on the two sides of the membrane, whereas the molecular motor motion is biased through a coupling to chemical energy. An externally applied force can oppose these biases. For both systems we solve lattice models exactly both with and without disorder. The models incorporate explicitly the coupling to the different chemical environments for polymer translocation and the coupling to the chemical energy (as well as nucleotide pairing energies) for molecular motors. Using the exact solutions and general arguments, we show that the heterogeneity leads to anomalous dynamics. Most notably, over a range of forces around the stall force (or stall tension for DNA polymerase/exonuclease systems) the displacement grows sublinearly as t(micro), with micro < 1. The range over which this behavior can be observed experimentally is estimated for several systems and argued to be detectable for appropriate forces and buffers. Similar sequence heterogeneity effects may arise in the packing of viral DNA.

Robust circadian clocks from coupled protein-modification and transcription–translation cycles

The cyanobacterium Synechococcus elongatus uses both a protein phosphorylation cycle and a transcription–translation cycle to generate circadian rhythms that are highly robust against biochemical noise. We use stochastic simulations to analyze how these cycles interact to generate stable rhythms in growing, dividing cells. We find that a protein phosphorylation cycle by itself is robust when protein turnover is low. For high decay or dilution rates (and compensating synthesis rates), however, the phosphorylation-based oscillator loses its integrity. Circadian rhythms thus cannot be generated with a phosphorylation cycle alone when the growth rate, and consequently the rate of protein dilution, is high enough; in practice, a purely posttranslational clock ceases to function well when the cell doubling time drops below the 24-h clock period. At higher growth rates, a transcription–translation cycle becomes essential for generating robust circadian rhythms. Interestingly, although a transcription–translation cycle is necessary to sustain a phosphorylation cycle at high growth rates, a phosphorylation cycle can dramatically enhance the robustness of a transcription–translation cycle at lower protein decay or dilution rates. In fact, the full oscillator built from these two tightly intertwined cycles far outperforms not just each of its two components individually, but also a hypothetical system in which the two parts are coupled as in textbook models of coupled phase oscillators. Our analysis thus predicts that both cycles are required to generate robust circadian rhythms over the full range of growth conditions.

Coupling mechanical deformations and planar cell polarity to create regular patterns in the zebrafish retina.

The orderly packing and precise arrangement of epithelial cells is essential to the functioning of many tissues, and refinement of this packing during development is a central theme in animal morphogenesis. The mechanisms that determine epithelial cell shape and position, however, remain incompletely understood. Here, we investigate these mechanisms in a striking example of planar order in a vertebrate epithelium: The periodic, almost crystalline distribution of cone photoreceptors in the adult teleost fish retina. Based on observations of the emergence of photoreceptor packing near the retinal margin, we propose a mathematical model in which ordered columns of cells form as a result of coupling between planar cell polarity (PCP) and anisotropic tissue-scale mechanical stresses. This model recapitulates many observed features of cone photoreceptor organization during retinal growth and regeneration. Consistent with the models predictions, we report a planar-polarized distribution of Crumbs2a protein in cone photoreceptors in both unperturbed and regenerated tissue. We further show that the pattern perturbations predicted by the model to occur if the imposed stresses become isotropic closely resemble defects in the cone pattern in zebrafish lrp2 mutants, in which intraocular pressure is increased, resulting in altered mechanical stress and ocular enlargement. Evidence of interactions linking PCP, cell shape, and mechanical stresses has recently emerged in a number of systems, several of which show signs of columnar cell packing akin to that described here. Our results may hence have broader relevance for the organization of cells in epithelia. Whereas earlier models have allowed only for unidirectional influences between PCP and cell mechanics, the simple, phenomenological framework that we introduce here can encompass a broad range of bidirectional feedback interactions among planar polarity, shape, and stresses; our model thus represents a conceptual framework that can address many questions of importance to morphogenesis.

Hydrodynamics of monolayer domains at the air–water interface

Molecules at the air–water interface often form inhomogeneous layers in which domains of different densities are separated by sharp interfaces. Complex interfacial pattern formation may occur through the competition of short‐ and long‐range forces acting within the monolayer. The overdamped hydrodynamics of such interfacial motion is treated here in a general manner that accounts for dissipation both within the monolayer and in the subfluid. Previous results on the linear stability of interfaces are recovered and extended, and a formulation applicable to the nonlinear regime is developed. A simplified dynamical law valid when dissipation in the monolayer itself is negligible is also proposed. Throughout the analysis, special attention is paid to the dependence of the dynamical behavior on a characteristic length scale set by the ratio of the viscosities in the monolayer and in the subphase.

A dynamical model of ommatidial crystal formation

The crystalline photoreceptor lattice in the Drosophila eye is a paradigm for pattern formation during development. During eye development, activation of proneural genes at a moving front adds new columns to a regular lattice of R8 photoreceptors. We present a mathematical model of the governing activator–inhibitor system, which indicates that the dynamics of positive induction play a central role in the selection of certain cells as R8s. The “switch and template” patterning mechanism we observe is mathematically very different from the well-known Turing instability. Unlike a standard lateral inhibition model, our picture implies that R8s are defined before the appearance of the complete group of proneural cells. The model reproduces the full time course of proneural gene expression and accounts for specific features of the refinement of proneural groups that had resisted explanation. It moreover predicts that perturbing the normal template can lead to eyes containing stripes of R8 cells. We observed these stripes experimentally after manipulation of the Notch and scabrous genes. Our results suggest an alternative to the generally assumed mode of operation for lateral inhibition during development; more generally, they hint at a broader role for bistable switches in the initial establishment of patterns as well as in their maintenance.