Predictions from a stochastic polymer model for the MinDE dynamics in E.coli

The spatiotemporal oscillations of the Min proteins in the bacterium Escherichia coli play an important role in cell division. A number of different models have been proposed to explain the dynamics from the underlying biochemistry. Here, we extend a previously described discrete polymer model from a deterministic to a stochastic formulation. We express the stochastic evolution of the oscillatory system as a map from the probability distribution of maximum polymer length in one period of the oscillation to the probability distribution of maximum polymer length half a period later and solve for the fixed point of the map with a combined analytical and numerical technique. This solution gives a theoretical prediction of the distributions of both lengths of the polar MinD zones and periods of oscillations – both of which are experimentally measurable. The model provides an interesting example of a stochastic hybrid system that is, in some limits, analytically tractable.

💡 Research Summary

The paper addresses the well‑known oscillatory behavior of the Min proteins (MinD and MinE) in Escherichia coli, which is essential for positioning the division septum at mid‑cell. While many deterministic reaction‑diffusion or polymerization models have been proposed, they generally predict only average quantities and fail to capture the experimentally observed variability in the length of the polar MinD zones and the period of oscillation. To overcome this limitation, the authors extend a previously introduced discrete polymer model into a fully stochastic framework, treating the system as a hybrid of continuous growth dynamics and discrete disassembly events.

In the stochastic formulation, each MinD polymer is characterized by its length (L). The key observable is the maximum polymer length reached during a half‑cycle, denoted (L_{\max}). The authors describe the evolution of the probability distribution (P(L)) of this quantity over one full oscillation as a mapping operator (\mathcal{M}). The operator incorporates (i) deterministic elongation at a constant rate (v_g), which translates into a continuous probability flux, and (ii) stochastic disassembly driven by MinE, modeled as a Poisson jump process with rate (k_{\text{off}}). By applying (\mathcal{M}) twice—once for each pole—the distribution after a full period is obtained: (P_{n+1}(L)=\mathcal{M}