Efficient Multiscale Models of Polymer Assembly

Protein polymerization and bundling play a central role in cell physiology. Predictive modeling of these processes remains an open challenge, especially when the proteins involved become large and their concentrations high. We present an effective kinetics model of filament formation, bundling, and depolymerization after GTP hydrolysis, which involves a relatively small number of species and reactions, and remains robust over a wide range of concentrations and timescales. We apply this general model to study assembly of FtsZ protein, a basic element in the division process of prokaryotic cells such as Escherichia coli, Bacillus subtilis, or Caulobacter crescentus. This analysis demonstrates that our model outperforms its counterparts in terms of both accuracy and computational efficiency. Because our model comprises only 17 ordinary differential equations, its computational cost is orders-of-magnitude smaller than the current alternatives consisting of up to 1000 ordinary differential equations. It also provides, to our knowledge, a new insight into the characteristics and functioning of FtsZ proteins at high concentrations. The simplicity and versatility of our model render it a powerful computational tool, which can be used either as a standalone descriptor of other biopolymers' assembly or as a component in more complete kinetic models.