Search

Standard pedagogy introduces optics as though it were a consequence of Maxwell's equations, and only grudgingly admits, usually in a rushed aside, that light has a particulate character that can somehow be reconciled with the wave picture. Recent revolutionary advances in optical imaging, however, make this approach more and more unhelpful: How are we to describe two-photon imaging, FRET, localization microscopy, and a host of related techniques to students who think of light primarily as a wave? I was surprised to find that everything I wanted my biophysics students to know about light, including image formation, x-ray diffraction, and even Bessel beams, could be expressed as well (or better) from the quantum viewpoint pioneered by Richard Feynman. Even my undergraduate students grasp this viewpoint as well as (or better than) the traditional one, and by mid-semester they are already well positioned to integrate the latest advances into their understanding. Moreover, I have found that this approach clarifies my own understanding of new techniques.

Resonance energy transfer has become an indispensable experimental tool for single-molecule and single-cell biophysics. Its physical underpinnings, however, are subtle: It involves a discrete jump of excitation from one molecule to another, and so we regard it as a strongly quantum-mechanical process. And yet, its kinetics differ from what many of us were taught about two-state quantum systems; quantum superpositions of the states do not seem to arise; and so on. Although J. R. Oppenheimer and T. F\"orster navigated these subtleties successfully, it remains hard to find an elementary derivation in modern language. The key step involves acknowledging quantum decoherence. Appreciating that aspect can be helpful when we attempt to extend our understanding to situations where F\"orster's original analysis is not applicable.

Stochastic simulation can make the molecular processes of cellular control more vivid than the traditional differential-equation approach by generating typical system histories instead of just statistical measures such as the mean and variance of a population. Simple simulations are now easy for students to construct from scratch, that is, without recourse to black-box packages. In some cases, their results can also be compared directly to single-molecule experimental data. After introducing the stochastic simulation algorithm, this article gives two case studies, involving gene expression and error correction, respectively. Code samples and resulting animations showing results are given in the online supplements.

Mitosis in the early syncytial Drosophila embryo is highly correlated in space and time, as manifested in mitotic wavefronts that propagate across the embryo. In this paper we investigate the idea that the embryo can be considered a mechanically-excitable medium, and that mitotic wavefronts can be understood as nonlinear wavefronts that propagate through this medium. We study the wavefronts via both image analysis of confocal microscopy videos and theoretical models. We find that the mitotic waves travel across the embryo at a well-defined speed that decreases with replication cycle. We find two markers of the wavefront in each cycle, corresponding to the onsets of metaphase and anaphase. Each of these onsets is followed by displacements of the nuclei that obey the same wavefront pattern. To understand the mitotic wavefronts theoretically we analyze wavefront propagation in excitable media. We study two classes of models, one with biochemical signaling and one with mechanical signaling. We find that the dependence of wavefront speed on cycle number is most naturally explained by mechanical signaling, and that the entire process suggests a scenario in which biochemical and mechanical signaling are coupled.

We present an algorithm to identify individual neural spikes observed on high-density multi-electrode arrays (MEAs). Our method can distinguish large numbers of distinct neural units, even when spikes overlap, and accounts for intrinsic variability of spikes from each unit. As MEAs grow larger, it is important to find spike-identification methods that are scalable, that is, the computational cost of spike fitting should scale well with the number of units observed. Our algorithm accomplishes this goal, and is fast, because it exploits the spatial locality of each unit and the basic biophysics of extracellular signal propagation. Human intervention is minimized and streamlined via a graphical interface. We illustrate our method on data from a mammalian retina preparation and document its performance on simulated data consisting of spikes added to experimentally measured background noise. The algorithm is highly accurate.

The force generated between actin and myosin acts predominantly along the direction of the actin filament, resulting in relative sliding of the thick and thin filaments in muscle or transport of myosin cargos along actin tracks. Previous studies have also detected lateral forces or torques that are generated between actin and myosin, but the origin and biological role of these sideways forces is not known. Here we adapt an actin gliding filament assay in order to measure the rotation of an actin filament about its axis (twirling) as it is translocated by myosin. We quantify the rotation by determining the orientation of sparsely incorporated rhodamine-labeled actin monomers, using polarized total internal reflection (polTIRF) microscopy. In order to determine the handedness of the filament rotation, linear incident polarizations in between the standard s- and p-polarizations were generated, decreasing the ambiguity of our probe orientation measurement four-fold. We found that whole myosin II and myosin V both twirl actin with a relatively long (micron), left-handed pitch that is insensitive to myosin concentration, filament length and filament velocity.

We study the role of microtubule movement in bidirectional organelle transport in Drosophila S2 cells and show that EGFP-tagged peroxisomes in cells serve as sensitive probes of motor induced, noisy cytoskeletal motions. Multiple peroxisomes move in unison over large time windows and show correlations with microtubule tip positions, indicating rapid microtubule fluctuations in the longitudinal direction. We report the first high-resolution measurement of longitudinal microtubule fluctuations performed by tracing such pairs of co-moving peroxisomes. The resulting picture shows that motor-dependent longitudinal microtubule oscillations contribute significantly to cargo movement along microtubules. Thus, contrary to the conventional view, organelle transport cannot be described solely in terms of cargo movement along stationary microtubule tracks, but instead includes a strong contribution from the movement of the tracks., Comment: 24 pages, 5 figures

In many cases, transcriptional regulation involves the binding of transcription factors at sites on the DNA that are not immediately adjacent to the promoter of interest. This action at a distance is often mediated by the formation of DNA loops: Binding at two or more sites on the DNA results in the formation of a loop, which can bring the transcription factor into the immediate neighborhood of the relevant promoter. Though there have been a variety of insights into the combinatorial aspects of transcriptional control, the mechanism of DNA looping as an agent of combinatorial control in both prokaryotes and eukaryotes remains unclear. We use single-molecule techniques to dissect DNA looping in the lac operon. In particular, we measure the propensity for DNA looping by the Lac repressor as a function of the concentration of repressor protein and as a function of the distance between repressor binding sites. As with earlier single-molecule studies, we find (at least) two distinct looped states and demonstrate that the presence of these two states depends both upon the concentration of repressor protein and the distance between the two repressor binding sites. We find that loops form even at interoperator spacings considerably shorter than the DNA persistence length, without the intervention of any other proteins to prebend the DNA. The concentration measurements also permit us to use a simple statistical mechanical model of DNA loop formation to determine the free energy of DNA looping, or equivalently, the J-factor for looping.

We calculate the probability of DNA loop formation mediated by regulatory proteins such as Lac repressor (LacI), using a mathematical model of DNA elasticity. Our model is adapted to calculating quantities directly observable in Tethered Particle Motion (TPM) experiments, and it accounts for all the entropic forces present in such experiments. Our model has no free parameters; it characterizes DNA elasticity using information obtained in other kinds of experiments. [...] We show how to compute both the "looping J factor" (or equivalently, the looping free energy) for various DNA construct geometries and LacI concentrations, as well as the detailed probability density function of bead excursions. We also show how to extract the same quantities from recent experimental data on tethered particle motion, and then compare to our model's predictions. [...] Our model successfully reproduces the detailed distributions of bead excursion, including their surprising three-peak structure, without any fit parameters and without invoking any alternative conformation of the LacI tetramer. Indeed, the model qualitatively reproduces the observed dependence of these distributions on tether length (e.g., phasing) and on LacI concentration (titration). However, for short DNA loops (around 95 basepairs) the experiments show more looping than is predicted by the harmonic-elasticity model, echoing other recent experimental results. Because the experiments we study are done in vitro, this anomalously high looping cannot be rationalized as resulting from the presence of DNA-bending proteins or other cellular machinery. We also show that it is unlikely to be the result of a hypothetical "open" conformation of the LacI tetramer., Comment: See the supplement at http://www.physics.upenn.edu/~pcn/Ms/TowlesEtalSuppl.pdf . This revised version accepted for publication at Physical Biology

Tethered particle motion is an experimental technique to monitor conformational changes in single molecules of DNA in real time, by observing the position fluctuations of a micrometer-size particle attached to the DNA. This article reviews some recent work on theoretical problems inherent in the interpretation of TPM experiments, both in equilibrium and dynamical aspects., Comment: 19pp. Accepted for publication in Curr Op Colloid Interf Science

In many biochemical processes, proteins bound to DNA at distant sites are brought into close proximity by loops in the underlying DNA. For example, the function of some gene-regulatory proteins depends on such DNA looping interactions. We present a new technique for characterizing the kinetics of loop formation in vitro, as observed using the tethered particle method, and apply it to experimental data on looping induced by lambda repressor. Our method uses a modified (diffusive) hidden Markov analysis that directly incorporates the Brownian motion of the observed tethered bead. We compare looping lifetimes found with our method (which we find are consistent over a range of sampling frequencies) to those obtained via the traditional threshold-crossing analysis (which can vary depending on how the raw data are filtered in the time domain). Our method does not involve any time filtering and can detect sudden changes in looping behavior. For example, we show how our method can identify transitions between long-lived, kinetically distinct states that would otherwise be difficult to discern.

Tethered particle experiments use light microscopy to measure the position of a micrometer-sized bead tethered to a microscope slide via a ~micrometer length polymer, in order to infer the behavior of the invisible polymer. Currently, this method is used to measure rate constants of DNA loop formation and breakdown mediated by repressor protein that binds to the DNA. We report a new technique for measuring these rates using a modified hidden Markov analysis that directly incorporates the diffusive motion of the bead, which is an inherent complication of tethered particle motion because it occurs on a time scale between the sampling frequency and the looping time. We compare looping lifetimes found with our method, which are consistent over a range of sampling frequencies, to those obtained via the traditional threshold-crossing analysis, which vary depending on how the raw data are filtered in the time domain. Our method does not involve such filtering, and so can detect short-lived looping events and sudden changes in looping behavior., Comment: 3 page pdf including 3 figures corrections: 2nd page, 1st column, values of diffusion coefficient, spring constant and the decay time were typed incorrectly. No conlcusions were affected

DNA loop formation is one of several mechanisms used by organisms to regulate genes. The free energy of forming a loop is an important factor in determining whether the associated gene is switched on or off. In this paper we use an elastic rod model of DNA to determine the free energy of forming short (50--100 basepair), protein mediated DNA loops. Superhelical stress in the DNA of living cells is a critical factor determining the energetics of loop formation, and we explicitly account for it in our calculations. The repressor protein itself is regarded as a rigid coupler; its geometry enters the problem through the boundary conditions it applies on the DNA. We show that a theory with these ingredients is sufficient to explain certain features observed in modulation of in vivo gene activity as a function of the distance between operator sites for the lac repressor. We also use our theory to make quantitative predictions for the dependence of looping on superhelical stress, which may be testable both in vivo and in single-molecule experiments such as the tethered particle assay and the magnetic bead assay., Comment: Accepted for publication in Phys Rev E. Contains the forbidden appendix they refused to publish

The mechanical properties of DNA play a critical role in many biological functions. For example, DNA packing in viruses involves confining the viral genome in a volume (the viral capsid) with dimensions that are comparable to the DNA persistence length. Similarly, eukaryotic DNA is packed in DNA-protein complexes (nucleosomes) in which DNA is tightly bent around protein spools. DNA is also tightly bent by many proteins that regulate transcription, resulting in a variation in gene expression that is amenable to quantitative analysis. In these cases, DNA loops are formed with lengths that are comparable to or smaller than the DNA persistence length. The aim of this review is to describe the physical forces associated with tightly bent DNA in all of these settings and to explore the biological consequences of such bending, as increasingly accessible by single-molecule techniques., Comment: 24 pages, 9 figures

We describe a simple numerical simulation, suitable for an undergraduate project (or graduate problem set), of the Brownian motion of a particle in a Hooke-law potential well. Understanding this physical situation is a practical necessity in many experimental contexts, for instance in single molecule biophysics; and its simulation helps the student to appreciate the dynamical character of thermal equilibrium. We show that the simulation succeeds in capturing behavior seen in experimental data on tethered particle motion., Comment: Submitted to American Journal of Physics

We propose an alternative mechanism for intracellular cargo transport which results from motor induced longitudinal fluctuations of cytoskeletal microtubules (MT). The longitudinal fluctuations combined with transient cargo binding to the MTs lead to long range transport even for cargos and vesicles having no molecular motors on them. The proposed transport mechanism, which we call ``hitchhiking'', provides a consistent explanation for the broadly observed yet still mysterious phenomenon of bidirectional transport along MTs. We show that cells exploiting the hitchhiking mechanism can effectively up- and down-regulate the transport of different vesicles by tuning their binding kinetics to characteristic MT oscillation frequencies.

Many proteins interact with and deform double-stranded DNA in cells. Single-molecule experiments have studied the elasticity of DNA with helix-deforming proteins, including proteins that bend DNA. These experiments increase the need for theories of DNA elasticity which include helix-deforming proteins. Previous theoretical work on bent DNA has examined a long DNA molecule with many nonspecifically binding proteins. However, recent experiments used relatively short DNA molecules with a single, well-defined bend site. Here we develop a simple theoretical description of the effect of a single bend. We then include the description of the bend in the finite worm like chain model (FWLC) of short DNA molecules attached to beads. We predict how the DNA force-extension relation changes due to formation of a single permanent kink, at all values of the applied stretching force. Our predictions show that high-resolution single-molecule experiments could determine the bend angle induced upon protein binding., Comment: 16 pages, 5 figures

Single-molecule experiments in which force is applied to DNA or RNA molecules have enabled important discoveries of nucleic acid properties and nucleic acid-enzyme interactions. These experiments rely on a model of the polymer force-extension behavior to calibrate the experiments; typically the experiments use the worm-like chain (WLC) theory for double-stranded DNA and RNA. This theory agrees well with experiments for long molecules. Recent single-molecule experiments have used shorter molecules, with contour lengths in the range of 1-10 persistence lengths. Most WLC theory calculations to date have assumed infinite molecule lengths, and do not agree well with experiments on shorter chains. Key physical effects that become important when shorter molecules are used include (i) boundary conditions which constrain the allowed fluctuations at the ends of the molecule and (ii) rotational fluctuations of the bead to which the polymer is attached, which change the apparent extension of the molecule. We describe the finite worm-like chain (FWLC) theory, which takes into account these effects. We show the FWLC predictions diverge from the classic WLC solution for molecules with contour lengths a few times the persistence length. Thus the FWLC will allow more accurate experimental calibration for relatively short molecules, facilitating future discoveries in single-molecule force microscopy., Comment: 30 pages, 10 figures

DNA bending on length scales shorter than a persistence length plays an integral role in the translation of genetic information from DNA to cellular function. Quantitative experimental studies of these biological systems have led to a renewed interest in the polymer mechanics relevant for describing the conformational free energy of DNA bending induced by protein-DNA complexes. Recent experimental results from DNA cyclization studies have cast doubt on the applicability of the canonical semiflexible polymer theory, the wormlike chain (WLC) model, to DNA bending on biological length scales. This paper develops a theory of the chain statistics of a class of generalized semiflexible polymer models. Our focus is on the theoretical development of these models and the calculation of experimental observables. To illustrate our methods, we focus on a specific toy model of DNA bending. We show that the WLC model generically describes the long-length-scale chain statistics of semiflexible polymers, as predicted by the Renormalization Group. In particular, we show that either the WLC or our new model adequate describes force-extension, solution scattering, and long-contour-length cyclization experiments, regardless of the details of DNA bend elasticity. In contrast, experiments sensitive to short-length-scale chain behavior can in principle reveal dramatic departures from the linear elastic behavior assumed in the WLC model. We demonstrate this explicitly by showing that our toy model can reproduce the anomalously large short-contour-length cyclization J factors observed by Cloutier and Widom. Finally, we discuss the applicability of these models to DNA chain statistics in the context of future experiments.

The importance of nonlinearities in material constitutive relations has long been appreciated in the continuum mechanics of macroscopic rods. Although the moment (torque) response to bending is almost universally linear for small deflection angles, many rod systems exhibit a high-curvature softening. The signature behavior of these rod systems is a kinking transition in which the bending is localized. Recent DNA cyclization experiments by Cloutier and Widom have offered evidence that the linear-elastic bending theory fails to describe the high-curvature mechanics of DNA. Motivated by this recent experimental work, we develop a simple and exact theory of the statistical mechanics of linear-elastic polymer chains that can undergo a kinking transition. We characterize the kinking behavior with a single parameter and show that the resulting theory reproduces both the low-curvature linear-elastic behavior which is already well described by the Wormlike Chain model, as well as the high-curvature softening observed in recent cyclization experiments., Comment: Revised for PRE. 40 pages, 12 figures