Survival of the aligned: ordering of the plant cortical microtubule array

The cortical array is a structure consisting of highly aligned microtubules which plays a crucial role in the characteristic uniaxial expansion of all growing plant cells. Recent experiments have shown polymerization-driven collisions between the membrane-bound cortical microtubules, suggesting a possible mechanism for their alignment. We present both a coarse-grained theoretical model and stochastic particle-based simulations of this mechanism, and compare the results from these complementary approaches. Our results indicate that collisions that induce depolymerization are sufficient to generate the alignment of microtubules in the cortical array.

💡 Research Summary

The cortical microtubule array is a hallmark of plant cell morphogenesis, providing a scaffold that directs the uniaxial expansion of the cell wall. Although it has long been known that cortical microtubules are membrane‑bound, grow by polymerization, and frequently collide, the physical basis for their striking alignment has remained elusive. In this study the authors combine a coarse‑grained continuum theory with stochastic particle‑based simulations to test the hypothesis that collision‑induced depolymerization (catastrophe) is sufficient to generate the observed order.

The continuum model treats the microtubule population as a density field ρ(x,θ) that depends on position x and orientation θ. Growth is modeled as a constant elongation speed v, while collisions occur with a rate proportional to the product of local densities and the sine of the intersecting angle. Upon collision three outcomes are possible: (i) zippering, where both filaments continue growing in a common direction; (ii) induced reorientation, where one filament rotates to align with the other; and (iii) catastrophe, where one filament undergoes rapid depolymerization and disappears. The key control parameter is the probability pc of catastrophe per collision. By inserting these rules into a kinetic equation for ρ, the authors derive a nonlinear partial differential equation that predicts a nematic order parameter S(t). Linear stability analysis shows that for pc above a critical value pc* the isotropic state becomes unstable and a nematically ordered state emerges. Numerical integration confirms a sharp transition in S as pc is increased, reproducing the high alignment measured in vivo.

To validate the continuum predictions, the authors develop an explicit Monte‑Carlo simulation in which individual microtubules are represented as line segments that polymerize at speed v, undergo stochastic nucleation, and interact via the same three‑outcome collision rules. Collision detection is performed using a distance‑and‑angle criterion; when a collision is detected, a random draw determines which of the three outcomes occurs. The simulation tracks filament length, orientation, and lifetime over many growth cycles. Parameter values (v≈1 µm min⁻¹, average initial length ≈5 µm, nucleation rate, and collision frequency) are chosen to match published experimental measurements. The results mirror the continuum theory: when pc ≥ 0.3 the system rapidly self‑organizes into a highly aligned array, with the nematic order parameter rising from near zero to >0.7. Moreover, the length distribution collapses toward longer, more persistent filaments, while short filaments are eliminated by catastrophe, providing a mechanistic explanation for the observed selection of long, parallel microtubules in the cortex.

Comparison with experimental data shows that the model reproduces both the average microtubule length and the degree of alignment reported for Arabidopsis hypocotyl cells. The best quantitative agreement is obtained for pc≈0.4, suggesting that in vivo collisions frequently trigger depolymerization of one partner. This supports the view that collision‑induced catastrophe is not a pathological side‑effect but a functional driver of cortical array ordering.

The authors conclude that (1) microtubule collisions, when coupled to a realistic probability of catastrophe, are sufficient to generate the high degree of alignment seen in plant cortical arrays; (2) the agreement between continuum and particle‑based approaches demonstrates that the essential physics can be captured at multiple scales; and (3) this mechanism provides a flexible means for cells to remodel their cortical array in response to developmental cues or environmental stresses, because altering the catastrophe probability (e.g., via microtubule‑associated proteins) can shift the system between isotropic and ordered states. Future work is suggested to integrate this ordering mechanism with cell‑wall synthesis dynamics and to explore how external signals modulate the catastrophe probability, thereby linking microtubule physics to the broader context of plant morphogenesis.