Publications

The organization of the cytokinetic ring at the cell equator of dividing animal and fungi cells depends crucially on the anillin scaffold proteins. In fission yeast, anillin related Mid1 binds to the plasma membrane and helps anchor and organize a medial broad band of cytokinetic nodes, which are the precursors of the contractile ring. Similar to other anillins, Mid1 contains a C terminal globular domain with two potential regions for membrane binding, the Pleckstrin Homology (PH) and C2 domains, and an N terminal intrinsically disordered region that is strongly regulated by phosphorylation. Previous studies have shown that both PH and C2 domains can associate with the membrane, preferring phosphatidylinositol-(4,5)-bisphosphate (PIP2) lipids. However, it is unclear if they can simultaneously bind to the membrane in a way that allows dimerization or oligomerization of Mid1, and if one domain plays a dominant role. To elucidate Mid1’s membrane binding mechanism, we used the available structural information of the C terminal region of Mid1 in all-atom molecular dynamics (MD) near a membrane with a lipid composition based on experimental measurements (including PIP2 lipids). The disordered L3 loop of C2, as well as the PH domain, separately bind the membrane through charged lipid contacts. In simulations with the full C terminal region started away from the membrane, Mid1 binds through the L3 loop and is stabilized in a vertical orientation with the PH domain away from the membrane. However, a configuration with both C2 and PH initially bound to the membrane remains associated with the membrane. These multiple modes of binding may reflect Mid1’s multiple interactions with membranes and other node proteins, and ability to sustain mechanical forces.

The F-BAR protein Cdc15 is essential for cytokinesis in Schizosaccharomyces pombe and plays a key role in attaching the cytokinetic ring (CR) to the plasma membrane (PM). Cdc15’s abilities to bind to the membrane and oligomerize via its F-BAR domain are inhibited by phosphorylation of its intrinsically disordered region (IDR). Multiple cell polarity kinases regulate Cdc15 IDR phosphostate, and of these the DYRK kinase Pom1 phosphorylation sites on Cdc15 have been shown in vivo to prevent CR formation at cell tips. Here, we compared the ability of Pom1 to control Cdc15 phosphostate and cortical localization to that of other Cdc15 kinases: Kin1, Pck1, and Shk1. We identified distinct but overlapping cohorts of Cdc15 phosphorylation sites targeted by each kinase, and the number of sites correlated with each kinases’ abilities to influence Cdc15 PM localization. Coarse-grained simulations predicted that cumulative IDR phosphorylation moves the IDRs of a dimer apart and toward the F-BAR tips. Further, simulations indicated that the overall negative charge of phosphorylation masks positively charged amino acids necessary for F-BAR oligomerization and membrane interaction. Finally, simulations suggested that dephosphorylated Cdc15 undergoes phase separation driven by IDR interactions. Indeed, dephosphorylated but not phosphorylated Cdc15 undergoes liquid–liquid phase separation to form droplets in vitro that recruit Cdc15 binding partners. In cells, Cdc15 phosphomutants also formed PM-bound condensates that recruit other CR components. Together, we propose that a threshold of Cdc15 phosphorylation by assorted kinases prevents Cdc15 condensation on the PM and antagonizes CR assembly.

The fluid–structure interactions between flexible fibres and viscous flows play an essential role in various biological phenomena, medical problems and industrial processes. Of particular interest is the case of particles transported freely in time-dependent flows. This work elucidates the dynamics and morphologies of actin filaments under oscillatory shear flows by combining microfluidic experiments, numerical simulations and theoretical modelling. Our work reveals that, in contrast to steady shear flows, in which small orientational fluctuations from a flow-aligned state initiate tumbling and deformations, the periodic flow reversal allows the filament to explore many different configurations at the beginning of each cycle. Investigation of filament motion during half time periods of oscillation highlights the critical role of the initial filament orientation on the emergent dynamics. This strong coupling between orientation and deformation results in new deformation regimes and novel higher-order buckling modes absent in steady shear flows. The primary outcome of our analysis is the possibility of suppression of buckling instabilities for certain combinations of the oscillation frequency and initial filament orientation, even in very strong flows. We explain this unusual behaviour through a weakly nonlinear Landau theory of buckling, in which we treat the filaments as inextensible Brownian Euler–Bernoulli rods whose hydrodynamics is described by local slender-body theory.

The nucleation of actin filament branches by the Arp2/3 complex involves activation through nucleation promotion factors (NPFs), recruitment of actin monomers, and binding of the complex to the side of actin filaments. Because of the large system size and processes that involve flexible regions and diffuse components, simulations of branch formation using all-atom molecular dynamics are challenging. We applied a coarse-grained model that retains amino-acid level information and allows molecular dynamics simulations in implicit solvent, with globular domains represented as rigid bodies and flexible regions allowed to fluctuate. We used recent electron microscopy structures of the inactive Arp2/3 complex bound to NPF domains and to mother actin filament for the activated Arp2/3 complex. We studied interactions of Arp2/3 complex with the activating VCA domain of the NPF Wiskott-Aldrich syndrome protein, actin monomers, and actin filament. We found stable configurations with one or two actin monomers bound along the branch filament direction and with CA domain of VCA associated to the strong and weak binding sites of the Arp2/3 complex, supporting prior structural studies and validating our approach. We reproduced delivery of actin monomers and CA to the Arp2/3 complex under different conditions, providing insight into mechanisms proposed in previous studies. Simulations of active Arp2/3 complex bound to a mother actin filament indicate the contribution of each subunit to the binding. Addition of the C-terminal tail of Arp2/3 complex subunit ArpC2, which is missing in the cryo-EM structure, increased binding affinity, indicating a possible stabilizing role of this tail.

Recent experiments have shown that the nematode {\it T. aceti} can assemble into collectively undulating groups at the edge of fluid drops. This coordinated state consists of metachronal waves and drives fluid circulation inside the drop. We find that the circulation velocity is about 2 mm/s and nearly half the speed of the metachronal wave. We develop a quasi two-dimensional hydrodynamics model using the Stokes flow approximation. The periodic motion of the nematodes constitute our moving boundary condition that drives the flow. Our model suggests that large amplitude excursions of the nematodes tails produce the fluid circulation. We discuss the constraints on containers that would enhance fluid motion, which could be used in the future design of on demand flow generating systems.

Inspired by the recent realization of a two-dimensional (2-D) chiral fluid as an active monolayer droplet moving atop a 3-D Stokesian fluid, we formulate mathematically its free-boundary dynamics. The surface droplet is described as a general 2-D linear, incompressible and isotropic fluid, having a viscous shear stress, an active chiral driving stress and a Hall stress allowed by the lack of time-reversal symmetry. The droplet interacts with itself through its driven internal mechanics and by driving flows in the underlying 3-D Stokes phase. We pose the dynamics as the solution to a singular integral–differential equation, over the droplet surface, using the mapping from surface stress to surface velocity for the 3-D Stokes equations. Specializing to the case of axisymmetric droplets, exact representations for the chiral surface flow are given in terms of solutions to a singular integral equation, solved using both analytical and numerical techniques. For a disc-shaped monolayer, we additionally employ a semi-analytical solution that hinges on an orthogonal basis of Bessel functions and allows for efficient computation of the monolayer velocity field, which ranges from a nearly solid-body rotation to a unidirectional edge current, depending on the subphase depth and the Saffman–Delbrück length. Except in the near-wall limit, these solutions have divergent surface shear stresses at droplet boundaries, a signature of systems with codimension-one domains embedded in a 3-D medium. We further investigate the effect of a Hall viscosity, which couples radial and transverse surface velocity components, on the dynamics of a closing cavity. Hall stresses are seen to drive inward radial motion, even in the absence of edge tension.

To prepare gametes with the appropriate number of chromosomes, mammalian oocytes undergo two sequential cell divisions. During each division, a large, long-lived, microtubule-based organelle called the meiotic spindle assembles around condensed chromosomes. Although meiotic spindles have been intensively studied for several decades, as force-generating mechanical objects, they remain very poorly understood. In materials physics, coarse-grained theories have been essential in understanding the large-scale behavior of systems composed of many interacting particles. It is unclear, however, if this approach can succeed in capturing the properties of active, biochemically complex, living materials like the spindle. Here, we show that a class of models based on nematic liquid crystal theory can describe important aspects of the organelle-scale structure and dynamics of spindles in living mouse oocytes. Using our models to interpret quantitative polarization microscopy data, we measure for the first time material properties relating to stress propagation in living oocytes, including the nematic diffusivities corresponding to splay and bend deformations. Unlike the reconstituted amphibian spindles that were previously studied in vitro, nematic elastic stress is exponentially screened in the microtubule network of living mammalian oocytes, with a screening length of order one micron. This observation can be explained by the relatively high volume fraction of embedded chromosomes in mammalian meiotic spindles, which cause long voids in the microtubule network and so disrupt orientational stress propagation.

We re-examine the model proposed by Alexander et al. (Phys. Fluids, vol. 18, 2006, 062103) for the closing of a circular hole in a molecularly thin incompressible Langmuir film situated on a Stokesian subfluid. For simplicity their model assumes that the surface phase is inviscid which leads to the result that the cavity area decreases at a constant rate determined by the ratio of edge tension to subfluid viscosity. We reformulate the problem, allowing for a regularising monolayer viscosity. The viscosity-dependent corrections to the hole dynamics are analysed and found to be non-trivial, even when the monolayer viscosity is small; these corrections may explain the departure of experimental data from the theoretical prediction when the hole radius becomes comparable to the Saffman–Delbrück length. Through fitting, under these relaxed assumptions, we find the edge tension could be as much as six times larger ( ∼
4.0 pN) than reported previously.

Cytosolic Ca2+ is a highly dynamic, tightly regulated and broadly conserved cellular signal. Ca2+ dynamics have been studied widely in cellular monocultures, yet organs in vivo comprise heterogeneous populations of stem and differentiated cells. Here, we examine Ca2+ dynamics in the adult Drosophila intestine, a self-renewing epithelial organ in which stem cells continuously produce daughters that differentiate into either enteroendocrine cells or enterocytes. Live imaging of whole organs ex vivo reveals that stem-cell daughters adopt strikingly distinct patterns of Ca2+ oscillations after differentiation: enteroendocrine cells exhibit single-cell Ca2+ oscillations, whereas enterocytes exhibit rhythmic, long-range Ca2+ waves. These multicellular waves do not propagate through immature progenitors (stem cells and enteroblasts), of which the oscillation frequency is approximately half that of enteroendocrine cells. Organ-scale inhibition of gap junctions eliminates Ca2+ oscillations in all cell types – even, intriguingly, in progenitor and enteroendocrine cells that are surrounded only by enterocytes. Our findings establish that cells adopt fate-specific modes of Ca2+ dynamics as they terminally differentiate and reveal that the oscillatory dynamics of different cell types in a single, coherent epithelium are paced independently.

Mixtures of filaments and molecular motors form active materials with diverse dynamical behaviors that vary based on their constituents’ molecular properties. To develop a multiscale of these materials, we map the nonequilibrium phase diagram of microtubules and tip-accumulating kinesin-4 molecular motors. We find that kinesin-4 can drive either global contractions or turbulent like extensile dynamics, depending on the concentrations of both microtubules and a bundling agent. We also observe a range of spatially heterogeneous nonequilibrium phases, including finite-sized radial asters, 1D wormlike chains, extended 2D bilayers, and system-spanning 3D active foams. Finally, we describe intricate kinetic pathways that yield microphase-separated structures and arise from the inherent frustration between the orientational order of filamentous microtubules and the positional order of tip-accumulating molecular motors. Our work reveals a range of novel active states. It also shows that the form of active stresses is not solely dictated by the properties of individual motors and filaments, but is also contingent on the constituent concentrations and spatial arrangement of motors on the filaments.