Cargo navigation across 3D microtubule intersections.

The eukaryotic cell's microtubule cytoskeleton is a complex 3D filament network. Microtubules cross at a wide variety of separation distances and angles. Prior studies in vivo and in vitro suggest that cargo transport is affected by intersection geometry. However, geometric complexity is not yet widely appreciated as a regulatory factor in its own right, and mechanisms that underlie this mode of regulation are not well understood. We have used our recently reported 3D microtubule manipulation system to build filament crossings de novo in a purified in vitro environment and used them to assay kinesin-1-driven model cargo navigation. We found that 3D microtubule network geometry indeed significantly influences cargo routing, and in particular that it is possible to bias a cargo to pass or switch just by changing either filament spacing or angle. Furthermore, we captured our experimental results in a model which accounts for full 3D geometry, stochastic motion of the cargo and associated motors, as well as motor force production and force-dependent behavior. We used a combination of experimental and theoretical analysis to establish the detailed mechanisms underlying cargo navigation at microtubule crossings.

Actin-myosin spatial patterns from a simplified isotropic viscoelastic model.

F-actin networks are involved in cell mechanical processes ranging from motility to endocytosis. The mesoscale architecture of assemblies of individual F-actin polymers that gives rise to micrometer-scale rheological properties is poorly understood, despite numerous in vivo and vitro studies. In vitro networks have been shown to organize into spatial patterns when spatially confined, including dense spherical shells inside spherical emulsion droplets. Here we develop a simplified model of an isotropic, compressible, viscoelastic material continually assembling and disassembling. We demonstrate that spherical shells emerge naturally when the strain relaxation rate (corresponding to internal network reorganization) is slower than the disassembly rate (corresponding to F-actin depolymerization). These patterns are consistent with recent experiments, including a collapse of shells to a central high-density focus of F-actin when either assembly or disassembly is reduced with drugs. Our results demonstrate how complex spatio-temporal patterns can emerge without spatially distributed force generation, polar alignment of F-actin polymers, or spatially nonuniform regulation of F-actin by upstream biochemical networks.

Estimating the bending modulus of a FtsZ bacterial-division protein filament.

FtsZ, a cytoskeletal protein homologous to tubulin, is the principle constituent of the division ring in bacterial cells. It is known to have force-generating capacity in vitro and has been conjectured to be the source of the constriction force in vivo. Several models have been proposed to explain the generation of force by the Z ring. Here we re-examine data from in vitro experiments in which Z rings formed and constricted inside tubular liposomes, and we carry out image analysis on previously published data with which to better estimate important model parameters that have proven difficult to measure by direct means. We introduce a membrane-energy-based model for the dynamics of multiple Z rings moving and colliding inside a tubular liposome and a fluid model for the drag of a Z ring as it moves through the tube. Using this model, we estimate an effective membrane bending modulus of 500-700 pN nm. If we assume that FtsZ force generation is driven by hydrolysis into a highly curved conformation, we estimate the FtsZ filament bending modulus to be 310-390 pN nm(2). If we assume instead that force is generated by the non-hydrolysis-dependent intermediate curvature conformation, we find that B(f)>1400 pN nm(2). The former value sits at the lower end of the range of previously estimated values and, if correct, may raise challenges for models that rely on filament bending to generate force.

Mechanical modulation of receptor-ligand interactions at cell-cell interfaces.

Cell surface receptors have been extensively studied because they initiate and regulate signal transduction cascades leading to a variety of functional cellular outcomes. An important class of immune receptors (e.g., T-cell antigen receptors) whose ligands are anchored to the surfaces of other cells remain poorly understood. The mechanism by which ligand binding initiates receptor phosphorylation, a process termed "receptor triggering", remains controversial. Recently, direct measurements of the (two-dimensional) receptor-ligand complex lifetimes at cell-cell interface were found to be smaller than (three-dimensional) lifetimes in solution but the underlying mechanism is unknown. At the cell-cell interface, the receptor-ligand complex spans a short intermembrane distance (15 nm) compared to long surface molecules (LSMs) whose ectodomains span >40 nm and these LSMs include phosphatases (e.g., CD45) that dephosphorylate the receptor. It has been proposed that size-based segregation of LSMs from a receptor-ligand complex is a mechanism of receptor triggering but it is unclear whether the mechanochemistry supports such small-scale segregation. Here we present a nanometer-scale mathematical model that couples membrane elasticity with the compressional stiffness and lateral mobility of LSMs. We find robust supradiffusive segregation of LSMs from a single receptor-ligand complex. The model predicts that LSM redistribution will result in a time-dependent tension on the complex leading to a decreased two-dimensional lifetime. Interestingly, the model predicts a nonlinear relationship between the three- and two-dimensional lifetimes, which can enhance the ability of receptors to discriminate between similar ligands.

A CLASP-modulated cell edge barrier mechanism drives cell-wide cortical microtubule organization in Arabidopsis.

It is well known that the parallel order of microtubules in the plant cell cortex defines the direction of cell expansion, yet it remains unclear how microtubule orientation is controlled, especially on a cell-wide basis. Here we show through 4D imaging and computational modelling that plant cell polyhedral geometry provides spatial input that determines array orientation and heterogeneity. Microtubules depolymerize when encountering sharp cell edges head-on, whereas those oriented parallel to those sharp edges remain. Edge-induced microtubule depolymerization, however, is overcome by the microtubule-associated protein CLASP, which accumulates at specific cell edges, enables microtubule growth around sharp edges and promotes formation of microtubule bundles that span adjacent cell faces. By computationally modelling dynamic 'microtubules on a cube' with edges differentially permissive to microtubule passage, we show that the CLASP-edge complex is a 'tuneable' microtubule organizer, with the inherent flexibility to generate the numerous cortical array patterns observed in nature.

A mechanochemical model explains interactions between cortical microtubules in plants.

Microtubules anchored to the two-dimensional cortex of plant cells collide through plus-end polymerization. Collisions can result in rapid depolymerization, directional plus-end entrainment, or crossover. These interactions are believed to give rise to cellwide self-organization of plant cortical microtubules arrays, which is required for proper cell wall growth. Although the cell-wide self-organization has been well studied, less emphasis has been placed on explaining the interactions mechanistically from the molecular scale. Here we present a model for microtubule-cortex anchoring and collision-based interactions between microtubules, based on a competition between cross-linker bonding, microtubule bending, and microtubule polymerization. Our model predicts a higher probability of entrainment at smaller collision angles and at longer unanchored lengths of plus-ends. This model addresses observed differences between collision resolutions in various cell types, including Arabidopsis cells and Tobacco cells.

Mechanisms of self-organization of cortical microtubules in plants revealed by computational simulations.

Microtubules confined to the two-dimensional cortex of elongating plant cells must form a parallel yet dispersed array transverse to the elongation axis for proper cell wall expansion. Some of these microtubules exhibit free minus-ends, leading to migration at the cortex by hybrid treadmilling. Collisions between microtubules can result in plus-end entrainment ("zippering") or rapid depolymerization. Here, we present a computational model of cortical microtubule organization. We find that plus-end entrainment leads to self-organization of microtubules into parallel arrays, whereas catastrophe-inducing collisions do not. Catastrophe-inducing boundaries (e.g., upper and lower cross-walls) can tune the orientation of an ordered array to a direction transverse to elongation. We also find that changes in dynamic instability parameters, such as in mor1-1 mutants, can impede self-organization, in agreement with experimental data. Increased entrainment, as seen in clasp-1 mutants, conserves self-organization, but delays its onset and fails to demonstrate increased ordering. We find that branched nucleation at acute angles off existing microtubules results in distinctive sparse arrays and infer either that microtubule-independent or coparallel nucleation must dominate. Our simulations lead to several testable predictions, including the effects of reduced microtubule severing in katanin mutants.

Pulling helices inside bacteria: imperfect helices and rings.

We study steady-state configurations of intrinsically-straight elastic filaments constrained within rod-shaped bacteria that have applied forces distributed along their length. Perfect steady-state helices result from axial or azimuthal forces applied at filament ends, however azimuthal forces are required for the small pitches observed for MreB filaments within bacteria. Helix-like configurations can result from distributed forces, including coexistence between rings and imperfect helices. Levels of expression and/or bundling of the polymeric protein could mediate this coexistence.

Force generation by a dynamic Z-ring in Escherichia coli cell division.

FtsZ, a bacterial homologue of tubulin, plays a central role in bacterial cell division. It is the first of many proteins recruited to the division site to form the Z-ring, a dynamic structure that recycles on the time scale of seconds and is required for division to proceed. FtsZ has been recently shown to form rings inside tubular liposomes and to constrict the liposome membrane without the presence of other proteins, particularly molecular motors that appear to be absent from the bacterial proteome. Here, we propose a mathematical model for the dynamic turnover of the Z-ring and for its ability to generate a constriction force. Force generation is assumed to derive from GTP hydrolysis, which is known to induce curvature in FtsZ filaments. We find that this transition to a curved state is capable of generating a sufficient force to drive cell-wall invagination in vivo and can also explain the constriction seen in the in vitro liposome experiments. Our observations resolve the question of how FtsZ might accomplish cell division despite the highly dynamic nature of the Z-ring and the lack of molecular motors.

Heterocyst patterns without patterning proteins in cyanobacterial filaments.

We have quantitatively modeled heterocyst differentiation after fixed nitrogen step-down in the filamentous cyanobacterium Anabaena sp. PCC 7120 without lateral inhibition due to the patterning proteins PatS or HetN. We use cell growth and division together with fixed-nitrogen dynamics and allow heterocysts to differentiate upon the local exhaustion of available fixed nitrogen. Slow transport of fixed nitrogen along a shared periplasmic space allows for fast growing cells to differentiate ahead of their neighbors. Cell-to-cell variability in growth rate determines the initial heterocyst pattern. Early release of fixed nitrogen from committed heterocysts allows a significant fraction of vegetative cells to be retained at later times. We recover the experimental heterocyst spacing distributions and cluster size distributions of Khudyakov and Golden [Khudyakov, I.Y., Golden, J.W., 2004. Different functions of HetR, a master regulator of heterocyst differentiation in Anabaena sp PCC 7120, can be separated by mutation. Proc. Natl. Acad. Sci. U. S. A. 101, 16040-16045].

Steady-state helices of the actin homolog MreB inside bacteria: dynamics without motors.

Within individual bacteria, we combine force-dependent polymerization dynamics of individual MreB protofilaments with an elastic model of protofilament bundles buckled into helical configurations. We use variational techniques and stochastic simulations to relate the pitch of the MreB helix, the total abundance of MreB, and the number of protofilaments. By comparing our simulations with mean-field calculations, we find that stochastic fluctuations are significant. We examine the quasistatic evolution of the helical pitch with cell growth, as well as time scales of helix turnover and de novo establishment. We find that while the body of a polarized MreB helix treadmills toward its slow-growing end, the fast-growing tips of laterally associated protofilaments move toward the opposite fast-growing end of the MreB helix. This offers a possible mechanism for targeted polar localization without cytoplasmic motor proteins.