Source author record

M. D. Betterton

M. D. Betterton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2016arXiv

Motor protein accumulation on antiparallel microtubule overlaps

Biopolymers serve as one-dimensional tracks on which motor proteins move to perform their biological roles. Motor protein phenomena have inspired theoretical models of one-dimensional transport, crowding, and jamming. Experiments studying the motion of Xklp1 motors on reconstituted antiparallel microtubule overlaps demonstrated that motors recruited to the overlap walk toward the plus end of individual microtubules and frequently switch between filaments. We study a model of this system that couples the totally asymmetric simple exclusion process (TASEP) for motor motion with switches between antiparallel filaments and binding kinetics. We determine steady-state motor density profiles for fixed-length overlaps using exact and approximate solutions of the continuum differential equations and compare to kinetic Monte Carlo simulations. Overlap motor density profiles and motor trajectories resemble experimental measurements. The phase diagram of the model is similar to the single-filament case for low switching rate, while for high switching rate we find a new low density-high density-low density-high density phase. The overlap center region, far from the overlap ends, has a constant motor density as one would naively expect. However, rather than following a simple binding equilibrium, the center motor density depends on total overlap length, motor speed, and motor switching rate. The size of the crowded boundary layer near the overlap ends is also dependent on the overlap length and switching rate in addition to the motor speed and bulk concentration. The antiparallel microtubule overlap geometry may offer a previously unrecognized mechanism for biological regulation of protein concentration and consequent activity.

preprint2015arXiv

Hysteresis, reentrance, and glassy dynamics in systems of self-propelled rods

Non-equilibrium active matter made up of self-driven particles with short-range repulsive interactions is a useful minimal system to study active matter as the system exhibits collective motion and nonequilibrium order-disorder transitions. We studied high-aspect-ratio self-propelled rods over a wide range of packing fraction and driving to determine the nonequilibrium state diagram and dynamic properties. Flocking and nematic-laning states occupy much of the parameter space. In the flocking state the average internal pressure is high and structural and mechanical relaxation times are long, suggesting that rods in flocks are in a translating glassy state despite overall flock motion. In contrast, the nematic-laning state shows fluid-like behavior. The flocking state occupies regions of the state diagram at both low and high packing fraction separated by nematic-laning at low driving and a history-dependent region at higher driving; the nematic-laning state transitions to the flocking state for both compression and expansion. We propose that the laning-flocking transitions are a type of glass transition which, in contrast to other glass-forming systems, can show fluidization as density increases. The fluid internal dynamics and ballistic transport of the nematic-laning state may promote collective dynamics of rod-shaped microorganisms.

preprint2015arXiv

Multiscale modeling and simulation of microtubule/motor protein assemblies

Microtubules and motor proteins self organize into biologically important assemblies including the mitotic spindle and the centrosomal microtubule array. Outside of cells, microtubule-motor mixtures can form novel active liquid-crystalline materials driven out of equilibrium by ATP-consuming motor proteins. Microscopic motor activity causes polarity-dependent interactions between motor proteins and microtubules, but how these interactions yield such larger-scale dynamical behavior such as complex flows and defect dynamics is not well understood. We develop a multiscale theory for microtubule-motor systems in which Brownian dynamics simulations of polar microtubules driven by motors are used to study microscopic organization and stresses created by motor-mediated microtubule interactions. We identify polarity-sorting and crosslink tether relaxation as two polar-specific sources of active destabilizing stress. We then develop a continuum Doi-Onsager model that captures polarity sorting and the hydrodynamic flows generated by these polar-specific active stresses. In simulations of active nematic flows on immersed surfaces, the active stresses drive turbulent flow dynamics and continuous generation and annihilation of disclination defects. The dynamics follow from two instabilities, and accounting for the immersed nature of the experiment yields unambiguous characteristic length and time scales. When turning off the hydrodynamics in the Doi-Onsager model, we capture formation of polar lanes as observed in the Brownian dynamics simulation.

preprint2015arXiv

Multiscale polar theory of microtubule and motor-protein assemblies

Microtubules and motor proteins are building blocks of self-organized subcellular biological structures such as the mitotic spindle and the centrosomal microtubule array. These same ingredients can form new "bioactive" liquid-crystalline fluids that are intrinsically out of equilibrium and which display complex flows and defect dynamics. It is not yet well understood how microscopic activity, which involves polarity-dependent interactions between motor proteins and microtubules, yields such larger scale dynamical structures. In our multiscale theory, Brownian dynamics simulations of polar microtubule ensembles driven by crosslinking motors allow us to study microscopic organization and stresses. Polarity sorting and crosslink relaxation emerge as two polar-specific sources of active destabilizing stress. On larger length scales, our continuum Doi-Onsager theory captures the hydrodynamic flows generated by polarity-dependent active stresses. The results connect local polar structure to flow structures and defect dynamics.

preprint2013arXiv

Biophysics of filament length regulation by molecular motors

Regulating physical size is an essential problem that biological organisms must solve from the subcellular to the organismal scales, but it is not well understood what physical principles and mechanisms organisms use to sense and regulate their size. Any biophysical size-regulation scheme operates in a noisy environment and must be robust to other cellular dynamics and fluctuations. This work develops theory of filament length regulation inspired by recent experiments on kinesin-8 motor proteins, which move with directional bias on microtubule filaments and alter microtubule dynamics. Purified kinesin-8 motors can depolymerize chemically-stabilized microtubules. In the length-dependent depolymerization model, the rate of depolymerization tends to increase with filament length, because long filaments accumulate more motors at their tips and therefore shorten more quickly. When balanced with a constant filament growth rate, this mechanism can lead to a fixed polymer length. However, the mechanism by which kinesin-8 motors affect the length of dynamic microtubules in cells is less clear. We study the more biologically realistic problem of microtubule dynamic instability modulated by a motor-dependent increase in the filament catastrophe frequency. This leads to a significant decrease in the mean filament length and a narrowing of the filament length distribution. The results improve our understanding of the biophysics of length regulation in cells.

preprint2010arXiv

Linear aggregation and liquid-crystalline order: comparison of Monte Carlo simulation and analytic theory

Many soft-matter and biophysical systems are composed of monomers which reversibly assemble into rod-like aggregates. The aggregates can then order into liquid-crystal phases if the density is high enough, and liquid-crystal ordering promotes increased growth of aggregates. Systems that display coupled aggregation and liquid-crystal ordering include wormlike micelles, chromonic liquid crystals, DNA and RNA, and protein polymers and fibrils. Coarse-grained molecular models that capture key features of coupled aggregation and liquid-crystal ordering common to many different systems are lacking; in particular, the role of monomer aspect ratio and aggregate flexibility in controlling the phase behavior are not well understood. Here we study a minimal system of sticky cylinders using Monte Carlo simulations and analytic theory. Cylindrical monomers interact primarily by hard-core interactions but can stack and bind end to end. We present results for several different cylinder aspect ratios and a range of end-to-end binding energies. The phase diagrams are qualitatively similar to those of chromonic liquid crystals, with an isotropic-nematic-columnar triple point. The location of the triple point is sensitive to the monomer aspect ratio.We find that the aggregate persistence length varies with temperature in a way that is controlled by the interaction potential; this suggests that the form of the interaction potential affects the phase behavior of the system. Our analytic theory shows improvement compared to previous theory in quantitatively predicting the I-N transition for relatively stiff aggregates, but requires a better treatment of aggregate flexibility.

preprint2008arXiv

Microtubule depolymerization by the kinesin-8 motor Kip3p: a mathematical model

Proteins from the kinesin-8 family promote microtubule (MT) depolymerization, a process thought to be important for the control of microtubule length in living cells. In addition to this MT shortening activity, kinesin 8s are motors that show plus-end directed motility on MTs. Here we describe a simple model that incorporates directional motion and destabilization of the MT plus end by kinesin 8. Our model quantitatively reproduces the key features of length-vs-time traces for stabilized MTs in the presence of purified kinesin 8, including length-dependent depolymerization. Comparison of model predictions with experiments suggests that kinesin 8 depolymerizes processively, i.e., one motor can remove multiple tubulin dimers from a stabilized MT. Fluctuations in MT length as a function of time are related to depolymerization processivity. We have also determined the parameter regime in which the rate of MT depolymerization is length dependent: length-dependent depolymerization occurs only when MTs are sufficiently short; this crossover is sensitive to the bulk motor concentration.

preprint2006arXiv

DNA entropic elasticity for short molecules attached to beads

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.

preprint2005arXiv

Incorporating expression data in metabolic modeling: a case study of lactate dehydrogenase

Integrating biological information from different sources to understand cellular processes is an important problem in systems biology. We use data from mRNA expression arrays and chemical kinetics to formulate a metabolic model relevant to K562 erythroleukemia cells. MAP kinase pathway activation alters the expression of metabolic enzymes in K562 cells. Our array data show changes in expression of lactate dehydrogenase (LDH) isoforms after treatment with phorbol 12-myristate 13-acetate (PMA), which activates MAP kinase signaling. We model the change in lactate production which occurs when the MAP kinase pathway is activated, using a non-equilibrium, chemical-kinetic model of homolactic fermentation. In particular, we examine the role of LDH isoforms, which catalyze the conversion of pyruvate to lactate. Changes in the isoform ratio are not the primary determinant of the production of lactate. Rather, the total concentration of LDH controls the lactate concentration.

preprint2002arXiv

Opening of DNA double strands by helicases. Active versus passive opening

Helicase opening of double-stranded nucleic acids may be "active" (the helicase directly destabilizes the dsNA to promote opening) or "passive" (the helicase binds ssNA available due to a thermal fluctuation which opens part of the dsNA). We describe helicase opening of dsNA, based on helicases which bind single NA strands and move towards the double-stranded region, using a discrete ``hopping'' model. The interaction between the helicase and the junction where the double strand opens is characterized by an interaction potential. The form of the potential determines whether the opening is active or passive. We calculate the rate of passive opening for the helicase PcrA, and show that the rate increases when the opening is active. Finally, we examine how to choose the interaction potential to optimize the rate of strand separation. One important result is our finding that active opening can increase the unwinding rate by 7 fold compared to passive opening.

preprint2000arXiv

Formation of Structure in Snowfields: Penitentes, Suncups, and Dirt Cones

Penitentes and suncups are structures formed as snow melts, typically high in the mountains. When the snow is dirty, dirt cones and other structures can form instead. Building on previous field observations and experiments, this work presents a theory of ablation morphologies, and the role of surface dirt in determining the structures formed. The glaciological literature indicates that sunlight, heating from air, and dirt all play a role in the formation of structure on an ablating snow surface. The present work formulates a mathematical model for the formation of ablation morphologies as a function of measurable parameters. The dependence of ablation morphologies on weather conditions and initial dirt thickness are studied, focusing on the initial growth of perturbations away from a flat surface. We derive a single-parameter expression for the melting rate as a function of dirt thickness, which agrees well with a set of measurements by Driedger. An interesting result is the prediction of a dirt-induced travelling instability for a range of parameters.

preprint1999arXiv

Collapsing Bacterial Cylinders

Under special conditions bacteria excrete an attractant and aggregate. The high density regions initially collapse into cylindrical structures, which subsequently destabilize and break up into spherical aggregates. This paper presents a theoretical description of the process, from the structure of the collapsing cylinder to the spacing of the final aggregates. We show that cylindrical collapse involves a delicate balance in which bacterial attraction and diffusion nearly cancel, leading to corrections to the collapse laws expected from dimensional analysis. The instability of a collapsing cylinder is composed of two distinct stages: Initially, slow modulations to the cylinder develop, which correspond to a variation of the collapse time along the cylinder axis. Ultimately, one point on the cylinder pinches off. At this final stage of the instability, a front propagates from the pinch into the remainder of the cylinder. The spacing of the resulting spherical aggregates is determined by the front propagation.

preprint1998arXiv

Electrostatic Edge Instability of Lipid Membranes

This paper considers the effect of electrostatics on the stability of a charged membrane. We show that at low ionic strength and high surface charge density, repulsion between charges on the membrane renders it unstable to the formation of holes. A free straight edge is unstable to modulations with wavelength longer than the Debye screening length. Hence at low ionic strength, membranes will disintegrate into vesicles. We use these results to interpret stable holes in red blood cell ghosts (Steck et. al., Science 168:255(1970)).