A PLETHORA-auxin transcription module controls cell division plane rotation through MAP65 and CLASP.

Despite their pivotal role in plant development, control mechanisms for oriented cell divisions have remained elusive. Here, we describe how a precisely regulated cell division orientation switch in an Arabidopsis stem cell is controlled by upstream patterning factors. We show that the stem cell regulatory PLETHORA transcription factors induce division plane reorientation by local activation of auxin signaling, culminating in enhanced expression of the microtubule-associated MAP65 proteins. MAP65 upregulation is sufficient to reorient the cortical microtubular array through a CLASP microtubule-cell cortex interaction mediator-dependent mechanism. CLASP differentially localizes to cell faces in a microtubule- and MAP65-dependent manner. Computational simulations clarify how precise 90° switches in cell division planes can follow self-organizing properties of the microtubule array in combination with biases in CLASP localization. Our work demonstrates how transcription factor-mediated processes regulate the cellular machinery to control orientation of formative cell divisions in plants.

Poroelasticity of (bio)polymer networks during compression: theory and experiment.

Soft living tissues like cartilage can be considered as biphasic materials comprising a fibrous complex biopolymer network and a viscous background liquid. Here, we show by a combination of experiment and theoretical analysis that both the hydraulic permeability and the elastic properties of (bio)polymer networks can be determined with simple ramp compression experiments in a commercial rheometer. In our approximate closed-form solution of the poroelastic equations of motion, we find the normal force response during compression as a combination of network stress and fluid pressure. Choosing fibrin as a biopolymer model system with controllable pore size, measurements of the full time-dependent normal force during compression are found to be in excellent agreement with the theoretical calculations. The inferred elastic response of large-pore (µm) fibrin networks depends on the strain rate, suggesting a strong interplay between network elasticity and fluid flow. Phenomenologically extending the calculated normal force into the regime of nonlinear elasticity, we find strain-stiffening of small-pore (sub-µm) fibrin networks to occur at an onset average tangential stress at the gel-plate interface that depends on the polymer concentration in a power-law fashion. The inferred permeability of small-pore fibrin networks scales approximately inverse squared with the fibrin concentration, implying with a microscopic cubic lattice model that the number of protofibrils per fibrin fiber cross-section decreases with protein concentration. Our theoretical model provides a new method to obtain the hydraulic permeability and the elastic properties of biopolymer networks and hydrogels with simple compression experiments, and paves the way to study the relation between fluid flow and elasticity in biopolymer networks during dynamical compression.

How selective severing by katanin promotes order in the plant cortical microtubule array.

Plant morphogenesis requires differential and often asymmetric growth. A key role in controlling anisotropic expansion of individual cells is played by the cortical microtubule array. Although highly organized, the array can nevertheless rapidly change in response to internal and external cues. Experiments have identified the microtubule-severing enzyme katanin as a central player in controlling the organizational state of the array. Katanin action is required both for normal alignment and the adaptation of array orientation to mechanical, environmental, and developmental stimuli. How katanin fulfills its controlling role, however, remains poorly understood. On the one hand, from a theoretical perspective, array ordering depends on the "weeding out" of discordant microtubules through frequent catastrophe-inducing collisions among microtubules. Severing would reduce average microtubule length and lifetime, and consequently weaken the driving force for alignment. On the other hand, it has been suggested that selective severing at microtubule crossovers could facilitate the removal of discordant microtubules. Here we show that this apparent conflict can be resolved by systematically dissecting the role of all of the relevant interactions in silico. This procedure allows the identification of the sufficient and necessary conditions for katanin to promote array alignment, stresses the critical importance of the experimentally observed selective severing of the "crossing" microtubule at crossovers, and reveals a hitherto not appreciated role for microtubule bundling. We show how understanding the underlying mechanism can aid with interpreting experimental results and designing future experiments.

A computational framework for cortical microtubule dynamics in realistically shaped plant cells.

Plant morphogenesis is strongly dependent on the directional growth and the subsequent oriented division of individual cells. It has been shown that the plant cortical microtubule array plays a key role in controlling both these processes. This ordered structure emerges as the collective result of stochastic interactions between large numbers of dynamic microtubules. To elucidate this complex self-organization process a number of analytical and computational approaches to study the dynamics of cortical microtubules have been proposed. To date, however, these models have been restricted to two dimensional planes or geometrically simple surfaces in three dimensions, which strongly limits their applicability as plant cells display a wide variety of shapes. This limitation is even more acute, as both local as well as global geometrical features of cells are expected to influence the overall organization of the array. Here we describe a framework for efficiently simulating microtubule dynamics on triangulated approximations of arbitrary three dimensional surfaces. This allows the study of microtubule array organization on realistic cell surfaces obtained by segmentation of microscopic images. We validate the framework against expected or known results for the spherical and cubical geometry. We then use it to systematically study the individual contributions of global geometry, cell-edge induced catastrophes and cell-face induced stability to array organization in a cuboidal geometry. Finally, we apply our framework to analyze the highly non-trivial geometry of leaf pavement cells of Arabidopsis thaliana, Nicotiana benthamiana and Hedera helix. We show that our simulations can predict multiple features of the microtubule array structure in these cells, revealing, among others, strong constraints on the orientation of division planes.

Confinement and crowding control the morphology and dynamics of a model bacterial chromosome.

Motivated by recent experiments probing the shape, size and dynamics of bacterial chromosomes in growing cells, we consider a polymer model consisting of a circular backbone to which side-loops are attached, confined to a cylindrical cell. Such a model chromosome spontaneously adopts a helical shape, which is further compacted by molecular crowders to occupy a nucleoid-like sub-volume of the cell. With increasing cell length, the longitudinal size of the chromosome increases in a non-linear fashion until finally saturating, its morphology gradually opening up while displaying a changing number of helical turns. For shorter cells, the chromosome extension varies non-monotonically with cell size, which we show is associated with a radial to longitudinal spatial reordering of the crowders. Confinement and crowders constrain chain dynamics leading to anomalous diffusion. While the scaling exponent for the mean squared displacement of center of mass grows and saturates with cell length, that of individual loci displays a broad distribution with a sharp maximum.

Gravity-driven syneresis in model low-fat mayonnaise.

Low-fat food products often contain natural, edible polymers to retain the desired mouth feel and elasticity of their full-fat counterparts. This type of product, however, can suffer from syneresis: densification due to the expulsion of fluid. Gaining insight into the physical principles governing syneresis in such soft hybrid dispersions remains a challenge from a theoretical perspective, as experimental data are needed to establish a basis. We record non-accelerated syneresis in a model system for low-fat mayonnaise: a colloid polymer mixture, consisting of oil in water emulsion with starch in the aqueous phase. We find the flow rate of expelled fluid to be proportional to the difference in hydrostatic pressure over the system. The osmotic pressure of the added starch, while being higher than the hydrostatic pressure, does not prevent syneresis because the soluble starch is lost to the expelled fluid. From these findings, we conclude that forced syneresis in these systems can be described as a gravity-driven porous flow through the densely packed emulsion, explainable with a model based on Darcy's law.

The Landau-de Gennes approach revisited: A minimal self-consistent microscopic theory for spatially inhomogeneous nematic liquid crystals.

We design a novel microscopic mean-field theory of inhomogeneous nematic liquid crystals formulated entirely in terms of the tensor order parameter field. It combines the virtues of the Landau-de Gennes approach in allowing both the direction and magnitude of the local order to vary, with a self-consistent treatment of the local free-energy valid beyond the small order parameter limit. As a proof of principle, we apply this theory to the well-studied problem of a colloid dispersed in a nematic liquid crystal by including a tunable wall coupling term. For the two-dimensional case, we investigate the organization of the liquid crystal and the position of the point defects as a function of the strength of the coupling constant.

Defect structures mediate the isotropic-nematic transition in strongly confined liquid crystals.

Using Monte Carlo simulations, we study rod-like lyotropic liquid crystals confined to a square slab-like geometry with lateral dimensions comparable to the length of the particles. We observe that this system develops linear defect structures upon entering the planar nematic phase. These defect structures flank a lens-shaped nematic region oriented along a diagonal of the square box. We interpret these structures as a compromise between the 2-fold order of the bulk nematic phase and the 4-fold order imposed by the lateral boundaries. A simple Onsager-type theory that effectively implements these competing tendencies is used to model the phase behavior in the center of the box and shows that the free-energy cost of forming the defect structures strongly offsets the transition-inducing effects of both the transverse and lateral confinement.

Cortical microtubule arrays are initiated from a nonrandom prepattern driven by atypical microtubule initiation.

The ordered arrangement of cortical microtubules in growing plant cells is essential for anisotropic cell expansion and, hence, for plant morphogenesis. These arrays are dismantled when the microtubule cytoskeleton is rearranged during mitosis and reassembled following completion of cytokinesis. The reassembly of the cortical array has often been considered as initiating from a state of randomness, from which order arises at least partly through self-organizing mechanisms. However, some studies have shown evidence for ordering at early stages of array assembly. To investigate how cortical arrays are initiated in higher plant cells, we performed live-cell imaging studies of cortical array assembly in tobacco (Nicotiana tabacum) Bright Yellow-2 cells after cytokinesis and drug-induced disassembly. We found that cortical arrays in both cases did not initiate randomly but with a significant overrepresentation of microtubules at diagonal angles with respect to the cell axis, which coincides with the predominant orientation of the microtubules before their disappearance from the cell cortex in preprophase. In Arabidopsis (Arabidopsis thaliana) root cells, recovery from drug-induced disassembly was also nonrandom and correlated with the organization of the previous array, although no diagonal bias was observed in these cells. Surprisingly, during initiation, only about one-half of the new microtubules were nucleated from locations marked by green fluorescent protein-γ-tubulin complex protein2-tagged γ-nucleation complexes (γ-tubulin ring complex), therefore indicating that a large proportion of early polymers was initiated by a noncanonical mechanism not involving γ-tubulin ring complex. Simulation studies indicate that the high rate of noncanonical initiation of new microtubules has the potential to accelerate the rate of array repopulation.

Alignment of nematic and bundled semiflexible polymers in cell-sized confinement.

The finite size of cells poses severe spatial constraints on the network of semiflexible filaments called the cytoskeleton, a main determinant of cell shape. At the same time, the high packing density of cytoskeletal filaments poses mutual packing constraints. Here we investigate the competition between excluded volume interactions in the bulk and surface packing constraints on the orientational ordering of confined actin filaments as a function of filament density and the presence of crosslinks. We grow fluorescently labeled actin filaments in shallow (thickness dz 3 µm), rectangular microchambers with a systematically varied length (dy between 5 and 100 µm) and in-plane aspect ratio (dx/dy between 1 and 10). We determine the nematic director field by image analysis of fluorescence confocal images. We find that high-density (nematic) solutions respond sensitively to changes in the size and aspect ratio of the chambers. In small chambers (dy ≤ 20 µm), filaments align parallel to the long walls as soon as the aspect ratio is ≥1.5, indicating that surface-induced ordering dominates. In larger chambers, the filaments instead align along the chamber diagonal, indicating that bulk packing constraints dominate. The nematic order parameter is maximal in small and highly anisometric chambers. In contrast to the nematic solutions, low-density (isotropic) solutions are rather insensitive to confinement. Bundled actin solutions behave similarly to nematic solutions, but are less well-ordered. Our observations imply that the orientational order of actin filaments in flat confining geometries is primarily determined by a balance between bulk and surface packing constraints with a minimal effect of the enthalpic cost of filament bending. Our assay provides an interesting platform for the future reconstitution of more complex, active cytoskeletal systems with actively treadmilling filaments or molecular motors.

Colloidal liquid crystals in rectangular confinement: theory and experiment.

We theoretically and experimentally study nematic liquid crystal equilibria within shallow rectangular wells. We model the wells within a two-dimensional Oseen-Frank framework, with strong tangent anchoring, and obtain explicit analytical expressions for the director fields and energies of the 'diagonal' and 'rotated' solutions reported in the literature. These expressions separate the leading-order defect energies from the bulk distortion energy for both families of solutions. The continuum Oseen-Frank study is complemented by a microscopic mean-field approach. We numerically minimize the mean-field functional, including the effects of weak anchoring, variable order and random initial conditions. In particular, these simulations suggest the existence of higher-energy metastable states with internal defects. We compare our theoretical results to experimental director profiles, obtained using two types of filamentous virus particles, wild-type fd-virus and a modified stiffer variant (Y21M), which display nematic ordering in rectangular chambers, as found by confocal scanning laser microscopy. We combine our analytical energy expressions with experimentally recorded frequencies of the different equilibrium states to obtain explicit estimates for the extrapolation length, defined to be the ratio of the nematic elastic constant to the anchoring coefficient, of the fd-virus.

The effect of anisotropic microtubule-bound nucleations on ordering in the plant cortical array.

The highly orientationally ordered cortical microtubule array in plant cells is a key component for cell growth and development. Recent experimental and computational work has shown that the anisotropic nucleation of new microtubules from pre-existing microtubules has a major effect on the alignment process. We formulate a theoretical model to investigate the role of the microtubule-bound nucleation on the self-organization of the dynamical cortical microtubules. A bifurcation analysis of the stability of the disordered phase of the model reveals that the effective degree of co-aligned nucleation is the main determinant of the location of the transition. Increased co-aligned nucleation creates a positive feedback effect on the ordering process that can significantly widen the ordered region. We validate these predictions by comparing to the results of particle-based simulations.

Molecular Dynamics Simulation of a Feather-Boa Model of a Bacterial Chromosome.

The chromosome of a bacterium consists of a mega-base pair-long circular DNA, which self-organizes within the micron-sized bacterial cell volume, compacting itself by three orders of magnitude. Unlike eukaryotic chromosomes, it lacks a nuclear membrane and freely floats in the cytosol confined by the cell membrane. It is believed that strong confinement, cross-linking by associated proteins, and molecular crowding all contribute to determine chromosome size and morphology. Modelling the chromosome simply as a circular polymer decorated with closed side loops in a cylindrical confining volume has been shown to already recapture some of the salient properties observed experimentally. Here we describe how a computer simulation can be set up to study structure and dynamics of bacterial chromosomes using this model.

A probabilistic algorithm for optimising the steady-state diffusional flux into a partially absorbing body.

Cells in an aqueous environment absorb diffusing nutrient molecules through nanoscale protein channels in their outer membranes. Assuming that there are constraints on the number of such channels a cell can produce, we ask the question: given a nondepleting source of nutrients, what is the optimal distribution of these channels over the cell surface? We coarse-grain this problem, phrasing it as a diffusion problem with position-dependent Robin boundary conditions on the surface. The aim is to maximize the steady-state total flux through the partially absorbing surface under an integral constraint on the local reactivities. We develop an algorithm to tackle this problem that uses the stored and processed results of a particle-based simulation with reflective boundary conditions to a posteriori estimate absorption flux at essentially negligible additional computational cost. We validate the algorithm against a few cases for which analytical or semi-analytical results are available. We apply it to two examples: a spherical cell in the presence of a point source and a spheroidal cell with an isotropic source at infinity. In the former case, there is a significant gain relative to the homogeneous case, while in the latter case the gain is only [Formula: see text].

Spontaneous helicity of a polymer with side loops confined to a cylinder.

Inspired by recent experiments on the spatial organization of bacterial chromosomes, we consider a type of "bottle-brush" polymer consisting of a flexible backbone chain, to which flexible side loops are connected. We show that such a model with an open linear backbone spontaneously adopts a helical structure with a well-defined pitch when confined to small cylindrical volume. This helicity persists over a range of sizes and aspect ratios of the cylinder, provided the packing fraction of the chain is suitably large. We analyze these results in terms of the interplay between the effective stiffness and actual intrachain packing effects caused by the side loops in response to the confinement.

Frustration-induced complexity in order-disorder transitions of the J_{1}-J_{2}-J_{3} Ising model on the square lattice.

We revisit the field-free Ising model on a square lattice with up to third-neighbor (NNNN) interactions, also known as the J_{1}-J_{2}-J_{3} model, in the mean-field approximation. Using a systematic enumeration procedure, we show that the region of phase space in which the high-temperature disordered phase is stable against all modes representing periodic magnetization patterns up to a given size is a convex polytope that can be obtained by solving a standard vertex enumeration problem. Each face of this polytope corresponds to a set of coupling constants for which a single set of modes, equivalent up to a symmetry of the lattice, bifurcates from the disordered solution. While the structure of this polytope is simple in the half-space J_{3}>0, where the NNNN interaction is ferromagnetic, it becomes increasingly complex in the half-space J_{3}<0, where the antiferromagnetic NNNN interaction induces strong frustration. We then pass to the limit N→∞ giving a closed-form description of the order-disorder surface in the thermodynamic limit, which shows that for J_{3}<0, the emergent ordered phases will have a "devil's surface"-like mode structure. Finally, using Monte Carlo simulations, we show that for small periodic systems, the mean-field analysis correctly predicts the dominant modes of the ordered phases that develop for coupling constants associated with the centroid of the faces of the disorder polytope.

Microtubule organization and cell geometry.

A characteristic feature of nondividing animal cells is the radial organization of microtubules (MTs), emanating from a single microtubule organizing center (MTOC). As generically these cells are not spherically symmetric, this raises the question of the influence of cell geometry on the orientational distribution of microtubules. We present a systematic study of this question in a simplified setting where MTs are nucleated from a single fixed MTOC in the center of an elliptical cell geometry. Within this context we introduce four models of increasing complexity, each one introducing additional mechanisms that govern the interaction of the MTs with the cell boundary. In order, we consider the cases: MTs that can bind to the boundary with a fixed mean residence time (M0), force-producing MTs that can slide on the boundary towards the cell poles (MS), MTs that interact with a generic polarity factor that is transported and deposited at the boundary, and which in turn stabilizes the MTs at the boundary (MP), and a final model in which both sliding and stabilization by polarity factors is taken into account (MSP). In the baseline model (M0), the exponential length distribution of MTs causes most of the interactions at the cell boundary to occur along the shorter transverse direction in the cell, leading to transverse biaxial order. MT sliding (MS) is able to reorient the main axis of this biaxial order along the longitudinal axis. The polarization mechanism introduced in MP and MSP overrules the geometric bias towards bipolar order observed in M0 and MS, and allows the establishment of unipolar order either along the short (MP) or the long cell axis (MSP). The behavior of the latter two models can be qualitatively reproduced by a very simple toy model with discrete MT orientations.

Modeling the molecular structures and dynamics responsible for the remarkable mechanical properties of a plant cell wall.

The primary plant cell wall is a hydrated meshwork of polysaccharides that is strong enough to withstand large mechanical stresses imposed by turgor while remaining pliant in ways that permit growth. To understand how its macromolecular architecture produces its complex mechanical properties, Zhang et al.1 computationally assembled a realistic network of cellulose microfibrils, hemicellulose, and pectin. The simulated wall responded to computationally applied stress like the real wall on which it was based. The model showed the location and chemical identity of stress-bearing components. It showed that cellulose microfibril interactions and movements dominated the wall's mechanical behavior, while hemicellulose and pectin had surprisingly minor effects.